[R] An error in fitting a non linear regression
Petr PIKAL
petr.pikal at precheza.cz
Tue Feb 24 11:48:06 CET 2009
Hi
r-help-bounces at r-project.org napsal dne 24.02.2009 11:31:22:
>
> Hi,
>
> Thank you for the reply and suggestions.
>
> I have two questions?
> 1) If I want to use log, it seems that I have to take log from both
sides of
> the model which will lead to lm(log(q)~log(-depth)). What is
tehdifference
> between this syntax and lm(log(q) ~ I(-depth))?
If you have
y = a*exp(-b*x) then log of this equation is
log(y) = log(a) - b * x
at least I was told that by my teacher back at school some decades ago
that log(exp(x)) = x.
You can prove it by
log(exp(whatever))
interactively in R
Regards
Petr
>
> 2) How can I calculate the R-squared of a fitted non linear model?
>
> Regards
> Saeed
>
>
> Christian Ritz-3 wrote:
> >
> > Hi Saeed,
> >
> > one approach is to try out several initial value combinations for a
and b.
> >
> > It often helps to find initial values of the same order of magnitude
and
> > of the same sign
> > as the final estimates.
> >
> > To get such initial values, you could linearize the model:
> >
> > lm(log(q) ~ I(-depth))
> >
> >
> > and supply the estimated coefficients from the linear regression as
> > starting values:
> >
> > nreg <- nls(q ~ a*exp(-b*depth), start = list(a = 0.76168, b =
-0.08484))
> > summary(nreg)
> >
> >
> > Christian
> >
> > ______________________________________________
> > R-help at r-project.org mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide
> > http://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
> >
> >
>
> --
> View this message in context:
http://www.nabble.com/An-error-in-fitting-a-non-
> linear-regression-tp22118160p22179525.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help
mailing list