input_data <- MetaHDInput(realdata)
Y <- input_data$Y
Slist <- input_data$Slist
## MULTIVARIATE RANDOM-EFFECTS META-ANALYSIS, ESTIMATED WITH REML
model <- MetaHD(Y, Slist, method = "reml", bscov = "unstructured")
model$estimate
model$pVal
Run the code above in your browser using DataLab