Displaying 1 result from an estimated 1 matches for "mu_23".
Did you mean:
mu_2
2007 Jan 10
2
problems with optim, "for"-loops and machine precision
...l should fit the data either equally
or better than the 6-parameter model. Therefore the difference of the
minimal NLLs should be 0 or larger.
For 500 data sets I use the following code:
require(odesolve)
res=matrix(nrow=500,ncol=18)
colnames(res)=c("a_23","beta_23","mu_23","d_23","I_23","M_23","NLL_23",
"a_21","beta_21","mu_21","e_21","d_21","I_21","M_21","NLL_21",
"NLL23_min_21","conv23","con...