projects
/
valse.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
données simulées
[valse.git]
/
R
/
generateIO.R
diff --git
a/R/generateIO.R
b/R/generateIO.R
index
0e776d0
..
4527f08
100644
(file)
--- a/
R/generateIO.R
+++ b/
R/generateIO.R
@@
-16,16
+16,18
@@
generateIO = function(covX, covY, pi, beta, n)
k = dim(covY)[3]
Y = matrix(0,n,m)
k = dim(covY)[3]
Y = matrix(0,n,m)
- BX = array(0, dim=c(n,m,k))
+ require(mvtnorm)
+ X = rmvnorm(n, mean = rep(0,p), sigma = covX)
require(MASS) #simulate from a multivariate normal distribution
for (i in 1:n)
{
require(MASS) #simulate from a multivariate normal distribution
for (i in 1:n)
{
+
for (r in 1:k)
{
BXir = rep(0,m)
for (mm in 1:m)
for (r in 1:k)
{
BXir = rep(0,m)
for (mm in 1:m)
- B
xir[[mm]
] = X[i,] %*% beta[,mm,r]
+ B
Xir[mm
] = X[i,] %*% beta[,mm,r]
Y[i,] = Y[i,] + pi[r] * mvrnorm(1,BXir, covY[,,r])
}
}
Y[i,] = Y[i,] + pi[r] * mvrnorm(1,BXir, covY[,,r])
}
}