版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、第7部分假设检验 第一节 二项分布 二项分布是从著名的贝努里试验中推导而来。所谓贝努里试验,是指只有两种可能结果的随机试验。在实际问题中,有许多随机现象只包含两个结果,如男与女,是与非,生与死,同意与不同意,赞成与反对等等。通常,我们把其中比较关注那个结果称为“成功”,另一个结果则称为“失败”。每当情况如同贝努里试验,是在相同的条件下重复n次,考虑的是“成功”的概率,且各次试验相互独立,就可利用与二项分布有关的统计检验。虽然许多分布较之二项分布更实用,但二项分布简单明了,况且其他概率分布的使用和计算逻辑与之相同。所以要理解统计检验以及它所涉及的许多新概念,人们几乎都乐意从二项分布的讨论入手。
2、7/20/202221. 二项分布的数学形式 从掷硬币的试验入手。假定二项试验由重复抛掷n次硬币组成,已知硬币面朝上(成功)的概率是p,面朝下(失败)的概率是q (显然有 q1p)。这样,对试验结果而言,成功的次数(即硬币面朝上的次数)X是一个离散型随机变量,它的可能取值是0,1,2,3,n。而对X的一个具体取值x而言,根据乘法规则,我们立刻可以就试验结果计算出一种特定排列方式(先x次面朝上,而后nx次面朝下)实现的概率,即 ppppqqqqpxqn-x7/20/20223 由于正确解决概率问题,光考虑乘法规则是不够的,还要考虑加法规则,于是就x次成功和(nx)次失败这个宏观结果而言所包含的所
3、有排列的方式数,用符号表示 这样,我们就得到了二项试验中随机变量X的概率分布,即 7/20/20224 譬如,二项试验是将一枚硬币重复做8次抛掷,假设这枚硬币是无偏的,即pq0.5,那么恰好得到5次面朝上的概率是 硬币面朝上数x 概率P(X=x) 012345678 1/256= .004 8/256= .031 28/256= .109 56/256= .219 70/256= .274 56/256= .219 28/256= .109 8/256= .031 1/256= .004合 计 1.000 同理,我们也可以求出这个二项试验中硬币刚好为0,1,2,8次面朝上的各种宏观结果的概率,
4、全部写出来就是右表。 7/20/202252. 二项分布讨论X 0 1 2 n合计P(X) 二项分布为离散型随机变量的分布。每当试验做的是在相同的条件下n次重复的贝努利试验时,随机变量X共有n+1个取值。二项分布可以用分布律(见上表)和折线图(见右图)来表示。 当P=0.5时二项分布的图形是对称的。7/20/20226 E(X)=np, D(X)= 2= npq 二项分布受 p 和 n 变化的影响,只要确定了 p和 n,成功次数 X 的分布也随之确定。因此,二项分布还可简写作 B(x;n,p)。 二项分布的概率值除了根据公式直接进行计算外,还可查表求得。二项分布表的编制方法有两种:一种依据概率
5、分布律 P(x) 编制(见附表2);另一种依据分布函数 F(x) 编制(见附表3)。 其中7/20/20227 例 某特定社区人口的10%是少数民族,现随机抽取6人,问其中恰好2人是少数民族的概率是多少? 解 解法一:根据(7.3)式直接计算 解法二:根据附表2中纵列n6和横行p0.1所对应x值,可直接查得B(x;6,0.1)的概率值 B (2;6,0.1)00984 解法三:根据附表3求得 B (2;6,0.1)F(2) F(3 ) 0.11430.01590.0984 7/20/20228第二节 统计检验的基本步骤 二项分布是用数学或演绎推理的方法求得的一种理论分布。认识到概率分布是先验的
6、理论分布这一点很重要,因为我们不禁要问,既然试验或抽样调查的结果仅与随机变量可能取值中的一个相联系,那么实际试验或样本调查对结果的概率分布及前提假设有没有一个检验的问题?具体来讲,对于一枚硬币被重复抛掷8次的二项试验,经验告诉我们,一共有9种可能的结果,而且实现这些结果的机会是大不相同的。研究者实际上从来不用经验的方法求得概率分布,因为通常我们只对一项试验进行一次或几次,抽取样本也是一个或至多不过几个。既然二项分布是按照数学规则得到的,那么对这9种结果的可能性我们应该作出何种评价呢?如果实际试验(或抽样)得到的结果偏巧就是先验概率预示的最不可能出现的结果,那么我们是认定纯属巧合,还是开始对用数
7、学或演绎推理方法求得的概率以及理想试验的种种前提假设产生怀疑?更准确地说,在一枚硬币被重复抛掷8次的这个二项试验中,究竟出现什么结果时,我们应该对二项分布及其前提假设产生怀疑呢?是不是只要不是得到4次成功4次失败这个最大可能性结果就开始怀疑,还是仅当出现8次成功或一次也不成功这两个极端情况时才产生怀疑呢?这就是统计检验的核心问题。 7/20/20229大数定理表明:就大量观察而言,事件的发生具有一定的规律性。根据概率的大小,人们处理的态度和方式很不一样。在日常生活中,人们往往习惯于把概率很小的事件,当作一次观察中是极不可能看到的事件。例如,人们出门做事就有可能遇到不测事故,但却很少人因此而不敢
8、出门。原因是:小概率事件极不可能发生。 统计检验是指先建立一个关于总体情况的假设,继而抽取一个随机样本,然后以样本的统计量或者统计性质来检定假设。 统计检验的依据是小概率原理:一是认为小概率事件在一次观察中是极少出现的;二是如果在一次观察中出现了小概率事件,那么应该否定原有事件具有小概率的说法或者假设。7/20/202210(1)建立假设(2)求抽样分布(4)计算检验统计量(3)选择显著性水平和否定域(5)判定所所包有含统的计步检骤验 根据以往多年的统计表明,上海财大英语的平均成绩为90分,随机抽取100个学生,其平均成绩为80分,问今年财大学生的英语成绩是否下降?7/20/202211 1建
9、立假设 统计检验是将抽样结果和抽样分布相对照而作出判断的工作。取得抽样结果,依据描述性统计的方法就足够了。抽样分布则不然,它无法从资料中得到,非利用概率论不可。而不对待概括的总体和使用的抽样程序做某种必要的假设,这项工作将无法进行。比如通过掷硬币的实验得到二项分布,必须假设:样本是随机的,试验中各次抛掷相互独立;硬币是无偏的(或称是诚实的),即pq05。概括地说,必须首先就研究总体和抽样方案都做出假设,再加上概率论,我们就可以对各种可能结果做具体的概率陈述了。 7/20/202212 2求抽样分布 在做了必要的假设之后,我们就能用数学推理过程来求抽样分布了。比如在这一章开头,在硬币重复抛掷n次
10、的理想实验中,我们计算了成功次数为x的宏观结果所具有的概率,得到二项分布。如果前提假设变动了,还可以求出其他形式的概率分布,如正态分布、泊松分布、卡方分布等等,它们都有特定的方程式。由于数学上已经取得的成果,实际上统计工作者要做的这项工作往往并不是真的去求抽样分布的数学形式,而是根据具体需要,确定特定问题的统计检验应该采用哪种分布的现成的数学用表。 7/20/202213 3选择显著性水平和否定域 在统计检验中,那些不大可能的结果称为否定域。如果这类结果真的发生了,我们将否定假设;反之就不否定假设。 在统计检验中,通常把被检验的那个假设称为零假设(用符号H0表示),并用它和其他备择假设(用符号
11、H1表示)相对比。零假设与备择假设否定域7/20/202214 在统计检验中,无论是拒绝或者接受原假设,都不可能做到百分之百的正确,都有一定的错误。第一类错误是,零假设H0实际上是正确的,却被否定了。第二类错误则是,H0实际上是错的,却没有被否定。 遗憾的是,不管我们如何选择否定域,都不可能完全避免第一类错误和第二类错误,也不可能同时把犯两类错误的危险压缩到最小。对任何一个给定的检验而言,第一类错误的危险越小,第二类错误的概率就越大;反之亦然。一般来讲,不可能具体估计出第二类错误的概率值。第一类错误则不然,犯第一类错误的概率是否定域内各种结果的概率之和。 两类错误及其关系7/20/202215
12、 被我们事先选定的可以犯第一类错误的概率,叫做检验的显著性水平(用表示),它决定了否定域的大小。因此,有人也把第一类错误称之错误 。相应地第二类错误被人称为 错误。 在原假设成立的条件下,统计检验中所规定的小概率标准一般取为=0.05或=0.01。 由所决定的否定域与接受域之间的分界值被称为临界值, 如Z 。 如果抽样分布是连续的,否定域可以建立在想要建立的任何水平上,否定域的大小可以和显著性水平的要求一致起来(后面的正态检验就如此)。如果抽样分布是非连续的,就要用累计概率的方法找出一组构成否定域的结果。显著性水平7/20/202216 根据否定域位置的不同,可以将假设检验分为双侧检验和单侧检
13、验。 在统计中,必须把否定域分配到抽样分布的两端的检验,被称为双侧检验。 在统计中,可以事先能预测偏差方向,因而可以把否定域集中到抽样分布更合适的一端的检验,被称为单侧检验。 双侧检验和单侧检验7/20/202217 奈曼皮尔逊 (NeymanPearson)提出了一个原则 “在控制犯第一类错误的概率不超过指定值的条件下, 尽量使犯第二类错误 小”按这种法则做出的检验称为“显著性检验”, 称为显著性水平或检验水平。 在同样显著性水平的条件下,单侧检验比双侧检验更合适。因为否定域被集中到抽样分布更合适的一侧,这样在犯第一类错误的危险不变的情况下,减少了犯第二类错误的危险。 7/20/202218
14、4计算检验统计量 在完成了上述工作之后,接下来就是做一次与理想试验尽量相同的实际抽样(比如实际做一次重复抛掷硬币的试验),并从获取的样本资料算出检验统计量。检验统计量是关于样本的一个综合指标,但与我们后面参数估计中将要讨论的统计量有所不同,它不用作估测,而只用作检验。 5判定 假设检验系指拒绝或保留零假设的判断,又称显著性检定。在选择否定域并计算检验统计量之后,我们完成最后一道手续,即根据试验或样本结果决定假设的取与舍。如果结果落在否定域内,我们将在已知犯第一类错误概率的条件下,否定零假设。反之,如果结果落在否定域外,则不否定零假设,与此同时,我们就有了犯第二类错误的危险。 7/20/2022
15、19 例 若想通过抛掷10次硬币的实验来检验这个硬币无偏的零假设,通过双侧检验0.10显著性水平,请指出否定域。如果单侧检验(p0.4)7/20/202220第三节 正态分布 如果说二项分布是离散型随机变量最具典型意义的概率分布,那么连续型随机变量最具典型意义的概率分布就是正态分布了。一般地讲,若影响某一变量的随机因素很多,而每个因素所起的作用不太大且相互独立,则这个变量服从正态分布。更为重要的是,不论总体是否服从正态分布,只要样本容量n足够大,样本平均数的抽样分布就趋于正态分布。 正态分布是最重要的概率分布:(1)许多自然现象和社会现象,都可用正态分布加以叙述;(2)当样本足够大时,都可用正
16、态近似法解决变量的概率分布问题;(3)许多统计量的抽样分布呈正态分布。7/20/2022211. 正态分布的数学形式正态分布性质: (1)正态曲线以x=呈钟型对称 均值=中位数=众数 (2)在x=处,概率密度最大;当区间离越远,x落在这个区间的概率越小。 7/20/202222 (3)正态曲线的外形由值确定。对于固定的值,不同均值的正态曲线的外形完全相同,差别只在于曲线在横轴方向上整体平移了一个位置 。 (5)E(X)= D(X)= 2 (4)对于固定的值,改变值,值越小,正态曲线越陡峭;值越大,正态曲线越低平。 (总之,正态分布曲线的位置是由决定的,而正态分布曲线的“高、矮、胖、瘦” 由决定
17、的。) 7/20/202223 2. 标准正态分布 Z分数(标准正态变量)用Z分数表达的标准正态分布,其概率密度为一般正态分布的表示标准正态分布的表示7/20/2022243. 正态曲线下的面积 但积分毕竟太麻烦了,更何况许多人对积分运算不熟悉,为此须计算出现成的数值表供使用者查找。由于正态曲线的优良性质,这项工作可以卓有成效地完成:经过X的标准分 ,可以将任何正态分布N(,2)转换成标准正态分布N(0,1);运用分布函数的定义,并利用正态曲线的对称性,通过下式(分布函数)可以计算编制出正态分布表(见附4)。 7/20/202225 采用标准正态变量表达正态分布,使标准差得到了进一步阐明。我们
18、看到,标准差是计算总体单位分布及其标志值变异范围的主要依据,下图说明了这一点。(1)变量值在【 -, + 】之间的概率为0.6826。(2)变量值在【 -2, +2 】之间的概率为0.9546。(3)变量值在【 -3, +3 】之间的概率为0.9973。7/20/202226 例 设随机变量X服从正态分布N(168,122),试求P(X143)。 总之,决定任意两点间的面积都完全是可能的。比如向均值两侧移1.96个标准差,曲线下方便包含了大约95的面积;如移动2.58个标准差,则面积几乎是99。附录4已编制了关于Z和标准正态曲线所含面积之间关系的精确数值表,即Z从0到+变化,相应区间含的面积从
19、0变至0.5。7/20/202227 解 已知168,12 z 是负值,表示X的取值处于均值左边。由于曲线完全对称,所以使用正态分布表时可以忽略 z 的正负号。查表可知,正态曲线在均值与z2.08之间所含面积是0.4812。由于总面积的一半是0.5,因P(X143)可以由下面计算求得 P(X143)0.5P(0Z2.08) 0.50.4812 1.88 这说明,X的取值小于或等于143的概率大约是2。由于即将讨论的正态检验几乎都要涉及概率分布的尾端,所以此例说明的是一个非常普遍的问题。 7/20/202228 4. 二项分布的正态近似法 通过前面的讨论,我们已经知道二项分布受成功事件概率p和重
20、复次数n两个参数的影响,只要确定了p和n,二项分布也随之确定了。 但是,二项分布的应用价值实际上受到了n的很大限制。也就是说,只有当n较小时,我们才能比较方便地计算二项分布。所幸的是,二项分布是以正态分布为极限的。所以当n很大时,只要p或q不近于零,我们就可以用正态近似来解决二项分布的计算问题。即以n p、n p q2,将B(x;n,p)视为N(n p,n p q)进行计算。在社会统计中,当n 30,n p、n q均不小于5时,对二项分布作正态近似是可靠的。 7/20/202229 第四节 中心极限定理 一旦统计的学习进入到推论统计,我们就必须同时与三种不同的分布概念打交道,即总体分布、样本分
21、布、抽样分布。为了不产生混淆,视分布不同,将统计指标的符号加以区别是完全必要的。对那些反映标志值集中趋势和离中趋势的综合指标,尤其对均值和标准差(或方差)。均值标准差总体分布样本分布抽样分布 抽样分布特指样本统计量作为随机变量的概率分布。用数学语言来说,抽样分布是运用数理统计的方法,把具体概率赋予样本的所有可能结果的一种理论分布。 在一个总体中可以产生无数个样本,所以样本统计量(比如均值 )必定是随机变量。 这样就提出一个问题:如果样本统计量作为随机变量,它的概率分布是什么样呢?7/20/202230 1中心极限定理 我们知道,概率论中用来阐明大量随机现象平均结果的稳定性的定理,是著名的大数定
22、理。其具体内容是:频率稳定于概率,平均值稳定于期望值。但是,大量随机现象的稳定性不仅表现在平均结果上,同时也表现在分布上,这就是中心极限定理所要阐明的内容。显然,推论统计需要有一座能够架通抽样调查和抽样分布的桥梁。中心极限定理告诉我们:如果从任何一个具有均值和方差2的总体(可以具有任何分布形式)中重复抽取容量为n的随机样本,那么当n变得很大时,样本均值的抽样分布接近正态,并具有均值和方差 。 7/20/202231 (2)由于抽样分布的标准差要比总体标准差小,并且 ,所以如右图所示,样本容量越大,抽样分布的峰态愈陡峭,由样本结果来推断总体参数的可靠性也随之提高。 无疑,中心极限定理大大拓展了正
23、态分布的适用面,同时我们得到了以下重要信息: (1)虽然样本的均值可能和总体均值有差别,但我们可期望这些将聚集在的周围。因此均值抽样分布的算术平均数能和总体的均值很好地重合,这就是为什么总体均值和抽样分布的均值用同一个来表示的缘故。7/20/202232 统计检验应用正态分布和二项分布有两点区别:抽样分布在这里是连续的而非离散的,否定域的大小可以和显著性水平的要求精确地一致起来。计算检验统计量不再像在应用二项分布时那样,可以不劳而获了。很显然,为了能使用现成的正态分布表,关键是要从样本资料中计算出在N(0,1)形式下的统计量Z,再根据Z是否落在否定城内而对被检验假设的取舍作出决定。 在上一节我
24、们曾引出 。Z 的这种形式适用于N(,2)的总体,但并不适用于取正态的抽样分布。正如我们反复强调的那样,统计检验单纯依靠样本自身是得不出结果的,必须首先在一系列假设的基础上求出抽样分布。如果这些假设实际上正确,那么抽样分布将告诉我们得到一个给定的的可能性是多少。在抽样分布中,随机变量的取值是每个 ,均值是,标准差是 。因此 Z 如果作为检验统计量,应该用替换X,用 替换,不动,因而有 。7/20/202233 一个完整的假设应该包括零假设和备择假设。在Z检验中,假设有三种写法:第一种:H0: =90 H1: 90第二种:H0: =90 H1: 90双侧检验单侧检验(右)单侧检验(左)7/20/
25、202234思考题1、如何理解两类错误及其关系?2、如何理解小概率原理在假设检验中的重要性?3、假设检验的基本原理和步骤是什么?4、什么是原假设域备择假设?5、如何理解双侧检验与单侧检验?7/20/202235 例 一位研究者试图检验某一社会调查所运用的抽样程序,该项调查是由一些缺乏经验的访问员进行的。研究者怀疑属于干部和知识分子的家庭抽得过多。过去的统计资料表明,该街区的家庭收入是7500元,标准差是1500元;此次调查共抽取100个家庭,样本平均收入是7900元。问:该研究人员是否有理由怀疑该样本有偏估?(选用=0.05)第五节 总体均值和成数的单样本检验1已知,对总体均值的检验实际上是要
26、检验“随机抽样”这个零假设 7/20/202236 解 根据题意,可做如下假设,并做单侧检验 因=0.05,查表得Z 0.05=1.65,故否定域为根据中心极限定理,检验统计量计算得 检验统计量Z的计算表明,样本均值比总体均值大267个标准差( ),超过了显著性水平规定的临界值,调查者应该否定“随机抽样”的零假设。也就是说,由于抽样在程序上不合要求,这项社会调查有必要重新组织。 7/20/202237 中心极限定理实际解决了大样本均值的检验问题。假定样本比较大(n50,这在社会调查中一般都能得到满足),样本均值的抽样分布就与总体分布无关,而服从正态分布。当H0成立时,样本均值的观察值比较集中地
27、分布在总体均值周围;当H0不成立时, 将对有明显偏离的趋势。因而,我们可以在选定的显著性水平上,通过计算检验统计量Z,对零假设进行检定。 注:当未知时,只要样本量很大,就可用S来代替 。但对于小样本,Z检验就要用 t 检验来替代了,而且还必须严格限于正态总体。7/20/202238 解 根据题意,可作如下的假设,并做双侧检验 H0:2330元 H1:2330元因0.05,查正态分布表得Z/21.96,故否定域|Z|1.96 计算检验统计量 Z 1.20196 所以,不能认为该单位人均月收入不是2330元,即不能认为该统计报表有误。 例 某单位统计报表显示,人均月收入为2330元,为了验证该统计
28、报表的正确性,作了共81人的抽样调查,样本人均月收入为2350元,标准差为150元,问能否说明该统计报表显示的人均 收入的数字有误(取显著性水平0.05)。 此乃“总体均值”零假设的检验7/20/202239提醒:单侧 Z 双侧 Z /20.050.010.0011.652.333.091.962.583.307/20/202240 2.小样本总体均值的检验(学生t分布) 中心极限定理解决了大样本均值的检验问题。但是当n较小时,用这种方法求出的概率可能是错误的,有必要做某种修正。于是有人设计了另一种检验统计量 这个统计量最初是由戈塞特(1876一1937)用笔名“学生”发表,所以这个统计量的抽
29、样分布称为学生t分布。比较t和Z,我们注意到它们的分子相同,而分母却稍有不同:为S所代替(这一点无须解释);根号下是n1。7/20/202241 当Z为t替代时,虽用因子nI所导致的修正看起来不大,但在样本容量较小时,这种修正就会起很大作用了。所以当不知道值、且样本容量较小时,我们应该考虑应用t分布而不是Z分布。 采用n1的原因:样本数据的离散程度小于总体数据的离散程度。 n1实际为自由度数k 。7/20/202242通过下式可以编制出t 分布表(见附录5) 注意t(k)写法的含义,它表示自由度数为k的t分布当分布函数 时随机变量t的临界值。当n变大,t 分布将越来越接近正态分布,应用t分布还
30、是Z分布就无所谓了。事实上随着n变大,S逐渐变成的精确估计量,因而分母项无论使用S还是,差别都非常小。但特别要留意,使用t 分布的条件比使用Z分布的条件苛刻,即必须假定总体为正态。 7/20/202243 例 已知初婚年龄服从正态分布。根据10人的调查有 = 23.5岁,S=3岁,问是否可以认为该地区的平均初婚年龄已超过20岁?(=0.01) 解 H0:=20;H1:20 因为n小,又不知值,因此用t检验 对自由度9来讲,单侧检验和显著性水平0.01,查表知否定域为t值等于或大于 2.821。再计算检验统计量 因此拒绝H0,即可以认为在显著性水平为0.01的条件下,该地区的初婚年龄已超过20岁
31、。7/20/202244 3.大样本成数的检验 有时,需要对总体中具有某种特征的单位在总体中所占的的比例 p(即总体成数)作显著性检验,如人口中的失业率、学龄儿童中的失学率等等。成数检验与二项检验的联系是不言而愈的。因为在二项检验中,随机变量是样本的“成功”次数x。而在成数检验中,随机变量是样本的“成功”比例 (即样本成数),这样在 n 一定的情况下,显然有 7/20/202245 既然 是一个随机变量,那么把具体概率赋予样本成数的每一个取值,我们就得到了样本成数的抽样分布。根据中心极限定理,我们不难想见,当n足够大时,样本成数的抽样分布也服从正态分布。由于数学上很容易证明 , ,这样一来,对于大样本(n30,np5),成数的检验统计量
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 幼儿园学生管理工作制度
- 幼儿园托班全天工作制度
- 幼儿园教师招生工作制度
- 幼儿园早中晚检工作制度
- 幼儿园活动监督工作制度
- 幼儿园督察工作制度汇编
- 幼儿园结核防治工作制度
- 幼儿园评价幼儿工作制度
- 幼儿园食堂卫生工作制度
- 基于遗传算法的XR公司的生产车间设施规划
- 《百年孤独(节选)》课件+2025-2026学年统编版高二语文选择性必修上册
- 青海招警考试真题及答案
- DB11∕T 2271-2024 村庄供水站建设导则
- 江苏省低空空域协同管理办法(试行)
- 肺癌营养支持治疗
- 施工协调费协议书
- 皮肤生理学试题及答案
- 《资治通鉴》与为将之道知到课后答案智慧树章节测试答案2025年春武警指挥学院
- 2018天成消防B-TG-TC5000火灾报警控制器消防联动控制器安装使用说明书
- 配电柜拆除施工方案
- 银行客户满意度调查手册
评论
0/150
提交评论