Description Usage Arguments Details Value Author(s) References See Also Examples

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

1 | ```
objfnA_CRD(des,ntmt,sige,means,probs=c(1))
``` |

`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. |

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`

.

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

Stephen Bush ([email protected])

Katya Ruggiero ([email protected])

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

`findOptimalExactDesign`

, `updateDesign_CRD`

1 2 3 4 5 | ```
# Finding the A-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)
objfnA_CRD(c(1,1,1,1,2,2,2,3,3,3),ntmt=3,sige=0,means=c(1,2,4))
``` |

designGLMM documentation built on May 30, 2017, 6:24 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.