- class = sample(1:k, 1, prob=pi)
- X[i,] = mvrnorm(1, meanX[,class], covX[,,class])
- Y[i,] = mvrnorm(1, X[i,] %*% beta[,,class], covY[,,class])
+ class[i] = sample(1:k, 1, prob=pi)
+ X[i,] = mvrnorm(1, meanX, covX)
+ print(X[i,])
+ print(beta[,,class[i]])
+ Y[i,] = mvrnorm(1, X[i,] %*% beta[,,class[i]], covY[,,class[i]])