R/PikPPS.r

In TeachingSampling: Selection of Samples and Parameter Estimation in Finite Population

Defines functions PikPPS

Documented in PikPPS

#' @export

PikPPS<-function(n,x){
pik<- n*x/sum(x)
while((sum(pik>1))!=0){
s<-which(pik>=1)
new=(1:length(pik))[-s]
pik[s]=1
txnew<-sum(x[s])
for(k in new){
pik[k]<- (n-length(s))*x[k]/(sum(x)-txnew)
}
}
pik
}

Try the TeachingSampling package in your browser