首 页IT知识库翔宇问吧收藏内容
当前位置:算搜IT乐园IT知识库R语言

使用R语言对两独立样本进行t检验

减小字体增大字体作者:suansuo.cn  来源:本站原创  发布时间:2019-03-06 18:25:54

1、什么是两独立样本t检验

根据样本数据对两个样本来自的两个独立总体的均值是否有显著差异进行判断。

2、使用的前提条件

(1)两个样本应该是相互独立的;

(2)样本来自的两个总体应该服从正态分布。

3、问题描述

设总体x1服从正态分布N(μ1,σ12),总体x2服从正态分布N(μ222),分别从两个总体中抽取两组样本:(X11,x12,...,x1m)和(x21,x22,...,x2n),且两个样本相互独立。则检验μ1和μ2是否有显著性差异。

检验方法:

(1)提出假设

H0:μ1 = μ2

H1:μ1 ≠ μ2

(2)建立检验统计量

一般分两种情况,方差齐而未知;另外一种情况两总体方差不齐而未知。这两种情况计算统计量t的方法不同。所以在实际中应判断一下两独立总体的样本方差是否相等,并选择不同的方法来判断。

(3)计算检验统计量的观测值和p值

(4)在给定的显著性水平上做出判断决策

4、R语言中进行两个独立样本t检验的方法

在R语言中使用t.test()函数可以实现对两个独立样本的t检验。

t.test()函数的使用形式可以参见本站中的文章《使用R语言进行单样本t检验》中关于t.test()函数的介绍。

5、R语言中进行两个独立样本t检验举例

(1)使用甲乙两台机床来加工同种零件,两种机床加工零件的尺寸服从正态分布,且方差相同,为检验甲乙两台机床加工的零件尺寸均值是否一致(α=0.05),从两种机床加工的零件中分别抽取若干零件测得的其尺寸如下:

甲:20.5,19.8,19.7,20.4,20.1,20.0,19.0,19.9

乙:20.7,19.8,19.5,20.8,20.4,19.6,20.2

在R中编写程序如下:

#甲机床加工数据
x1<-c(20.5,19.8,19.7,20.4,20.1,20.0,19.0,19.9)
#乙机床加工数据
x2<-c(20.7,19.8,19.5,20.8,20.4,19.6,20.2)
#进行两样本t检验,指定方差相等
t.test(x1,x2,var.equal=T)

 计算结果如下图所示:

R语言进行双样本T检验

 从计算结果来看,t=-0.85485,p-value=0.4081,由于计算的p值>0.05,则不能拒绝原假设,即认为两台机床加工零件的尺寸均值没有显著性区别,也即认为均值相同。同时结果也给出了两台机床加工零件尺寸样本数据的均值分别为:19.925和20.14286。

由于题目中给出了两台机床的方差相同,所以这里在使用t.test()函数进行检验时,指定了参数var.equal=T.

(2)为了检验两种供热设备的能耗是否相同。测量两种设备(电动气闸和热活化气闸)的能耗数据如下 。您要确定是否有证据证明这两种设备之间的差值不为零,以比较出这两种设备的功效。
数据如下:

热活化气闸数据:

13.07,7.6,3.2,9.28,9.73,11.73,9.67,10.76,11.05,17.63,15.58,12.53,11.87,14.19,6.84,11.89,7.41,7.42,
10.83,9.44,12.94,13.15,11.69,7.73,11.94,13.62,17.07,14.66,9.56,12.37,8.33,8.67,11.27,11.67,9.37,8.93,
8.41,12.85,5.27,10.02,7.87,11.82,14.42,13.69,10.77,15.26,14.53,6.84,10.92,13.05

电动气闸数据:

8.25,9.66,8.33,8.82,12.06,9.67,17.51,10.79,13.59,7.99,12.64,14.42,9.25,7.79,
11.29,10.26,9.46,14.77,7.21,4.29,9.81,8.41,6.78,16.3,9.01,11.41,12.37,13.28,7.24,
10.55,13.89,10.72,9.22,10.61,10.04,10.2,20.55,11.75,7.08,5.5

在R中编写程序如下:

#热活化气闸数据
x1<-c(13.07,7.6,3.2,9.28,9.73,11.73,9.67,10.76,11.05,17.63,15.58,12.53,11.87,
14.19,6.84,11.89,7.41,7.42,10.83,9.44,12.94,13.15,11.69,7.73,11.94,13.62,17.07,
14.66,9.56,12.37,8.33,8.67,11.27,11.67,9.37,8.93,8.41,12.85,5.27,10.02,7.87,
11.82,14.42,13.69,10.77,15.26,14.53,6.84,10.92,13.05)
#电动气闸数据
x2<-c(8.25,9.66,8.33,8.82,12.06,9.67,17.51,10.79,13.59,7.99,12.64,14.42,9.25,
7.79,11.29,10.26,9.46,14.77,7.21,4.29,9.81,8.41,6.78,16.3,9.01,11.41,12.37,
13.28,7.24,10.55,13.89,10.72,9.22,10.61,10.04,10.2,20.55,11.75,7.08,5.5)

#若方差相等
t.test(x1,x2,var.equal=T,conf.level=0.95)

#若假定方差不等
t.test(x1,x2,var.equal=F)

假定方差相等时的运行结果如下:

 R语言中进行双样本t检验的例子

假定方差不等时的运行结果如下:

R语言中进行双样本t检验举例说明

从上面两个运行结果来看,p值都大于0.05,则可以认为两种供暖设备的能耗在统计上没有显著区别。

在实际中可以使用R语言提供的相关函数(如var.test())对两组数据的方差齐次性进行检验,以便选择使用不同的方法进行检验。

微信搜索“优雅的代码”关注本站的公众号,或直接使用微信扫描下面二维码关注本站公众号,以获取最新内容。

个人成长离不开各位的关注,你的关注就是我继续前行的动力。

知识评论评论内容只代表网友观点,与本站立场无关!

   评论摘要(共 0 条,得分 0 分,平均 0 分) 查看完整评论
愿您的精彩评论引起共鸣,带来思考和价值。
用户名:查看更多评论
分 值:100分 90分 80分 70分 60分 40分 20分
内 容:
验证码:
关于本站 | 网站帮助 | 广告合作 | 网站声明 | 友情连接 | 网站地图
本站部分内容来自互联网,如有侵权,请来信告之,谢谢!
Copyright © 2007-2021 suansuo.cn. All Rights Reserved.
大盘指数基金买哪个好,最近买的基金一直在跌。就几天之内跌了快10%了,应该止损吗?_册肮股票网b股是什么意思,股票中的A.B股是什么意思?_厕仄股票网华宝证券科创板知识问卷答案,如何开通低佣金股票证券账户?_侧根股票网资本市场主要包括哪些市场,投行部与资本市场部哪个部门更好?_侧给股票网速动比率的计算方法充分考虑了,小白头一回参加券商研究所面试,我该从哪些方面准备?PS:面试时候要求带上写的小型行研报告。_侧谮股票网双龙战法指标下载,一个投资的故事告诉你:死啃七句口诀,穷人都能变富翁!_策髦股票网什么属于资本市场产品,资本市场有什么职能?中国资本市场学院顾丽诗,热点文章推荐邵和顾txt的小说《月亮进入我的相思局》分享全文阅读流动负债怎么计算,估算负荷电流准确吗?该如何估算?如何准确计算?中级财务管理视频课程,中级财务管理:第二章 财务管理基础(4)长期资产负债率计算,资产负债率高如何纳税筹划?现金债务总额比计算公式,企业财务绩效定量评价指标计算公式个人负债率是什么意思,负债率过高?别慌,玩卡区教你负债高如何办理大额信用卡?光大证券科创板开户答题,新手炒股入门怎么开户?科创板开通试题答案,2019年一季度,VC/PE成绩单出炉:募资低潮持续,科创板落地迎来新生科创板交易特别规定由谁制定,在科创板交易制度新规下,我们如何做好二级市场的营销工作?开通科创板知识测试题答案,聚源科创板课堂(三):正式开板!已经拥有姓名的企业都是哪路“豪杰”?东方财富科创板知识测评的答案,【干货】科创板首批9家受理企业基本面深度解读征信一直有逾期会有什么后果?,蚂蚁借呗逾期一天会影响征信记录吗?负债率百分之90,如何评价2021年3月19日A股大跌?计算流动比率,资本金比例的计算基数,究竟是总投资(包含建设期利息)还是建设投资(不包含建设期利息)?个人负债率怎么计算,001:如何区分有息负债和无息负债?负债率低的公司,怎么求医药行业平均资产负债率?控股和股份的区别,如何认识控股股东与实际控制人的不同?股东频繁减持,频繁减持,股价较高点狂跌近八成,赛诺医疗“两条腿”走路能行吗?新股上市流通股比例,属于「知乎官方指南」的问答有哪些?把上市公司股票全买了,买两只股票好还是八只以上好?上市公司流通股和非流通股的比例,有哪些必须掌握的股市术语?股票今日行情分析最新,20180521股市分析及20180522股市预测分析?为什么股票发行价格不能低于面值,企业为什么不调整可转债的发行条款从而使可转债的理论价值低于面值100呢?