Nothing
## ExpSims <-
## expand.grid(n = c(10, 20, 40, 80), rep = 1:2000) %>%
## group_by(n, rep) %>%
## mutate(
## pval = pval(t.test(rexp(n), mu = 1)),
## dist = paste0("Exp(1); n=", n))
##
## TSims <-
## expand.grid(n = c(10, 20, 40, 80), rep = 1:2000) %>%
## group_by(n, rep) %>%
## mutate(
## pval = pval(t.test(rt(n, df = 3), mu = 0)),
## dist = paste0("t(3); n=", n))
##
## gf_qq( ~ pval, data = bind_rows(ExpSims, TSims),
## distribution = qunif, geom = "line") %>%
## gf_abline(slope = 1, intercept = 0, color = "red",
## linetype = "dashed", alpha = 0.6) %>%
## gf_facet_wrap( ~ dist, nrow = 2)
## gf_qq( ~ pval, data = bind_rows(ExpSims, TSims), na.rm = TRUE,
## distribution = qunif, geom = "line") %>%
## gf_abline(slope = 1, intercept = 0, color = "red",
## linetype = "dashed", alpha = 0.6) %>%
## gf_lims(x = c(0, 0.2), y = c(0, 0.2)) %>%
## gf_facet_wrap( ~ dist, nrow = 2)
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.