projects
/
valse.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
prepare EMGLLF / EMGrank wrappers, simplify folder generateTestData
[valse.git]
/
src
/
test
/
generate_test_data
/
generateRunSaveTest_EMGLLF.R
diff --git
a/src/test/generate_test_data/generateRunSaveTest_EMGLLF.R
b/src/test/generate_test_data/generateRunSaveTest_EMGLLF.R
index
1c16318
..
fe8f847
100644
(file)
--- a/
src/test/generate_test_data/generateRunSaveTest_EMGLLF.R
+++ b/
src/test/generate_test_data/generateRunSaveTest_EMGLLF.R
@@
-1,12
+1,14
@@
+source("EMGLLF.R")
+
generateRunSaveTest_EMGLLF = function(n=200, p=15, m=10, k=3, mini=5, maxi=10,
gamma=1., lambda=0.5, tau=1e-6)
{
generateRunSaveTest_EMGLLF = function(n=200, p=15, m=10, k=3, mini=5, maxi=10,
gamma=1., lambda=0.5, tau=1e-6)
{
- testFolder = "data/"
+ testFolder = "
../
data/"
dir.create(testFolder, showWarnings=FALSE, mode="0755")
require(valse)
dir.create(testFolder, showWarnings=FALSE, mode="0755")
require(valse)
- params = valse:::basic
_Init_
Parameters(n, p, m, k)
-
io = generateIO
default(n, p, m, k)
+ params = valse:::basic
Init
Parameters(n, p, m, k)
+
xy = valse:::generateXY
default(n, p, m, k)
#save inputs
write.table(as.double(params$phiInit), paste(testFolder,"phiInit",sep=""),
#save inputs
write.table(as.double(params$phiInit), paste(testFolder,"phiInit",sep=""),
@@
-25,9
+27,9
@@
generateRunSaveTest_EMGLLF = function(n=200, p=15, m=10, k=3, mini=5, maxi=10,
row.names=F, col.names=F)
write.table(as.double(lambda), paste(testFolder,"lambda",sep=""),
row.names=F, col.names=F)
row.names=F, col.names=F)
write.table(as.double(lambda), paste(testFolder,"lambda",sep=""),
row.names=F, col.names=F)
- write.table(as.double(
io
$X), paste(testFolder,"X",sep=""),
+ write.table(as.double(
xy
$X), paste(testFolder,"X",sep=""),
row.names=F, col.names=F)
row.names=F, col.names=F)
- write.table(as.double(
io
$Y), paste(testFolder,"Y",sep=""),
+ write.table(as.double(
xy
$Y), paste(testFolder,"Y",sep=""),
row.names=F, col.names=F)
write.table(as.double(tau), paste(testFolder,"tau",sep=""),
row.names=F, col.names=F)
row.names=F, col.names=F)
write.table(as.double(tau), paste(testFolder,"tau",sep=""),
row.names=F, col.names=F)
@@
-35,7
+37,7
@@
generateRunSaveTest_EMGLLF = function(n=200, p=15, m=10, k=3, mini=5, maxi=10,
row.names=F, col.names=F)
res = EMGLLF(params$phiInit,params$rhoInit,params$piInit,params$gamInit,mini,maxi,
row.names=F, col.names=F)
res = EMGLLF(params$phiInit,params$rhoInit,params$piInit,params$gamInit,mini,maxi,
- gamma,lambda,
io$X,io
$Y,tau)
+ gamma,lambda,
xy$X,xy
$Y,tau)
#save outputs
write.table(as.double(res$phi), paste(testFolder,"phi",sep=""), row.names=F, col.names=F)
#save outputs
write.table(as.double(res$phi), paste(testFolder,"phi",sep=""), row.names=F, col.names=F)