Skip to content
Prev 5002 / 5636 Next

[R-meta] PET-PEESE Question

Tori Peña 
Hello,

Thank you again for all your guidance.

I have been trying to conduct PET-PEESE analyses on several rma.mv objects
by adding variance as a moderator.  The following code works well:

#pet-peese test
m_1_pet <- rma.mv(yi = yi,
                     V = vi,
                     mods = ~ vi,
                     slab = STUDY.ID,
                     data = Relationship_TE_Occurrence,
                     random = ~ 1 | STUDY.ID/EFFECT_SIZE_ID,
                     test = "t",
                     method = "REML")

summary(m_1_pet, digits=3)

However, when I try to conduct the PEESE part I receive this error message:
"Model matrix not of full rank. Cannot fit model" for the following code:

Relationship$SE_Mod <- Relationship$vi^2

m_1_peese <- rma.mv(yi = yi,
                     V = vi,
                     mods = ~ SE_Mod,
                     slab = STUDY.ID,
                     data = Relationship_TE_Occurrence,
                     random = ~ 1 | STUDY.ID/EFFECT_SIZE_ID,
                     test = "t",
                     method = "REML")

summary(m_1_peese, digits=3)

Does anyone have any insight into why the model would fit for the PET part
but not the PEESE? I'm not sure how to fix this issue.

Thanks,
Tori