Assaf oron
2009-Jan-15 05:59 UTC
[R] Question about correct syntax to use with nlme / nlmer
Hi, I turn to the help list after spending hours on the package help, looking online, asking colleagues and trying various versions on R itself. The work I need help on involves a confidentiality contract, so I am limited in details on the specific model I used. However, my problem is really simple in essence, and is described below. The theoretical model for response is roughly : y = (intercept+linear effects of certain covariates) + alpha * exp (beta*t). (with beta generally believed to be positive, i.e., y begins at some horizontal asymptote, and increases unboundedly as time passes) There is one grouping variable (let's call it "ID"). Most of the fixed effects need to have a random component, i.e. for different ID's they would vary around the fixed mean. Here are my questions: 1. Do nlme/nlmer accept on the RHS of the main formula, only SelfStart-type objects? (I assume yes, because I tried to code the RHS directly in the formula with no success). 2. If so, is there any SelfStart object I can take "off the shelf" and use for this - or must I create my own? (would need guidance to do the latter) 3. In any case, can you instruct me about the proper syntax for the main, fixed and random formulae - or for the seemingly-niftier nlmer interface? Here are some of the errors I get: Error in nlmeCall[[i]] <- NULL : subscript out of bounds Error in nlme.formula(): element 1 is empty; the part of the args list of 'is.numeric' being evaluated was: (start) By the way, the data integrity is fine. I have no problem running lmer on the same dataset. I use R 2.8.1 for Windows and the latest package versions. Thanks in advance, Assaf [[alternative HTML version deleted]]