S <- simple_rs(N = 100)
table(S)
S <- simple_rs(N = 100, prob = 0.3)
table(S)
# A probability of inclusion that varies unit by unit
S <- simple_rs(N = 100, prob_unit = seq(0.1, 0.9, length.out = 100))
table(S)
Run the code above in your browser using DataLab