R help - May 2007 - A simple question regarding plot of nls objects

Hi,

I was trying to run the example of Indomethacin kinetics from the book:

## From Pinheiro/Bates, Mixed-Effects-Models in S and S-Plus,
## Springer, Second Printing 2001, Section 6.2
library(nlme)
plot(Indometh)
fm1Indom.nls <- nls(conc~SSbiexp(time,A1,lrc1,A2,lrc2), data=Indometh)
summary(fm1Indom.nls)
plot(fm1Indom.nls,Subject~resid(.),abline=0)
## ....

the last plot command gives me the error message: Subject not found in data

What point am I missing?

Thanks in advance,
Hans
-- 

Dr. Hans Mielke
Federal Institute for Risk Assessment
Thielallee 88-92
14195 Berlin
Germany
Phone: +49 30 8412-3969

Hi
r-help-bounces at stat.math.ethz.ch napsal dne 11.05.2007 11:12:21:
> Hi,
> 
> I was trying to run the example of Indomethacin kinetics from the book:
> 
> ## From Pinheiro/Bates, Mixed-Effects-Models in S and S-Plus,
> ## Springer, Second Printing 2001, Section 6.2
> library(nlme)
> plot(Indometh)
> fm1Indom.nls <- nls(conc~SSbiexp(time,A1,lrc1,A2,lrc2), data=Indometh)
> summary(fm1Indom.nls)
> plot(fm1Indom.nls,Subject~resid(.),abline=0)
> ## ....
> 
> the last plot command gives me the error message: Subject not found in 
data
> 
> What point am I missing?
Works for me in R 2.6.0devel. So check your versions of R and nlme. And 
look also to your objects and search path 
ls(), search()

if there is not some naming conflict.

Regards
Petr

> 
> Thanks in advance,
> Hans
> -- 
> 
> Dr. Hans Mielke
> Federal Institute for Risk Assessment
> Thielallee 88-92
> 14195 Berlin
> Germany
> Phone: +49 30 8412-3969
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide 
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

Mielke Hans wrote:
> Hi,
>
> I was trying to run the example of Indomethacin kinetics from the book:
>
> ## From Pinheiro/Bates, Mixed-Effects-Models in S and S-Plus,
> ## Springer, Second Printing 2001, Section 6.2
> library(nlme)
> plot(Indometh)
> fm1Indom.nls <- nls(conc~SSbiexp(time,A1,lrc1,A2,lrc2), data=Indometh)
> summary(fm1Indom.nls)
> plot(fm1Indom.nls,Subject~resid(.),abline=0)
> ## ....
>
> the last plot command gives me the error message: Subject not found in data
>
> What point am I missing?
>
> Thanks in advance,
> Hans
>   
Point-5 maybe? This works nicely for me in R 2.5.0. (On openSuSE 10.2 if
that matters)

-- 
   O__  ---- Peter Dalgaard             ?ster Farimagsgade 5, Entr.B
  c/ /'_ --- Dept. of Biostatistics     PO Box 2099, 1014 Cph. K
 (*) \(*) -- University of Copenhagen   Denmark          Ph:  (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)                  FAX: (+45) 35327907

The example works fine for me on R 2.5.0 and WinXP...

Mielke Hans wrote:
> Hi,
>
> I was trying to run the example of Indomethacin kinetics from the book:
>
> ## From Pinheiro/Bates, Mixed-Effects-Models in S and S-Plus,
> ## Springer, Second Printing 2001, Section 6.2
> library(nlme)
> plot(Indometh)
> fm1Indom.nls <- nls(conc~SSbiexp(time,A1,lrc1,A2,lrc2), data=Indometh)
> summary(fm1Indom.nls)
> plot(fm1Indom.nls,Subject~resid(.),abline=0)
> ## ....
>
> the last plot command gives me the error message: Subject not found in data
>
> What point am I missing?
>
> Thanks in advance,
> Hans
>