R/sigmaup.R

sigmaup <-
function(p,G,lg,yg,sigma){
	for(g in 1:G){
		sigma[,,g] <- lg[,,g] %*% t(lg[,,g]) + yg[,,g]
	}
	sigma
}

Try the mmtfa package in your browser

Any scripts or data that you put into this service are public.

mmtfa documentation built on May 2, 2019, 1:45 p.m.