search for: nlenv

...= 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) , c...

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