CMMM 计算和数学方法在医学 1748 - 6718 1748 - 670 x Hindawi出版公司 10.1155 / 2015/172918 172918年 研究文章 敏感的比例估计连续抽样的随机响应数据 http://orcid.org/0000 - 0002 - 9849 - 1762 1、2 http://orcid.org/0000 - 0003 - 2984 - 6192 Zongda 3 Jiayong 4 http://orcid.org/0000 - 0002 - 3846 - 7672 通用电气 1 1 公共卫生学院,医学院 苏州大学 苏州215123 中国 suda.edu.cn 2 数学科学学院 德州大学 德州253023 中国 dzu.edu.cn 3 公共卫生部门 浙江医科大学 杭州310053 中国 zjmc.net.cn 4 危重病医学,Linshu县人民医院 山东临沂,276700 中国 2015年 18 5 2015年 2015年 31日 10 2014年 22 11 2014年 05年 12 2014年 18 5 2015年 2015年 版权©2015 Bo Yu et al。 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。

本文考虑了二项比例的估计问题的敏感或羞辱属性感兴趣的人群中。随机响应技术提出了保护被调查者的隐私和减少反应偏差而引起敏感属性的信息。在许多敏感问题调查,同样的人口通常是在每个场合反复采样。在本文中,我们应用连续抽样方案改善敏感比例的估计当前的场合。

1。介绍

社会调查有时包括侮辱或敏感问题的询问,如习惯性逃税、性行为、物质滥用和过度赌博,很难获得有效和值得信赖的信息。如果受访者被要求直接对有争议的问题,它往往导致拒绝或不真实的答案,特别是当他们有侮辱的行为。为了克服这个困难,华纳( 1]介绍了随机响应技术来估计的人口比例轴承侮辱或敏感等特性在一个特定的社区。这种技术可以让被申请人如实回答敏感问题没有暴露的令人尴尬的行为。在华纳的开创性工作 1),一些研究人员在这个领域作出了重要贡献,比如Christofides [ 2, 3辛格],[ 4金),和以拦 5)、黄( 6, 7],辛格和Sedory [ 8),常和郭 9),Arnab et al。 10]。所有这些结果都是基于样本有一次,在目前的研究并非如此。

在许多敏感问题调查,同样的人口经常反复每一次采样,以便发展随着时间的推移。在这种情况下,使用连续抽样方案可以有吸引力的替代提高水平的估计时间点或测量两个时间点之间的变化。在连续两次采样,先前的理论 11, 12)旨在提供最优估计量的意思是对当前(第二)。连续采样也被Narain[讨论一些细节 13],Raj [ 14辛格],[ 15),Ghangurde和饶 16],奥卡福[ 17),Arnab和奥卡福 18),Biradar和辛格 19辛格),g . n和v . k .辛格( 20.),Artes et al。 21),等等,辛格et al。 22]。但是没有努力已经取得了估计当前时刻敏感的无限的人口的比例。这动机导致作者考虑的问题估计敏感或羞辱属性的二项比例的人口连续两次采样的兴趣。此外,集群抽样时通常是首选的目标人群是地理多样化。在本文中,我们利用旋转集群样品设计构建一个类估计的随机响应的情况下调查。剩下的纸是组织如下。节 2,我们提出了一种新的科学调查方法与集群使用西蒙斯模型旋转抽样。节 3,相应的公式方法调查发现之后,上述提到的方法和相应的公式成功地设计和应用于婚前性行为的调查在学生中东吴大学的部分 4。部分 5包含的结论。

2。该调查方法 2.1。西蒙斯模型

西蒙斯模型基于华纳霍维茨提出的随机响应技术et al。 23]。基本思想是建立一个随机关系之间的个人和两个不相关的问题。西蒙斯设计由两个不相关的问题,A和B,在概率基础上回答,一个是“你拥有敏感特征”和B是一个非敏感的问题,如“你的生日号码是奇怪的。“两个问题提出了A和B与预设受访者概率 P 1 - - - - - - P ,分别。简单随机抽样的替换(SRSWR)。选定的被调查者被要求选择一个问题a或B和报告“是”如果他/她的实际地位匹配与选择的问题,否则“不”。

2.2。西蒙斯模型集群旋转抽样

在以下采样两次被认为是估计人口比例与敏感特征二次旋转时抽样单位是集群。西蒙斯模型下的抽样步骤部分集群旋转如下。

首先,人口分为初级抽样单位(或集群)和集群内的单位二级抽样单位(人)。

其次,在第一次的随机样本 n 集群替换是从人口。绘制集群内的人们被要求选择一个问题a或B和报告“是”如果他/她的实际地位匹配与选择的问题,“不”否则,使用西蒙斯模型。

第三,在第二次 n 集群随机选择的第一次被保留,其余 u = n - - - - - - 集群取而代之的是一个全新的选择。所有的人总在第二次集群内调查使用西蒙斯模型。

3所示。公式推导 3.1。估计量的人口比例在第二次及其方差

考虑的随机样本 n 集群替代来自人口由 N 集群和 th群 单位 ( = 1、2 , , N )

在第二个(当前)场合 n 集群随机选择的第一次被保留,其余 u = n - - - - - - 集群取而代之的是一个全新的选择。让 一个 1 , l 的数量 l 集群(包括保留 l 单位)与研究第一次下的敏感特性 ( l = 1、2 , , ) ,让 一个 1 , + k 的数量 k 集群(包括旋转 + k 单位)的研究第一次下敏感特性,分别 ( k = 1、2 , , u ) 一个 2 , l 的数量吗 l 集群(包括保留 l 单位)与敏感特征研究第二(当前)的场合 ( l = 1、2 , , ) 一个 2 , + j 的数量吗 j 新的集群(包括 + j 单位)与敏感特征研究第二(当前)的场合 ( j = 1、2 , , u ) 。同样,让 π 1 , l 的比例 l th留存的集群研究第一次下的敏感特性 ( l = 1、2 , , ) ,让 π 1 , + k 的比例 k th旋转集群与研究第一次下的敏感特性 ( k = 1、2 , , u ) ,分别。 π 2 , l 的比例吗 l th保留第二集群与敏感特征研究(当前)场合 ( l = 1、2 , , ) π 2 , + j 的比例吗 j th新鲜集群下的敏感特性研究第二(当前)场合 ( j = 1、2 , , u ) 。假设方差和第一次和第二次之间的相关系数是常数和整体校正系数将被忽略。

定义如下:

π 1 :敏感特征的人口比例第一次;

π 2 :敏感特征的人口比例在第二次;

π 1 的比例: 保留第一次集群敏感特征;

π 2 的比例: 保留与敏感特点第二次集群;

π 1 u 的比例: u 旋转集群第一次敏感特征;

π 2 u 的比例: u 新的集群与第二次敏感的特征。

以下是根据公式和结果由科克伦( 24]。

的估计量 π 1 (1) π ^ 1 = = 1 n 一个 1 , = 1 n

的估计量 π 1 (2) π ^ 1 = = 1 一个 1 , = 1

的估计量 π 1 u (3) π ^ 1 u = j = 1 u 一个 1 , + j j = 1 u + j

的估计量 π 2 (4) π ^ 2 = j = 1 一个 2 , j j = 1 j

的估计量 π 2 u (5) π ^ 2 u = j = 1 u 一个 2 , + j j = 1 u + j 考虑一个广义估计量 π ^ 2 的人口比例的敏感特征第二次或电流场合 (6) π ^ 2 = 一个 π ^ 1 u + b π ^ 1 + c π ^ 2 u + d π ^ 2 , 在哪里 一个 , b , c , d 是合适的常数。

我们有 (7) E π ^ 2 = 一个 + b π 1 + c + d π 2 因为估计量 π ^ 2 π 2 是一个无偏估计量的 π 2 ,我们有 一个 + b = 0 c + d = 0

因此,估计量 6的形式 (8) π ^ 2 = 一个 π ^ 1 u - - - - - - π ^ 1 + c π ^ 2 u + 1 - - - - - - c π ^ 2 估计量的方差 π ^ 2 (9) Var π ^ 2 = 一个 2 Var ( π ^ 1 u ) + 一个 2 Var ( π ^ 1 ) + c 2 Var ( π ^ 2 u ) + 1 - - - - - - c 2 Var ( π ^ 2 ) - - - - - - 2 一个 ( 1 - - - - - - c ) ( π ^ 1 , π ^ 2 ) = 一个 2 年代 1 2 u + 一个 2 年代 1 2 + c 2 年代 2 2 u + 1 - - - - - - c 2 年代 2 2 - - - - - - 2 一个 1 - - - - - - c ρ 年代 1 年代 2 = 一个 2 年代 1 2 u + 年代 1 2 + c 2 年代 2 2 u + 1 - - - - - - c 2 年代 2 2 - - - - - - 2 一个 1 - - - - - - c ρ b 年代 1 年代 2 其他协方差为零。

估计量的方差最小化 π ^ 2 关于 一个 c N 是足够大, (10) Var ( π ^ 2 ) 一个 = 2 一个 1 u + 1 年代 1 2 - - - - - - 2 ( 1 - - - - - - c ) ρ b 年代 1 年代 2 = 0 , Var ( π ^ 2 ) c = 2 c 年代 2 2 u - - - - - - 2 1 - - - - - - c 年代 2 2 + 2 一个 ρ b 年代 1 年代 2 = 0 然后我们得到 (11) 一个 1 u + 1 年代 1 2 = 1 - - - - - - c ρ b 年代 1 年代 2 , 一个 ρ b 年代 1 年代 2 = c 年代 2 2 u - - - - - - 1 - - - - - - c 年代 2 2 我们得到 (12) 1 / u + 1 / ρ b / = ( 1 - - - - - - c ) ρ b / c / u - - - - - - 1 - - - - - - c / 一个人 (13) u + / u ρ b / = ( 1 - - - - - - c ) ρ b / c ( u + ) - - - - - - u / u , u + = n

我们有 (14) n / u ρ b = ( 1 - - - - - - c ) ρ b c n - - - - - - u / u

因此, (15) c = u / n - - - - - - u 2 / n 2 ρ b 2 1 - - - - - - u 2 / n 2 ρ b 定义 μ = u / n

我们得到了 (16) c = μ - - - - - - μ 2 ρ b 2 1 - - - - - - μ 2 ρ b 2 通过 16,我们得到 (17) 一个 = ( 1 - - - - - - c ) ( ρ b 年代 2 ) ( / u + 1 ) 年代 1 , u + = n μ = u / n

一个人 (18) 一个 = ( 1 - - - - - - c ) μ ( ρ b 年代 2 ) 年代 1 通过 16 18,我们得到 (19) 一个 = 年代 2 年代 1 · μ 1 - - - - - - μ ρ b 1 - - - - - - μ 2 ρ b 2 , 在哪里 (20) 年代 h 2 = 1 N - - - - - - 1 = 1 N π h , - - - - - - π h 2 , h = 1、2 ; μ = u n , ρ b = = 1 N ( π 2 , - - - - - - π 2 ) ( π 1 , - - - - - - π 1 ) = 1 N π 2 , - - - - - - π 2 2 = 1 N π 1 , - - - - - - π 1 2

定理1。

西蒙斯模型部分集群旋转,下一个 (21) π ^ 2 = 年代 2 年代 1 · μ 1 - - - - - - μ ρ b 1 - - - - - - μ 2 ρ b 2 ( π ^ 1 u - - - - - - π ^ 1 ) + μ - - - - - - μ 2 ρ b 2 1 - - - - - - μ 2 ρ b 2 π ^ 2 u + ( 1 - - - - - - μ 1 - - - - - - μ 2 ρ b 2 ) π ^ 2 , 和估计量的方差 π ^ 2 (22) Var 最小值 P ^ 2 = 1 - - - - - - μ ρ b 2 1 - - - - - - μ 2 ρ b 2 年代 2 2 n

备注2。

在实践中, ρ b 年代 h 2 是未知的。的估计量 ρ b (23) ρ ^ b = k = 1 ( π 2 , k - - - - - - π ^ 2 ) ( π 1 , k - - - - - - π ^ 1 ) k = 1 π 2 , k - - - - - - π ^ 2 2 k = 1 π 1 , k - - - - - - π ^ 1 2 的估计量 年代 h 2 (24) 年代 h 2 = 1 - - - - - - 1 = 1 π h , - - - - - - π ^ h 2 , h = 1、2

定理3。

西蒙斯模型部分集群旋转,下一个最优转速 (25) μ = 1 + 1 - - - - - - ρ b 2 - - - - - - 1 和最优估计量的方差 π ^ 2 (26) Var 最小值 选择 π ^ 2 = 1 + 1 - - - - - - ρ b 2 年代 2 2 2 n 实际上,抽样调查的成本通常代表下面的简单函数,根据科克伦( 24]: (27) C T = c 0 + c 1 + c 2 u , 在哪里 C T 是抽样的总成本, c 0 调查的基本成本, c 1 是调查一个保留集群的平均基本成本第二次,然后呢 c 2 是调查一个新鲜的平均基本成本集群第二次。

定理4。

在给定的抽样调查的成本 C T ,一个 (28) Var 选择 π ^ 2 = c 1 1 - - - - - - ρ b 2 + c 2 年代 2 2 2 C T - - - - - - c 0 和样本大小的估计部分集群旋转 (29) n = C T - - - - - - c 0 1 + 1 - - - - - - ρ b 2 c 1 1 - - - - - - ρ b 2 + c 2 , 在哪里 (30) 年代 2 2 = 1 N - - - - - - 1 = 1 N π 2 , - - - - - - π 2 2

3.2。估计量的< inline-formula > < mml:数学xmlns: mml = " http://www.w3.org/1998/Math/MathML " id = " M118 " > < mml: mrow > < mml: msub > < mml: mrow > < mml: mi >π< / mml: mi > < / mml: mrow > < mml: mrow > < mml: mi > h < / mml: mi > < mml:莫>,< / mml:莫> < mml: mi >我< / mml: mi > < / mml: mrow > < / mml: msub > < / mml: mrow > < / mml:数学> < / inline-formula >

R h , 被选择的比例 th集群(包括 单位)与非敏感特性研究 h th场合; λ h , h , 表示的数量和比例的“是”的答案 th集群,分别 λ ^ h , = h , / , h = 1、2 , ( = 1、2 , , n )

从总概率公式(参见[ 25]),我们可以得到的 (31) π h , = λ h , - - - - - - ( 1 - - - - - - P ) R h , P , h = 1、2 , = 1、2 , , n 因此 (32) 一个 h , = π h , , h = 1、2 , = 1、2 , , n

4所示。应用程序 4.1。调查设计

这项调查是关于婚前性行为在苏州大学独墅湖校区的学生。我们把每个类的集群平均45人/班。在第一次(2011)、12类随机从所有类。所有的人选择调查的12类西蒙斯模型敏感问题。第二次(2013),8的12类随机选择的第一次被保留,其余4类替换为一个新的选择。然后选中的类中的所有人包括8保留类和4新鲜类调查的西蒙斯模型敏感问题。

在我们的设计中,每个人被要求画一个球放回随机从包含6个红球和一袋4白色球与已知的概率(红球的比例是0.6)。如果一个红球是由被申请人选择,那么他或她就会被问到那个敏感的问题,一个是“你是一个组的成员有婚前性行为。“如果一个白色的球被选中,他或她会回答的非敏感问题B, B是“你的学生号码是奇怪的。“被申请人报告“是”如果他/她的实际地位匹配与选择的问题,“不”。

两次的所有问卷都进行了检查,确保它们是独立完成的,没有遗漏的问题。调查的回收率为100%,没有失败的问卷。所有的数据处理和分析通过Excel 2003和SAS 9.13。

4.2。结果 4.2.1。准备调查的结果

在我们的设计中,每个人被要求画一个球放回随机从包含6个红球和一袋4白色球与已知的概率(红球的比例是0.6)。如果一个红球是由被申请人选择,那么他或她就会被问到那个敏感的问题,一个是“你是一个组的成员有婚前性行为。“如果一个白色的球被选中,他或她会回答的非敏感问题B, B是“你的学生号码是奇怪的。“被申请人报告” 是的“如果他/她的实际状态匹配与选择的问题,“ 没有“否则。根据 31日,我们得到样本比例的大学生婚前性行为 π h , , h = 1、2 ,如表所示 1

每个类的比例与婚前性行为在第一和第二次。

类数 π 1 , π 2 ,
1 0.2624 0.2348
2 0.1631 0.1945
3 0.2101 0.2264
4 0.2063
5 0.1556 0.1986
6 0.2390
7 0.1783
8 0.1970 0.1550
9 0.0123 0.0114
10 0.0476 0.0738
11 0.0455
12 0.1185 0.1187
13 0.2035
14 0.1587
15 0.1926
16 0.1583
4.2.2。估计量的人口比例在第二次及其方差

通过 1的估计人口比例与婚前性行为在第一次如下: π ^ 1 = 0.142933

根据 24, 2, 3,我们有 (33) 年代 1 2 = 0.005853 , π ^ 1 = 0.1458 , π ^ 1 u = 0.1372 , 分别。

根据调查的结果,婚前性行为在苏州大学独墅湖校区学生在第二次,从公式 4 5, (34) π ^ 2 = 0.1562 , π ^ 2 u = 0.1783 通过 23 24,我们获得 年代 2 2 = 0.00426 2 ρ ^ b = 0.936 ,分别。

从公式 25, μ = 0.7985 ;然后根据公式 21,我们得到 (35) π ^ 2 = 0.2912 ( 0.1372 - - - - - - 0.1458 ) + 0.5453 × 0.1783 + ( 1 - - - - - - 0.5453 ) × 0.1562 = 0.1657 使用 22,我们得到 弗吉尼亚州 r 最小值 ( 选择 ) ( π ^ 2 ) = 0.00024 。因此,标准偏差如下: (36) 弗吉尼亚州 r 最小值 ( 选择 ) ( π ^ 2 ) = 0.01549 所以,95%置信区间的人口比例与婚前性 ( 0.1353,0.1961 ]

5。讨论和结论

总之,在这项研究中,我们提出了一种新的抽样方法解决问题的敏感问题调查重复随着时间的推移,这是第一次尝试由作者在这个方向。那么相应的公式估计的人口比例与敏感特性及其方差提出的抽样方法。此外,公式给出的最优转速和样本容量成本的抽样调查。

上述方法和相应的公式成功地设计和应用于婚前性行为调查在苏州大学独墅湖校区。总之,设计采样方法和相应的公式具有重要的理论和应用价值实现连续测量的敏感问题。

6。定理的证明 定理的证明< xref ref-type =“声明”掉= " thm1 " > < / xref > 1。

使用的最优值 一个 c 给出的 16 19,估计量 π ^ 2 减少到 21

通过 9, 16, 19,我们有 (37) 弗吉尼亚州 r 最小值 P ^ 2 = 1 - - - - - - μ ρ b 2 1 - - - - - - μ 2 ρ b 2 年代 2 2 n

定理的证明< xref ref-type =“声明”掉= " thm2 " > 3 < / xref >。

的最优值 μ 给出了进一步最小化 22关于 μ , (38) 弗吉尼亚州 r 最小值 ( π ^ 2 ) μ = 0 所以 (39) μ = 1 + 1 - - - - - - ρ b 2 - - - - - - 1 替换 39 22我们有最优估计量的方差 π ^ 2 作为 (40) 弗吉尼亚州 r 最小值 选择 π ^ 2 = 1 + 1 - - - - - - ρ b 2 年代 2 2 2 n

定理的证明< xref ref-type =“声明”掉= " thm3 " > < / xref >。

ByTheorem 3, (41) μ = 1 + 1 - - - - - - ρ b 2 - - - - - - 1 , μ = u n , u = n - - - - - - 替换 41 27,我们获得 (42) n = C T - - - - - - c 0 1 + 1 - - - - - - ρ b 2 c 1 1 - - - - - - ρ b 2 + c 2 假设平均集群组成 - - - - - - 单位;然后 (43) n = - - - - - - C T - - - - - - c 0 1 + 1 - - - - - - ρ b 2 c 1 1 - - - - - - ρ b 2 + c 2 替换 42 26,我们有 (44) 弗吉尼亚州 r 选择 π ^ 2 = c 1 1 - - - - - - ρ b 2 + c 2 年代 2 2 2 C T - - - - - - c 0

利益冲突

作者宣称没有利益冲突有关的出版。

确认

作者想表达深深的感谢相关裁判,仔细阅读论文和评论大大改善。本文由国家自然科学基金的资助支持中国。高g . 81273188)。作者感谢高g(通讯作者)他宝贵的帮助。

华纳 s . L。 随机反应:一项调查技术消除逃避回答偏见 美国统计协会杂志》上 1965年 60 309年 63年 69年 10.1080 / 01621459.1965.10480775 2 - s2.0 - 0013776710 Christofides t . C。 广义随机响应技术 Metrika 2003年 57 2 195年 200年 10.1007 / s001840200216 MR1969252 2 - s2.0 - 10044243431 Christofides t . C。 在分层抽样的随机响应 杂志的统计规划和推理 2005年 128年 1 303年 310年 10.1016 / j.jspi.2003.11.001 MR2110192 2 - s2.0 - 15344343319 辛格 g . N。 使用链式的比率差异估计量在连续采样 应用数学和统计的国际期刊 2006年 6 41 49 人类。 以拦 m E。 一个分层随机响应模型无关的问题 统计文件 2007年 48 2 215年 233年 10.1007 / s00362 - 006 - 0327 - 6 MR2295812 2 - s2.0 - 33847395801 K.-C。 使用可选的随机响应技术估计为敏感特征 质量和数量 2008年 42 5 679年 686年 10.1007 / s11135 - 007 - 9082 - 6 2 - s2.0 - 50249155326 K.-C。 无偏估计量的意思是,方差和敏感水平有限总体抽样量化特征 Metrika 2010年 71年 3 341年 352年 2 - s2.0 - 77952097928 10.1007 / s00184 - 009 - 0234 - 7 MR2610002 辛格 年代。 Sedory 美国一个。 模拟一个真正的平等保护的三个估计研究的受访者在随机响应抽样 Statistica Neerlandica 2012年 66年 4 442年 451年 10.1111 / j.1467-9574.2012.00524.x MR2983304 2 - s2.0 - 84867403774 周宏儒。 M.-P。 人口比例随机响应采样使用加权估计置信区间施工 Metrika 2012年 75年 5 655年 672年 10.1007 / s00184 - 011 - 0346 - 8 MR2946653 2 - s2.0 - 84862592642 Arnab R。 辛格 年代。 D。 使用两个扑克牌在随机响应技术复杂的调查设计 通信数据。理论和方法 2012年 41 16 - 17 3198年 3210年 10.1080 / 03610926.2012.682634 MR3003862 2 - s2.0 - 84866050916 安杰森 r . J。 统计调查的调查获取农场的事实 爱荷华州农业实验台研究公告 1942年 304年 1 104年 耶茨 F。 人口普查和调查抽样方法 1949年 英国伦敦 查尔斯·格里芬 Narain r D。 在连续的场合抽样的递推公式 农业统计的印度社会杂志》上 1953年 5 96年 99年 MR0062385 拉吉 D。 采样两次与概率比例的大小 《数理统计 1965年 36 327年 330年 10.1214 /中耳炎/ 1177700297 MR0171365 辛格 D。 估计在连续使用多级抽样设计 美国统计协会杂志》上 1968年 63年 99年 112年 MR0228142 Ghangurde p D。 j . N。 一些结果在两次取样 数论派 1969年 31日 463年 472年 MR0270531 奥卡福 f . C。 抽样的理论和应用两次估计当前的人口比例 Statistica 1992年 42 137年 147年 MR1212969 Arnab R。 奥卡福 f . C。 双采样两次 巴基斯坦《统计 1992年 8 3 9 18 MR1225063 Biradar r S。 辛格 h·P。 连续采样的场合上使用辅助信息 加尔各答统计协会公告 2001年 51 203 - 204 243年 251年 MR1900665 辛格 g . N。 辛格 诉K。 使用连续采样的辅助信息 农业统计的印度社会杂志》上 2001年 54 1 1 12 MR2059533 Artes R。 伊娃 M。 加西亚 l 阿米莉亚 V。 估计目前的人口比例在连续采样 农业统计的印度社会杂志》上 2001年 54 3 342年 354年 MR1914137 辛格 h·P。 裁缝 R。 辛格 年代。 人类。 在连续抽样总体方差估计 质量和数量 2011年 45 3 477年 494年 10.1007 / s11135 - 009 - 9309 - 9 2 - s2.0 - 79952191687 霍维茨 d·G。 沙阿 b . V。 席梦思床品公司 w·R。 不相关的随机响应模型的问题 《社会统计学部分:美国统计协会 1967年 326年 65年 72年 科克伦 w·G。 抽样技术 1977年 3日 纽约,纽约,美国 约翰威利& Sons MR0474575 Z。 抽样技术及其应用 2005年 1日 中国,北京 清华大学出版社