R/winit.R

In mmtfa: Model-Based Clustering and Classification with Mixtures of Modified t Factor Analyzers

winit <-
function(x,n,G,mug,sigmainv,vg,p,sg,zmat){
	w <- matrix(0,n,G)
	delta <- matrix(0,n,G)
	for(g in 1:G){
		delta[,g] <- mahalanobis(x, mug[g,], sigmainv[,,g], inverted=TRUE)
		w[,g] <- (vg[g]+p)/(vg[g]+delta[,g])
	}
	w
}