1 basic_Init_Parameters = function(n,p,m,k){
2 phiInit = array(0, dim=c(p,m,k))
4 piInit = (1.0/k)*rep.int(1,k)
6 rhoInit = array(0, dim=c(m,m,k))
12 gamInit = 0.1*array(1, dim=c(n,k))
14 R = sample(1:k,n, replace= TRUE)
19 gamInit = gamInit/sum(gamInit[1,])
22 return(list(phiInit, rhoInit, piInit, gamInit))