Learn R Programming

BayesfMRI (version 0.3.11)

GLMEM_fixptseparate: Fixed point function for the joint BayesGLMEM update algorithm

Description

Fixed point function for the joint BayesGLMEM update algorithm

Usage

GLMEM_fixptseparate(theta, spde, model_data, Psi, K, A, cl, Ns = 50)

Value

a vector with the same length as theta, the EM updates

Arguments

theta

a list containing kappa2, phi, and sigma2, in that order

spde

the spde object

model_data

the model_data object containing y and X

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.