| 1 | m=11 |
| 2 | p=10 |
| 3 | |
| 4 | covY = array(0,dim = c(m,m,2)) |
| 5 | covY[,,1] = diag(m) |
| 6 | covY[,,2] = diag(m) |
| 7 | |
| 8 | Beta = array(0, dim = c(p, p, 2)) |
| 9 | Beta[1:4,1:4,1] = 3*diag(4) |
| 10 | Beta[1:4,1:4,2] = -2*diag(4) |
| 11 | |
| 12 | Data = generateXY(100, c(0.5,0.5), rep(0,p), Beta, diag(m), covY) |