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

De: Thiago De paula protásio <depaulaprotasio@yahoo.com.br>
Assunto: [R-br] R² ajustado de Modelos Não Lineares!
Para: "r-br" <r-br@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@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br