Rename main function into runValse, remove testthat folder since nobody's gonna write...
[valse.git] / reports / essai16mars.R
... / ...
CommitLineData
1p = 10
2q = 8
3k = 2
4D = 20
5
6meanX = matrix(nrow=p,ncol=k)
7meanX[,1] = rep(0,p)
8meanX[,2] = rep(1,p)
9
10covX = array(dim=c(p,p,k))
11covX[,,1] = 0.1*diag(p)
12covX[,,2] = 0.5*diag(p)
13
14covY = array(dim = c(q,q,k))
15covY[,,1] = 0.1*diag(q)
16covY[,,2] = 0.2*diag(q)
17
18beta = array(dim = c(p,q,2))
19beta[,,2] = matrix(c(rep(2,(D)),rep(0, p*q-D)))
20beta[,,1] = matrix(c(rep(1,D),rep(0, p*q-D)))
21
22n = 100
23
24pi = c(0.4,0.6)
25
26data = generateXY(meanX,covX,covY, pi, beta, n)
27
28X = data$X
29Y = data$Y
30
31res_valse = valse(X,Y)