[R] nlmeODE

Leonid Gibiansky LGibiansky at quantpharm.com
Fri Nov 20 20:37:56 CET 2009

Dear All,

One more nlmeODE question:

I have a differential equation where one of the parameters depends on 
time. Values of this parameter at fixed time points are known but need 
to be interpolated between these known points. Can it be implemented in 
nlmeODE solver?


Leonid Gibiansky, Ph.D.
President, QuantPharm LLC
web:    www.quantpharm.com
e-mail: LGibiansky at quantpharm.com
tel:    (301) 767 5566

David Winsemius wrote:
> On Nov 20, 2009, at 10:52 AM, gcheer3 wrote:
>> Thank you again, Joh
>> Based on your code seems I solved the problem.
>> new code as follows:
>> cs<-expression(delta[1], delta[2])
>> plot(sigma, delta1, ylim=range(-0.5, 2), xlab='sigma', ylab='delta', 
>> pch=22,
>> type='o')
>> points(sigma, delta2, col='red', axes=FALSE, type='o')
>> legend("topleft", cs,    pch = c(21, 22), lty=1, col=2:1)
>> only concern is there is a waring message in R: Warning message:
>> In plot.xy(xy.coords(x, y), type = type, ...) :
>>  "axes" is not a graphical parameter
> You should not be specifying axis types or flags to a function, [ 
> points() in this case] that is not intended to do anything with the 
> axes. To be Biblical: If thine warning offend thee, pluck it[the 
> argument] out.

More information about the R-help mailing list