练习R：qqnorm+t.test函数配对t检验

动手练习学得更快！

为比较两种方法对乳酸饮料中脂肪含悬测定结果是否不同，某人随机抽取了10份乳酸饮料制品。分别用脂肪酸水解法和哥特里一罗紫法测定其结果(单位:%)如下表，试问两法测定结果是否有显著差异？(a=0.01)。

同一受试对象分别接受两种不同的处理，这是配对数据，差异检验我们用配对t检验。

先来看一下glf法和zfsf法分别测定的脂肪含量平均水平。

library(psych)describe(mydata[,2:3])

gfl法均值水平0.8±0.18

zfsf法均值水平0.52±0.19

画个正态qq图看一下是否满足正态分布的要求：

diff <- mydata$glf-mydata$zfsfqqnorm(diff)qqline(diff,col="red")

嗯基本满足。继续t检验。

t.test(mydata$glf,mydata$zfsf,paired=T)

t检验结果表明，t=7.926，P<005，按显著性a=0.05水平，可认为两种方法对脂肪含量测定结果不同，glf法测定的脂肪含量较高（0.8±0.18）。

本文完

文/图=数据小兵