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

View source: R/EM_utils.R

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


mandymejia/BayesfMRI documentation built on April 12, 2025, 3:44 p.m.