1 generateIO = function(covX, covY, pi, beta, n)
11 BX = array(0, dim=c(n,m,k))
13 require(MASS) #simulate from a multivariate normal distribution
20 Bxir[[mm]] = X[i,] %*% beta[,mm,r]
21 Y[i,] = Y[i,] + pi[r] * mvrnorm(1,BXir, covY[,,r])
25 return (list(X=X,Y=Y))