In [1]:
require(ggplot2)
require(reshape)
In [2]:
n <- 15
xs <- 0:n
ps <- seq(0.1, 0.9, 0.1)
In [3]:
df <- data.frame(x=xs, sapply(ps, dbinom, x=xs, size=n))
df <- melt(df, id.vars='x', variable_name='distr')
In [4]:
ggplot(df, aes(x, value)) + geom_point() + facet_grid(distr ~ .)
In [5]:
ns <- c(10, 20, 50, 200)
p <- 0.05
In [6]:
df <- data.frame(x=0:200, sapply(ns, dbinom, x=0:200, prob=p))
df <- melt(df, id.vars='x', variable_name='distr')
In [7]:
ggplot(df, aes(x, value)) + geom_point() + facet_grid(distr ~ .)
In [ ]: