[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

-- 
*Tori Peña, Ph.D. *(she/her/ella)
Cognitive Psychology
Dept. of Psychology
Stony Brook University
Stony Brook, NY 11790-2500
[image: Stony Brook University logo]

	[[alternative HTML version deleted]]