From: Benjamin Goehry Date: Mon, 9 Jan 2017 13:43:16 +0000 (+0100) Subject: typo X-Git-Url: https://git.auder.net/%7B%7B%20asset%28%27mixstore/css/user/doc/pieces/mini-custom.min.css?a=commitdiff_plain;h=71a8ee557968aa1130b9c5e47690cf73631474a4;p=valse.git typo --- 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,])