[R] nlmeODE

Leonid Gibiansky LGibiansky at quantpharm.com
Fri Nov 20 21:31:51 CET 2009


Please, ignore this question (see below): I found the way how to do it

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



Leonid Gibiansky wrote:
> 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?
> 
> Thanks
> Leonid
> 
> --------------------------------------
> 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