[R] Non-Linear Regression Problem
sundar.dorai-raj at PDF.COM
Fri Apr 16 15:40:59 CEST 2004
> Dear all,
> First of all i like to thank specially Spencer Graves and Douglas Bates
> for their kind help and suggestions. I estimated the parameter nicely.
> But, for my previous posted problem, the rational of the following line
> R> fit0 <- lm(log(Y)~t-1, draper.data)
> is not clear to me. How come this gives the initial value?
> Any response / help / comment / suggestion / idea / web-link / replies
> will be greatly appreciated.
> Thanks in advance for your time.
I don't recall your previous post, but my guess is that you're trying to
Y = exp(A * t)
Taking the log of both sides gives
log(Y) = A * t
This a linear model with no intercept. To fit this in R, use -1 in the
formula to drop the intercept.
Ahat <- coef(lm(log(Y) ~ t - 1))
Note, I would not use "t" as a variable name, since it is also a
function in the base package (transpose of a matrix). Most of the time R
can make the distinction, but why chance it?
More information about the R-help