Inputparametters.R + SelecModel

[valse.git] / InputParameters / generateIO.R

diff --git a/InputParameters/generateIO.R b/InputParameters/generateIO.R
new file mode 100644 (file)
index 0000000..9e84af5
--- /dev/null
+++ b/InputParameters/generateIO.R
@@ -0,0 +1,25 @@
+library(MASS) #simulate from a multivariate normal distribution
+
+generateIO = function(meanX, covX, covY, pi, beta, n){ #don't need meanX
+  size_covX = dim(covX)
+  p = size_covX[1]
+  k = size_covX[3]
+  
+  size_covY = dim(covY)
+  m = size_covY[1]
+  
+  Y = matrix(0,n,m)
+  BX = array(0, dim=c(n,m,k))
+  
+  for(i in 1:n){
+    for(r in 1:k){
+      BXir = rep(0,m)
+      for(mm in 1:m){
+        Bxir[[mm]] = X[i,] %*% beta[,mm,r]
+      }
+      Y[i,]=Y[i,] + pi[[r]] * mvrnorm(1,BXir, covY[,,r])
+    }
+  }
+  
+  return(list(X,Y))
+}
\ No newline at end of file