R语言实验七.docx
《R语言实验七.docx》由会员分享，可在线阅读，更多相关《R语言实验七.docx（16页珍藏版）》请在冰豆网上搜索。
R语言实验七
实验7假设检验
（一）
、实验目的：
1.掌握重要的参数检验方法（单个总体的均值检验，两个总体的均值检验，成对样本
的均值的检验，两个总体方差的检验，二项分布总体的检验）；
2.掌握若干重要的非参数检验方法（Pearson拟合优度2检验，Kolmogorov-Smirnov单样本和双样本检验）。
、实验内容：
练习：
要求：
①完成练习并粘贴运行截图到文档相应位置（截图方法见下），并将所有自己输入文字
的字体颜色设为红色（包括后面的思考及小结），②回答思考题，③简要书写实验小结。
④修改本文档名为“本人完整学号姓名1”，其中1表示第1次实验，以后更改为2,3,…。
如文件名为“1305543109张立1”，表示学号为1305543109的张立同学的第1次实验，注意文件名中没有空格及任何其它字符。
最后连同数据文件、源程序文件等（如果有的话，本次实验没有），一起压缩打包发给课代表，压缩包的文件名同上。
截图方法：
法1:
调整需要截图的窗口至合适的大小，并使该窗口为当前激活窗口（即该窗口在屏幕最前方），按住键盘Alt键（空格键两侧各有一个）不放，再按键盘右上角的截图键（通常印有“印屏幕”或“PrScrn”等字符），即完成截图。
再粘贴到word文档的相应位置即可。
法2：
利用QQ输入法的截屏工具。
点击QQ输入法工具条最右边的“扳手”图标[ZEI，选择其中的“截屏”工具。
）
1.自行完成教材第五章的例题。
2.（习题5.1）正常男子血小板计数均值为225X109/L,今测得20名男性油漆作业
工人的血小板计数值（单位：
109/L）
220188162230145160238188247113
126245164231256183190158224175
问油漆工人的血小板计数与正常成年男子有无差异？
解：
提出假设：
H0：
油漆工人的血小板计数与正常成年男子无差异
H1：
油漆工人的血小板计数与正常成年男子有差异源代码及运行结果：
（复制到此处，不需要截图）
>x<-c（220,188,162,230,145,160,238,188,247,113,126,245,164,231,256,183,190,158,224,175）
>t.test（x,mu=225）
结论：
OneSamplet-test
data:
x
t=-3.4783,df=19,p-value=0.002516
alternativehypothesis:
truemeanisnotequalto225
95percentconfideneeinterval:
172.3827211.9173
sampleestimates:
meanofx
192.15
P=0.002516<0.05,拒绝原假设，认为油漆工人的血小板计数与正常成年男子有差异
3.（习题5.2）已知某种灯泡寿命服从正态分布，在某星期所生产的该灯泡中随机抽取10只，测得其寿命（单位：
小时）为
1067919119678511269369181156920948
求这个星期生产出的灯泡能使用1000小时以上的概率。
解：
源代码及运行结果：
（复制到此处，不需要截图）
>x<-c（1067,919,1196,785,1126,936,918,1156,920,948）
>p<-pnorm（1000,mean（x）,sd（x））
>1-p
[1]0.4912059
结论：
这个星期生产出的灯泡能使用1000小时以上的概率为0.4912059
4.（习题5.3）为研究某铁剂治疗和饮食治疗營养性缺铁性贫血的效果，将16名患者
按年龄、体重、病程和病情相近的原则配成8对，分别使用饮食疗法和补充铁剂治
疗的方法，3个月后测得两种患者血红资白如下表所示，问两种方法治疗后的患者血红蛋白有无差异？
铁剂和饮食两种方法治疗后患者血红蛋白值（g/L）
铁剂治疗组
113120138120100118138123
饮食治疗组
138116125136110132130110
解：
提出假设：
H0：
两种方法治疗后的患者血红蛋白无差异
H1：
两种方法治疗后的患者血红蛋白有差异源代码及运行结果：
（复制到此处，不需要截图）
>x<-c（113,120,138,120,100,118,138,123）
>y<-c（138,116,125,136,110,132,130,110）
>t.test（x,y,paired=T）
Pairedt-test
data:
xandy
t=-0.65127,df=7,p-value=0.5357
alternativehypothesis:
truedifferenceinmeansisnotequalto095percentconfidenceinterval:
-15.6288918.878891
sampleestimates:
meanofthedifferences
-3.375
结论：
p=0.5357>0.05,不拒绝原假设，两种方法治疗后的患者血红蛋白无差异
5.（习题5.4）为研究国产四类新药阿卡波糖股囊效果，某医院用40名H型糖尿病
病人进行同期随机对照实验。
试验者将这些病人随机等分到试验组（阿卡波糖股囊
组）和对照组（拜唐苹股囊组），分别测得试验开始前和8周后空腹血糖，算得空腹血糖下降值，如下所示。
能否认为国产四类新药阿卡波糖股囊与拜唐苹股囊对空腹血糖的降糖效果不同？
试验组与对照组空腹腔血糖下降值（mmol/L）
试验组
-0.70
-5.60
2.00
2.80
0.70
3.50
4.00
5.80
7.10
-0.50
（n1=20）:
>.50
-1.60
1.70
3.00
0.40
4.50
4.60
2.50
6.00
-1.40
对照组
3.70
6.50
5.00
5.20
0.80
0.20
0.60
3.40
6.60
-1.10
（n2=20）6
）.00
3.80
2.00
1.60
2.00
2.20
1.20
3.10
1.70
-2.00
（1）检验试验组和对照组的的数据是否来自正态分布，采用正态性W检验方法（见第
3章）、Kolmogorov-Smirnov检验方法和Pearson拟合优度2检验；
解：
提出假设：
H0：
认为国产四类新药阿卡波糖股囊与拜唐苹股囊对空腹血糖的降糖效果不同
H1:
认为国产四类新药阿卡波糖股囊与拜唐苹股囊对空腹血糖的降糖效果相同
1正态性W检验方法源代码及运行结果：
（复制到此处，不需要截图）
>x<-c（-0.70,-5.60,2.00,2.80,0.70,3.50,4.00,5.80,7.10,-0.50,2.50,-1.60,1.70,3.00,0.40,4.50,4.60,2.50,6.00,-1.4）
>shapiro.test（x）
Shapiro-Wilknormalitytest
data:
x
W=0.9699,p-value=0.7527
>y<-c（3.70,6.50,5.00,5.20,0.80,0.20,0.60,3.40,6.60,-1.10,6.00,3.80,2.00,1.60,2.00,2.20,1.20,
3.10,1.70,-2.00）
>shapiro.test（y）
Shapiro-Wilknormalitytest
data:
y
W=0.97098,p-value=0.7754
结论：
试验组p=0.7527>0.05,对照组p=0.7754>0.05,所以检验试验组和对照组的的数据是来自
正态分布
2Kolmogorov-Smirnov检验方法
源代码及运行结果：
（复制到此处，不需要截图）
>ks.test（x,"pnorm",mean（x）,sd（x））
One-sampleKolmogorov-Smirnovtest
data:
x
D=0.10652,p-value=0.9771
alternativehypothesis:
two-sided
Warningmessage:
Inks.test（x,"pnorm",mean（x）,sd（x））:
Kolmogorov-Smirnov检验里不应该有连结
>ks.test（y,"pnorm",mean（y）,sd（y））
One-sampleKolmogorov-Smirnovtest
data:
y
D=0.11969,p-value=0.9368
alternativehypothesis:
two-sided
Warningmessage:
Inks.test（y,"pnorm",mean（y）,sd（y））:
Kolmogorov-Smirnov检验里不应该有连结
结论：
试验组p=0.9771>0.05,对照组p=0.9368>0.05,所以检验试验组和对照组的的数据是来自
正态分布
3Pearson拟合优度2检验
源代码及运行结果：
（复制到此处，不需要截图）
>x<-c（-0.70,-5.60,2.00,2.80,0.70,3.50,4.00,5.80,7.10,-0.50,2.50,-1.60,1.70,3.00,0.40,4.50,4.6
0,2.50,6.00,-1.4）
>A<-table（cut（x,br=c（-6,-3,0,3,6,9）））
>p<-pnorm（c（-3,0,3,6,9）,mean（x）,sd（x））
>p
[1]0.048947120.249900090.620022880.900758560.98828138
>p<-c（p[1],p[2]-p[1],p[3]-p[2],p[4]-p[3],1-p[4]）
>p
[1]0.048947120.200952980.370122780.280735680.09924144
>chisq.test（A,p=p）
Chi-squaredtestforgivenprobabilities
data:
A
X-squared=0.56387,df=4,p-value=0.967
Warningmessage:
Inchisq.test（A,p=p）:
Chi-squared近似算法有可能不准
>y<-c（3.70,6.50,5.00,5.20,0.80,0.20,0.60,3.40,6.60,-1.10,6.00,3.80,2.00,1.60,2.00,2.20,1.20,
3.10,1.70,-2.00）
>B<-table（cut（y,br=c（-2,1,2,4,7）））
>p<-pnorm（c（-2,1,2,4,7）,mean（y）,sd（y））
>p
[1]0.028017520.250999070.398121950.715004470.96465615
>p<-
升级会员 