[R] nls, nlrq, and box-cox transformation
Philippe Grosjean
phgrosjean at sciviews.org
Fri Nov 21 08:36:32 CET 2003
>> You suggest the solution yourself: transform the equation to have all
>> parameters at the right, thus:
>>
>> y ~ ((b0 + b1 * x) * t + 1) ^ 1/t
>>
>Bit this is still not correct, since the transformation changes
>the scale of the variance, and lesat squares will not be correct.
>There is needed a factor (the jacobian) to compensate for this,
>Kjetil Halvorsen
OK, sorry you are correct: one would need also to calculate residuals as
(y - ymodel)^2*t instead of (y - ymodel)^2 in the case of nls. This effects
also nlrq, although in a somewhat reduced manner.
Best,
Philippe Grosjean
More information about the R-help
mailing list