-
- betaInit1[,,r,repet] = ginv(crossprod(X[Z_indice,])) %*% crossprod(X[Z_indice,], Y[Z_indice,])
+ if (length(Z_indice) == 1) {
+ betaInit1[,,r,repet] = ginv(crossprod(t(X[Z_indice,]))) %*%
+ crossprod(t(X[Z_indice,]), Y[Z_indice,])
+ } else {
+ betaInit1[,,r,repet] = ginv(crossprod(X[Z_indice,])) %*%
+ crossprod(X[Z_indice,], Y[Z_indice,])
+ }