Displaying 1 result from an estimated 1 matches for "mu_mod".
Did you mean:
md_mod
2006 Aug 26
1
problems with loop
...ngth(sim[,1]))
{
if (tmp[k,1]==sim[n,1] && v[k]<=sim[n,2]){tmp[k,2]=1}
if(tmp[k,1]==sim[n,1] && v[k]>sim[n,2]) {tmp[k,2]=0}
}
}
return(tmp)
}
# The model to be fitted
Model_1=function(beta0_mod,mu_mod,k_mod,I_mod)
{
parms = c(beta0=beta0_mod, mu=mu_mod, k=k_mod)
my.atol = 1e-6
times1=c(0,0.002884615,0.003846154,0.005769231,0.009615385,0.01,0.019230769)
times2=c(0.028846154,0.038461538,0.057692308,
0.076923077,0.083333333,0.096153846)
times3=c(0.115384615,0.125,0.1346...