SEXP lambda_,
SEXP X_,
SEXP Y_,
- SEXP tau_
+ SEXP eps_
) {
// Get matrices dimensions
int n = INTEGER(getAttrib(X_, R_DimSymbol))[0];
int maxi = INTEGER_VALUE(maxi_);
double gamma = NUMERIC_VALUE(gamma_);
double lambda = NUMERIC_VALUE(lambda_);
- double tau = NUMERIC_VALUE(tau_);
+ double eps = NUMERIC_VALUE(eps_);
// Get pointers from SEXP arrays ; WARNING: by columns !
double* phiInit = REAL(phiInit_);
PROTECT(phi = allocArray(REALSXP, dimPhiS));
PROTECT(rho = allocArray(REALSXP, dimRho));
PROTECT(pi = allocVector(REALSXP, k));
- PROTECT(LLF = allocVector(REALSXP, maxi-mini+1));
+ PROTECT(LLF = allocVector(REALSXP, maxi));
PROTECT(S = allocArray(REALSXP, dimPhiS));
PROTECT(affec = allocVector(INTSXP, n));
double *pPhi=REAL(phi), *pRho=REAL(rho), *pPi=REAL(pi), *pLLF=REAL(LLF), *pS=REAL(S);
// Call to EMGLLF //
////////////////////
- EMGLLF_core(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,lambda,X,Y,tau,
+ EMGLLF_core(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,lambda,X,Y,eps,
pPhi,pRho,pPi,pLLF,pS,pAffec,
n,p,m,k);