[R-br] R² ajustado de Modelos Não Lineares!

[Emmanuel Arnhold]

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 <depaulaprotasio em yahoo.com.br> escreveu:


De: Thiago De paula protásio <depaulaprotasio em yahoo.com.br>
Assunto: [R-br] R² ajustado de Modelos Não Lineares!
Para: "r-br" <r-br em listas.c3sl.ufpr.br>
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
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110426/454f19a7/attachment-0001.html>