pps (version 1.0)

sampfordpi: Joint inclusion probabilities for Sampford's PPS sampling method

Description

Compute joint inclusion probabilities for Sampford's method of PPS sampling

Usage

sampfordpi(sizes,n)

Arguments

sizes

A vector of the sizes of the units in the population

n

The sample size

Value

Returns a matrix with the inclusion probability pi(i) for each unit i in the population and with the joint inclusion probability pi(i,j) of units i and j in position (i,j) in the matrix, where i and j are not equal. Note that the size of the matrix is NxN, where N is the population size.

Examples

# NOT RUN {
sizes <- c(9,2,5,17,4,21,15,7,4,11,23,23,14)
piij <- sampfordpi(sizes,4)
weights <- 1/diag(piij)  # the weights one would use for estimation
# }