# Compare binomial S4 with normal n <- 4 x <- 0:n plot(x, xlim=c(-1,5), ylim=c(0,0.5), type="n", xlab="x and n", ylab="dens(X=x) and P(S4=n)", main="Comparison of S4 with N(2,1)") points(x, dbinom(x, n, 0.5), pch=16) x <- (-1000:5000)/1000 points(x, dnorm(x, mean=2), type="l")