Displaying 2 results from an estimated 2 matches for "fm2indom".
Did you mean:
fm1indom
2007 Aug 23
0
Lost in substitute: nlsList and nlme with dynamic formula
...e.
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.lis,random=pdDiag(A1+lrc1+A2~1))
#Error in switch(mode(x), "NULL" = structure(NULL, class = "formula"), :
# invalid formula
# B.D. Ripley val(subst) workaround
# http://finzi.psych.upenn.edu/R/Rhelp02a/ar...
2001 Dec 03
0
problems with nmle
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:...