R语言实验七Word格式.docx

1、油漆工人的血小板计数与正常成年男子有差异源代码及运行结果：（复制到此处，不需要截图）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.002516alternative hypothesis: true mean is not equal to 22595 percent confidenc

2、e interval: 172.3827 211.9173sample estimates:mean of x 192.15P=0.0025160.05,拒绝原假设，认为油漆工人的血小板计数与正常成年男子有差异1.（习题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-pnorm（1

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

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

5、.601.703.000.404.504.606.00-1.40对照组3.706.505.005.200.800.200.603.406.60-1.10（n2 = 20）3.801.602.201.203.10-2.00（1）检验试验组和对照组的的数据是否来自正态分布，采用正态性W检验方法（见第3章）、Kolmogorov-Smirnov检验方法和Pearson拟合优度2检验；认为国产四类新药阿卡波糖股嚢与拜唐苹股嚢对空腹血糖的降糖效果不同认为国产四类新药阿卡波糖股嚢与拜唐苹股嚢对空腹血糖的降糖效果相同正态性W检验方法-c（-0.70,-5.60,2.00,2.80,0.70,3.50,4.0

6、0,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-Wilk normality testW = 0.9699, p-value = 0.7527y0.05,对照组p=0.77540.05,所以检验试验组和对照组的的数据是来自正态分布Kolmogorov-Smirnov检验方法 ks.test（x,pnorm,mean（x）,sd（x） One-sample Kolmogorov-Smirnov testD = 0.10652, p-value = 0.977

7、1 two-sidedWarning message:In ks.test（x, , mean（x）, sd（x） : Kolmogorov - Smirnov检验里不应该有连结 ks.test（y,mean（y）,sd（y）D = 0.11969, p-value = 0.9368In ks.test（y, , mean（y）, sd（y） :试验组p=0.97710.05,对照组p=0.9368Pearson拟合优度2检验 A-table（cut（x,br=c（-6,-3,0,3,6,9）-pnorm（c（-3,0,3,6,9）,mean（x）,sd（x） p1 0.04894712 0.

8、24990009 0.62002288 0.90075856 0.98828138-c（p1,p2-p1,p3-p2,p4-p3,1-p4）1 0.04894712 0.20095298 0.37012278 0.28073568 0.09924144 chisq.test（A,p=p） Chi-squared test for given probabilities AX-squared = 0.56387, df = 4, p-value = 0.967In chisq.test（A, p = p） : Chi-squared近似算法有可能不准 B0.05,对照组的p= 3.483e-060.05,因此试验组和对照组都服从正态分布（2）用 t检验两组数据均值是否有差异，分别用方差相同模型、方差不同模型和成对t检验模型；两组数据均值没有差异两组数据均值是有差异方差相同模型（复制到此