## (standardized) mean difference to sensitivity/specificity
## equal variances
library(ggplot2)
delta <- seq(from = 2, to = 18, by = 0.05)
res <- sapply(delta, md2zfactor)
DF <- data.frame(SMD = delta, zfactor = res)
ggplot(DF, aes(x = SMD, y = zfactor)) +
geom_line() + xlab("(standardized) mean difference") +
ylab("z-factor") + ggtitle("SD1 = SD2 = 1") +
geom_hline(yintercept = 1, linetype = "dotted")
## unequal variances
delta <- seq(from = 2.5, to = 20, by = 0.05)
res <- sapply(delta, md2zfactor, sd1 = 1, sd2 = 2)
DF <- data.frame(MD = delta, zfactor = res)
ggplot(DF, aes(x = MD, y = zfactor)) +
geom_line() + xlab("mean difference") +
ylab("z-factor") + ggtitle("SD1 = 1, SD2 = 2") +
geom_hline(yintercept = 1, linetype = "dotted")
Run the code above in your browser using DataLab