Pessoal, enviai um trabalho pra um Congresso e voltou pedindo algumas correções, dentre elas pedindo Anova dos modelos que eu utilizei. Como estou descobrindo o R agora, preciso da ajuda de voces. A partir do modelo: Y=a+b+c+ab+ac+bc, eu consigo recuperar uma anova? Tenho os valores de a, b e c e tb tenho os Y originais e os preditos. Bastaria escrever o modelo e pedir a anova?
A principio sim, mas se coloque no lugar de quem leu sua mensgem.. dificil ajudar mais objetivamente Tente e envie um codigo minimo reproduzivel Leia tb o "posting guide"desta lista On Wed, 11 Jan 2012, Julia Maria wrote:
Pessoal, enviai um trabalho pra um Congresso e voltou pedindo algumas correções, dentre elas pedindo Anova dos modelos que eu utilizei. Como estou descobrindo o R agora, preciso da ajuda de voces. A partir do modelo: Y=a+b+c+ab+ac+bc, eu consigo recuperar uma anova? Tenho os valores de a, b e c e tb tenho os Y originais e os preditos. Bastaria escrever o modelo e pedir a anova?
Se você ajustou o modelo no R usando funções como lm(), aov(), é só pedir m0 <- lm(y~ bla bla bla seu modelo aqui) anova(m0) Se você ajustou em outro programa, este deve ter meios de fornecer isso para você. Como o Prof Paulo disse, envie uma mensagem mais completa. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Pois bem Walmes, o que eu quero saber na verdade é o seguinte: Dá pra fazer: m0 <- lm(y~ 1a+2b+3c+4ab+5ac+6bc) anova(m0) ? Eu acho que realmenteo não consigo explicar o problema. É que eu ja tenho os coeficientes. Acho que deve ser até um erro absurdo esta minha pergunta. Por que até onde eu sei, o correto seria pedir algo do tipo: m0 <- lm(y~ a+b+c+ab+ac+bc) anova(m0) Sem os coeficientes, dái eu obteria os coeficientes e depois a análise certo? parece que estou querendo fazer o contrário.
Você obteve os coeficientes como? Usando a lm() do R? Se você pediu summary() para ver os coeficientes, peça anova() para ver a análise de variância. Até onde sei, não existe função que retorne a anova recebendo apenas os coeficientes do modelo, são necessários o vetor de dados e a matriz de delineamento para isso. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Não sei se deu pra entender, mas eu já tenho no caso um modelo ajustado. ________________________________ De: Julia Maria <juliam_est@yahoo.com.br> Para: "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br> Enviadas: Quarta-feira, 11 de Janeiro de 2012 10:04 Assunto: [R-br] Anova de um modelo. Pessoal, enviai um trabalho pra um Congresso e voltou pedindo algumas correções, dentre elas pedindo Anova dos modelos que eu utilizei. Como estou descobrindo o R agora, preciso da ajuda de voces. A partir do modelo: Y=a+b+c+ab+ac+bc, eu consigo recuperar uma anova? Tenho os valores de a, b e c e tb tenho os Y originais e os preditos. Bastaria escrever o modelo e pedir a anova? _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
participantes (3)
-
Julia Maria -
Paulo Justiniano -
Walmes Zeviani