objfnD_CRD: Calculate D-optimality criterion for Completely Randomised...

In sabush/designGLMM: Finding Optimal Block Designs for a Generalised Linear Mixed Model

Description

This function calculates the determinant of the Fisher information matrix for a completely randomised design with a Poisson response.

Usage

objfnD_CRD(des, ntmt, sige, means, probs=c(1))

Arguments

des	The design under consideration, a list of length b x k containing treatment indices (1, 2, ..., v)
ntmt	The number of distinct treatments (v)
sige	Standard deviation of excess error
means	A list of length v containing conditional means for each treatment e.g. c(1,1,2) for three treatments with means 1, 1, and 2 respectively
probs	a list of probabilities specifying the probability that each step of the simulated annealing substitutes a certain number of design points. The first entry corresponds to the probability that only one substitution is made in a simulated annealing step, the second is the probability that two substitutions are made and so on. By default this is set to c(1) which means that only one substitution is made in each simulated annealing step.

Details

This function is designed to work with findOptimalExactDesign, and as such shares the arguments of updateDesign_CRD. It can, however, be used on its own. The probs argument is not used in this function, but is in updateDesign_CRD.

Value

Returns the negative of the determinant of the Fisher information matrix for the proveded design.

Author(s)

Stephen Bush ([e-mail])

Katya Ruggiero ([e-mail])

References

Bush, S., and Ruggiero, K. (2016) Optimal block designs for experiments with responses drawn from a Poisson distribution, Under Review, preprint available at http://arxiv.org/abs/1601.00477

See Also

findOptimalExactDesign, updateDesign_CRD

Examples

# Finding the D-optimality objective value for the design (1,1,1,1,2,2,2,3,3,3)
# where there are three treatments, the treatment means are 1, 2, and 4, and
# there is no overdispersion (sige=0)

objfnD_CRD(c(1,1,1,1,2,2,2,3,3,3),ntmt=3,sige=0,means=c(1,2,4))

sabush/designGLMM documentation built on May 29, 2019, 12:21 p.m.