Displaying 1 result from an estimated 1 matches for "beta_21".
Did you mean:
beta_1
2007 Jan 10
2
problems with optim, "for"-loops and machine precision
...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","conv21")
for (s in 1:500)
{
assign("data",read.table(paste("populations/TE23simset_",s,"....