Commit | Line | Data |
---|---|---|
39046da6 BA |
1 | vec_bin = function(X,r) |
2 | { | |
3 | Z = c() | |
4 | indice = c() | |
5 | j = 1 | |
6 | for (i in 1:length(X)) | |
7 | { | |
8 | if(X[i] == r) | |
9 | { | |
10 | Z[i] = 1 | |
11 | indice[j] = i | |
12 | j=j+1 | |
13 | } else | |
14 | Z[i] = 0 | |
15 | } | |
16 | return (list(Z=Z,indice=indice)) | |
35b42a4b B |
17 | } |
18 | ||
39046da6 BA |
19 | initSmallEM = function(k,X,Y,tau) |
20 | { | |
21 | n = nrow(Y) | |
22 | m = ncol(Y) | |
23 | p = ncol(X) | |
35b42a4b | 24 | |
39046da6 BA |
25 | betaInit1 = array(0, dim=c(p,m,k,20)) |
26 | sigmaInit1 = array(0, dim = c(m,m,k,20)) | |
27 | phiInit1 = array(0, dim = c(p,m,k,20)) | |
28 | rhoInit1 = array(0, dim = c(m,m,k,20)) | |
29 | piInit1 = matrix(0,20,k) | |
30 | gamInit1 = array(0, dim=c(n,k,20)) | |
31 | LLFinit1 = list() | |
32 | ||
33 | require(MASS) #Moore-Penrose generalized inverse of matrix | |
34 | for(repet in 1:20) | |
35 | { | |
36 | clusters = hclust(dist(y)) #default distance : euclidean | |
37 |