Displaying 2 results from an estimated 2 matches for "modelresid".
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: ww <- eval(modelExpression[[2]], envir = nlEnv)
Browse[1]> n
Error: subscript out of bounds
Here the fm1Indom.nlme object is calculated with a fo...
2001 May 30
2
environments
I would like to be able, inside a function, to create a new function, and
use it as part of a formula as an argument to, say, gnls or nlme. for
example:
MyTop <- function(data=dta) {
Cexp <- function(dose,A,B,m){...}
Model <- as.formula(paste("y","~ Cexp(",paste(formals(Cexp),collapse
=", "),")"))
MyCall <-