# Statistics Lab mathematics M3; hypothetical models only
n <- 10; p <- 0.30
c(expected=n*p, variance=n*p*(1-p))
dbinom(2, size=n, prob=p)
set.seed(2026)
sim <- rbinom(10000, size=n, prob=p)
c(mean=mean(sim), variance=var(sim))