+++ /dev/null
-testEMGLLF = function()
-{
- testFolder = "data/"
- array_delimiter = " "
-
- #get dimensions
- dimensions = read.table(paste(testFolder,"dimensions",sep=""), header=FALSE,
- sep=array_delimiter)
- n = dimensions[1]
- p = dimensions[2]
- m = dimensions[3]
- k = dimensions[4]
-
- #get all input arrays
- phiInit = read.table(paste(testFolder,"phiInit",sep=""), header=FALSE, sep=array_delimiter)
- rhoInit = read.table(paste(testFolder,"rhoInit",sep=""), header=FALSE, sep=array_delimiter)
- piInit = read.table(paste(testFolder,"piInit",sep=""), header=FALSE, sep=array_delimiter)
- gamInit = read.table(paste(testFolder,"gamInit",sep=""), header=FALSE, sep=array_delimiter)
- mini = read.table(paste(testFolder,"mini",sep=""), header=FALSE, sep=array_delimiter)
- maxi = read.table(paste(testFolder,"maxi",sep=""), header=FALSE, sep=array_delimiter)
- gamma = read.table(paste(testFolder,"gamma",sep=""), header=FALSE, sep=array_delimiter)
- lambda = read.table(paste(testFolder,"lambda",sep=""), header=FALSE, sep=array_delimiter)
- X = read.table(paste(testFolder,"X",sep=""), header=FALSE, sep=array_delimiter)
- Y = read.table(paste(testFolder,"Y",sep=""), header=FALSE, sep=array_delimiter)
- tau = read.table(paste(testFolder,"tau",sep=""), header=FALSE, sep=array_delimiter)
-
- #run EMGLLF.c
- EMG = .Call("EMGLLF_core",phiInit,rhoInit,piInit1,gamInit,mini,maxi,gamma,lambda,X,Y,tau,
- PACKAGE="valse")
- phi = EMG$phi
- rho = EMG$rho
- pi = EMG$pi
- LLF = EMG$LLF
- S = EMG$S
-
- #get all stored outputs
- ref_phi = read.table(paste(testFolder,"phi",sep=""), header=FALSE, sep=array_delimiter)
- ref_rho = read.table(paste(testFolder,"rho",sep=""), header=FALSE, sep=array_delimiter)
- ref_pi = read.table(paste(testFolder,"pi",sep=""), header=FALSE, sep=array_delimiter)
- ref_LLF = read.table(paste(testFolder,"LLF",sep=""), header=FALSE, sep=array_delimiter)
- ref_S = read.table(paste(testFolder,"S",sep=""), header=FALSE, sep=array_delimiter)
-
- #check that output correspond to stored output
- tol = 1e-5;
- checkOutput("phi",phi,ref_phi,tol);
- checkOutput("rho",rho,ref_rho,tol);
- checkOutput("pi",pi,ref_pi,tol);
- checkOutput("LLF",LLF,ref_LLF,tol);
- checkOutput("S",S,ref_S,tol);
-}