[R] problems with nmle

[kjetil halvorsen]

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(modelExpression[[2]], envir = nlEnv)
Browse[1]> n
Error: subscript out of bounds


Here the fm1Indom.nlme object is calculated with a formula somewhat
different from the one given
in the text, because that gives convergence problems: 


 fm1Indom.nlme <- nlme(fm1Indom.lis, random=pdDiag(A1+lrc1+A2+lrc2 ~ 1)
,
                 control=list(msMaxIter=1000,
msVerbose=TRUE,niterEM=100, 
                       returnObject=TRUE))

The algorithm says it has not converged, although the tracing output
multiple
times indicate that the current solution is probably a solution. The
summary
of the obtained object is very close to that given in the book.

Kjetil Halvorsen
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._