ANOVA de duas vias no R (Parte 1)

[@cienciaabessa]

Aê! O rstatix saiu!

[@lucimariacostapereira7108]

Guriaaaa to tão felizzzz que tu voltou a fazer videos, pelo amorrrrr de DEUS não pare maissss

[@FernandaPeres]

Hahaha! Pode deixar! Eu voltei para ficar, rs. Estou também produzindo conteúdo no Instagram, se interessar: instagram.com/estatisticaaplicada Abraços!

[@andrepegoraro]

Fernanda, muito obrigado. Você realmente tem didática e oratória muito boas. Estou cursando a disciplina de Estatística e Modelos de Previsão do curso de Engenharia de Produção Mecânica da UFSC e, desde o começo do semestre eu tenho assistido seus vídeos como complementos e tira-dúvidas. Certamente recomendarei seu canal para todos os próximos alunos da disciplina, assim mais pessoas passam bem capacitadas pela mesma. As vezes, como produtora de conteúdo, deve ser difícil perceber o tamanho do impacto das suas atitudes e seus vídeos. Venho aqui te dizer que você faz um belo serviço para a educação brasileira com isso. Sucesso!

[@higoralmeida7441]

Parabéns pelo trabalho! Conteúdo incrível!

[@alexdonascimento3312]

Estava esperando o vídeo exatamente de Anova de duas vias. E muito obrigado por sanar a duvida, sobre os pressupostos, sempre vejo por ai falando, Dados brutos vs Resíduos. Como sempre, excelente vídeo aula.

[@rafaelscari8449]

Didática e conteúdo de muita qualidade. Obrigado! Sucesso!

[@Joao99444276]

Fernanda Peres, muito bom! Amo suas aula...!

[@carolinabender]

Esse é o melhor canal de estatística, sem dúvidas!!! 😁 Tem algum vídeo (ou dica) para fazer ANOVA de duas vias para within-subjects design?

[@FernandaPeres]

Eu tenho vídeo explicando isso, chama ANOVA de medidas repetidas :)

[@anap.6299]

Eu estou fazendo mestrado e preciso rodar minhas análises no R. Não sei o que faria sem seus vídeos pra me ajudar. Vira e mexe assisto outros, mas só consigo aprender contigo explicando. Agora já vou direto na sua playlist... só vou pra outro canal se não tiver o assunto lá. Muuuuuuito obg por compartilhar.

[@FernandaPeres]

Ahhh, que feedback gostoso de ler! Obrigada!!! Fico feliz de te ajudar ❤️

[@vicenteneto2862]

Muito obrigado pela aula, muito boa e compartilhando nos grupos. Muito obrigado!!!

[@pedroabneiva]

Que meu professor ruim n me veja aqui :D Muito bom!!!!

[@eltoncorrea8319]

Muito obrigado pela aula!! 💻📊

[@didimarinhoo]

parabéns ! ótima explicação

[@samarabarbosa5335]

Estou adorando seu canal Fernanda. Tem me ajudado bastante a iniciar no R. Suas explicações são claras e objetivas e as informações muito completas. Muito Obrigada e Parabéns pelo trabalho. Seria possível gravar um vídeo de como trabalhar com a análise fatorial com adicional no R? Sou da área de ciências agrárias e muitos de meus experimentos requerem este tipo de análise por trabalhar com tratamento "controle".

[@FernandaPeres]

Oi, Samara, está nos planos, mas vou demorar um pouco para eu gravar esses vídeos...

[@julianavieira5891]

Obrigada pelo conteúdo! Seus vídeos de ajudam muito! Qual teste não parâmetro que substitui a ANOVA de duas vias?

[@brunaramosdesouza5987]

Tem algum vídeo que ensina ANOVA de 3 fatores?

[@biancagomes9516]

Adoro suas aulas! Tem como fazer um vídeo ensinando Anova de três vias no 3?

[@ayrtonsennan]

Ou queria saber se posso usar dessa forma para analisar um teste sensorial, cuja as fontes de variações são os provadores e as amostras?

[@ritag.m9208]

Parabéns Patrícia excelente explicação. Tenho uma dúvida, e quando os dados não atendem os pressupostos de normalidade e homogeneidade dos dados mesmo com transformação, o que você sugere que seja feito? Um modelo linear generalizado? Ou uma regressão logística utilizando o pacote lme4? O que fazer?

[@FernandaPeres]

Se a variável é quantitativa, modelo linear generalizado com alguma distribuição que faça sentido para ela. Logística só para variável dependente dicotômica. E a que o lme4 faz é para medidas repetidas.

[@ritag.m9208]

​@@FernandaPeresEntendi, obrigada!!!!!!! ❤❤❤

[@milenaoliveira2626]

Oi Fernanda, tudo bem?! Amo suas aulas!! Por favor me tira essa dúvida: Caso a interação entre dois fatores (Tratamento *Combinações) não tenha sido significativa (p>0.05), posso utilizar o resultado do Tukey HSD (letras) mesmo assim, num gráfico de barras, por exemplo? Houve interação para tratamento. Estou utilizando o Programa JMP.

[@caiopluvier8181]

Que incrível! Parabéns pelo trabalho! Será que poderia ensinar a fazer curva ROC no R, para os próximos vídeos? :D

[@FernandaPeres]

Obrigada, Caio! Vou gravar, sim. Mas no vídeo de regressão logística, então vai demorar um pouquinho... Abraços!

[@jorge3633]

Tenho uma dúvida. Sou muito leigo com estatistica. Tenho uns dados que são 3 grupos e 1 de controle, cada grupo possui dados de 3 órgãos diferentes. Fiz no R aov(valores ~ grupo * orgaos) mas não atenderam os pressupostos de normalidade. Qual análise eu preciso fazer? Um Kruskal wallis analisando a valores ~ interation(grupo:orgaos) seria o certo? Desde já, agradeço

[@verenacarmelo8086]

Fernanda, tem video sobre anova de MR?

[@FernandaPeres]

Medidas repetidas? Tem, sim. Tem de medidas repetidas e de delineamento misto.

[@fredericobaptista7103]

Olá, Fernanda. Muito obrigado pela partilha deste conteúdo. Entretanto, tenho uma dúvida. O banco de dados que você utiliza diz respeito a um estudo transversal, ou seja, os dados da memória foram todos recolhidos num único momento do tempo e as variáveis independentes são o gênero e a quantidade de álcool consumida neste mesmo momento. Se estivéssemos a analisar os dados de um estudo longitudinal onde as variáveis independentes fossem o grupo (experimental versus controlo) e o tempo (T0 e T1), como que a minha base de dados deveria estar organizada? Digo isso, pois na minha base de dados tenho uma coluna para o grupo, uma coluna para os valores da intensidade da dor no momento T0 e uma outra coluna para os valores da intensidade da dor no momento T1. Neste caso, penso que a verificação dos pressupostos deve levar em consideração apenas os valores da intensidade da dor para o momento T0 (baseline). Mas na hora de construir o modelo da ANOVA propriamente dito, como faço para organizar essas colunas, considerando que terei de usar os valores de T0 e de T1? Não sei se me fiz entender 🙈

[@FernandaPeres]

Oi, recomendo um vídeo que eu tenho sobre organização da base de dados. É um dos últimos publicados. E, você vai fazer uma ANOVA com medidas repetidas (ou ANOVA mista, a depender do delineamento). Assistindo aos vídeos desses testes você vai entender como organizar a base e avaliar os pressupostos.

[@fredericobapt]

Perfeito 😊 Obrigado!

[@influence8811]

O teste de Levene é o mais indicado para ANOVA em esquema fatorial?

[@jessica123456679]

Boa tarde Fernanda, muito obrigada. Uma dúvida, alguns dos meus dados não atenderam as premissas da ANOVA Fatorial, mesmo após a transformação em log. Optei por usar o GLM para esses dados, você sabe me informar se posso incluir a interação no modelo do glm (ex. y ~ x1 +x2 + x1*x2)? Desde já agradeço!!

[@FernandaPeres]

Pode, sim!

[@jessica123456679]

Muito obrigada!!@@FernandaPeres

[@brunasa9747]

Olá meu bem, eu estou amando estudar por seus vídeos. São muito bem explicativos, principalmente no R, que é o programa que muita gente ama, mas tem dificuldade em mexer (eeeeuuuuu rsrsr). Queria saber se você pretende lançar algum vídeo ou se você tem alguma material que possa me ajudar com a PERMANOVA. Os materiais que eu estou tendo acesso não estão me ajudando muito :/ obrigada e mais uma vez, parabéns!

[@rafaelarauber2400]

Adorei o vídeo, mas não consegue encontrar o pacote do Shapiro_test, poderia escrever eles no script porque pelo sua fala eu não consegui identificar o nome com precisão. Ou colocar na descrição os pacotes que usou.

[@FernandaPeres]

Tá no script, como todos os outros pacotes usados. É o rstatix.

[@alexdonascimento3312]

Fernanda !! sempre tenho que atender as duas condições, a normalidade e homogeneidade de variâncias. porem caso minha homogeneidade seja diferente, mas a normalidade passou, poderia fazer o Anova do mesmo jeito ? *Essa parte sempre tenho duvidadas hahahaha Obrigado ^^ Você salva vidas!!! com seus videos :)

[@FernandaPeres]

Oi, Alex. Então, já ouvi professores dizerem que precisa da normalidade ou da homogeneidade, mas nunca achei uma referência que corroborasse isso. Tem referência que diz que a ANOVA é robusta a violações da normalidade quando o n é grande e a distribuição não é muito assimétrica. Mas, não é robusta a violações de homogeneidade. Nesses casos, uma opção é transformar os dados (por exemplo, logaritmizar) e verificar se isso resolve o problema da homogeneidade. Se resolver, então você roda a ANOVA com os dados transformados. Tem também a opção de fazer testes que não tenham esse pressuposto. O R deixa fazer ANOVAs de duas vias robustas, as quais que eu saiba não têm esse pressuposto (quero estudar um pouco mais antes de gravar vídeo sobre elas). Se for o seu caso, talvez compense dar uma estudada nisso. O livro do Andy Field (Descobrindo a estatística usando o R) ensina a fazer. Que ótimo que o canal te ajuda!! Abraços!

[@alexdonascimento3312]

@@FernandaPeres Nossa!!! Muito obrigado pela atenção. Irei procurar esse livro com certeza. Esse canal (você) salva a vida nas horas de compreender o "Mundo da Estatística"

[@vinnyciusnunes1019]

Fernanda, aditividade do modelo é um pressuposto? Algumas referências dizem que sim.

[@FernandaPeres]

Olha, as referências que eu uso (Livros: Descobrindo a estatística usando o SPSS, Statistical reasoning in the behavioral sciences) não falam de aditividade ser um pressuposto. Eu dei uma pesquisada, e tem fontes que falam isso, mas pelo que entendi, esse pressuposto só vale para modelos sem interação. Os que eu construo são sempre com interação.

[@jacksonlima6057]

Fernanda, você poderia fazer ANOVA em DBC com comparação de médias?

[@gabrielachavezcossio8027]

Oi Fernanda parabéns pelos vídeos. Um favor você poderia me informar quais seriam esses pacotes que fazem ANOVA de duas vias robustas? não consegui transforma-los em paramétricos, Obrigada.

[@FernandaPeres]

Oi, Gabriela, o pacote é o WRS2. Mas recomendo ir seguindo o tutorial do livro do Andy Field para R que é tranquilo de fazer, ele explica inclusive o post hoc. O livro é o "Discovering Statistics using R". Também tem outros tutoriais desse pacote na internet.

[@jorgemmmmteixeira]

Alô mestre. Pergunta de noob, mas neste caso o factor "tempo" tem características de variável contínua ou é uma variável discreta? Grato.

[@FernandaPeres]

Oi, Jorge, tudo bem? Eu não entendi muito bem a qual variável tempo você se refere... Você diz a variável latência (que eu não usei nesse exemplo)?

[@jorgemmmmteixeira]

@@FernandaPeres tudo bem e consigo, obrigado? Refiro-me ao "tempo" que passou entre a primeiro e segunda medição. Apesar de podermos quantificar esse período (horas, dias, semanas), devemos encarar esse factor como sendo categórico (simplesmente são períodos distintos)?

[@jorgemmmmteixeira]

Ou então minha pergunta não faz sentido nenhum e é melhor ignorar. Também é possível 😁

[@eduardopereirad]

Boa noite! Qual o método mais adequado quando eu tenho duas variáveis independentes e uma dependente? Ex.: 4 tipos de exposição, 4 tipos de inseticida e variável resposta consumo de presa.

[@FernandaPeres]

ANOVA de duas vias mesmo.

[@eduardopereirad]

@@FernandaPeres e se nesse mesmo exemplo eu tivesse ao invés de uma, 5 variáveis resposta?

[@FernandaPeres]

Aí ou você analisa separadamente, uma ANOVA para cada, ou, se fizer sentido teórico, faz uma MANOVA.

[@eduardopereirad]

@@FernandaPeres Certo. Obrigado!

[@kateyaas]

Oii, Fernanda! O que eu faço se o R acaba reconhecendo minha variável numérica que seria dbl como um fator? Tem algo que eu possa fazer pra melhorar esse reconhecimento? Acabo não conseguindo fazer os testes porque ele pede variável numérica

[@FernandaPeres]

Você tem como converter um fator em variável numérica. Dá uma pesquisada, acho que com as.numeric(). Mas é estranho o R reconhecer um número como fator. Verifique se você está declarando o separador de decimal correto, se em algum dos casos não tem uma letra...

[@FernandaPeres]

Ah, essa parte de definir separador de decimal eu ensino nos primeiros vídeos da playlist.

[@juliacordeiro1212]

Como fazer quando houver mais de 2 variáveis independentes?

[@FernandaPeres]

É só inserir mais variáveis no modelo. Eu fiz de duas vias, mas não há limites.

[@juliacordeiro1212]

@@FernandaPeres Ao rodar o shapiro_test, aparece o seguinte erro: Erro: Problem with `mutate()` input `data`. x Problem with `mutate()` input `data`. x sample size must be between 3 and 5000 i Input `data` is `map(.data$data, .f, ...)`. i Input `data` is `map(.data$data, .f, ...)`. Run `rlang::last_error()` to see where the error occurred. O que fazer? Existe algum comando que permita verificar a normalidade por grupo quando as amostras são muito grandes?

[@FernandaPeres]

@@juliacordeiro1212 Isso é mesmo uma limitação do teste de Shapiro-Wilk, cujo cálculo do p não é validado para amostras maiores que 5000. Nesses casos, a sugestão é avaliar a normalidade por gráficos: histograma ou qqplot são boas opções.

[@joabmartins7095]

Quando os dados são não normais, qual o teste não paramétrico equivalente à ANOVA de duas vias?

[@FernandaPeres]

Então, não existe um teste não-paramétrico. Mas há no R pacotes que fazem ANOVA de duas vias robustas (usando a mediana ou a média aparada). Ainda vou gravar vídeos com eles. Fora isso, você tem a opção de transformar os dados ou usar um modelo que não tenha normalidade como pressuposto (como os modelos lineares generalizados).

[@tribology]

@@FernandaPeres boa tarde Fernanda parabéns pelos seu vídeos. Você poderia me informar quais seriam esses pacotes que fazem ANOVA de duas vias robustas? Estou com um questionário Likert de 2 fatores que não consegui transforma-los em paramétricos. Abraços e parabéns mais uma vez.

[@potiramelo9193]

@@FernandaPeres o software past utiliza o permanova two-way

[@caduguimaraes]

Crie um github....