Fixed point function for the joint BayesGLMEM update algorithm
GLMEM_fixptseparate(theta, spde, model_data, Psi, K, A, cl, Ns = 50)
a vector with the same length as theta
, the EM updates
a list containing kappa2, phi, and sigma2, in that order
the spde object
the model_data object containing y
and X
a conversion matrix (N by V) (or N by n)
number of covariates
The value for Matrix::crossprod(X%*%Psi) (saves time on computation)
parallelization cluster
The number of samples used to approximate traces using the Hutchinson estimator. If set to 0, the exact trace is found.