[R] Ajusts in nls
Ronaldo Reis Jr.
chrysopa at insecta.ufv.br
Fri Mar 28 19:00:51 CET 2003
Hi,
I make an ajust using NLS.
The ajust is:
> library(nls)
> modelo <- nls(y~a-b*exp(-c*x),start=list(a=100,b=90,c=0.1609))
> summary(modelo)
Formula: y ~ a - b * exp(-c * x)
Parameters:
Estimate Std. Error t value Pr(>|t|)
a 115.2527 2.9139 39.55 < 2e-16 ***
b 118.6876 7.8925 15.04 < 2e-16 ***
c 0.1235 0.0171 7.22 2.44e-09 ***
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1
Residual standard error: 13.21 on 51 degrees of freedom
Correlation of Parameter Estimates:
a b
b 0.0784
c -0.6587 0.4612
My doubt is about the model significance, I read in list some questions about
the nonsense of this, but I need of the model significance and the R^2.
I make this:
modelonulo <- lm(y~1)
> anova(modelo,modelonulo,test="F")
Analysis of Variance Table
Model 1: y ~ a - b * exp(-c * x)
Model 2: y ~ 1
Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
1 51 8897
2 53 59008 -2 -50111 143.62 < 2.2e-16 ***
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1
And the R^2 is:
> 50111/59008
[1] 0.8492238
It is """correct"""?
Thanks
Ronaldo
--
Sutil como gato que vai pegar passarinho.
--
| // | \\ [*****************************][*******************]
|| ( õ õ ) [Ronaldo Reis Júnior ][PentiumIII-600 ]
| V [UFV/DBA-Entomologia ][HD: 30 + 10 Gb ]
|| / \ [36571-000 Viçosa - MG ][RAM: 128 Mb ]
| /(.''`.)\ [Fone: 31-3899-2532 ][Video: SiS620-8Mb ]
||/(: :' :)\ [chrysopa at insecta.ufv.br ][Modem: Pctel-onboar]
|/ (`. `'` ) \[ICQ#: 5692561 ][Kernel: 2.4.18 ]
|| ( `- ) [*****************************][*******************]
||| _/ \_Powered by GNU/Debian W/Sarge D+ || Lxuser#: 205366
More information about the R-help
mailing list