Estou usando o pacote metafor em R. Eu ajustei um modelo de efeitos aleatórios com um preditor contínuo da seguinte maneira
SIZE=rma(yi=Ds,sei=SE,data=VPPOOLed,mods=~SIZE)
Qual produz a saída:
R^2 (amount of heterogeneity accounted for): 63.62%
Test of Moderators (coefficient(s) 2):
QM(df = 1) = 9.3255, p-val = 0.0023
Model Results:
se zval pval ci.lb ci.ub
intrcpt 0.3266 0.1030 3.1721 0.0015 0.1248 0.5285 **
SIZE 0.0481 0.0157 3.0538 0.0023 0.0172 0.0790 **
Abaixo, plotei a regressão. Os tamanhos dos efeitos são plotados proporcionalmente ao inverso do erro padrão. Percebo que essa é uma afirmação subjetiva, mas o valor de R2 (63% de variação explicada) parece muito maior do que o refletido pelo modesto relacionamento mostrado no gráfico (mesmo levando em consideração os pesos).
Para mostrar o que quero dizer, se eu fizer a mesma regressão com a função lm (especificando os pesos do estudo da mesma maneira):
lmod=lm(Ds~SIZE,weights=1/SE,data=VPPOOLed)
Então o R2 cai para 28% de variação explicada. Isso parece mais próximo do que as coisas são (ou pelo menos, minha impressão de que tipo de R2 deve corresponder ao enredo).
Sei que, depois de ler este artigo (incluindo a seção de meta-regressão): ( http://www.metafor-project.org/doku.php/tips:rma_vs_lm_and_lme ), essas diferenças na maneira como as funções lm e rma se aplicam pesos podem influenciar os coeficientes do modelo. No entanto, ainda não está claro para mim por que os valores de R2 são muito maiores no caso de meta-regressão. Por que um modelo que parece ter um ajuste modesto é responsável por mais da metade da heterogeneidade dos efeitos?
O valor R2 maior é porque a variação é particionada de maneira diferente no caso meta-analítico? (variabilidade da amostragem versus outras fontes) Especificamente, o R2 reflete a porcentagem de heterogeneidade representada na parte que não pode ser atribuída à variabilidade da amostragem ? Talvez haja uma diferença entre "variância" em uma regressão não meta-analítica e "heterogeneidade" em uma regressão meta-analítica que eu não estou apreciando.
Receio que declarações subjetivas como "Não parece certo" sejam tudo o que tenho para continuar aqui. Qualquer ajuda com a interpretação de R2 no caso de meta-regressão seria muito apreciada.