R语言实验七Word格式.docx

1、利用QQ输入法的截屏工具。点击 QQ输入法工具条最右边的“扳手”图标 ZEI，选择 其中的“截屏”工具。）1. 自行完成教材第五章的例题。2. （习题5.1）正常男子血小板计数均值为 225 X 109/L,今测得20名男性油漆作业工人的血小板计数值（单位： 109/L）220 188 162 230 145 160 238 188 247 113126 245 164 231 256 183 190 158 224 175问油漆工人的血小板计数与正常成年男子有无差异？解：提出假设：H0：油漆工人的血小板计数与正常成年男子无差异H1：油漆工人的血小板计数与正常成年男子有差异 源代码及运行结果：

2、（复制到此处，不需要截图）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）结论：One Sample t-testdata: xt = -3.4783, df = 19, p-value = 0.002516alter native hypothesis: true mean is not equal to 22595 perce nt con fide nee in terval:172.3827 211.9

3、173sample estimates:mean of x192.15P=0.0025160.05,拒绝原假设，认为油漆工人的血小板计数与正常成年男子有差异3. （习题5.2）已知某种灯泡寿命服从正态分布，在某星期所生产的该灯泡中随机抽 取10只，测得其寿命（单位：小时）为1067 919 1196 785 1126 936 918 1156 920 948求这个星期生产出的灯泡能使用 1000小时以上的概率。源代码及运行结果： x-c（1067, 919, 1196, 785, 1126, 936, 918, 1156, 920, 948） p-p no rm（1000,mea n（ x）,

4、sd（x） 1-p1 0.4912059这个星期生产出的灯泡能使用 1000小时以上的概率为 0.49120594. （习题5.3）为研究某铁剂治疗和饮食治疗營养性缺铁性贫血的效果， 将16名患者按年龄、体重、病程和病情相近的原则配成 8对，分别使用饮食疗法和补充铁剂治疗的方法，3个月后测得两种患者血红资白如下表所示，问两种方法治疗后的患者 血红蛋白有无差异？铁剂和饮食两种方法治疗后患者血红蛋白值 （g/L）铁剂治疗组113 120 138 120 100 118 138 123饮食治疗组138 116 125 136 110 132 130 110两种方法治疗后的患者血红蛋白无差异两种方法治

5、疗后的患者血红蛋白有差异 源代码及运行结果：-c（113,120,138,120,100,118,138,123） y0.05,不拒绝原假设，两种方法治疗后的患者血红蛋白无差异5.（习题5.4）为研究国产四类新药阿卡波糖股囊效果，某医院用 40名H型糖尿病病人进行同期随机对照实验。试验者将这些病人随机等分到试验组 （阿卡波糖股囊组）和对照组（拜唐苹股囊组），分别测得试验开始前和 8周后空腹血糖，算得空 腹血糖下降值，如下所示。能否认为国产四类新药阿卡波糖股囊与拜唐苹股囊对空 腹血糖的降糖效果不同？试验组与对照组空腹腔血糖下降值（ mmol/L）试验组-0.70-5.602.002.800.70

6、3.504.005.807.10-0.50（n1 = 20）:.50-1.601.703.000.404.504.602.506.00-1.40对照组3.706.505.005.200.800.200.603.406.60-1.10（n2 = 20） 6）.003.801.602.201.203.10-2.00（1）检验试验组和对照组的的数据是否来自正态分布， 采用正态性 W检验方法（见第3章）、Kolmogorov-Smirnov 检验方法和 Pearson拟合优度 2检验；认为国产四类新药阿卡波糖股囊与拜唐苹股囊对空腹血糖的降糖效果不同H1:认为国产四类新药阿卡波糖股囊与拜唐苹股囊对空腹血

7、糖的降糖效果相同1 正态性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.6 0,2.50,6.00,-1.4） shapiro.test（x）Shapiro-Wilk n ormality testW = 0.9699, p-value = 0.7527y0.05,对照组p=0.77540.05,所以检验试验组和对照组的的数据是来自正态分布2 Kolmogorov-Smirnov 检验方法 ks.test（x,p norm,mea n

8、（ x）,sd（x）On e-sample Kolmogorov-Smir nov testD = 0.10652, p-value = 0.9771 two-sidedWarning message:In ks.test（x, pno rm, mean（ x）, sd（x）:Kolmogorov - Smirnov 检验里不应该有连结 ks.test（y,p no rm,mea n（ y）,sd（y）D = 0.11969, p-value = 0.9368In ks.test（y, , mean（ y）, sd（y）:试验组p=0.97710.05,对照组p=0.93683 Pearson

9、拟合优度2检验-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 norm（c（-3,0,3,6,9）,mea n（x）,sd（x） p1 0.04894712 0.24990009 0.62002288 0.90075856 0.98828138-c（p1,p2-p1,p3-p2,p4-p3,1-p4）1 0.04894712 0.20095298 0.370

10、12278 0.28073568 0.09924144 chisq.test（A,p=p）Chi-squared test for give n probabilities AX-squared = 0.56387, df = 4, p-value = 0.967In chisq.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）1 0.02801752 0.25099907 0.39812195 0.71500447 0.96465615-