Nothing
Plot_data <-
expand.grid(
pi = c(0.5, 0.3, 0.1, 0.05),
n = c(10, 40, 80, 800),
p = ppoints(50)
) %>%
mutate(
label = paste("n=", n, "; pi=", pi, sep = ""),
group = paste("pi=", pi, "; n=", n, sep = ""),
y = qbinom(p, n, pi),
x = qnorm(p, n * pi, sqrt(n * pi * (1-pi)))
)
gf_qq( ~ y, data = Plot_data, color = ~factor(pi)) %>%
gf_qqline(color = "gray50") %>%
gf_facet_wrap( ~ label, scales = "free", ncol = 4) %>%
gf_labs(
y = expression(qbinom(p, n, pi)),
x = expression(qnorm(p, n * pi, sqrt(n * pi * (1-pi))))
) %>%
gf_theme(legend.position = "top") %>%
gf_refine(guides(color = guide_legend(title = "pi:")))
# p <- xyplot(
# y ~ x | paste("n=", n, sep = "") * paste("pi=", pi, sep = ""),
# data = Plot_data,
# scales = list(relation = "free"),
# cex = 0.6,
# ylab = expression(qbinom(p, n, pi)),
# xlab = expression(qnorm(p, n * pi, sqrt(n * pi * (1-pi)))),
# panel = function(x, y, ...){
# panel.abline(0, 1, ...)
# panel.xyplot(x, y, ...)
# })
# latticeExtra::useOuterStrips(p)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.