1 testEMGLLF = function()
7 dimensions = read.table(paste(testFolder,"dimensions",sep=""), header=FALSE,
15 phiInit = read.table(paste(testFolder,"phiInit",sep=""), header=FALSE, sep=array_delimiter)
16 rhoInit = read.table(paste(testFolder,"rhoInit",sep=""), header=FALSE, sep=array_delimiter)
17 piInit = read.table(paste(testFolder,"piInit",sep=""), header=FALSE, sep=array_delimiter)
18 gamInit = read.table(paste(testFolder,"gamInit",sep=""), header=FALSE, sep=array_delimiter)
19 mini = read.table(paste(testFolder,"mini",sep=""), header=FALSE, sep=array_delimiter)
20 maxi = read.table(paste(testFolder,"maxi",sep=""), header=FALSE, sep=array_delimiter)
21 gamma = read.table(paste(testFolder,"gamma",sep=""), header=FALSE, sep=array_delimiter)
22 lambda = read.table(paste(testFolder,"lambda",sep=""), header=FALSE, sep=array_delimiter)
23 X = read.table(paste(testFolder,"X",sep=""), header=FALSE, sep=array_delimiter)
24 Y = read.table(paste(testFolder,"Y",sep=""), header=FALSE, sep=array_delimiter)
25 tau = read.table(paste(testFolder,"tau",sep=""), header=FALSE, sep=array_delimiter)
28 EMG = .Call("EMGLLF_core",phiInit,rhoInit,piInit1,gamInit,mini,maxi,gamma,lambda,X,Y,tau,
36 #get all stored outputs
37 ref_phi = read.table(paste(testFolder,"phi",sep=""), header=FALSE, sep=array_delimiter)
38 ref_rho = read.table(paste(testFolder,"rho",sep=""), header=FALSE, sep=array_delimiter)
39 ref_pi = read.table(paste(testFolder,"pi",sep=""), header=FALSE, sep=array_delimiter)
40 ref_LLF = read.table(paste(testFolder,"LLF",sep=""), header=FALSE, sep=array_delimiter)
41 ref_S = read.table(paste(testFolder,"S",sep=""), header=FALSE, sep=array_delimiter)
43 #check that output correspond to stored output
45 checkOutput("phi",phi,ref_phi,tol);
46 checkOutput("rho",rho,ref_rho,tol);
47 checkOutput("pi",pi,ref_pi,tol);
48 checkOutput("LLF",LLF,ref_LLF,tol);
49 checkOutput("S",S,ref_S,tol);