3f62d540 |
1 | p = 10 |
2 | q = 15 |
3 | k = 2 |
4 | D = 20 |
5 | |
6 | meanX = rep(0,p) |
7 | covX = 0.1*diag(p) |
f227455a |
8 | |
3f62d540 |
9 | covY = array(dim = c(q,q,k)) |
10 | covY[,,1] = 0.1*diag(q) |
11 | covY[,,2] = 0.2*diag(q) |
f227455a |
12 | |
3f62d540 |
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))) |
f227455a |
16 | |
3f62d540 |
17 | n = 100 |
f227455a |
18 | |
19 | pi = c(0.4,0.6) |
20 | |
f227455a |
21 | data = generateXY(meanX,covX,covY, pi, beta, n) |
22 | |
23 | X = data$X |
24 | Y = data$Y |
25 | |
3f62d540 |
26 | res_valse = valse(X,Y) |