统计学-R语言实例（部分

最新推荐文章于 2022-07-06 23:51:29 发布

原创

最新推荐文章于 2022-07-06 23:51:29 发布 · 3.3k 阅读

28 ·

CC 4.0 BY-SA版权

文章标签：

#数据分析

1t检验
提前，满足正态性
1. 1单样本t检验
已知某水样中含碳酸钙的真值为20.7mg/L，现用某法重复测定该水样12次，碳酸钙的含量分别为…问该法测定碳酸钙含量所得的均值与诊治有无显著差异？
12次测定中，碳酸钙的含量分别为20.99,20.41,20.10,20.00,20.91,22.60,20.99,20.42,20.90,22.99,23.12,20.89
输入：

> x <- c(20.99,20.41,20.10,20.00,20.91,22.60,20.99,20.42,20.90,22.99,23.12,20.89)
> t.test(x,alternative = 'two.sided',mu = 20.7)

输出：

	One Sample t-test

data:  x
t = 1.5665, df = 11, p-value = 0.1455
alternative hypothesis: true mean is not equal to 20.7
95 percent confidence interval:
 20.50020 21.88646
sample estimates:
mean of x 
 21.19333

结果解读：p-value = 0.1455>0.05，所以该法测定碳酸钙含量所得的均值与诊治没有显著差异
在这里插入图片描述

1.2 配对样本t检验

判断简便法和常规法测定尿铅含量的差别有无统计意义，对12份人尿同时用两种方法进行测定，所得结果如下表所示，请分析两种测定方法的测量结果是否不同？
法一的测量结果：2.41,2.90,2.75,2.23,3.67,4.49,5.16,5.45,2.06,1.64,1.06,0.77
法二的测量结果：2.80,3.04,1.88,3.43,3.81,4.00,4.44,5.41,1.24,1.83,1.45,0.92
输入：

> x <- c(2.41,2.90,2.75,2.23,3.67,4.49,5.16,5.45,2.06,1.64,1.06,0.77)
> y <- c(2.80,3.04,1.88,3.43,3.81,4.00,4.44,5.41,1.24,1.83,1.45,0.92)
> t.test(x,y,paired=T)

或

> data=data.frame(x=c(2.41,2.90,2.75,2.23,3.67,4.49,5.16,5.45,2.06,1.64,1.06,0.77,2.80,3.04,1.88,3.43,3.81,4.00,
4.44,5.41,1.24,1.83,1.45,0.92),g=factor(rep(1:2,c(12,12))))
>  t.