A function that calculates the Monte Carlo standard error for the Monte Carlo maximum likelihood estimates returned from `glmm`

.

1 | ```
mcse(object)
``` |

`object` |
An object of class |

With maximum likelihood performed by Monte Carlo likelihood approximation, there are two sources of variability: there is variability from sample to sample and from Monte Carlo sample (of generated random effects) to Monte Carlo sample. The first source of variability (from sample to sample) is measured using standard error, which appears with the point estimates in the `summary`

tables. The second source of variability is due to the Monte Carlo randomness, and this is measured by the Monte Carlo standard error.

A large Monte Carlo standard error indicates the Monte Carlo sample size `m`

is too small.

`mcse ` |
The Monte Carlo standard errors for the Monte Carlo maximum likelihood estimates returned from |

Christina Knudson

Geyer, C. J. (1994). On the convergence of Monte Carlo maximum likelihood calculations. Journal of the Royal Statistical Society, Series B, 61:261–274.

Knudson, C. (2016). Monte Carlo likelihood approximation for generalized linear mixed models. PhD thesis, University of Minnesota.

`glmm`

for model fitting.

1 2 3 4 5 6 7 8 9 | ```
library(glmm)
set.seed(1234)
data(salamander)
#To get more accurate answers for this model, use m=10^4 or 10^5
# and doPQL=TRUE.
m<-10
sal<-glmm(Mate~0+Cross,random=list(~0+Female,~0+Male),varcomps.names=c("F","M"),
data=salamander,family.glmm=bernoulli.glmm,m=m,debug=TRUE,doPQL=FALSE)
mcse(sal)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.

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.