Antecedentes:
eu gostaria de realizar uma meta-regressão usando estudos que tenham (1) vários resultados / construções (= multivariados) e (2) tamanhos de efeitos múltiplos para cada um desses resultados devido a diferentes medidas. Aqui está um esquema que, com sorte, explica melhor:
- Estudo 1, Resultado A, Tamanho do Efeito 1
- Estudo 1, Resultado A, Tamanho do Efeito 2
- Estudo 1, Resultado B, Tamanho do Efeito 3
- Estudo 2, Resultado A, Tamanho do Efeito 4
- Estudo 2, Resultado C, Tamanho do Efeito 5
- Estudo 2, Resultado C, Tamanho do Efeito 6
- ...
Os estudos comparam as médias de dois grupos com diferentes resultados e os tamanhos dos efeitos são g de Hedge.
Um exemplo prático seria "Memória de Trabalho", que pode ser dividida em diferentes resultados (Baddeley, 1974), por exemplo, "Laço Fonológico", "Bloco de Notas Visoespacial" ou "Executivo Central".
Por exemplo, o Estudo 1 avalia "Loop fonológico" (resultado A) com duas medidas diferentes (= tamanho do efeito 1 e 2) e "executivo central" (resultado B) com uma medida (= tamanho do efeito 3).
Problema:
Uma abordagem multivariada adequada requer conhecer todas as correlações entre tamanhos de efeito e resultados para estimar as covariâncias. No entanto, não conheço (1) as correlações entre diferentes tamanhos de efeito no mesmo estudo e (2) as correlações entre os resultados de diferentes estudos. Posso calculá-los ou tentar encontrar pelo menos algumas correlações para trabalhar, mas isso significaria muitas pesquisas adicionais na literatura que eu gostaria de evitar.
Solução (até agora):
deparei-me com alguns métodos que lidam com problemas semelhantes.
A estimativa de variância robusta (Hedges, 2010) é uma boa abordagem para lidar com vários tamanhos de efeito. No entanto, ainda tenho que adivinhar uma correlação e realizar uma análise de sensibilidade e também não parece possível comparar vários resultados (ou seja, apenas meta-regressão univariada).
A abordagem multinível de Van den Noorgate (2014) é promissora, uma vez que não requer estimativa de nenhuma correlação, permitindo variação entre os tamanhos dos efeitos e entre os tamanhos dos estudos. Uma meta-análise multivariada multinível (= resultados diferentes e tamanhos de efeitos múltiplos, como no esquema acima) e uma meta-regressão univariada multinível (= tamanhos de efeitos múltiplos, mas nenhuma diferenciação entre resultados) é descrita.
Usando o pacote metafor em R, estou pensando se posso combinar as duas abordagens multiníveis e executar uma meta-regressão multivariada multinível. Exemplos para uma meta-análise multinível e uma meta-regressão multivariada usando metafor são dados aqui http://www.metafor-project.org/doku.php/analyses:konstantopoulos2011 (multinível) e aqui http: //www.metafor- project.org/doku.php/analyses:berkey1998 (multivariado). (Observe que o exemplo multinível vinculado acima na verdade descreve uma abordagem para lidar com dependências hierárquicas (por exemplo, estudos conduzidos pelo mesmo laboratório de pesquisa). Em vez disso, uso a abordagem multinível descrita por Van den Noorgate.)
Variáveis:
ES: Tamanhos de efeito (g de Hedge)
VI: Variação dos tamanhos de efeito
Pub_Year: Ano de publicação como um preditor na meta-regressão
ES_ID: Todo tamanho de efeito tem um ID exclusivo, independentemente de qual estudo ou resultado eles pertencem.
ID do resultado : Os mesmos resultados têm o mesmo ID (por exemplo, "Loop fonológico" = 1, "Executivo central" = 2), independentemente de qual estudo eles pertencem.
ID do estudo: Os tamanhos de efeito do mesmo estudo têm o mesmo ID (por exemplo, tamanhos de efeito do estudo 1 = 1, tamanhos de efeito do estudo 2 = 2), independentemente de qual resultado eles pertencem.
Código R na metafor para a meta-análise multivariada multinível:
rma.mv (ES, VI, mods = ~ Outcome_ID -1, aleatório = lista (~ 1 | Study_ID, ~ 1 | ES_ID), data = data.set)
- mods = ~ Outcome_ID -1 exige uma abordagem multivariada e lista os tamanhos médios de efeito para cada resultado.
- random = list (~ 1 | Study_ID, ~ 1 | ES_ID) é a abordagem multinível descrita por Van den Noorgate. Permite variação aleatória entre os tamanhos dos efeitos nos estudos (~ 1 | Study_ID) e entre os tamanhos dos efeitos (~ 1 | ES_ID). Você também pode realizar essa análise usando o pacote metaSEM. Os resultados são idênticos.
Código R na metafora para a meta-regressão multivariada multinível:
rma.mv (ES, VI, mods = ~ ID_de resultado + Resultado: I (Média de anos_de_ Pub (média de anos de Pub)) -1, aleatória = lista (~ 1 | ID de estudo, ~ 1 | ES_ID), data = data.set)
- mods = ~ Resultado_ID + Resultado: I (Pub_Year-mean (Pub_Year)) -1 agora pede uma meta-regressão multivariada com o Publication Year centrado em torno da média como preditor.
Usando a opção profile () na metaforça, os gráficos de probabilidade do perfil parecem bons. No entanto, ainda estou me perguntando se não superparâmetro o modelo ou se há algo errado ao combinar os argumentos mods e aleatórios dessa maneira.
Aguardamos a sua opinião, sugestões, idéias, outras abordagens, tudo ;-) Obrigado!
Atualização, resposta à resposta de Wolfgang:
Primeiro de tudo: muito obrigado pela resposta detalhada e pelos links adicionais que você forneceu. Eu não sabia sobre a lista de discussão R-sig-mixed-models. Então obrigado! Eu aprecio muito.
Deixe-me tentar resumir tudo e adaptá-lo à minha situação para ver se entendo as coisas aqui. Eu posso fazer o seguinte:
Obtenção de correlações: Infelizmente, as correlações não são relatadas. Inicialmente, a meta-análise consistiu em mais de 50 estudos. Quase metade dos estudos tinha dados ausentes ou não relatados. Todos os autores desses estudos foram contatados e recebi 4 respostas de 26 solicitações (após 2 meses de espera). Mas esse é um problema geral de relatório que não deve ser discutido aqui.
Se eu fizer um palpite aproximado de todas as correlações, eu posso:
Realizar uma meta-análise multivariada e meta-regressão, como no Berkey et al. (1998) exemplo e faça uma análise de sensibilidade.
Use este modelo de metanálise multivariada e trabalhe com a função robust (). No entanto, nenhuma meta-regressão baseada na função robust () parece possível no metafor. E a função robust () descrita no blog de James Pustejovsky funciona apenas com meta-regressões univariadas. Portanto, se eu entendi direito, as estimativas da função robust () são mais ou menos para confirmar as estimativas do meu modelo já ajustado (?).
Vá diretamente para métodos robustos e use o pacote robumeta. No entanto, nenhuma meta-análise multivariada é possível. Encontrei um código SASpara lidar com esse problema. Mas o código foi desenvolvido há 3 anos e parece que nunca foi realmente discutido. No final, ao usar o robumeta, tenho que resumir muitos resultados diferentes em uma enorme meta-análise ou preciso realizar várias meta-análises univariadas para cada resultado que eu gostaria de evitar.Se não quiser adivinhar nenhuma correlação, posso seguir a abordagem multinível descrita por Van den Noorgate usando metafor, metaSEM ou SAS. No entanto, existem algumas restrições ao usar essa abordagem em comparação com uma abordagem multivariada com base em correlações. Também não tenho certeza se é possível uma meta-regressão multivariada multinível. O pacote metaSEM descreve apenas uma meta-análise multivariada multinível ou uma meta-regressão univariada multinível.
Infelizmente, não estou familiarizado com o uso de métodos de reamostragem na meta-análise. Estudei seus exemplos, mas não tenho certeza de como isso pode me ajudar a resolver o problema "correlação / multivariada". Você quer dizer que eu deveria tentar estimar as correlações usando o bootstrap? E, nesse caso, não tenho certeza de quais valores devem se correlacionar, pois o número de médias ou tamanhos de efeito dentro e entre os estudos diferem.
A simplificação do modelo descrito por Riley e colegas parece interessante. Eu tenho isso em mente, embora eu gostaria de trabalhar com um dos métodos descritos acima.
robust()
nos rma.mv
modelos que você escreveu na sua pergunta (com Study_ID
a variável de cluster). Como esses modelos assumem erros de amostragem independentes, eles são mal especificados. A abordagem robusta fornecerá resultados consistentes mesmo quando a parte da matriz de variância-covariância do modelo for especificada incorretamente.