Displaying 1 result from an estimated 1 matches for "wvar0".
Did you mean:
var0
2013 May 02
0
modMCMC runs in FME package
...tout <- seq(0, 50, by = 1)
out <- model(pars, tout)
return(modCost(obs = observed, model = out))
}
prior <- function(p)
return( sum(((p-c(0.06,18))/c(0.04,6))^2 ))
final <- modMCMC(p = c(k=0.06,v=18), f = obj, prior=prior,lower=c(0.0001,0.1),jump = NULL, niter = 1000,updatecov=100, wvar0 = 1,var0=NULL,burninlength = 50)
summary(final)
plot(final)
the values for "observed" were generated with k=0.05, v=20, which are the retunrs I expect from the modMCMC run. Please comment on:
1. is there any way to get modMCMC to run multiple times faster with this data?
2. my accepted...