今天有朋友咨询我怎么写PPS抽样的代码,试着找了下,找到一个实现PPS抽样的R包。
百度百科:
PPS 抽样是指按概率比例抽样,属于概率抽样中的一种。是指在多阶段抽样中,尤其是二阶段抽样中,初级抽样单位被抽中的机率取决于其初级抽样单位的规模大小,初级抽样单位规模越大,被抽中的机会就越大,初级抽样单位规模越小,被抽中的机率就越小。就是将总体按一种准确的标准划分出容量不等的具有相同标志的单位在总体中不同比率分配的样本量进行的抽样。
R语言源码:
function (m, x)
{
N <- length(x)
pk <- x/sum(x)
cumpk <- cumsum(pk)
U <- runif(m)
ints <- cbind(c(0, cumpk[-N]), cumpk)
sam <- rep(