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

Mauro Sznelwar sznelwar em uol.com.br
Quarta Abril 27 22:57:58 BRT 2011


Obrigado pelo retorno, mas originou uma só coluna em vez de duas, veja um trecho:
Dias
1  Tambaqui
2         0
3      6.09
4         0
5      6.09
6         0
7      6.09
8        30
9     65.73
10       30
11    52.00
12       30
13    53.13
14       60
15   105.88
16       60
17   115.80
18       60

        Olá Mauro. O comando "clipboard" extrai da área de trasnferência. 

        Faça o seguinte.

        digite o comando abaixo e deixe o cursor no final do texto, após )
        x=read.table('clipboard', h=T) 



        Vá para a tabela em um arquivo de texto ou planilha e copie para a área de transferencia



        Volte para a área de trabalho do R e digite enter



        Pronto, o R recupera da área de transfêrência e vc cria o objeto x com a tabela que vc quiser.



        Valeu.



        --- Em qua, 27/4/11, Mauro Sznelwar <sznelwar em uol.com.br> escreveu:


          De: Mauro Sznelwar <sznelwar em uol.com.br>
          Assunto: Re: [R-br] R² ajustado de Modelos Não Lineares!
          Para: r-br em listas.c3sl.ufpr.br
          Data: Quarta-feira, 27 de Abril de 2011, 1:58


          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/20110427/ca144843/attachment-0001.html>


Mais detalhes sobre a lista de discussão R-br