ANOVA de uma via com medidas repetidas no R
Vložit
- čas přidán 7. 09. 2024
- O teste ANOVA de uma via com medidas repetidas permite que você teste o efeito de uma variável independente intra-sujeito sobre uma variável dependente. É um teste muito usado para avaliar alterações ao longo do tempo.
Nesse vídeo, eu explico como realizar no R esse teste e as comparações entre pares apropriadas a ele (com correção de Bonferroni). Além disso, explico uma parte fundamental para realizar o teste: como transpor o banco de dados do formato "wide" ("amplo") para o "long" ("longo"). Ao final, como sempre, te explico como reportar os resultados.
Download do script:
drive.google.c...
Download do banco de dados:
drive.google.c...
Quer me contratar para uma consultoria ou para analisar os seus dados? No link abaixo tem a explicação desses serviços e um formulário para você entrar em contato comigo: fernandafperes...
Precisa fazer esse teste, mas não é fã do R? Esse vídeo ensina a fazê-lo no SPSS:
• ANOVA 1 via com medida...
Acompanhe o conteúdo de estatística também pelo Instagram:
/ estatisticaaplicada
Acompanhe os cursos de estatística aplicada:
No R: • Curso de Estatística A...
No SPSS: • Curso de Estatística A...
Aula extremamente clara! Sua didática é difícil de ver por aqui! Obrigado pelo trabalho, Fernanda!
Fernanda! Você é maravilhosa... suas aulas são muito esclarecedoras e muito, muito didáticas. Tenho aprendido demais com elas. Parabéns!!!
Que legal, no aguardo do método de duas variáveis independentes em medidas repetidas ....... parabéns pela excelente aula...
Obrigada, Manoel! Vai ter vídeo com duas VI, sim. Abraços!
Didatica incrivel! Vai direto ao ponto!
Fernanda, excelente trabalho! Muito obrigado por compartilhar este conteúdo!
suas aulas são incríveis não pare por favor !
Parabéns pelo seu incrível trabalho!!!!!!
Ótimo conteúdo!!
Para a instalação e library uso essa função, ajuda bastante pois já faz algumas validações:
ipak
Excelente explicação! Lembrando que dá pra customizar o boxplot com a função do r base deixando ele publicável.
Ah, até dá, mas com ggplot2 fica tão mais bonito, né 🥰
MARAVILHOSA!!! Muito obrigada pela partilha de conhecimentos
Descobri seu canal e estou amando os vídeos e aprendendo bastante sobre algumas análises que eu não sabia como realizar no R. Seus vídeos são ótimos!! E obrigado por compartilhar essas informações!
Excelente Fernanda! Super didática 👏🏼👏🏼👏🏼👏🏼👏🏼
Ótimo vídeo, muito explicativo. E ótimo humor. Parabéns!
Hahaha! Obrigada!! :)
Aula muito boa!Parabéns.
Muchas gracias
Excelente vídeo! Tem como fazer uma aula teórica sobre os principais características do teste de Multicomparação de Tukey?
O professor 1 é o cara, todos os profs. tinham que ser como ele. hahahhaa
Hahahahah! Não acho que eu sou a professora 1, não, tem que ser um pouco rígido 💁🏻♀️
Excelente!!!
Parabéns pelo conteúdo, adoro os seus vídeos. No caso se as premissas de normalidade não fossem atendidas, teria alguma análise semelhante não paramétrica ou como você seguiria??
Para a ANOVA de medidas repetidas com uma via, tem uma alternativa não paramétrica que é o teste de Friedman.
Mas, para delineamentos mais complexos, as alternativas são GEE ou modelos lineares generalizados mistos.
Olá Fernanda! Seus vídeos são ótimos. Você recomenda algum livro que aborde esses assuntos? Mais uma vez, parabéns pelo vídeos!
Parabéns pela aula!!! faz um tempo que estava e busca dessa temática na plataforma R.
Tenho mais umas dúvidas, não sei se posso fazer a pergunta por aqui.
Como efetuar o Intervalo de Confiança para a diferença na comparação entre os pares, abordado do Teste Post Hoc?
Obrigado e mais uma vez, parabéns!!
gêniaaaa
Parabéns pelos vídeos Fernanda, são bem legais! Uma dúvida, neste tipo de anova existe algum teste post-hoc em que seja possível ver quais alunos tiveram notas diferentes de acordo com o professor?
Excelente! quando vai fazer vídeos de tamanho do efeito? tem com você fazer alfo de diferença mínima detectável?
Professora, quando há necessidade de remoção dos outliers, eu opto por remover os valores ou o individuo?
Parabéns pela clareza! Qual o software você utiliza para criar os vídeos?
Obrigada, Fernanda. Gravo com o OBS Studio e edito com o HitFilm.
Fernanda uma dúvida.. Neste caso com os professores significativos.. A efeito de comparação poderia se aplicar teste tukey como pos-hoc?
Acredito que sim. Eu nunca usei, mas deve ter função. Só é importante garantir que o post-hoc está considerando que é medidas repetidas ;)
Oieee professora....amei sua aula!!
Uma dúvida, se os meus dados não atenderem as premissas do teste, como ocorreu no seu exemplo, não é possivel fazer uma correção, por exemplo transformar os valores em log?
No meu caso, só um fator não apresentou normalidade, por isso pensei em fazer a transformação ou utilizar o GLMM, caso não seja possivel.
Desde já agradeço!!
Oi! Nesse caso você tem três opções:
- Usar um teste não-paramétrico (o equivalente seria o Friedman)
- Transformar os dados (por exemplo, por logaritmização)
- Usar um modelo linear generalizado
Todas essas opções são válidas e têm suas vantagens e desvantagens. Eu tendo a preferir não transformar os dados, porque isso dificulta a interpretação dos resultados.
@@FernandaPeres Muito obrigadaa!!!
Olá Fernanda, você acha que seria correto analisar o efeito do professor utilizando modelos mistos (professor fixo e ID aleatório)?
Se sim, está correto dizer que não precisaria testar a esfericidade pois em modelos mistos há correção por indivíduo
Fernanda, na minha base de dados eu tenho variáveis que são contínuas e variáveis que são categóricas. Quando eu reestruturo o banco de dados para o formato "long", todas as variáveis (categóricas e contínuas) ficam listadas na segunda coluna do banco de dados e os respetivos valores/categorias listados na terceira coluna do banco de dados. Portanto, quando tento verificar os pressupostos de normalidade e outliers, aparece uma mensagem de erro e penso que é porque existem variáveis contínuas e categóricas todas misturadas. Qual sua sugestão para resolver esse problema? Devo ter uma base de dados só com as variáveis contínuas? Ou é possível dizer ao R quais são as variáveis que me interessam?
Já consegui encontrar a função que resolve isso no pacote tidyverse. Se alguém precisar, aqui está a função:
# Remover do banco de dados todas as colunas com variáveis categóricas (pacote: tidyverse)
dados_numeric % select(where(is.numeric))
Fernanda, me surgiu uma dúvida.. Posso usar esse teste para avaliar se uma pessoa que irá marcar os pontos de uma morfometria geométrica está calibrado? No caso seria, o mesmo avaliador marcando os pontos em um mesmo conjunto de dados após um intervalo de tempo.. reunir esses dados e fazer a análise para ver se há diferença entre os pontos que foram marcados.
como faço uma anova de duas vias com medidas repetidas no R?
Fernanda você ministra aulas particulares? Confesso que fiquei em dúvida em aplicar quando tem um ensaio comum DBC que tenha estas medidas repetidas no tempo. Se sim. como poderia entrar em contato?
Fernanda, fiz ANOVA de medidas repetidas para uma distribuição normal.
Onde a minha variável de agrupamento tem 4 categorias.
E dai do Mauchly's Test for Sphericity` para frente não apareceu nada.
Será que é porque os dados são esféricos?
Olá, Fernanda! Na função melt é preciso mesmo escrever o nome de todas as colunas? É que a minha base de dados tem 219 colunas (metade das colunas para os dados pré-intervenção e a outra metade para os dados do pós-intervenção). Existe alguma forma de automaticamente selecionar todas as colunas ou tenho mesmo de escrever uma a uma?
Faz um tempo que eu parei de usar a melt. Usa a pivot_longer do tidyr. Lá você coloca a numeração das colunas. Tem bastante documentação sobre ela na internet.
Acho que consegui. Escrevi o código apenas assim:
dados_long
@@FernandaPeres Obrigado :)
Como sempre, uma ótima aula vídeo de como usar o R.
Meus experimentos, basicamente são anovas de amostras repetiras, ao longo do tempo, (24, 48 e 72 horas, após aplicação dos produtos).
Gostaria de fazer uma pergunta, é possível realizar uma análise de amostras repetidas (anovas de amostras repetiras), porem com mais níveis? Exemplo.
Tenho 7 produtos que tenho que avaliar suas respostas (são 5) ao longo de 3 dias. Porem ao invés de analisar cada resposta separadamente, seria possível fazer uma anova de amostras repetidas com todas essas variáveis respostas? Com todas essas informações a qual você apresentou no vídeo (como falei, são ótimas suas explicações).
Exemplo dos dados:
24 horas 48 horas 72 horas
1a 2b 3c 4d 5e 1a 2b 3c 4d 5e 1a 2b 3c 4d 5e
Produto1
Produto2
Produto3
Produto4
Produto5
Produto6
Produto7
Oi, Alex, dá, sim. Mas aí você precisa fazer ANOVAs com mais de uma via. Nos vídeos de ANOVA com medidas repetidas no SPSS eu explico isso. O vídeo explicando isso em R vai demorar um pouquinho para sair...
Abraços!
@@FernandaPeres
Olá, Fernanda. Muito obrigado esclarecer essa dúvida. Ficarei aguardando o vídeo explicativo no ambiente R.
Parabéns, mais uma vez pelo seu trabalho.
Olá Fernanda, excelente material, acho que o arquivo do download está errado, grato
Não está, não. Eu atualizei o script para corrigir códigos que pararam de funcionar com o tempo.
Baita video. Quando tento modificar o formato do banco, sempre aparece um erro (Error in match.names(clabs, names(xi)): names not match previus names). Só que já conferir os nomes das colunas e está tudo certo, porém o erro segue aparecendo. Tenho uma dúvida. Estou trabalhando com um banco que contém NA, pode ser isso que esteja ocasionando esse erro? Não consigo encontrar o real motivo pra isso.
Olha, essa função de reestruturar o banco é meio chatinha, vive dando erro. É possível que sejam os NA, sim. Recomendo jogar o erro no Google, em geral alguém já levou o problema para o StackOverflow e lá tem a solução ;)
Oi Fernanda, e quando não atende ao pressuposto da normalidade?
Uma alternativa é fazer o teste não-paramétrico de Friedman
Ótimo vídeo, muito obrigado por deixar mais claro.
Tenho uma dúvida em relação com dados da minha análise. Na minha coluna "ID" eu tenho 3 níveis, por exemplo, alunos do meio da sala (MS), alunos no fundo da sala (BS) e alunos na frente da sala (FS). Em cada nível há 5 alunos que as notas foram analisadas por 4 períodos durante o ano. Como eu posso avaliar se a posição da sala interfere na nota ao longo do ano? Sigo o mesmo caminho do vídeo?
Nesse caso, você precisa fazer uma ANOVA mista. Ainda não tem vídeo sobre ela no canal, mas é o próximo que vai sair!
(Mas tem esse tipo de ANOVA no SPSS, se interessar)
Mas, basicamente, você vai ter que adicionar uma variável between subjects no modelo, na função da ezANOVA.
@@FernandaPeres muitooo obrigado pelo feedback professora Fernanda. Infelizmente eu não tenho acesso ao SPSS. Só ao R e o SYSTAT (versão 12.0)...vou esperar pelo seu vídeo! Muito obrigado pela luz!!!
Professora Fernanda, o teste de esfericidade também é um pressuposto para medidas repetidas onde temos apenas duas categorias?
É que nos outros vídeos não foi feito (Teste-T pareado e Wilcoxon)
Não é. Ele só é um pressuposto para quando há mais de duas medidas repetidas. Tanto que se você fizer uma ANOVA só com dois tempos, o SPSS nem calcula a esfericidade :)
@@FernandaPeres Perfeito, obrigado pela rápida reposta. Aprecio muito sua didática, está ótimo assim.
prezada, o arquivo csv não abriu e acusa que o formato e sua extensão não são iguais. CSV vs. SYLK. Pode upa-lo novamente em formato e extensão csv? obrigado e parabéns pelo trabalho.
Oi, Otávio, baixei aqui o banco e o script e estão funcionando normalmente. Verifique se o problema não é que está no formato em português, logo o separador de coluna é o ponto vírgula.
@@FernandaPeres era isso! resolvi importanto para txt e tratando os dados no excel, em dados> textos para colunas. obrigado.
Professora, a ANOVA possui alguma limitação com relação ao tamanho do banco de dados (número de observações)? Posso usar com bancos grandes e/ou pequenos?
Oi, Larissa, não tem uma limitação. Mas, você precisa atender aos pressupostos de normalidade, ausência de outliers e homogeneidade de variâncias. Em geral, esses pressupostos não são atendidos em amostras pequenas. Tem também uma discussão se o teste que verifica a normalidade é confiável em amostras pequenas, então tem pessoas que sugerem que ela seja usada com n > 30.
Oiee professora...espero que esteja bem!!
Estou com um problema, os meus dados atenderam as premissas, e o resultado da ANOVA deu uma diferença significativa mas o test post-hoc não mostrou essa diferença, como devo proceder?
Desde já agradeço!!!
Oii, isso pode acontecer mesmo. Nesse caso, a ANOVA deve ser considerada um "falso positivo". Esse artigo aqui discute isso: pubmed.ncbi.nlm.nih.gov/29719361/
@@FernandaPeres Muitoo obrigadaa...vou lerr...
Eu tive outro problema, ao pedir o resultado da ANOVA o teste de esfericidade não esta aparecendo, seria possível fazer ele separado (ao rodar o comando da anova apareceu um aviso de produção de NA, mas nos meus dados brutos não tem nenhum valor faltando)?
@@jessica123456679 Provavelmente porque você só tem dois momentos. Aí não tem como rodar teste de esfericidade, isso não é um pressuposto quando há apenas dois momentos.
@@FernandaPeres Entendiiiii...otimooo!!!
Fernanda, desculpa, mas no minuto 29:32 vc afirma que o p-valor entre professor 3 e 1 não é menor que 0,05.
No caso, o p-valor entre o professor 3 e o professor 1 é 0.01943, e nesse caso é menor que 0,05 e portanto tem uma diferença entre o professor 3 e 1 certo? Na primeira vez eu não enxerguei isso, mas refazendo eu observei isso, mas posso estar enganado.
Realmente, eu vi errado. É menor que 0,05, sim!
Olá professora! não consigo abrir o arquivo do excel, a senhora podia me dar uma orientação?
obrigado.
Ele vem no formato SKL, alguma coisa do tipo...
@@domingos21855 Ele está no formato .csv, como os outros bancos. Acabei de checar aqui, tá tudo certo com ele.
@@FernandaPeres Está bem, obrigado prof.
o banco de dados nao é o mesmo usado, profª
Atualizei o link. Obrigada pelo aviso!
Os arquivos do dropbox foram excluídos =(
Atualizei o link