[R] nls, nlrq, and box-cox transformation
Philippe Grosjean
phgrosjean at sciviews.org
Thu Nov 20 15:24:19 CET 2003
>Dear r-help members
>I posted this message already yesterday, but don't know whether it
>reached you since I joined the group only yesterday.
>I would like to estimate the boxcox transformed model
> (y^t - 1)/t ~ b0 + b1 * x.
>Unfortunately, R returns with an error message when I try to
>perform this with the call
>nls( I((y^t - 1)/t) ~ I(b0 + b1*x),
> start = c(t=1,b0=0,b1=0), data = mydataframe)
>The error message is: Object "t" not found
>Apparently R seems not to accept parameters on the left hand
>side of a regression model. I know that my do-it-yourself
>strategy is not necessary, since the package box-cox is
>available. Unfortunately, I want the use the box-cox
>transformation in a quantile regression, i.e. I have to replace
>nls by nlrq in the call above.
>Any suggestions?
>Thanks and best regards,
> Johannes Ludsteck
You suggest the solution yourself: transform the equation to have all
parameters at the right, thus:
y ~ ((b0 + b1 * x) * t + 1) ^ 1/t
(double check if this is correct)
Best,
Philippe Grosjean
...........]<(({?<...............<?}))><...............................
) ) ) ) )
( ( ( ( ( Dr. Philippe Grosjean
) ) ) ) )
( ( ( ( ( Numerical Ecology Laboratory
) ) ) ) ) Mons-Hainaut University
( ( ( ( ( 8, Av. du Champ de Mars, 7000 Mons
) ) ) ) ) Belgium
( ( ( ( (
) ) ) ) ) e-mail: phgrosjean at sciviews.org
( ( ( ( ( SciViews project coordinator (http://www.sciviews.org)
) ) ) ) )
.......................................................................
More information about the R-help
mailing list