X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=test%2Ftest_EMGrank.R;fp=test%2Ftest_EMGrank.R;h=b140e5ee36cfb1fd39cb460d005126580c11ca7d;hp=9aec8bb9000186159be226d9c53ac4e1b94b3645;hb=f32535f2bc8d50470aa87204bbd7971805dbc9ef;hpb=7fd371e5317f9c61fe5a32daadbbac1c64b2dd31

diff --git a/test/test_EMGrank.R b/test/test_EMGrank.R
index 9aec8bb..b140e5e 100644
--- a/test/test_EMGrank.R
+++ b/test/test_EMGrank.R
@@ -4,13 +4,24 @@ testFolder = "data/"
 # NOTE: R typing is really terrible. as.double as.matrix ...and so on; don't remove.
 
 #get inputs
-#TODO
+npmk = as.matrix(read.table(paste(testFolder,"dimensions",sep="")))
+n = npmk[1]; p=npmk[2]; m=npmk[3]; k=npmk[4]
+Pi = as.double(as.matrix(read.table(paste(testFolder,"Pi",sep="")))[,])
+Rho = array(as.double(as.matrix(read.table(paste(testFolder,"Rho",sep="")))), dim=c(m,m,k))
+mini = as.integer(as.matrix(read.table(paste(testFolder,"mini",sep="")))[1])
+maxi = as.integer(as.matrix(read.table(paste(testFolder,"maxi",sep="")))[1])
+X = matrix(as.double(as.matrix(read.table(paste(testFolder,"X",sep="")))), n,p)
+Y = matrix(as.double(as.matrix(read.table(paste(testFolder,"Y",sep="")))), n,m)
+eps = as.double(as.matrix(read.table(paste(testFolder,"eps",sep="")))[1])
+rank = as.double(as.matrix(read.table(paste(testFolder,"rank",sep="")))[,])
 
 #get outputs
-#TODO
+phi = array(as.double(as.matrix(read.table(paste(testFolder,"phi",sep="")))), dim=c(p,m,k))
+LLF = as.double(as.matrix(read.table(paste(testFolder,"LLF",sep="")))[1])
 
-res = valse::EMGrank(...)
+res = valse::EMGrank(Pi,Rho,mini,maxi,X,Y,eps,rank,fast=TRUE)
 
 #compare outputs
 nd=7 #number of digits
-#TODO
+print( all(round(phi,nd) == round(res$phi,nd)) )
+print( all(round(LLF,nd) == round(res$LLF,nd)) )