GLMEM_fixptseparate: Fixed point function for the joint BayesGLM EM update...

In BayesfMRI: Spatial Bayesian Methods for Task Functional MRI Studies

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 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.

Value

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

BayesfMRI documentation built on April 4, 2025, 1:58 a.m.