search for: fm1indom

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 In...

Following the Indomethicin example in Pinheiro & Bates, chapter 6, page 277 etc, coming to the following comand: fm2Indom.nlme <- update( fm1Indom.nlme, random = pdDiag(A1 + lrc1 + A2 ~ 1) ) debugging nlme gives the following output: Browse[1]> n debug: modelResid <- ~eval(model, data.frame(data, getParsNlme(plist, fmap, rmapRel, bmap, groups, beta, bvec, b, level, N)))[naPat] Browse[1]> n debug: ww <- eval(modelExpressio...

DeaR I am trying to use a dynamically create formula with nlsList and nlme, but I cannot get the environment of the string-generated formal to work similarly to the manually entered one. Any idea? Dieter #----- library(nlme) # Pinheiro/Bates p 280 fm1Indom.lis = nlsList(conc~SSbiexp(time,A1,lrc1,A2,lrc2), data=Indometh) nlme(fm1Indom.lis,random=pdDiag(A1+lrc1+A2~1)) # works... # Simulating a dynamic formula form = as.formula("conc~SSbiexp(time,A1,lrc1,A2,lrc2)") fm2Indom.lis = nlsList(form,data=Indometh) # still correct nlme(fm2Indom.l...