From 0b216f854a21821f9be375d07c2932b31e227e78 Mon Sep 17 00:00:00 2001
From: emilie <emilie@devijver.org>
Date: Tue, 6 Dec 2016 13:29:09 +0100
Subject: [PATCH] =?utf8?q?donn=C3=A9es=20simul=C3=A9es?=
MIME-Version: 1.0
Content-Type: text/plain; charset=utf8
Content-Transfer-Encoding: 8bit

---
 R/generateIO.R        |   6 ++++--
 R/generateIOdefault.R |   3 +--
 data/TODO             |   2 ++
 data/data.RData       | Bin 0 -> 1010 bytes
 4 files changed, 7 insertions(+), 4 deletions(-)
 create mode 100644 data/data.RData

diff --git a/R/generateIO.R b/R/generateIO.R
index 0e776d0..4527f08 100644
--- 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)
-  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)
   {
+    
     for (r in 1:k)
     {
       BXir = rep(0,m)
       for (mm in 1:m)
-        Bxir[[mm]] = X[i,] %*% beta[,mm,r]
+        BXir[mm] = X[i,] %*% beta[,mm,r]
       Y[i,] = Y[i,] + pi[r] * mvrnorm(1,BXir, covY[,,r])
     }
   }
diff --git a/R/generateIOdefault.R b/R/generateIOdefault.R
index 85213cc..3613f2b 100644
--- a/R/generateIOdefault.R
+++ b/R/generateIOdefault.R
@@ -8,11 +8,10 @@
 #-----------------------------------------------------------------------
 generateIOdefault = function(n, p, m, k)
 {
-  covX = array(0, dim=c(p,p,k))
+  covX = diag(p)
   covY = array(0, dim=c(m,m,k))
   for(r in 1:k)
   {
-    covX[,,r] = diag(p)
     covY[,,r] = diag(m)
   }
   
diff --git a/data/TODO b/data/TODO
index c0603b4..a3bb58d 100644
--- a/data/TODO
+++ b/data/TODO
@@ -1,2 +1,4 @@
 Trouver un jeu de données (+) intéressant (que les autres) ?
 Ajouter toy datasets pour les tests (ou piocher dans MASS ?)
+
+ED : j'ai simulé un truc basique via 'generateIOdefault(10,5,6,2)'
diff --git a/data/data.RData b/data/data.RData
new file mode 100644
index 0000000000000000000000000000000000000000..a9f09e143adfa04f595a40b993efa688ebd05a4b
GIT binary patch
literal 1010
zcmb2|=3oE==GaM@0asWhj{V=w{n|mQeJZ0?yPJwYdQQuSf71okLr$bBM0OY(I0_^N
zB&Bp3`1J}H1sw5cv5*q#c4oZ4)NO)xLc}`Vgvf{|MK_sSQ!S%5``1+5dp@($-u83#
zp1vQqvVLU=u{f$`m}ogJUAagqv{a)-YT?IAwW0g3haMC=(tYo)pZ}B}Im)Wfqw37M
zWHq>#seesqJ-Xrh&8h>cU7Ms{Pc;s>{_>Gnk{ma8v(K(=EA2c!ovIAGWBg-p<(4N+
z&Qf|3)n>DAWjo2gp5^VBVNq72`zi0NR`Z-o!5<GYS_c@sSohoUh~xxAxuEnd9ywpv
zYhOr9D_qHXb{D(F@<dCu$W=XF?eA_Tt`=J=^ZZZF4ZnFe4YrpTW$b%XJ7N2sujjpa
z6*I$=euV!$z44m$htFFq4Et(czyGZ9U;GKXl!$lh-Wadu-8mt7&O(0!eU%iYTb9;y
z3hlYQK>byB8^_~ri)SIJYwqq~KmPsi#A&B5ADFr8v6P)q`5X_~_rCt?jv79f`*w)G
z^~d>#Epo3dY$v?g-tGNLhv{eJL0|XZrXd%S-{miP-@f(Tru8%Ga(_JP{naxec-6fg
z??v0TzxCL4Ro6B?QCp{S*4?ukGMZ;KcfL`3TNq)fRNftF-`6!IzfUsV{`l0Z{3nI0
zn;Znb?9q={z^uJdZ&OP8JF|c-8#*$kzWvp&>bNIm)w)a1%I3ZKUnzHNT1QY~;(_A^
z3wCUr-8p;1#yw?w_i9ZK=}t*u7o8k4OZES!Yi)LCtpC*hzwkg;Z?1Vj`L{S5k6XWk
zC4H|KuF42|kr=|idZFUyShetkbxWd7NpT!f%@*!TT0L!vgyk%ebrp>Nmlrjcm9bA3
z+L5_!&iz*9jn|gzKb1<H|15xY<Kq4$$DX}iH&ynK>1k#svud#=*2Qvlv(7nfTqt>2
zZ7X}*q`4Jq{)GgcHM8hAxUpqvL&kEojbf*=`aa6OWn3P)x=ZAGRLp+C+=r#_VvE9`
zS--Mk^IWl}$87Z<L&+mY&72wYs~R!{zJ=}l*Lmr4)Vgvb_CTdTmU$XKMVCeNADSP!
zxo~HyOG^C1rp3n3?!WZE<S6!erp&QJ3E#|j&kj0QeL}lZNG_&Qh~@cYx2dMvCq7RQ
z{kgl8FHes5xk1^igD!TrXU4p8TBvBTO|9+VBkt|T-A=K+v9UVhsk57dTaDd1O@E5R
z^zyWvCF>c~vdlKEiqmeXyMDXKe0}YSd?~-5mc49Bwf&ZN*Y<C)*LjuMy-8YKTKeiW
zW1a^q7oSUgIX7UA+`Qe}BR-1x6xwCh#f0}L*?+yq!tFdmu-tY9p9TA6iFJZf|DS4q
ZxUx2P@q+zQtQ&von>I_e{9$Hb008yy=cxby

literal 0
HcmV?d00001

-- 
2.53.0