#' Demonstration for the systematic sampling
#'
#' The whole sample frame is denoted by a matrix (\code{nrow * ncol}) in the
#' plane, and the sample points with equal intervals are drawn out according to
#' a random starting point. The points being sampled are marked by red circles.
#' @param nrow the desired number of rows of the sample frame.
#' @param ncol the desired number of columns of the sample frame.
#' @param size the sample size.
#' @param p.col,p.cex different colors / magnification rate to annotate the
#' population and the sample
#' @return None (invisible \code{NULL}).
#' @author Yihui Xie
#' @references Examples at \url{https://yihui.org/animation/example/sample-system/}
#' @seealso \code{\link{sample}}, \code{\link{sample.simple}},
#' \code{\link{sample.cluster}}, \code{\link{sample.ratio}},
#' \code{\link{sample.strat}}
#' @export
sample.system = function(
nrow = 10, ncol = 10, size = 15, p.col = c('blue', 'red'), p.cex = c(1, 3)
) {
n = nrow * ncol
if (size > n)
stop('sample size must be smaller than the population')
x = cbind(rep(1:ncol, nrow), gl(nrow, ncol))
nmax = ani.options('nmax')
for (i in 1:nmax) {
dev.hold()
plot(x, pch = 19, col = p.col[1], cex = p.cex[1], axes = FALSE,
ann = FALSE, xlab = '', ylab = '')
points(x[seq(sample(n, 1), by = n%/%size, length = size)%%n, ],
col = p.col[2], cex = p.cex[2])
box(lwd = 1)
ani.pause()
}
invisible(NULL)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.