R语言实验8.docx
《R语言实验8.docx》由会员分享，可在线阅读，更多相关《R语言实验8.docx（13页珍藏版）》请在冰豆网上搜索。
R语言实验8
实验8假设检验
（二）
一、实验目的：
1.掌握若干重要的非参数检验方法（2检验——列联表独立性检验，Mcnemar检验——对一个样本两种研究方法是否有差异的检验，符号检验，Wilcoxon符号秩检验，Wilcoxon秩和检验）；
2.掌握另外两个相关检验：
Spearman秩相关检验，Kendall秩相关检验。
二、实验内容：
练习：
要求：
①完成练习并粘贴运行截图到文档相应位置（截图方法见下），并将所有自己输入文字的字体颜色设为红色（包括后面的思考及小结），②回答思考题，③简要书写实验小结。
④修改本文档名为“本人完整学号姓名1”，其中1表示第1次实验，以后更改为2,3,...。
如文件名为“09张立1”，表示学号为09的张立同学的第1次实验，注意文件名中没有空格及任何其它字符。
最后连同数据文件、源程序文件等（如果有的话，本次实验没有），一起压缩打包发给课代表，压缩包的文件名同上。
截图方法：
法1：
调整需要截图的窗口至合适的大小，并使该窗口为当前激活窗口（即该窗口在屏幕最前方），按住键盘Alt键（空格键两侧各有一个）不放，再按键盘右上角的截图键（通常印有“印屏幕”或“PrScrn”等字符），即完成截图。
再粘贴到word文档的相应位置即可。
法2：
利用QQ输入法的截屏工具。
点击QQ输入法工具条最右边的“扳手”图标
，选择其中的“截屏”工具。
）
1.自行完成教材第五章的例题。
2.（习题）为研究分娩过程中使用胎儿电子监测仪对剖腹产率有无影响，对5824例分娩的经产妇进行回顾性调查，结果如下表所示，试进行分析。
5824例经产妇回顾性调查结果
剖腹产
胎儿电子监测仪
合计
使用
未使用
是
358
229
587
否
2492
2745
5237
合计
2850
2974
5824
解：
提出假设：
H0：
分娩过程中使用胎儿电子监测仪对剖腹产率无影响
H1：
分娩过程中使用胎儿电子监测仪对剖腹产率有影响
源代码及运行结果：
（复制到此处，不需要截图）
>x<-c（358,2492,229,2745）
>dim（x）<-c（2,2）
>（x）
Pearson'sChi-squaredtestwithYates'continuitycorrection
data:
x
X-squared=,df=1,p-value=
结论：
P=<,拒绝原假设，分娩过程中使用胎儿电子监测仪对剖腹产率有影响
3.（习题）在高中一年级男生中抽取300名考察其两个属性：
B是1500米长跑，C是每天平均锻炼时间，得到4×3列联表，如下表所示。
试对=，检验B与C是否独立。
300名高中学生体育锻炼的考察结果
1500米
长跑记录
锻炼时间
合计
2小时以上
1〜2小时
1小时以下
5"01'〜5"30'
45
12
10
67
5"31'〜6"00'
46
20
28
94
6"01'〜6"30'
28
23
30
81
6"31'〜7"00'
11
12
35
58
合计
130
67
103
300
解：
提出假设：
H0：
B与C独立。
H1：
B与C不独立。
源代码及运行结果：
（复制到此处，不需要截图）
>x<-c（45,46,28,11,12,20,23,12,10,28,30,35）
>dim（x）<-c（4,3）
>（x）
Pearson'sChi-squaredtest
data:
x
X-squared=,df=6,p-value=
结论：
P=<,拒绝原假设，B与C不独立，有关系。
4.（习题）为比较两种工艺对产品的质量是否有影响，对其产品进行抽样检查，其结果如下表所示。
试进行分析。
两种工艺下产品质量的抽查结果
合格
不合格
合计
工艺一
3
4
7
工艺二
6
4
10
合计
9
8
17
解：
提出假设：
H0：
两种工艺对产品的质量没影响
H1：
两种工艺对产品的质量有影响
源代码及运行结果：
（复制到此处，不需要截图）
>x<-c（3,6,4,4）
>dim（x）<-c（2,2）
>（x）
Fisher'sExactTestforCountData
data:
x
p-value=
alternativehypothesis:
trueoddsratioisnotequalto1
95percentconfidenceinterval:
5.
sampleestimates:
oddsratio
结论：
P=>，接受原假设，两种工艺对产品的质量没影响
5.（习题）应用核素法和对比法检测147例冠心病患者心脏收缩运动的符合情况，其结果如下表所示。
试分析这两种方法测定结果是否相同。
两法检查室壁收缩运动的符合情况
对比法
核素法
合计
正常
减弱
异常
正常
58
2
3
63
减弱
1
42
7
50
异常
8
9
17
34
合计
67
53
27
147
解：
提出假设：
H0：
这两种方法测定结果不相同
H1：
这两种方法测定结果相同
源代码及运行结果：
（复制到此处，不需要截图）
>x<-c（58,1,8,2,42,9,3,7,17）
>dim（x）<-c（3,3）
>（x）
McNemar'sChi-squaredtest
data:
x
McNemar'schi-squared=,df=3,p-value=
结论：
P=>，因此，不能认为这两种方法测定结果不相同
6.（习题）在某养鱼塘中，根据过去经验，鱼的长度的中位数为，现对鱼塘中鱼的长度进行一次估测，随机地从鱼塘中取出10条鱼长度如下：
将它们作为一个样本进行检验。
试分析，该鱼塘中鱼的长度是在中位数之上，还是在中位数之下。
（1）用符号检验分析；
（2）用Wilcoxon符号秩检验分析。
解：
（1）用符号检验分析
提出假设：
H0：
M>=
H1：
M<
源代码及运行结果：
（复制到此处，不需要截图）
>x<-c,,,,,,,,,
>（sum（x>,length（x）,al="l"）
Exactbinomialtest
data:
sum（x>andlength（x）
numberofsuccesses=1,numberoftrials=10,p-value=
alternativehypothesis:
trueprobabilityofsuccessislessthan
95percentconfidenceinterval:
sampleestimates:
probabilityofsuccess
（2）用Wilcoxon符号秩检验分析
提出假设：
H0：
M>=
H1：
M<
源代码及运行结果：
（复制到此处，不需要截图）
x<-c,,,,,,,,,
>（x,mu=,al="l",exact=F）
Wilcoxonsignedranktestwithcontinuitycorrection
data:
x
V=,p-value=
alternativehypothesis:
truelocationislessthan
结论：
两种方法的P都是小于，拒绝原假设，中位数小于
7.（习题）用两种不同的测定方法，测定同一种中草药的有效成分，共重复20次，得到实验结果如下表所示。
两种不同的测定方法得到的结果
方法
A
方法
B
（1）试用符号检验法来测定有无显著性差异；
（2）试用Wilcoxon符号秩检验法检验两测定有无显著差异；
（3）试用Wilcoxon秩和检验法检验两测定有无显著差异；
（4）对数据作正态性和方差齐性检验，该数据是否能作t检验，如果能，请作t检验；
（5）分析各种的检验方法，试说明哪种检验法效果最好。
解：
（1）符号检验法
提出假设：
H0：
两测定无显著差异
H1：
两测定有显著差异
源代码及运行结果：
（复制到此处，不需要截图）
>x<-c,,,,,,,,,,,,,,,,,,,
>y<-c,,,,,,,,,,,,,,,,,,,
>（sum（x>y）,length（x））
Exactbinomialtest
data:
sum（x>y）andlength（x）
numberofsuccesses=14,numberoftrials=20,p-value=
alternativehypothesis:
trueprobabilityofsuccessisnotequalto
95percentconfidenceinterval:
sampleestimates:
probabilityofsuccess
结论：
P=>,接受原假设，两测定无显著差异
（2）Wilcoxon符号秩检验法
提出假设：
H0：
两测定无显著差异
H1：
两测定有显著差异
源代码及运行结果：
（复制到此处，不需要截图）
>x<-c,,,,,,,,,,,,,,,,,,,
>y<-c,,,,,,,,,,,,,,,,,,,
>（x,y,paired=T,exact=F）
Wilcoxonsignedranktestwithcontinuitycorrection
data:
xandy
V=136,p-value=
alternativehypothesis:
truelocationshiftisnotequalto0
结论：
P=<,拒绝原假设，两测定有显著差异
（3）Wilcoxon秩和检验法
提出假设：
H0：
两测定无显著差异
H1：
两测定有显著差异
源代码及运行结果：
（复制到此处，不需要截图）
>x<-c,,,,,,,,,,,,,,,,,,,
>y<-c,,,,,,,,,,,,,,,,,,,
>（x,y,exact=F）
Wilcoxonranksumtestwithcontinuitycorrection
data:
xandy
W=,p-value=
alternativehypothesis:
truelocationshiftisnotequalto0
结论：
P=<,拒绝原假设，两测定有显著差异
（4）两独立样本t检验，要求判断独立性、正态性和方差齐性
①正态性检验
提出假设：
H0：
两组数据服从正态分布
H1：
两组数据不服从正态分布
源代码及运行结果：
（复制到此处，不需要截图）
>x<-c,,,,,,,,,,,,,,,,,,,
>y<-c,,,,,,,,,,,,,,,,,,,
>（x,pnorm,mean（x）,sd（x））
One-sampleKolmogorov-Smirnovtest
data:
x
D=
升级会员 