1 #' Generate a sample of (X,Y) of size n
2 #' @param covX covariance for covariates
3 #' @param covY covariance for the response vector
4 #' @param pi proportion for each cluster
5 #' @param beta regression matrix
6 #' @param n sample size
7 #' @return list with X and Y
9 #-----------------------------------------------------------------------
10 generateIO = function(covX, covY, pi, beta, n)
20 BX = array(0, dim=c(n,m,k))
22 require(MASS) #simulate from a multivariate normal distribution
29 Bxir[[mm]] = X[i,] %*% beta[,mm,r]
30 Y[i,] = Y[i,] + pi[r] * mvrnorm(1,BXir, covY[,,r])
34 return (list(X=X,Y=Y))