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