如何解决Python中使用“samplics”时单个PSU的问题？

Question

我正在使用“samples”包来进行复杂的调查。一开始我使用以下顺序：

import numpy as np
import pandas as pd
import pyreadr
import samplics
from samplics.estimation import TaylorEstimator
from samplics.estimation import ReplicateEstimator
from samplics.utils import SinglePSUEst

我特别有兴趣获得

`total_1 =    TaylorEstimator('mean').estimate(y=tabla['TOT_1'],samp_weight=tabla['FAC_EXP'],stratum=tabla['EST_DIS'], psu=tabla['UPM_DIS'])`, 

但是输出是：

单个 PSU 的问题

我理解这个问题，但我无法找到在命令中包含 single_psu =“singleton”或“skip”选项的路径。

我希望找到一种方法来包含解决单个初级抽样单位问题的选项。你有什么主意吗？ ：）

非常感谢！

Answer 1

您可以在 Samplics 的 GitHub 页面上找到此问题的答案： https://github.com/samplics-org/samplics/issues/33

简而言之，您可以在 TaylorEstimator('mean').estimate 的调用中添加参数“single_psu”，例如

`total_1 =    TaylorEstimator('mean').estimate(y=tabla['TOT_1'],samp_weight=tabla['FAC_EXP'],stratum=tabla['EST_DIS'], psu=tabla['UPM_DIS'], sigle_psu= SinglePSUEst.skip, )

希望有帮助。