projects
/
valse.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
merge with remote
[valse.git]
/
test
/
generateRunSaveTest_EMGrank.R
diff --git
a/test/generateRunSaveTest_EMGrank.R
b/test/generateRunSaveTest_EMGrank.R
index
9969620
..
f6a3314
100644
(file)
--- a/
test/generateRunSaveTest_EMGrank.R
+++ b/
test/generateRunSaveTest_EMGrank.R
@@
-1,22
+1,21
@@
source("helper.R")
library(valse)
source("helper.R")
library(valse)
-generateRunSaveTest_EMGrank = function(n=200, p=15, m=10, k=3, mini=5, maxi=10, gamma=1.0,
- rank = c(1,2,4))
+generateRunSaveTest_EMGrank = function(n=200, p=15, m=10, k=3, mini=5, maxi=10, gamma=1.0, rank = c(1,2,4))
{
eps = 1e-6
{
eps = 1e-6
-
p
i = rep(1.0/k, k)
-
r
ho = array(dim=c(m,m,k))
+
P
i = rep(1.0/k, k)
+
R
ho = array(dim=c(m,m,k))
for(i in 1:k)
for(i in 1:k)
-
r
ho[,,i] = diag(1,m)
+
R
ho[,,i] = diag(1,m)
xy = generateXYdefault(n, p, m, k)
xy = generateXYdefault(n, p, m, k)
- testFolder = ".
.
/data/"
+ testFolder = "./data/"
dir.create(testFolder, showWarnings=FALSE, mode="0755")
#save inputs
dir.create(testFolder, showWarnings=FALSE, mode="0755")
#save inputs
-
write.table(as.double(rho), paste(testFolder,"rho
",sep=""),
+
write.table(as.double(Pi), paste(testFolder,"Pi
",sep=""),
row.names=F, col.names=F)
row.names=F, col.names=F)
-
write.table(as.double(pi), paste(testFolder,"pi
",sep=""),
+
write.table(as.double(Rho), paste(testFolder,"Rho
",sep=""),
row.names=F, col.names=F)
write.table(as.integer(mini), paste(testFolder,"mini",sep=""),
row.names=F, col.names=F)
row.names=F, col.names=F)
write.table(as.integer(mini), paste(testFolder,"mini",sep=""),
row.names=F, col.names=F)
@@
-33,7
+32,7
@@
generateRunSaveTest_EMGrank = function(n=200, p=15, m=10, k=3, mini=5, maxi=10,
write.table(as.integer(c(n,p,m,k)), paste(testFolder,"dimensions",sep=""),
row.names=F, col.names=F)
write.table(as.integer(c(n,p,m,k)), paste(testFolder,"dimensions",sep=""),
row.names=F, col.names=F)
- res = valse::EMGrank(
pi,r
ho,mini,maxi,xy$X,xy$Y,eps,rank,fast=FALSE)
+ res = valse::EMGrank(
Pi,R
ho,mini,maxi,xy$X,xy$Y,eps,rank,fast=FALSE)
#save output
write.table(as.double(res$phi),paste(testFolder,"phi",sep=""),row.names=F,col.names=F)
#save output
write.table(as.double(res$phi),paste(testFolder,"phi",sep=""),row.names=F,col.names=F)