Оценщик Хайека всегда жалуется
С использованием sampling пакет, если я запускаю следующий пример:
library(sampling)
# Belgian municipalities data base
data(belgianmunicipalities)
# Computes the inclusion probabilities
pik=inclusionprobabilities(belgianmunicipalities$Tot04,200)
N=length(pik)
n=sum(pik)
# Defines the variable of interest
y=belgianmunicipalities$TaxableIncome
# Draws a Poisson sample of expected size 200
s=UPpoisson(pik)
# Computes the Hajek estimator of the population mean
Hajekestimator(y[s==1],pik[s==1],type="mean")
R печатает следующее предупреждение:
Warning messages:
1: In Hajekestimator(y[s == 1], pik[s == 1]) :
the type estimator is missing
2: In Hajekestimator(y[s == 1], pik[s == 1]) :
by default the mean estimator is computed
Однако, если я укажу тип как предложено:
Hajekestimator(y[s==1], pik[s==1], type="mean")
R все еще печатает предупреждение
Warning message:
In Hajekestimator(y[s == 1], pik[s == 1], type = "mean") :
by default the mean estimator is computed
Что здесь происходит? Это ошибка? Как я могу получить sampling::Hajekestimator перестать жаловаться?