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

[Mauro Sznelwar]

Belíssima constribuição! Estava olhando vosso exemplo e rodei. Só uma pequena dúvida: como uso este comando do clipboard, pois para rodar eu copiei e colei num outro bloco de textos e salvei de novo.
        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 

          

       




-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110426/a4e79773/attachment-0001.html>