src/test/TEST R/testEMGLLF.R

   1 testEMGLLF = function(n,p,m,k, phiInit, rhoInit, gamInit, mini, maxi, gamma, lambda, X, Y, tau, ref_phi, ref_rho, ref_pi, ref_LLf, ref_S){
   2   testFolder = 'data/'
   3   delimiter = '\n'
   4
   5   EMG = .Call("EMGLLF_core",phiInit,rhoInit,piInit1,gamInit,mini,maxi,gamma,lambda,X,Y,tau)
   6   phi = EMG$phi
   7   rho = EMG$rho
   8   pi = EMG$pi
   9   LLF = EMG$LLF
  10   S = EMG$S
  11
  12
  13   tol = 1e-5;
  14   checkOutput('phi',phi,ref_phi,tol);
  15   checkOutput('rho',rho,ref_rho,tol);
  16   checkOutput('pi',pi,ref_pi,tol);
  17   checkOutput('LLF',LLF,ref_LLF,tol);
  18   checkOutput('S',S,ref_S,tol);
  19 }