13 nstart <- 3 #nstart-1 random starting points for each MC run
15 cmd_args <- commandArgs()
18 if (substr(arg,1,1)!='-')
20 spl <- strsplit(arg,'=')[[1]]
22 ncores <- as.integer(spl[2])
23 } else if (spl[1] == "N") {
24 N <- as.integer(spl[2])
25 } else if (spl[1] == "n") {
26 n <- as.integer(spl[2])
27 } else if (spl[1] == "d") {
28 d <- as.integer(spl[2])
29 } else if (spl[1] == "link") {
31 } else if (spl[1] == "nstart") {
37 matrix( c(1,-2, 3,1), ncol=K ), #d=2
38 matrix( c(1,2,-1,0,3, 2,-3,0,1,0), ncol=K ), #d=5
39 matrix( c(1,2,-1,0,3,4,-1,-3,0,2, 2,-3,0,1,0,-1,-4,3,2,0), ncol=K ) ) #d=10
40 beta <- betas[[ ifelse( d==2, 1, ifelse(d==5,2,3) ) ]]
43 list(n=n,p=p,beta=beta,b=b,optargs=list(K=K,link=link,nstart=nstart)), list(
48 op <- optimParams(K, fargs$optargs$link, fargs$optargs)
49 x_init <- list(p=rep(1/K,K-1), beta=fargs$mu, b=rep(0,K))
50 do.call(rbind,op$run(x_init))
56 op <- optimParams(K, fargs$optargs$link, fargs$optargs)
59 for (i in 1:fargs$optargs$nstart)
61 x_init <- list(p=rep(1/K,K-1), beta=i*fargs$mu, b=rep(0,K))
63 val <- op$f( op$linArgs(par) )
70 do.call(rbind,best_par)
72 prepareArgs = function(fargs) {
74 io = generateSampleIO(fargs$n, fargs$p, fargs$beta, fargs$b, fargs$optargs$link)
75 fargs$optargs$M <- computeMoments(io$X, io$Y)
76 mu <- computeMu(io$X, io$Y, fargs$optargs)
78 }, N=N, ncores=ncores, verbose=TRUE)
80 ms[[i]] <- alignMatrices(ms[[i]], ref=rbind(p,beta,b), ls_mode="exact")
82 ms_params <- list("N"=N, "nc"=ncores, "n"=n, "K"=K, "link"=link,
83 "p"=p, "beta"=beta, "b"=b, "nstart"=nstart)
85 save(ms, ms_params, file="multistart.RData")