R语言 实验七Word文档下载推荐.docx
《R语言 实验七Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《R语言 实验七Word文档下载推荐.docx(15页珍藏版)》请在冰豆网上搜索。
![R语言 实验七Word文档下载推荐.docx](https://file1.bdocx.com/fileroot1/2022-10/14/07bf00c9-4b31-43c8-a16f-678d2ae6820d/07bf00c9-4b31-43c8-a16f-678d2ae6820d1.gif)
再粘贴到word文档的相应位置即可。
法2:
利用QQ输入法的截屏工具。
点击QQ输入法工具条最右边的“扳手〞图标,选择其中的“截屏〞工具。
〕
1.自行完成教材第五章的例题。
2.〔习题5.1〕正常男子血小板计数均值为225×
109/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=
alternativehypothesis:
truemeanisnotequalto225
95percentconfidenceinterval:
sampleestimates:
meanofx
P=<
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
这个星期生产出的灯泡能使用1000小时以上的概率为
4.〔习题5.3〕为研究某铁剂治疗和饮食治疗營养性缺铁性贫血的效果,将16名患者按年龄、体重、病程和病情相近的原那么配成8对,分别使用饮食疗法和补充铁剂治疗的方法,3个月后测得两种患者血红资白如下表所示,问两种方法治疗后的患者血红蛋白有无差异?
铁剂和饮食两种方法治疗后患者血红蛋白值(g/L)
铁剂治疗组
113
120
138
100
118
123
饮食治疗组
116
125
136
110
132
130
两种方法治疗后的患者血红蛋白无差异
两种方法治疗后的患者血红蛋白有差异
-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
xandy
7
truedifferenceinmeansisnotequalto0
meanofthedifferences
-3.375
p=>
0.05,不拒绝原假设,两种方法治疗后的患者血红蛋白无差异
5.〔习题5.4〕为研究国产四类新药阿卡波糖股嚢效果,某医院用40名Ⅱ型糖尿病病人进行同期随机对照实验。
试验者将这些病人随机等分到试验组〔阿卡波糖股嚢组〕和对照组〔拜唐苹股嚢组〕,分别测得试验开始前和8周后空腹血糖,算得空腹血糖下降值,如下所示。
能否认为国产四类新药阿卡波糖股嚢与拜唐苹股嚢对空腹血糖的降糖效果不同?
试验组与对照组空腹腔血糖下降值〔mmol/L)
试验组
(n1=20)
对照组
(n2=20)
〔1〕检验试验组和对照组的的数据是否来自正态分布,采用正态性W检验方法〔见第3章〕、Kolmogorov-Smirnov检验方法和Pearson拟合优度χ2检验;
认为国产四类新药阿卡波糖股嚢与拜唐苹股嚢对空腹血糖的降糖效果不同
认为国产四类新药阿卡波糖股嚢与拜唐苹股嚢对空腹血糖的降糖效果相同
①正态性W检验方法
-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
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)
y
试验组p=>
0.05,对照组p=>
0.05,所以检验试验组和对照组的的数据是来自正态分布
②Kolmogorov-Smirnov检验方法
ks.test(x,"
pnorm"
mean(x),sd(x))
One-sampleKolmogorov-Smirnovtest
two-sided
Warningmessage:
Inks.test(x,"
mean(x),sd(x)):
Kolmogorov-Smirnov检验里不应该有连结
ks.test(y,"
mean(y),sd(y))
Inks.test(y,"
mean(y),sd(y)):
③Pearson拟合优度χ2检验
0,4.50,4.60,2.50,6.00,-1.4)
A<
-table(cut(x,br=c(-6,-3,0,3,6,9)))
-pnorm(c(-3,0,3,6,9),mean(x),sd(x))
p
-c(p[1],p[2]-p[1],p[3]-p[2],p[4]-p[3],1-p[4])
chisq.test(A,p=p)
Chi-squaredtestforgivenprobabilities
A
Inchisq.test(A,p=p):
Chi-squared近似算法有可能不准
B<
-table(cut(y,br=c(-2,1,2,4,7)))
-pnorm(c(-2,1,2,4,7),mean(y),sd(y))
-c(p[1],p[2]-p[1],p[3]-p[2],1-p[3])
chisq.test(B,p=p)
B
Inchisq.test(B,p=p):
试验组的p=>
0.05,对照组的p=>
0.05,因此试验组和对照组都服从正态分布
〔2〕用t检验两组数据均值是否有差异,分别用方差相同模型、方差不同模型和成对t检验模型;
两组数据均值没有差异
两组数据均值是有差异
①方差相同模型
t.test(x,y,var.equal=TRUE)
TwoSamplet-test
meanofxmeanofy
0.05,不拒绝原假设,两组数据均值没有差异