Skip to contents

OP on Phi and each Lambda_s, form SigmaPhi, SigmaLambdaList, estimate Psi_s by averaging fit$psi[[s]] over iterations, and build marginal covariance as SigmaPhi + SigmaLambda_s + diag(Psi_s).

Usage

postprocess_bmsfa(fit)

Arguments

fit

A bifa_fit from fit_bmsfa() or a raw MSFA::sp_msfa() fit.

Value

A list with Phi, SigmaPhi, LambdaList, SigmaLambdaList, PsiList, SigmaMarginal.