Dá uma olhada no exemplo abaixo.
Dias
Tambaqui
0
6.09
0
6.09
0
6.09
30
65.73
30
52.00
30
53.13
60
105.88
60
115.80
60
119.31
90
161.06
90
174.94
90
213.63
240
233.67
240
291.93
240
301.33
270
258.47
270
326.84
270
337.04
x=read.table('clipboard', h=T)
attach(x)
mt = nls(Tambaqui~a*(1+b*(exp(-c*Dias)))^-1 , start=c(a=300,b=4, c=0.01))
a=summary(mt)
attach(a)
sqet= sigma^2*df[2]
sqet # Soma de quadrados do erro
sqtt = var(Tambaqui)*(length(Tambaqui)-1)
R2t = 1- sqet/sqtt
R2t #coeficiente de determinação
R2tajustado=1-( ((length(Tambaqui)-1)/ (length(Tambaqui)-3))*(1-R2t))
R2tajustado #coeficiente de determinação ajustado
--- Em ter, 26/4/11, Thiago De paula protásio <[email protected]>
escreveu:
De: Thiago De paula protásio <[email protected]>
Assunto: [R-br] R² ajustado de Modelos Não Lineares!
Para: "r-br" <[email protected]>
Data: Terça-feira, 26 de Abril de 2011, 17:32
Pessoal;
Gostaria de perguntar-lhe uma coisa:
Como faço para extrair o R² ajustado dos modelos não lineares?
Agradeço-lhes qualquer ajuda!
Att.,
Thiago de Paula Protásio
Acadêmico de Engenharia Florestal
Universidade Federal de Lavras
Laboratório de Energia da Biomassa Florestal
Laboratório de Biomateriais
(035) 9183-2246
-----Anexo incorporado-----
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
