R语言实验8.docx

1、R语言实验8实验8假设检验（二）一、 实验目的：1. 掌握若干重要的非参数检验方法（?2检验列联表独立性检验，Mcnemar检验对一个样本两种研究方法是否有差异的检验，符号检验，Wilcoxon符号秩检验，Wilcoxon秩和检验）；2. 掌握另外两个相关检验：Spearman秩相关检验，Kendall秩相关检验。二、 实验内容：练习：要求：完成练习并粘贴运行截图到文档相应位置(截图方法见下)，并将所有自己输入文字的字体颜色设为红色（包括后面的思考及小结），回答思考题，简要书写实验小结。修改本文档名为“本人完整学号姓名1”，其中1表示第1次实验，以后更改为2,3,.。如文件名为“09张立1”，

2、表示学号为09的张立同学的第1次实验，注意文件名中没有空格及任何其它字符。最后连同数据文件、源程序文件等（如果有的话，本次实验没有），一起压缩打包发给课代表，压缩包的文件名同上。截图方法：法1：调整需要截图的窗口至合适的大小，并使该窗口为当前激活窗口（即该窗口在屏幕最前方），按住键盘Alt键（空格键两侧各有一个）不放，再按键盘右上角的截图键（通常印有“印屏幕”或“PrScrn”等字符），即完成截图。再粘贴到word文档的相应位置即可。法2：利用QQ输入法的截屏工具。点击QQ输入法工具条最右边的“扳手”图标，选择其中的“截屏”工具。）1. 自行完成教材第五章的例题。2. （习题5.11）为研究分

3、娩过程中使用胎儿电子监测仪对剖腹产率有无影响，对5824例分娩的经产妇进行回顾性调查，结果如下表所示，试进行分析。5824例经产妇回顾性调查结果剖腹产胎儿电子监测仪合计使用未使用是358229587否249227455237合计285029745824解：提出假设：H0：分娩过程中使用胎儿电子监测仪对剖腹产率无影响H1：分娩过程中使用胎儿电子监测仪对剖腹产率有影响源代码及运行结果：（复制到此处，不需要截图）xdim(x)chisq.test(x)PearsonsChi-squaredtestwithYatescontinuitycorrectiondata:xX-squared=37.414,

4、df=1,p-value=9.552e-10结论：P=9.552e-10xdim(x)chisq.test(x)PearsonsChi-squaredtestdata:xX-squared=40.401,df=6,p-value=3.799e-07结论：P=3.799e-07xdim(x)fisher.test(x)FishersExactTestforCountDatadata:xp-value=0.6372alternativehypothesis:trueoddsratioisnotequalto195percentconfidenceinterval:sampleestimates:o

5、ddsratio0.521271结论：P=0.63720.05，接受原假设，两种工艺对产品的质量没影响5. （习题5.14）应用核素法和对比法检测147例冠心病患者心脏收缩运动的符合情况，其结果如下表所示。试分析这两种方法测定结果是否相同。两法检查室壁收缩运动的符合情况对比法核素法合计正常减弱异常正常582363减弱142750异常891734合计675327147解：提出假设：H0：这两种方法测定结果不相同H1：这两种方法测定结果相同源代码及运行结果：（复制到此处，不需要截图）xdim(x)mcnemar.test(x)McNemarsChi-squaredtestdata:xMcNemar

6、schi-squared=2.8561,df=3,p-value=0.4144结论：P=0.41440.05，因此，不能认为这两种方法测定结果不相同6. （习题5.15）在某养鱼塘中，根据过去经验，鱼的长度的中位数为14.6cm，现对鱼塘中鱼的长度进行一次估测，随机地从鱼塘中取出10条鱼长度如下：13.3213.0614.0211.8613.5813.7713.5114.4214.4415.43将它们作为一个样本进行检验。试分析，该鱼塘中鱼的长度是在中位数之上，还是在中位数之下。（1）用符号检验分析；（2）用Wilcoxon符号秩检验分析。解：（1）用符号检验分析提出假设：H0：M=14.6H

7、1：Mxbinom.test(sum(x14.6),length(x),al=l)Exactbinomialtestdata:sum(x14.6)andlength(x)numberofsuccesses=1,numberoftrials=10,p-value=0.01074alternativehypothesis:trueprobabilityofsuccessislessthan0.595percentconfidenceinterval:0.00000000.3941633sampleestimates:probabilityofsuccess0.1（2）用Wilcoxon符号秩检验分

8、析提出假设：H0：M=14.6H1：M14.6源代码及运行结果：（复制到此处，不需要截图）xwilcox.test(x,mu=14.6,al=l,exact=F)Wilcoxonsignedranktestwithcontinuitycorrectiondata:xV=4.5,p-value=0.01087alternativehypothesis:truelocationislessthan14.6结论：两种方法的P都是小于0.05，拒绝原假设，中位数小于14.67. （习题5.16）用两种不同的测定方法，测定同一种中草药的有效成分，共重复20次，得到实验结果如下表所示。两种不同的测定方法得

9、到的结果方法48.033.037.548.042.540.042.036.011.322.0A36.027.314.232.152.038.017.320.021.046.1方法37.041.023.417.031.540.031.036.05.711.5B21.06.126.521.344.528.022.620.011.022.3（1）试用符号检验法来测定有无显着性差异；（2）试用Wilcoxon符号秩检验法检验两测定有无显着差异；（3）试用Wilcoxon秩和检验法检验两测定有无显着差异；（4）对数据作正态性和方差齐性检验，该数据是否能作t检验，如果能，请作t检验；（5）分析各种的检验方

10、法，试说明哪种检验法效果最好。解：（1）符号检验法提出假设：H0：两测定无显着差异H1：两测定有显着差异源代码及运行结果：（复制到此处，不需要截图）xybinom.test(sum(xy),length(x)Exactbinomialtestdata:sum(xy)andlength(x)numberofsuccesses=14,numberoftrials=20,p-value=0.1153alternativehypothesis:trueprobabilityofsuccessisnotequalto0.595percentconfidenceinterval:0.45721080.88

11、10684sampleestimates:probabilityofsuccess0.7结论：P=0.11530.05,接受原假设，两测定无显着差异（2）Wilcoxon符号秩检验法提出假设：H0：两测定无显着差异H1：两测定有显着差异源代码及运行结果：（复制到此处，不需要截图）xywilcox.test(x,y,paired=T,exact=F)Wilcoxonsignedranktestwithcontinuitycorrectiondata:xandyV=136,p-value=0.005191alternativehypothesis:truelocationshiftisnotequ

12、alto0结论：P=0.005191xywilcox.test(x,y,exact=F)Wilcoxonranksumtestwithcontinuitycorrectiondata:xandyW=274.5,p-value=0.04524alternativehypothesis:truelocationshiftisnotequalto0结论：P=0.04524xyks.test(x,pnorm,mean(x),sd(x)One-sampleKolmogorov-Smirnovtestdata:xD=0.14067,p-value=0.8235alternativehypothesis:t

13、wo-sidedWarningmessage:Inks.test(x,pnorm,mean(x),sd(x):Kolmogorov-Smirnov检验里不应该有连结ks.test(y,pnorm,mean(y),sd(y)One-sampleKolmogorov-Smirnovtestdata:yD=0.10142,p-value=0.973alternativehypothesis:two-sided结论：两组数据P值都大于0.05，接受原假设，两组数据服从正态分布方差齐性检验提出假设：H0：两组数据方差相同H1：两组数据方差不相同源代码及运行结果：（复制到此处，不需要截图）xyvar.te

14、st(x,y)Ftesttocomparetwovariancesdata:xandyF=1.1406,numdf=19,denomdf=19,p-value=0.7772alternativehypothesis:trueratioofvariancesisnotequalto195percentconfidenceinterval:0.45147882.8817689sampleestimates:ratioofvariances1.140639结论：P=0.77720.05,接受原假设，两组数据方差相同能否作t检验？如果能，请按照前面各题的格式、内容写出作t检验的过程综上，能作t检验提出

15、假设：H0：两测定无显着差异H1：两测定有显着差异源代码及运行结果：（复制到此处，不需要截图）xyt.test(x,y,var.equal=T)TwoSamplet-testdata:xandyt=2.2428,df=38,p-value=0.03082alternativehypothesis:truedifferenceinmeansisnotequalto095percentconfidenceinterval:0.812552915.8774471sampleestimates:meanofxmeanofy33.21524.870结论：P=0.03082xycor.test(x,y,m

16、ethod=spearman,exact=F)Spearmansrankcorrelationrhodata:xandyS=10,p-value=5.484e-05alternativehypothesis:truerhoisnotequalto0sampleestimates:rho0.9393939结论：P=5.484e-05xycor.test(x,y,method=kendall,exact=F)Kendallsrankcorrelationtaudata:xandyz=3.3094,p-value=0.000935alternativehypothesis:truetauisnote

17、qualto0sampleestimates:tau0.8222222结论：P=0.000935xywilcox.test(x,y,exact=F)Wilcoxonranksumtestwithcontinuitycorrectiondata:xandyW=266,p-value=0.05509alternativehypothesis:truelocationshiftisnotequalto0结论：P=0.055090.05,接受原假设，即认为新方法的疗效不显着地优于原疗法思考：1. ?2检验针对的哪一类型的数据？A.连续型变量；B.无序分类变量；C.有序分类变量B.无序分类变量2. 符号

18、检验、Wilcoxon符号秩检验和Wilcoxon秩和检验针对的哪一类型的数据？A.连续型变量；B.无序分类变量；C.有序分类变量C.有序分类变量3. 顺序统计量与秩统计量是一样的吗？在R语言中分别利用什么函数可以求得？秩统计计量是按顺序排列还是按降序排列？不一样。顺序统计量是将样本观测值从小到大排顺序。在R中，可以由sort()函数求出一个顺序统计量。而秩统计量是在非参数检验中广泛运用的统计量，它的一个重要的特性是分布无关性。其样本不必取自同一总体。在R中，可以由rank()函数求出一个秩统计量。秩统计计量是按顺序排列的。4. 非参数检验中在检验一个样本是否来自某个总体，或者检验两个总体是否

19、有显着性差异时，为什么不能用均值这个统计量，而是用中位数这个统计量？由于不再涉及到分布类型，因此不能使用均值这一与总体分布有关的参数加以检验。所以，改用另一个同样描述集中趋势的统计量中位数M来代替均值U。5. 为什么符号秩检验的功效比符号检验好？（P239）因为符号检验利用了观测值跟原假设的中心位置之差的符号来进行检验，但是它只计算符号的个数，而不考虑每个符号差中所包含的绝对值的大小，而符号秩检验还利用了绝对值的大小，考虑了绝对值距中心位置的远近，自然比仅仅利用符号更有效。6. 符号检验和Wilcoxon符号秩检验都可以用于单个样本的检验，来判断一个样本是否来自某个总体；也可以用于两配对（成对）样本的检验，来判断两个相关总体是否存在显着性差异。Wilcoxon秩和检验主要是用于哪类问题的检验？两个配对样本的检验、两个独立样本的检验。7. 参数检验中的均值检验，两配对样本t检验可以转化为单样本t检验，即检验两个样本的差值序列的均值是否为0。非参数检验中的中位数检验，两个