update the model selection step. Beginning for the plots
[valse.git] / reports / essai16mars.R
1 p = 10
2 q = 15
3 k = 2
4 D = 20
5
6 meanX = rep(0,p)
7 covX = 0.1*diag(p)
8
9 covY = array(dim = c(q,q,k))
10 covY[,,1] = 0.1*diag(q)
11 covY[,,2] = 0.2*diag(q)
12
13 beta = array(dim = c(p,q,2))
14 beta[,,2] = matrix(c(rep(2,(D)),rep(0, p*q-D)))
15 beta[,,1] = matrix(c(rep(1,D),rep(0, p*q-D)))
16
17 n = 100
18
19 pi = c(0.4,0.6)
20
21 data = generateXY(meanX,covX,covY, pi, beta, n)
22
23 X = data$X
24 Y = data$Y
25
26 res_valse = valse(X,Y)