## Not run:
# 1. Import a dataset of study by Fleiss (1993)
library(meta)
data("Fleiss1993bin")
data <- Fleiss1993bin
# 2. Calculate total sample size and standard error of each study
data$n <- data$n.asp + data$n.plac
data$se <- sqrt((1 / data$d.asp) - (1 / data$n.asp) + (1 / data$d.plac) - (1 / data$n.plac))
# 3. Test discordance in ranks between theoretical and observed study size.
output <- TestDiscordance(n = n, se = se, study = study, data = data)
# 4. Illustrate discordance plot
TestDiscordance(n = n, se = se, study = study, data = data, plot = TRUE)
## End(Not run)
Run the code above in your browser using DataLab