Rename EMGLLF and EMGrank with suffix _R in generate_tests/.
[valse.git] / reports / essai16mars.R
... / ...
CommitLineData
1p = 10
2q = 8
3k = 2
4D = 20
5
6meanX = rep(0,p)
7covX = 0.1*diag(p)
8
9covY = array(dim = c(q,q,k))
10covY[,,1] = 0.1*diag(q)
11covY[,,2] = 0.2*diag(q)
12
13beta = array(dim = c(p,q,2))
14beta[,,2] = matrix(c(rep(2,(D)),rep(0, p*q-D)))
15beta[,,1] = matrix(c(rep(1,D),rep(0, p*q-D)))
16
17n = 100
18
19pi = c(0.4,0.6)
20
21data = generateXY(meanX,covX,covY, pi, beta, n)
22
23X = data$X
24Y = data$Y
25
26res_valse = valse(X,Y)