| 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 | } |