R² ajustado de Modelos Não Lineares!
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
Thiago, Percebo que vc está querendo utilizar as mesmas técnicas para ajuste de modelos lineares em modelos não-lineares. Recomendo o estudo de modelos não-lineares antes de se aventurar neste mundo! Ontem foi o CV e agora o R2. Consulte algum livro sobre modelos não-lineares, a própria dissertação do Walmes que está no site de teses da UFLA vai dar uma idéia sobre o ajuste de modelos nao-lineares! Valeu!!! Fábio Mathias Corrêa Departamento de Estatística Universidade Estadual de Santa Cruz Tel.: 73-3680-5076 ________________________________ De: Thiago De paula protásio <depaulaprotasio@yahoo.com.br> Para: r-br <r-br@listas.c3sl.ufpr.br> Enviadas: Terça-feira, 26 de Abril de 2011 14:32:00 Assunto: [R-br] R² ajustado de Modelos Não Lineares! 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
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
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
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@uol.com.br> escreveu: De: Mauro Sznelwar <sznelwar@uol.com.br> Assunto: Re: [R-br] R² ajustado de Modelos Não Lineares! Para: r-br@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 -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
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@uol.com.br> escreveu: De: Mauro Sznelwar <sznelwar@uol.com.br> Assunto: Re: [R-br] R² ajustado de Modelos Não Lineares! Para: r-br@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
participantes (4)
-
Emmanuel Arnhold -
Fabio Mathias Corrêa -
Mauro Sznelwar -
Thiago De paula protásio