rm(list=ls())
set.seed(1)
library(couplingsmontecarlo)
graphsettings <- set_theme_chapter2()
## visualization of TV bounds
p <- function(x) dbeta(x, 5, 5)
q <- function(x) dbeta(x, 2, 1.5)
xseq <- seq(from=0, to=1, length.out=1e3)
pseq <- p(xseq)
qseq <- q(xseq)
g1 <- ggplot(data=data.frame(x=xseq, y=pseq), aes(x=x, y=y)) +
geom_polygon(fill=graphsettings$colors[1], alpha=0.75) +
geom_polygon(aes(y=qseq),
fill=graphsettings$colors[2], alpha=0.75) +
geom_polygon(aes(y=pmin(pseq,qseq)),
col='black', alpha=0.5, linetype=2) +
xlim(0,1)
g1
g2 <- ggplot(data=data.frame(x=xseq, y=pseq-qseq), aes(x=x, y=y)) +
geom_polygon(fill=rgb(0.9,0.9,0.9), alpha=1) +
stat_function(fun=function(x){ 1*(p(x)>q(x)) - 1*(p(x)<q(x)) }) +
xlim(0,1)
g2
ggsave(filename = "../tvdistance1.pdf", plot = gridExtra::grid.arrange(g1, g2, nrow = 2, ncol = 1))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.