X-Git-Url: https://git.auder.net/%3C?a=blobdiff_plain;ds=sidebyside;f=reports%2Fmultistart.R;h=9971cd2f7a243f72fcb7a46e3d2a8ad8a7d609d8;hb=b46623addc9e63019aa1df1dd2de800b48fdd609;hp=b3d1e100e4e4e2cba822281b860da9ac8bdd5304;hpb=476a79d4b332d378716e320fdfd95fa3e1d7be9f;p=morpheus.git diff --git a/reports/multistart.R b/reports/multistart.R index b3d1e10..9971cd2 100644 --- a/reports/multistart.R +++ b/reports/multistart.R @@ -46,7 +46,7 @@ ms <- multiRun( library(morpheus) K <- fargs$optargs$K op <- optimParams(K, fargs$optargs$link, fargs$optargs) - x_init <- c(rep(1/K,K-1), as.double(fargs$mu), rep(0,K)) + x_init <- list(p=rep(1/K,K-1), beta=fargs$mu, b=rep(0,K)) do.call(rbind,op$run(x_init)) }, function(fargs) { @@ -58,7 +58,7 @@ ms <- multiRun( best_par <- list() for (i in 1:fargs$optargs$nstart) { - x_init <- c(rep(1/K,K-1), as.double(i*fargs$mu), rep(0,K)) + x_init <- list(p=rep(1/K,K-1), beta=i*fargs$mu, b=rep(0,K)) par <- op$run(x_init) val <- op$f( op$linArgs(par) ) if (val < best_val)