+ #get dimensions
+ dimensions = read.table(strcat(testFolder,'dimensions'), header=FALSE)
+ n = dimensions[1]
+ p = dimensions[2]
+ m = dimensions[3]
+ k = dimensions[4]
+
+ #get all input arrays
+ phiInit = read.table(strcat(testFolder,'phiInit'), header=FALSE)
+ rhoInit = read.table(strcat(testFolder,'rhoInit'), header=FALSE)
+ piInit = t(read.table(strcat(testFolder,'piInit'), header=FALSE))
+ gamInit = read.table(strcat(testFolder,'gamInit'), header=FALSE)
+ mini = read.table(strcat(testFolder,'mini'), header=FALSE)
+ maxi = read.table(strcat(testFolder,'maxi'), header=FALSE)
+ gamma = read.table(strcat(testFolder,'gamma'), header=FALSE)
+ lambda = read.table(strcat(testFolder,'lambda'), header=FALSE)
+ X = rread.table(strcat(testFolder,'X'), header=FALSE)
+ Y = read.table(strcat(testFolder,'Y'), header=FALSE)
+ tau = read.table(strcat(testFolder,'tau'), header=FALSE)
+
+ #run EMGLLF.c