Commit | Line | Data |
---|---|---|
c3bc4705 BA |
1 | generateRunSaveTest_constructionModelesLassoRank = function(n=200, p=15, m=10, L=12, mini=5, |
2 | maxi=10, gamma=1.0, rangmin=3, rangmax=6) | |
3 | { | |
f1b906da BG |
4 | testFolder = "data/" |
5 | dir.create(testFolder, showWarnings=FALSE, mode="0755") | |
c3bc4705 | 6 | |
f1b906da | 7 | tau = 1e-6 |
f1b906da BG |
8 | pi = matrix(0, k,L) |
9 | for(i in 1:L){ | |
10 | pi[,i] = rep(1.0/k, k) | |
11 | } | |
12 | rho = array(0, dim=c(m,m,k,L)) | |
13 | for(l in 1:L){ | |
14 | for(r in 1:k){ | |
15 | rho[,,r,l] = diag(1,m) | |
16 | } | |
17 | } | |
9ade3f1b | 18 | require(valse) |
c3bc4705 | 19 | io = valse:::generateIOdefault(n, p, m, k) |
f1b906da BG |
20 | A1 = matrix(0,p,L) |
21 | for(i in 1:L){ | |
22 | A1[,i] = seq(1,p) | |
23 | } | |
c3bc4705 | 24 | |
f1b906da | 25 | #save inputs |
c3bc4705 BA |
26 | write.table(as.double(rho),paste(testFolder,"rho",sep=""), |
27 | row.names=F, col.names=F) | |
28 | write.table(as.double(pi), paste(testFolder,"pi",sep=""), | |
29 | row.names=F, col.names=F) | |
30 | write.table(as.integer(mini), paste(testFolder,"mini",sep=""), | |
31 | row.names=F, col.names=F) | |
32 | write.table(as.integer(maxi), paste(testFolder,"maxi",sep=""), | |
33 | row.names=F, col.names=F) | |
34 | write.table(as.double(io$X), paste(testFolder,"X",sep=""), | |
35 | row.names=F, col.names=F) | |
36 | write.table(as.double(io$Y), paste(testFolder,"Y",sep=""), | |
37 | row.names=F, col.names=F) | |
38 | write.table(as.double(tau),paste(testFolder,"tau",sep=""), | |
39 | row.names=F, col.names=F) | |
40 | write.table(as.double(A1),paste(testFolder,"A1",sep=""), | |
41 | row.names=F, col.names=F) | |
42 | write.table(as.integer(rangmin),paste(testFolder,"rangmin",sep=""), | |
43 | row.names=F, col.names=F) | |
44 | write.table(as.integer(rangmax),paste(testFolder,"rangmax",sep=""), | |
45 | row.names=F, col.names=F) | |
46 | write.table(as.integer(c(n,p,m,k)),paste(testFolder,"dimensions",sep=""), | |
47 | row.names=F, col.names=F) | |
48 | ||
49 | res = constructionModelesLassoRank(pi,rho,mini,maxi,X,Y,tau,A1,rangmin,rangmax) | |
50 | ||
f1b906da | 51 | #save output |
c3bc4705 BA |
52 | write.table(as.double(res$phi), paste(testFolder,"phi",sep=""), row.names=F, col.names=F) |
53 | write.table(as.double(res$llh), paste(testFolder,"llh",sep=""), row.names=F, col.names=F) | |
9ade3f1b | 54 | } |