[R] How to find best parameter values using deSolve n optim() ?
mhimanshu
bioinfo.himanshu at gmail.com
Thu Apr 5 18:32:32 CEST 2012
Hi Thomas,
Thank you so much for your suggestion.
I tried your code and it is working fine. Now when I change the values of Y
in yobs I am getting so many warnings.
say,
yobs <- data.frame(
time = 0:7,
Y = c(0.00, 3.40, 4.60 ,5.80, 5.80, 6.00, 6.00 ),
Z = c(0.1, 0.11, 0.119, 0.128, 0.136, 0.145, 0.153, 0.16)
)
So when i fit the model with the same code that you have written, i got the
following warnings:
DLSODA- Warning..Internal T (=R1) and H (=R2) are
such that in the machine, T + H = T on the next step
(H = step size). Solver will continue anyway.
In above, R1 = 0.1484502806322D+01 R2 = 0.2264549048113D-16
and I have got so many such warnings.
Can you explain me why this is happening?? and Secondly, I dont understand
why i am getting parameters values in negatives after fitting?? Can you
please help me out in this... :)
Thanks
Himanshu
--
View this message in context: http://r.789695.n4.nabble.com/How-to-find-best-parameter-values-using-deSolve-n-optim-tp4506042p4535368.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list