R语言实验七Word文件下载.docx
- 文档编号:19083532
- 上传时间:2023-01-03
- 格式:DOCX
- 页数:7
- 大小:28.78KB
R语言实验七Word文件下载.docx
《R语言实验七Word文件下载.docx》由会员分享,可在线阅读,更多相关《R语言实验七Word文件下载.docx(7页珍藏版)》请在冰豆网上搜索。
)1.自行完成教材第五章得例题。
2.(习题5、1)正常男子血小板计数均值为225109/L,今测得20名男性油漆作业工人得血小板计数值(单位:
109/L)220188162230145160238188247113126245164231256183190158224175问油漆工人得血小板计数与正常成年男子有无差异?
解:
提出假设:
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)结论:
OneSamplettestdata:
xt=3、4783,df=19,p-value=0、002516alternativehypothesis:
truemeanisnotequalto22595percentconfidenceinterval:
172、3827211、9173sampleestimates:
meanofx192、15P=0、0025160、05,拒绝原假设,认为油漆工人得血小板计数与正常成年男子有
差异3.(习题5、2)已知某种灯泡寿命服从正态分布,在某星期所生产得该灯泡中随机抽取10只,测得其寿命(单位:
小时)为1067919119678511269369181156920948求这个星期生产出得灯泡能使用1000小时以上得概率。
源代码及运行结果:
(复制到此处,不需要截图)〉x〈c(1067,919,1196,785,1126,936,918,1156,920,948)ppnorm(1000,mean(x),sd(x))1p[1]0、4912059结论:
这个星期生产出得灯泡能使用1000小时以上得概率为0、49120594.(习题5、3)为研究某铁剂治疗与饮食治疗營养性缺铁性贫血得效果,将16名患者按年龄、体重、病程与病情相近得原则配成8对,分别使用饮食疗法与补充铁剂治疗得方法,3个月后测得两种患者血红资白如下表所示,问两种方法治疗后得患者血红蛋白有无差异?
铁剂与饮食两种方法治疗后患者血红蛋白值(g/L)铁剂治疗组11饮食治疗组6110132130110解:
提出假设:
H0:
两种方法治疗后得患者血红蛋白无差异H1:
两种方法治疗后得患者血红蛋白有差异源代码及运行结果:
(复制到此处,不需要截图)〉x〈c(113,120,138,120,100,118,138,123)yc(138,116,125,136,110,132,130,110)t、test(x,y,paired=T)Pairedt-testdata:
xandyt=-0、65127,df=7,pvalue=0、5357alternativehypothesis:
truedifferenceinmeansisnotequalto095percentconfidenceinterval:
-15、6288918、878891sampleestimates:
meanofthedifferences-3、375结论:
p=0、53570、05,不拒绝原假设,两种方法治疗后得患者血红蛋白无差异5.(习题5、4)为研究国产四类新药阿卡波糖股嚢效果,某医院用40名Ⅱ型糖尿病病人进行同期随机对照实验.试验者将这些病人随机等分到试验组(阿卡波糖股嚢组)与对照组(拜唐苹股嚢组),分别测得试验开始前与8周后空腹血糖,算得空腹血糖下降值,如下所示。
能否认为国产四类新药阿卡波糖股嚢与拜唐苹股嚢对空腹血糖得降糖效果不同?
试验组与对照组空腹腔血糖下降值(mmol/L)
试验组-0、705、602、002、800、703、504、005、807、100、50(n1=20)2、501、601、703、000、404、504、602、506、001、40对照组3、706、505、005、200、800、200、603、406、60-1、10(n2=20)6、003、802、001、602、002、201、203、101、702、00
(1)检验试验组与对照组得得数据就是否来自正态分布,采用正态性W检验方法(见第3章)、Kolmogorov-Smirnov检验方法与Pearson拟合优度2检验;
H0:
认为国产四类新药阿卡波糖股嚢与拜唐苹股嚢对空腹血糖得降糖效果不同H1:
认为国产四类新药阿卡波糖股嚢与拜唐苹股嚢对空腹血糖得降糖效果相同①正态性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-Wilknormalitytestdata:
xW=0、9699,pvalue=0、7527>yc(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)ShapiroWilknormalitytestdata:
yW=0、97098,pvalue=0、7754结论:
试验组p=0、7527>0、05,对照组p=0、7754〉0、05,所以检验试验组与对照组得得数据就是来自正态分布②Kolmogorov-Smirnov检验方法源代码及运行结果:
(复制到此处,不需要截图)>ks、test(x,pnorm,mean(x),sd(x))One-sampleKolmogorov-Smirnovtestdata:
xD=0、10652,p-value=0、9771alternativehypothesis:
two-sidedWarningmessage:
Inks、test(x,pnorm,mean(x),sd(x)):
Kolmogorov-Smirnov检验里不应该有连结〉ks、test(y,pnorm,mean(y),sd(y))OnesampleKolmogorov-Smirnovtestdata:
yD=0、11969,p-value=0、9368alternativehypothesis:
two-sided
Warningmessage:
Inks、test(y,"pnorm,mean(y),sd(y)):
Kolmogorov-Smirnov检验里不应该有连结结论:
试验组p=0、97710、05,对照组p=0、9368〉0、05,所以检验试验组与对照组得得数据就是来自正态分布③Pearson拟合优度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、60,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、98828138p<-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)Chisquaredtestforgivenprobabilitiesdata:
AX-squared=0、56387,df=4,p-value=0、967Warningmessage:
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、96465615p〈c(p[1],p[2]-p[1],p[3]-p[2],1p[3])〉p[1]0、028017520、222981540、147122890、60187805chisq、test(B,p=p)Chisquaredtestforgivenprobabilitiesdata:
BX-squared=28、087,df=3,p-value=3、483e06Warningmessage:
Inchisq、test(B,p=p):
Chi-squared近似算法有可能不准结论:
试验组得p=0、967〉0、05,对照组得p=3、483e060、05,因此试验组与对照组都服从正态分布(2)用t检验两组数据均值就是否有差异,分别用方差相同模型、方差不同模型与成对
t检验模型;
解:
两组数据均值没有差异H1:
两组数据均值就是有差异①方差相同模型源代码及运行结果:
(复制到此处,不需要截图)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)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)〉t、test(x,y,var、equal=TRUE)TwoSamplettestdata:
xandyt=-0、64187,df=38,pvalue=0、5248alternativehypothesis:
truedifferenceinmeansisnotequalto095percentconfidenceinterval:
-2、3261791、206179sampleestimates:
meanofxmeanofy2、0652、625结论:
p=0、52480、05,不拒绝原假设,两组数据均值没有差异②方差不同模型源代码及运行结果:
(复制到此处,不需要截图)>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)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)t、test(x,y)WelchTwoSamplet-testdata:
xandyt=-0、64187,df=36、086,pvalue=0、525alternativehypothesis:
truedifferenceinmeansisnotequalto095percentconfidenceinterval:
-2、329261、20926sampleestimates:
meanofxmeanofy2、0652、625结论:
p=0、5250、05,不拒绝原假设,两组数据均值没有差异③成对t检验模型源代码及运行结果:
(复制到此处,不需要截图)
xc(-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)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)〉t、test(x,y,paired=T)Pairedttestdata:
xandyt=-0、64644,df=19,pvalue=0、5257alternativehypothesis:
truedifferenceinmeansisnotequalto095percentconfidenceinterval:
-2、3731461、253146sampleestimates:
meanofthedifferences-0、56结论:
p=0、5257〉0、05,不拒绝原假设,两组数据均值没有差异(3)检验试验组与对照组得方差就是否相同。
H0:
试验组与对照组得方差相同H1:
试验组与对照组得方差不相同源代码及运行结果:
(复制到此处,不需要截图)〉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)〉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)>var、test(x,y)Ftesttoparetwovariancesdata:
xandyF=1、5984,numdf=19,denomdf=19,p-value=0、3153alternativehypothesis:
trueratioofvariancesisnotequalto195percentconfidenceinterval:
0、63265054、0381795sampleestimates:
ratioofvariances1、598361结论:
p=0、31530、05,不拒绝原假设,试验组与对照组得方差相同6.(习题5、5)为研究某种新药对抗凝血酶活力得影响,随机安排新药组病人12例,对照组病人10例,分别测定其抗凝血酶活力(单位:
mm3),其结果如下:
新药组:
126125136128123138142116110108115140对照组:
162172177170175152157159160162
试分析新药组与对照组病人得抗凝血酶活力有无差别(=0、05)。
(1)检验两组数据就是否服从正态分布;
(2)检验两组样本方差就是否相同;
(3)选择最合适得检验方法检验新药组与对照组病人得抗凝血酶活力有无差别。
(1)检验两组数据就是否服从正态分布①新药组数据提出假设:
H0:
新药组数据服从正态分布H1:
新药组数据不服从正态分布源代码及运行结果:
(复制到此处,不需要截图)x<-c(126,125,136,128,123,138,142,116,110,108,115,140)〉ks、test(x,pnorm,mean(x),sd(x))One-sampleKolmogorov-Smirnovtestdata:
xD=0、14644,pvalue=0、9266alternativehypothesis:
two-sided②对照组数据提出假设:
对照组数据服从正态分布H1:
对照组数据不服从正态分布源代码及运行结果:
(复制到此处,不需要截图)>y-c(162,172,177,170,175,152,157,159,160,162)〉ks、test(y,pnorm,mean(y),sd(y))One-sampleKolmogorovSmirnovtestdata:
yD=0、22216,pvalue=0、707alternativehypothesis:
twosidedWarningmessage:
Inks、test(y,"pnorm,mean(y),sd(y)):
KolmogorovSmirnov检验里不应该有连结
(2)检验两组样本方差就是否相同;
两组样本方差相同H1:
两组样本方差不相同源代码及运行结果:
(复制到此处,不需要截图)>x<c(126,125,136,128,123,138,142,116,110,108,115,140)〉y-c(162,172,177,170,175,152,157,159,160,162)〉var、test(x,y)Ftesttoparetwovariancesdata:
xandyF=1、9646,numdf=11,denomdf=9,p-value=0、32alternativehypothesis:
trueratioofvariancesisnotequalto195percentconfidenceinterval:
0、50219437、0488630sampleestimates:
ratioofvariances1、964622(3)选择最合适得检验方法检验新药组与对照组病人得抗凝血酶活力有无差别。
H0:
新药组与对照组病人得抗凝血酶活力无差别H1:
新药组与对照组病人得抗凝血酶活力有差别源代码及运行结果:
(复制到此处,不需要截图)〉x〈-c(126,125,136,128,123,138,142,116,110,108,115,140)yc(162,172,177,170,175,152,157,159,160,162)>t、test(x,y,var、equal=TRUE)TwoSamplet-testdata:
xandyt=8、8148,df=20,pvalue=2、524e08alternativehypothesis:
truedifferenceinmeansisnotequalto095percentconfidenceinterval:
-48、24975-29、78358sampleestimates:
meanofxmeanofy125、5833164、6000结论:
p=2、524e080、05,拒绝原假设,新药组与对照组病人得抗凝血酶活力有差别7.(习题5、6)项调查显示某城市老年人口比重为14、7%。
该市老年研究协会为了检验该项调查就是否可靠,随机抽选了400名居民,发现其中有57人就是老年人。
问调查结果就是否支持该市老年人口比重为14、7%得瞧法(=0、05).(提示,此题就是二项分布总体得检验)解:
p=p0=0、147H1:
pp0源代码及运行结果:
(复制到此处,不需要截图)>binom、test(57,400,p=0、147)Exactbinomialtestdata:
57and400numberofsuccesses=57,numberoftrials=400,pvalue=0、8876alternativehypothesis:
trueprobabilityofsuccessisnotequalto0、14795percentconfidenceinterval:
0、10974770、1806511sampleestimates:
probabilityofsuccess0、1425
结论:
P值=0、8876>0、05,不拒绝原假设,调查结果支持该市老年人口比重为14、7%得瞧法8.(习题5、7)作性别控制试验,经某种处理后,共有雏鸡328只,其中公雏150只,母雏178只,试问这种处理能否增加母雏得比例?
(性别比应为1:
1).解:
这种处理不能增加母雏得比例H1:
这种处理能增加母雏得比例源代码及运行结果:
(复制到此处,不需要截图)>binom、test(178,328,p=0、5,alternative=greater")Exactbinomialtestdata:
178and328numberofsuccesses=178,numberoftrials=328,p-value=0、06794alternativehypothesis:
trueprobabilityofsuccessisgreaterthan0、595percentconfidenceinterval:
0、49576161、0000000sampleestimates:
probabilityofsuccess0、5426829结论:
P值=0、06794〉0、05,不拒绝原假设,这种处理不能增加母雏得比例9.(习题5、8)Mendel用豌豆得两对相对性状进行杂交实验,黄色圆滑种子与绿色皱缩种子得豌豆杂交后,第二代根据自由组合规律,理论分离比为黄圆:
黄皱:
绿圆:
绿皱=(9/16):
(3/16):
(3/16):
(1/16)实际实验值为:
黄圆315粒、黄皱101粒、绿圆108粒、绿皱32粒,共556粒。
问此结果就是否符合自由组合规律得理论分离比?
H0:
符合自由组合规律得理论分离比H1:
不符合自由组合规律得理论分离比源代码及运行结果:
(复制到此处,不需要截图)chisq、test(c(315,101,108,32),p=c(9,3,3,1)/16)Chi-squaredtestforgivenprobabilitiesdata:
c(315,101,108,32)X-squared=0、47002,df=3,pvalue=0、9254结论:
P值=0、92540、05,接受原假设,符合自由组合规律得理论分离比10.(习题5、9)观察每分钟进入某商店得人数X,任取200分钟,所得数据如下:
顾客人数012345频数9268281110试分析,能否认为每分钟顾客数X服从Poisson分布(=0、1).解:
能认为每分钟顾客数X服从Poisson分布
H1:
不能认为每分钟顾客数X服从Poisson分布源代码及运行结果:
(复制到此处,不需要截图)X<-0:
5;
Y〈c(92,68,28,11,1,0)>q〈ppois(X,mean(rep(X,Y)));
n<-length(Y)>p-numeric(
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 语言 实验