sunterpi2: Joint inclusion probabilities for Sunter sampling

In SunterSampling: Sunter's sampling design

Description

The function computes the joint (second-order) inclusion probabilities for Sunter sampling.

Usage

sunterpi2(x, n)

Arguments

x	the vector of values of the auxiliary variable (measure of size)
n	the sample size

Value

Returns a NxN matrix of the following form: the main diagonal contains the first-order inclusion probabilities for each unit k in the population; elements (k,l) are the joint inclusion probabilities of units k and l, with k not equal to l. N is the population size.

Author(s)

Alessandro Barbiero, Giancarlo Manzi

References

C-E. Sarndal, B. Swensson, J. Wretman: Model Assisted Survey Sampling, Springer, New York, 1992

See Also

sunter, sunterpi

Examples

# from Sarndal et al. (1992)
x<-c(40,25,20,10,5)
n<-2
sunterpi2(x,n)

Example output

          [,1]       [,2]       [,3]       [,4]       [,5]
[1,] 0.8000000 0.33333333 0.15555556 0.15555556 0.15555556
[2,] 0.3333333 0.50000000 0.05555556 0.05555556 0.05555556
[3,] 0.1555556 0.05555556 0.23333333 0.01111111 0.01111111
[4,] 0.1555556 0.05555556 0.01111111 0.23333333 0.01111111
[5,] 0.1555556 0.05555556 0.01111111 0.01111111 0.23333333

SunterSampling documentation built on May 1, 2019, 8:47 p.m.