This is the hyper-likelihood for updating the covariances. It is always based on mvtnorm::dmvnorm. The par is a matrix of parameters, one row per random effect, one column the set of parameters. It allows the Gibbs sampler to work by passing a single scalar parameter as the first argument.
1 2 | lmerBayes.hyperllike.sigma(testcov, fullcov, hypermean, modelpar, whichrow,
whichcol)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.