From 71a8ee557968aa1130b9c5e47690cf73631474a4 Mon Sep 17 00:00:00 2001 From: Benjamin Goehry <benjamin.goehry@math.u-psud.fr> Date: Mon, 9 Jan 2017 14:43:16 +0100 Subject: [PATCH] typo --- R/initSmallEM.R | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/R/initSmallEM.R b/R/initSmallEM.R index c24fca9..b70eea9 100644 --- a/R/initSmallEM.R +++ b/R/initSmallEM.R @@ -37,7 +37,7 @@ initSmallEM = function(k,X,Y,tau) Z_vec = Z_bin$vec #vecteur 0 et 1 aux endroits o? Z==r Z_indice = Z_bin$indice #renvoit les indices o? Z==r - betaInit1[,,r,repet] = ginv(t(X[Z_indice,])*X[Z_indice,])%*%t(X[Z_indice,])%*%Y[Z_indice,] + betaInit1[,,r,repet] = ginv( crossprod(X[Z_indice,]) ) %*% crossprod(X[Z_indice,], Y[Z_indice,]) sigmaInit1[,,r,repet] = diag(m) phiInit1[,,r,repet] = betaInit1[,,r,repet]/sigmaInit1[,,r,repet] rhoInit1[,,r,repet] = solve(sigmaInit1[,,r,repet]) @@ -48,7 +48,7 @@ initSmallEM = function(k,X,Y,tau) { for(r in 1:k) { - dotProduct = 3 * (Y[i,]%*%rhoInit1[,,r,repet]-X[i,]%*%phiInit1[,,r,repet]) %*% (Y[i,]%*%rhoInit1[,,r,repet]-X[i,]%*%phiInit1[,,r,repet]) + dotProduct = (Y[i,]%*%rhoInit1[,,r,repet]-X[i,]%*%phiInit1[,,r,repet]) %*% (Y[i,]%*%rhoInit1[,,r,repet]-X[i,]%*%phiInit1[,,r,repet]) Gam[i,r] = piInit1[repet,r]*det(rhoInit1[,,r,repet])*exp(-0.5*dotProduct) } sumGamI = sum(Gam[i,]) -- 2.44.0