Après mon merge héroïque
[valse.git] / reports / essai16mars.R
CommitLineData
3f62d540 1p = 10
64b28e3e 2q = 8
3f62d540 3k = 2
4D = 20
5
f72201cc
BA
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)
f227455a 13
3f62d540 14covY = array(dim = c(q,q,k))
15covY[,,1] = 0.1*diag(q)
16covY[,,2] = 0.2*diag(q)
f227455a 17
3f62d540 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)))
f227455a 21
3f62d540 22n = 100
f227455a 23
24pi = c(0.4,0.6)
25
f227455a 26data = generateXY(meanX,covX,covY, pi, beta, n)
27
28X = data$X
29Y = data$Y
30
3f62d540 31res_valse = valse(X,Y)