res <- rma.mv(yi,vi,data=mydata,mods=~factor(A)+factor(B)+C,
>                           random = ~1|studyID/combination,
>                           method="REML",tdist=TRUE))

robust(res, cluster=mydata$studyID)