X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=R%2FbasicInitParameters.R;h=6090d0aa775d393e85d9856602742f826710ae1c;hp=bc7e88bd2b3f0431ce2b827d252ea844fabc1a91;hb=39046da6016f15d625bd99cf0303ea8beb838c79;hpb=8266149c7d93aa0543cee2a1b22e1233e7b82617 diff --git a/R/basicInitParameters.R b/R/basicInitParameters.R index bc7e88b..6090d0a 100644 --- a/R/basicInitParameters.R +++ b/R/basicInitParameters.R @@ -1,23 +1,18 @@ -basic_Init_Parameters = function(n,p,m,k){ - phiInit = array(0, dim=c(p,m,k)) - - piInit = (1.0/k)*rep.int(1,k) - - rhoInit = array(0, dim=c(m,m,k)) - - for(i in 1:k){ - rhoInit[,,i] = diag(m) - } - - gamInit = 0.1*array(1, dim=c(n,k)) - - R = sample(1:k,n, replace= TRUE) - - for(i in 1:n){ - gamInit[i,R[i]] = 0.9 - } - gamInit = gamInit/sum(gamInit[1,]) - - - return(list(phiInit, rhoInit, piInit, gamInit)) +basic_Init_Parameters = function(n,p,m,k) +{ + phiInit = array(0, dim=c(p,m,k)) + + piInit = (1./k)*rep.int(1,k) + + rhoInit = array(0, dim=c(m,m,k)) + for(i in 1:k) + rhoInit[,,i] = diag(m) + + gamInit = 0.1*array(1, dim=c(n,k)) + R = sample(1:k,n, replace=TRUE) + for(i in 1:n) + gamInit[i,R[i]] = 0.9 + gamInit = gamInit/sum(gamInit[1,]) + + return (list(phiInit, rhoInit, piInit, gamInit)) }