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).
Arguments
- fit
A
bifa_fitfromfit_bmsfa()or a rawMSFA::sp_msfa()fit.