Fatorial com parcelas perdidas
Prezados, gostaria de saber se no caso de parcelas perdidas é só substituir o valor perdido par NA e usar a função aov() como no exemplo abaixo. fator1<-c(rep("A",12),rep("B",12)) fator2<-c(rep("C",4),rep("D",4),rep("E",4),rep("C",4),rep("D",4),rep("E",4)) resposta<-c(39.03794, 30.19444, 32.57167, 34.87484, NA ,37.15130, 39.67066, 32.34696, 19.96688, 63.33362, 43.74447, 47.39147, 37.37704, NA, 29.35858, 28.76528, 16.092950, 21.080503, 16.906132, 12.345620, 26.475927, 23.639545, 23.890539, NA) dados<-data.frame(fator1,fator2,resposta) anova<-aov(resposta~fator1*fator2) summary(anova) Existe algum pacote que realiza testes de comparação múltipla para dados desbalanciados. Att. Tiago.
O pacote TukeyC realiza um teste de tukey para dados desbalanceado Em 25 de novembro de 2012 20:07, tiago souza marçal < tiagosouzamarcal@hotmail.com> escreveu:
Prezados,
gostaria de saber se no caso de parcelas perdidas é só substituir o valor perdido par NA e usar a função aov() como no exemplo abaixo.
fator1<-c(rep("A",12),rep("B",12))
fator2<-c(rep("C",4),rep("D",4),rep("E",4), rep("C",4),rep("D",4),rep("E",4))
resposta<-c(39.03794, 30.19444, 32.57167, 34.87484, NA ,37.15130, 39.67066, 32.34696, 19.96688, 63.33362, 43.74447,
47.39147, 37.37704, NA, 29.35858, 28.76528, 16.092950, 21.080503, 16.906132, 12.345620, 26.475927, 23.639545, 23.890539, NA)
dados<-data.frame(fator1,fator2,resposta)
anova<-aov(resposta~fator1*fator2)
summary(anova)
Existe algum pacote que realiza testes de comparação múltipla para dados desbalanciados.
Att.
Tiago.
_______________________________________________ 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.
há estava esquecendo o pacote mutoss também trabalha para dados desbalanceado contendo um fator de variação somente (one way) Em 25 de novembro de 2012 20:33, Fernando Antonio de souza < nandodesouza@gmail.com> escreveu:
O pacote TukeyC realiza um teste de tukey para dados desbalanceado
Em 25 de novembro de 2012 20:07, tiago souza marçal < tiagosouzamarcal@hotmail.com> escreveu:
Prezados,
gostaria de saber se no caso de parcelas perdidas é só substituir o valor perdido par NA e usar a função aov() como no exemplo abaixo.
fator1<-c(rep("A",12),rep("B",12))
fator2<-c(rep("C",4),rep("D",4),rep("E",4), rep("C",4),rep("D",4),rep("E",4))
resposta<-c(39.03794, 30.19444, 32.57167, 34.87484, NA ,37.15130, 39.67066, 32.34696, 19.96688, 63.33362, 43.74447,
47.39147, 37.37704, NA, 29.35858, 28.76528, 16.092950, 21.080503, 16.906132, 12.345620, 26.475927, 23.639545, 23.890539, NA)
dados<-data.frame(fator1,fator2,resposta)
anova<-aov(resposta~fator1*fator2)
summary(anova)
Existe algum pacote que realiza testes de comparação múltipla para dados desbalanciados.
Att.
Tiago.
_______________________________________________ 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.
Obrigado pelas dicas Fernando. Att. Tiago. Date: Sun, 25 Nov 2012 20:38:23 -0200 From: nandodesouza@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Fatorial com parcelas perdidas há estava esquecendo o pacote mutoss também trabalha para dados desbalanceado contendo um fator de variação somente (one way) Em 25 de novembro de 2012 20:33, Fernando Antonio de souza <nandodesouza@gmail.com> escreveu: O pacote TukeyC realiza um teste de tukey para dados desbalanceado Em 25 de novembro de 2012 20:07, tiago souza marçal <tiagosouzamarcal@hotmail.com> escreveu: Prezados, gostaria de saber se no caso de parcelas perdidas é só substituir o valor perdido par NA e usar a função aov() como no exemplo abaixo. fator1<-c(rep("A",12),rep("B",12)) fator2<-c(rep("C",4),rep("D",4),rep("E",4),rep("C",4),rep("D",4),rep("E",4)) resposta<-c(39.03794, 30.19444, 32.57167, 34.87484, NA ,37.15130, 39.67066, 32.34696, 19.96688, 63.33362, 43.74447, 47.39147, 37.37704, NA, 29.35858, 28.76528, 16.092950, 21.080503, 16.906132, 12.345620, 26.475927, 23.639545, 23.890539, NA) dados<-data.frame(fator1,fator2,resposta) anova<-aov(resposta~fator1*fator2) summary(anova) Existe algum pacote que realiza testes de comparação múltipla para dados desbalanciados. Att. Tiago. _______________________________________________ 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. _______________________________________________ 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.
# install.pakages("easyanova") require(easyanova) resultado=ea2(dados) resultado names(resultado) ________________________________ De: tiago souza marçal <tiagosouzamarcal@hotmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Domingo, 25 de Novembro de 2012 20:07 Assunto: [R-br] Fatorial com parcelas perdidas Prezados, gostaria de saber se no caso de parcelas perdidas é só substituir o valor perdido par NA e usar a função aov() como no exemplo abaixo. fator1<-c(rep("A",12),rep("B",12)) fator2<-c(rep("C",4),rep("D",4),rep("E",4),rep("C",4),rep("D",4),rep("E",4)) resposta<-c(39.03794, 30.19444, 32.57167, 34.87484, NA ,37.15130, 39.67066, 32.34696, 19.96688, 63.33362, 43.74447, 47.39147, 37.37704, NA, 29.35858, 28.76528, 16.092950, 21.080503, 16.906132, 12.345620, 26.475927, 23.639545, 23.890539, NA) dados<-data.frame(fator1,fator2,resposta) anova<-aov(resposta~fator1*fator2) summary(anova) Existe algum pacote que realiza testes de comparação múltipla para dados desbalanciados. Att. Tiago. _______________________________________________ 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.
Obrigado pela ajuda Emmanuel, realmente as somas de quadrados foram bem diferentes das que eu obtive no CMR que enviei. Então a função aov() não deve ser utilizada quando temos dados desbalanceados? Att. Tiago. Date: Mon, 26 Nov 2012 04:59:32 -0800 From: emmanuelarnhold@yahoo.com.br To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Fatorial com parcelas perdidas # install.pakages("easyanova") require(easyanova) resultado=ea2(dados) resultado names(resultado) De: tiago souza marçal <tiagosouzamarcal@hotmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Domingo, 25 de Novembro de 2012 20:07 Assunto: [R-br] Fatorial com parcelas perdidas Prezados, gostaria de saber se no caso de parcelas perdidas é só substituir o valor perdido par NA e usar a função aov() como no exemplo abaixo. fator1<-c(rep("A",12),rep("B",12)) fator2<-c(rep("C",4),rep("D",4),rep("E",4),rep("C",4),rep("D",4),rep("E",4)) resposta<-c(39.03794, 30.19444, 32.57167, 34.87484, NA ,37.15130, 39.67066, 32.34696, 19.96688, 63.33362, 43.74447, 47.39147, 37.37704, NA, 29.35858, 28.76528, 16.092950, 21.080503, 16.906132, 12.345620, 26.475927, 23.639545, 23.890539, NA) dados<-data.frame(fator1,fator2,resposta) anova<-aov(resposta~fator1*fator2) summary(anova) Existe algum pacote que realiza testes de comparação múltipla para dados desbalanciados. Att. Tiago. _______________________________________________ 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. _______________________________________________ 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.
Tanto lm() quanto aov() podem ser usadas com dados desbalanceados. Acontece que o usuário deve estar ciente das conseguências do desbalanceamento: (1) diferença na precisão dos efeitos e (2) não ortogonalidade entre as fontes de variação. Ou seja, não ortogonalidade implica que as somas de quadrado sequencial (padrão do R) mudam de acordo com a ordem dos termos no modelo e que as médias amostrais não são estimadores livres de efeitos (deve-se usar médias ajustadas) e a diferença de precisão implica que você não deve usar métodos baseados em diferença mínima significativa única, com LSD, HSD, etc para comparar tais médias. Nestes casos é recomendado fazer as comparações de interesse e corrigir o p-valor devido à multiplicidade adotando por exemplo a correção de Bonferroni ou similar. Muitas opções estão disponíveis no pacote multcomp. Para médias ajustadas pode-se usar a doBy::popMeans(). À 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 skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Obrigado pela ajuda Walmes, vou ler a documentação do pacote. Att. Tiago. Date: Mon, 26 Nov 2012 13:21:03 -0200 From: walmeszeviani@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Fatorial com parcelas perdidas Tanto lm() quanto aov() podem ser usadas com dados desbalanceados. Acontece que o usuário deve estar ciente das conseguências do desbalanceamento: (1) diferença na precisão dos efeitos e (2) não ortogonalidade entre as fontes de variação. Ou seja, não ortogonalidade implica que as somas de quadrado sequencial (padrão do R) mudam de acordo com a ordem dos termos no modelo e que as médias amostrais não são estimadores livres de efeitos (deve-se usar médias ajustadas) e a diferença de precisão implica que você não deve usar métodos baseados em diferença mínima significativa única, com LSD, HSD, etc para comparar tais médias. Nestes casos é recomendado fazer as comparações de interesse e corrigir o p-valor devido à multiplicidade adotando por exemplo a correção de Bonferroni ou similar. Muitas opções estão disponíveis no pacote multcomp. Para médias ajustadas pode-se usar a doBy::popMeans(). À disposição. Walmes. ========================================================================== Walmes Marques ZevianiLEG (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 1173e-mail: walmes@ufpr.br skype: walmeszevianitwitter: @walmeszevianihomepage: http://www.leg.ufpr.br/~walmes linux user number: 531218========================================================================== _______________________________________________ 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.
Você pode usar aov. A diferença esta na função para obter a análise de variância. Você usaou anova(). Veja a função Anova() do pacote car. ?Anova Você também pode alterar o tipo de contraste na aov() ou na lm(). Observe na CMR abaixo as diferenças. fator1<-as.factor(c(rep("A",12),rep("B",12))) fator2<-as.factor(c(rep("C",4),rep("D",4),rep("E",4),rep("C",4),rep("D",4),rep("E",4))) resposta<-c(39.03794, 30.19444, 32.57167, 34.87484, NA ,37.15130, 39.67066, 32.34696, 19.96688, 63.33362, 43.74447,47.39147, 37.37704, NA, 29.35858, 28.76528, 16.092950, 21.080503, 16.906132, 12.345620, 26.475927, 23.639545, 23.890539, NA) dados<-data.frame(fator1,fator2,resposta) require(easyanova) m1<-aov(resposta~fator1*fator2, data=dados) m2<-aov(resposta~fator1*fator2, data=dados, contrasts=list(fator1=contr.sum, fator2=contr.sum)) resultado<-ea2(dados) anova(m1) anova(m2) Anova(m1, type=3) Anova(m2, type=3) resultado[1] ________________________________ De: tiago souza marçal <tiagosouzamarcal@hotmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Segunda-feira, 26 de Novembro de 2012 11:38 Assunto: Re: [R-br] Fatorial com parcelas perdidas Obrigado pela ajuda Emmanuel, realmente as somas de quadrados foram bem diferentes das que eu obtive no CMR que enviei. Então a função aov() não deve ser utilizada quando temos dados desbalanceados? Att. Tiago. ________________________________ Date: Mon, 26 Nov 2012 04:59:32 -0800 From: emmanuelarnhold@yahoo.com.br To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Fatorial com parcelas perdidas # install.pakages("easyanova") require(easyanova) resultado=ea2(dados) resultado names(resultado) ________________________________ De: tiago souza marçal <tiagosouzamarcal@hotmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Domingo, 25 de Novembro de 2012 20:07 Assunto: [R-br] Fatorial com parcelas perdidas Prezados, gostaria de saber se no caso de parcelas perdidas é só substituir o valor perdido par NA e usar a função aov() como no exemplo abaixo. fator1<-c(rep("A",12),rep("B",12)) fator2<-c(rep("C",4),rep("D",4),rep("E",4),rep("C",4),rep("D",4),rep("E",4)) resposta<-c(39.03794, 30.19444, 32.57167, 34.87484, NA ,37.15130, 39.67066, 32.34696, 19.96688, 63.33362, 43.74447, 47.39147, 37.37704, NA, 29.35858, 28.76528, 16.092950, 21.080503, 16.906132, 12.345620, 26.475927, 23.639545, 23.890539, NA) dados<-data.frame(fator1,fator2,resposta) anova<-aov(resposta~fator1*fator2) summary(anova) Existe algum pacote que realiza testes de comparação múltipla para dados desbalanciados. Att. Tiago. _______________________________________________ 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. _______________________________________________ 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. _______________________________________________ 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.
Agradeço novamente por sua ajuda Emmanuel. Att. Tiago. Date: Mon, 26 Nov 2012 11:32:08 -0800 From: emmanuelarnhold@yahoo.com.br To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Fatorial com parcelas perdidas Você pode usar aov. A diferença esta na função para obter a análise de variância. Você usaou anova(). Veja a função Anova() do pacote car. ?Anova Você também pode alterar o tipo de contraste na aov() ou na lm(). Observe na CMR abaixo as diferenças. fator1<-as.factor(c(rep("A",12),rep("B",12))) fator2<-as.factor(c(rep("C",4),rep("D",4),rep("E",4),rep("C",4),rep("D",4),rep("E",4))) resposta<-c(39.03794, 30.19444, 32.57167, 34.87484, NA ,37.15130, 39.67066, 32.34696, 19.96688, 63.33362, 43.74447,47.39147, 37.37704, NA, 29.35858, 28.76528, 16.092950, 21.080503, 16.906132, 12.345620, 26.475927, 23.639545, 23.890539, NA) dados<-data.frame(fator1,fator2,resposta) require(easyanova) m1<-aov(resposta~fator1*fator2, data=dados) m2<-aov(resposta~fator1*fator2, data=dados, contrasts=list(fator1=contr.sum, fator2=contr.sum)) resultado<-ea2(dados) anova(m1) anova(m2) Anova(m1, type=3) Anova(m2, type=3) resultado[1] De: tiago souza marçal <tiagosouzamarcal@hotmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Segunda-feira, 26 de Novembro de 2012 11:38 Assunto: Re: [R-br] Fatorial com parcelas perdidas Obrigado pela ajuda Emmanuel, realmente as somas de quadrados foram bem diferentes das que eu obtive no CMR que enviei. Então a função aov() não deve ser utilizada quando temos dados desbalanceados? Att. Tiago. Date: Mon, 26 Nov 2012 04:59:32 -0800 From: emmanuelarnhold@yahoo.com.br To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Fatorial com parcelas perdidas # install.pakages("easyanova") require(easyanova) resultado=ea2(dados) resultado names(resultado) De: tiago souza marçal <tiagosouzamarcal@hotmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Domingo, 25 de Novembro de 2012 20:07 Assunto: [R-br] Fatorial com parcelas perdidas Prezados, gostaria de saber se no caso de parcelas perdidas é só substituir o valor perdido par NA e usar a função aov() como no exemplo abaixo. fator1<-c(rep("A",12),rep("B",12)) fator2<-c(rep("C",4),rep("D",4),rep("E",4),rep("C",4),rep("D",4),rep("E",4)) resposta<-c(39.03794, 30.19444, 32.57167, 34.87484, NA ,37.15130, 39.67066, 32.34696, 19.96688, 63.33362, 43.74447, 47.39147, 37.37704, NA, 29.35858, 28.76528, 16.092950, 21.080503, 16.906132, 12.345620, 26.475927, 23.639545, 23.890539, NA) dados<-data.frame(fator1,fator2,resposta) anova<-aov(resposta~fator1*fator2) summary(anova) Existe algum pacote que realiza testes de comparação múltipla para dados desbalanciados. Att. Tiago. _______________________________________________ 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. _______________________________________________ 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. _______________________________________________ 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. _______________________________________________ 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.
A ordem com que você declara não muda a precisão dos efeitos estimados.
Ajuste com inversão de ordem e veja o erro padrão summary() do modelo. É o
mesmo.
> d1 <- read.table("http://www.leg.ufpr.br/~walmes/data/montgomery_14-1.txt
",
+ header=TRUE, sep="\t",
colClasses=c("factor","factor","integer"))
> str(d1)
'data.frame': 18 obs. of 3 variables:
$ vidro : Factor w/ 2 levels "1","2": 1 1 1 2 2 2 1 1 1 2 ...
$ fosforo : Factor w/ 3 levels "800","825","850": 1 1 1 1 1 1 2 2 2 2 ...
$ corrente: int 280 290 285 230 235 240 300 310 295 260 ...
> d1$corrente[1] <- NA
> m0 <- lm(corrente~vidro*fosforo, d1)
> m1 <- lm(corrente~fosforo*vidro, d1)
> anova(m0)
Df Sum Sq Mean Sq F value Pr(>F)
vidro 1 14309.0 14309.0 264.1663 4.884e-09 ***
fosforo 2 *869.5* 434.7 8.0261 0.007088 **
vidro:fosforo 2 81.5 40.8 0.7527 0.493861
Residuals 11 595.8 54.2
> anova(m1)
Df Sum Sq Mean Sq F value Pr(>F)
fosforo 2 *1119.2* 559.6 10.3312 0.002983 **
vidro 1 14059.3 14059.3 259.5560 5.361e-09 ***
fosforo:vidro 2 81.5 40.8 0.7527 0.493861
Residuals 11 595.8 54.2
> summary(m0)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 287.5000 5.2042 55.244 8.43e-15 ***
vidro2 -52.5000 *6.7185* -7.814 8.16e-06 ***
fosforo825 14.1667 6.7185 2.109 0.0587 .
fosforo850 0.8333 6.7185 0.124 0.9035
vidro2:fosforo825 -4.1667 9.0139 -0.462 0.6529
vidro2:fosforo850 -10.8333 9.0139 -1.202 0.2547
> summary(m1)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 287.5000 5.2042 55.244 8.43e-15 ***
fosforo825 14.1667 6.7185 2.109 0.0587 .
fosforo850 0.8333 6.7185 0.124 0.9035
vidro2 -52.5000 *6.7185* -7.814 8.16e-06 ***
fosforo825:vidro2 -4.1667 9.0139 -0.462 0.6529
fosforo850:vidro2 -10.8333 9.0139 -1.202 0.2547
À 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
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
participantes (5)
-
andrebvs -
Emmanuel Arnhold -
Fernando Antonio de souza -
tiago souza marçal -
Walmes Zeviani