Re: [R] scope, lme, ns, nlme, splines

2012-12-19 Thread Duncan Murdoch
On 12-12-19 5:02 PM, Jacob Wegelin wrote: Your solution works beautifully for predict.lme in the same data that were used to compute the model. But what if I want to compute the population fitted values in newdata? I didn't offer a solution, I offered a workaround for a bug. The solution is

Re: [R] scope, lme, ns, nlme, splines

2012-12-19 Thread Jacob Wegelin
Your solution works beautifully for predict.lme in the same data that were used to compute the model. But what if I want to compute the population fitted values in newdata? In the expanded example below, predict.lm is able to find the object called fixed. But predict.lme is unable to find it

[R] scope, lme, ns, nlme, splines

2012-12-06 Thread Jacob Wegelin
I want to fit a series of lme() regression models that differ only in the degrees of freedom of a ns() spline. I want to use a wrapper function to do this. The models will be of the form y ~ ns(x, df=splineDF) where splineDF is passed as an argument to a wrapper function. This works fine if

Re: [R] scope, lme, ns, nlme, splines

2012-12-06 Thread Duncan Murdoch
On 06/12/2012 2:28 PM, Jacob Wegelin wrote: I want to fit a series of lme() regression models that differ only in the degrees of freedom of a ns() spline. I want to use a wrapper function to do this. The models will be of the form y ~ ns(x, df=splineDF) where splineDF is passed as an argument