R语言实验8文档格式.docx
《R语言实验8文档格式.docx》由会员分享，可在线阅读，更多相关《R语言实验8文档格式.docx（14页珍藏版）》请在冰豆网上搜索。
2850
2974
5824
解：
提出假设：
H0：
分娩过程中使用胎儿电子监测仪对剖腹产率无影响
H1：
分娩过程中使用胎儿电子监测仪对剖腹产率有影响
源代码及运行结果：
（复制到此处，不需要截图）
>
x<
-c（358,2492,229,2745）
dim（x）<
-c（2,2）
chisq.test（x）
Pearson'
sChi-squaredtestwithYates'
continuitycorrection
data:
x
X-squared=37.414,df=1,p-value=9.552e-10
结论：
P=9.552e-10<
0.05,拒绝原假设，分娩过程中使用胎儿电子监测仪对剖腹产率有影响
3.（习题5.12）在高中一年级男生中抽取300名考察其两个属性：
B是1500米长跑，C是每天平均锻炼时间，得到4×
3列联表，如下表所示。
试对?
?
=0.05，检验B与C是否独立。
300名高中学生体育锻炼的考察结果
1500米
长跑记录
锻炼时间
2小时以上
1?
2小时
1小时以下
5"
01'
30'
45
12
10
67
31'
6"
00'
46
20
28
94
23
30
81
7"
11
35
58
130
103
300
B与C独立。
B与C不独立。
-c（45,46,28,11,12,20,23,12,10,28,30,35）
-c（4,3）
sChi-squaredtest
X-squared=40.401,df=6,p-value=3.799e-07
P=3.799e-07<
0.05,拒绝原假设，B与C不独立，有关系。
4.（习题5.13）为比较两种工艺对产品的质量是否有影响，对其产品进行抽样检查，其结果如下表所示。
试进行分析。
两种工艺下产品质量的抽查结果
合格
不合格
工艺一
3
4
7
工艺二
6
9
8
17
两种工艺对产品的质量没影响
两种工艺对产品的质量有影响
-c（3,6,4,4）
fisher.test（x）
Fisher'
sExactTestforCountData
p-value=0.6372
alternativehypothesis:
trueoddsratioisnotequalto1
95percentconfidenceinterval:
sampleestimates:
oddsratio
0.521271
P=0.6372>
0.05，接受原假设，两种工艺对产品的质量没影响
5.（习题5.14）应用核素法和对比法检测147例冠心病患者心脏收缩运动的符合情况，其结果如下表所示。
试分析这两种方法测定结果是否相同。
两法检查室壁收缩运动的符合情况
对比法
核素法
正常
减弱
异常
2
63
1
42
50
34
53
27
147
这两种方法测定结果不相同
这两种方法测定结果相同
-c（58,1,8,2,42,9,3,7,17）
-c（3,3）
mcnemar.test（x）
McNemar'
schi-squared=2.8561,df=3,p-value=0.4144
P=0.4144>
0.05，因此，不能认为这两种方法测定结果不相同
6.（习题5.15）在某养鱼塘中，根据过去经验，鱼的长度的中位数为14.6cm，现对鱼塘中鱼的长度进行一次估测，随机地从鱼塘中取出10条鱼长度如下：
13.3213.0614.0211.8613.5813.7713.5114.4214.4415.43
将它们作为一个样本进行检验。
试分析，该鱼塘中鱼的长度是在中位数之上，还是在中位数之下。
（1）用符号检验分析；
（2）用Wilcoxon符号秩检验分析。
（1）用符号检验分析
M>
=14.6
M<
14.6
-c（13.32,13.06,14.02,11.86,13.58,13.77,13.51,14.42,14.44,15.43）
binom.test（sum（x>
14.6）,length（x）,al="
l"
Exactbinomialtest
sum（x>
14.6）andlength（x）
numberofsuccesses=1,numberoftrials=10,p-value=0.01074
trueprobabilityofsuccessislessthan0.5
0.00000000.3941633
probabilityofsuccess
0.1
（2）用Wilcoxon符号秩检验分析
wilcox.test（x,mu=14.6,al="
exact=F）
Wilcoxonsignedranktestwithcontinuitycorrection
V=4.5,p-value=0.01087
truelocationislessthan14.6
两种方法的P都是小于0.05，拒绝原假设，中位数小于14.6
7.（习题5.16）用两种不同的测定方法，测定同一种中草药的有效成分，共重复20次，得到实验结果如下表所示。
两种不同的测定方法得到的结果
方法
48.0
33.0
37.5
42.5
40.0
42.0
36.0
11.3
22.0
A
27.3
14.2
32.1
52.0
38.0
17.3
20.0
21.0
46.1
37.0
41.0
23.4
17.0
31.5
31.0
5.7
11.5
B
6.1
26.5
21.3
44.5
28.0
22.6
11.0
22.3
（1）试用符号检验法来测定有无显着性差异；
（2）试用Wilcoxon符号秩检验法检验两测定有无显着差异；
（3）试用Wilcoxon秩和检验法检验两测定有无显着差异；
（4）对数据作正态性和方差齐性检验，该数据是否能作t检验，如果能，请作t检验；
（5）分析各种的检验方法，试说明哪种检验法效果最好。
（1）符号检验法
两测定无显着差异
两测定有显着差异
-c（48.0,33.0,37.5,48.0,42.5,40.0,42.0,36.0,11.3,22.0,36.0,27.3,14.2,32.1,52.0,38.0,17.3,20.0,21.0,46.1）
y<
-c（37.0,41.0,23.4,17.0,31.5,40.0,31.0,36.0,5.7,11.5,21.0,6.1,26.5,21.3,44.5,28.0,22.6,20.0,11.0,22.3）
y）,length（x））
y）andlength（x）
numberofsuccesses=14,numberoftrials=20,p-value=0.1153
trueprobabilityofsuccessisnotequalto0.5
0.45721080.8810684
0.7
P=0.1153>
0.05,接受原假设，两测定无显着差异
（2）Wilcoxon符号秩检验法
wilcox.test（x,y,paired=T,exact=F）
xandy
V=136,p-value=0.005191
truelocationshiftisnotequalto0
P=0.005191<
0.05,拒绝原假设，两测定有显着差异
（3）Wilcoxon秩和检验法
wilcox.test（x,y,exact=F）
Wilcoxonranksumtestwithcontinuitycorrection
W=274.5,p-value=0.04524
P=0.04524<
（4）两独立样本t检验，要求判断独立性、正态性和方差齐性
①正态性检验
两组数据服从正态分布
两组数据不服从正态分布
ks.test（x,pnorm,mean（x）,sd（x））
One-sampleKolmogorov-Smirnovtest
D=0.14067,p-value=0.8235
two-sided
Warningmessage:
Inks.test（x,pnorm,mean（x）,sd（x））:
Kolmogorov-Smirnov检验里不应该有连结
ks.test（y,pnorm,mean（y）,sd（y））
y
D=0.10142,p-value=0.973
两组数据P值都大于0.05，接受原假设，两组数据服从正态分布
②方差齐性检验
两组数据方差相同
两组数据方差不相同
var.test（x,y）
Ftesttocomparetwovariances
F=1.1406,numdf=19,denomdf=19,p-value=0.7772
trueratioofvariancesisnotequalto1
0.45147882.8817689
ratioofvariances
1.140639
P=0.7772>
0.05,接受原假设，两组数据方差相同
③能否作t检验？
如果能，请按照前面各题的格式、内容写出作t检验的过程
综上，能作t检验
t.test（x,y,var.equal=T）
TwoSamplet-test
t=2.2428,df=38,p-value=0.03082
truedifferenceinmeansisnotequalto0
0.812552915.8774471
meanofxmeanofy
33.21524.870
P=0.03082<
综上所述，Wilcoxon秩和检验法的差异性检出能力最强，符号检验的差异性检出最弱。
8.（习题5.17）调查某大学学生每周学习时间与得分的平均等级之间的关系，现抽查10个学生的资料如下：
学习时间
24
41
52
18
15
29
学习等级
5
其中等级10表示最好，1表示最差。
试用秩相关检验（Spearman检验和Kendall检验）分别分析学习等级与学习成绩有无关系。
（1）用Spearman秩检验方法
学习等级与学习成绩无关系
学习等级与学习成绩有关系
-c（24,17,20,41,52,23,46,18,15,29）
-c（8,1,4,7,9,5,10,3,2,6）
cor.test（x,y,method="
spearman"
Spearman'
srankcorrelationrho
S=10,p-value=5.484e-05
truerhoisnotequalto0
rho
0.9393939
P=5.484e-05<
0.05，拒绝原假设，学习等级与学习成绩有关系
（2）用Kendall秩检验方法
kendall"
Kendall'
srankcorrelationtau
z=3.3094,p-value=0.000935
truetauisnotequalto0
tau
0.8222222
P=0.000935<
9.（习题5.18）为比较一种新疗法对某种疾病的治疗效果，将40名患者随机地分为两组，每组20人，一组采用新疗法，另一组用原标准疗法。
经过一段时间的治疗后，对每个患者的疗效作仔细的评估，并划分为差、较差、一般、较好和好5个等级。
两组处于不同等级的患者人数下表所示。
试分析，由此结果能否认为新方法的疗效显着地优于原疗法。
（?
=0.05）
不同方法治疗后的结果
等级
差
较差
一般
较好
好
新疗法组
原疗法组
新方法的疗效不显着地优于原疗法
新方法的疗效显着地优于原疗法
-rep（1:
5,c（0,1,9,7,3））
5,c（2,2,11,4,1））
W=266,p-value=0.05509
P=0.05509>
0.05,接受原假设，即认为新方法的疗效不显着地优于原疗法
思考：
1.?
2检验针对的哪一类型的数据？
A.连续型变量；
B.无序分类变量；
C.有序分类变量
B.无序分类变量
2.符号检验、Wilcoxon符号秩检验和Wilcoxon秩和检验针对的哪一类型的数据？
3.顺序统计量与秩统计量是一样的吗？
在R语言中分别利用什么函数可以求得？
秩统计计量是按顺序排列还是按降序排列？
不一样。
顺序统计量是将样本观测值从小到大排顺序。
在R中，可以由sort（）函数求出一个顺序统计量。
而秩统计量是在非参数检验中广泛运用的统计量，它的一个重要的特性是分布无关性。
其样本不必取自同一总体。
在R中，可以由rank（）函数求出一个秩统计量。
秩统计计量是按顺序排列的。
4.非参数检验中在检验一个样本是否来自某个总体，或者检验两个总体是否有显着性差异时，为什么不能用均值这个统计量，而是用中位数这个统计量？
由于不再涉及到分布类型，因此不能使用均值这一与总体分布有关的参数加以检验。
所以，改用另一个同样描述集中趋势的统计量中位数M来代替均值U。
5.为什么符号秩检验的功效比符号检验好？
（P239）
因为符号检验利用了观测值跟原假设的中心位置之差的符号来进行检验，但是它只计算符号的个数，而不考虑每个符号差中所包含的绝对值的大小，而符号秩检验还利用了绝对值的大小，考虑了绝对值距中心位置的远近，自然比仅仅利用符号更有效。
6.符号检验和Wilcoxon符号秩检验都可以用于单个样本的检验，来判断一个样本是否来自某个总体；
也可以用于两配对（成对）样本的检验，来判断两个相关总体是否存在显着性差异。
Wilcoxon秩和检验主要是用于哪类问题的检验？
两个配对样本的检验、两个独立样本的检验。
7.参数检验中的均值检验，两配对样本t检验可以转化为单样本t检验，即检验两个样本的差值序列的均值是否为0。
非参数检验中的中位数检验，两个
升级会员 