【R语言】t检验

一、基本介绍

t检验(t-test)是用于比较两个样本均值是否存在显著差异的一种统计方法。

t.test()函数的调用格式:

t.test(x, y=NULL, alternative=c("two.sided", "less", "greater"), mu=0, pair=FALSE, var.equal=FALSE, conf.level=0.95, ...)

参数详解:

  1. x:用于均值比较的数值向量。
  2. y:用于均值比较的数值向量。当y=NULL时,表示为单样本t检验。
  3. alternative:表示指定假设检验的备择检验类型。two.sided:双侧检验(默认);less:左侧检验(备择假设为样本均值小于mu);greater:右侧检验(备择假设为样本均值大于mu)。
  4. mu:一个对比用的数值,默认为0。当检验为单样本t检验时,表示x所代表的总体均值是否与mu中指定的值相等;当检验为双样本t检验时,表示x和y代表的总体均值之差是否与mu中指定的值相等。
  5. pair:逻辑值,指定是否为配对样本,默认为FALSE,即独立样本。
  6. var.equal:逻辑值,指定两个样本的方差是否相等(默认为TRUE,即假设方差相等;如果方差可能不相等,应设置为FALSE以执行Welch t检验)
  7. conf.level:置信区间的置信水平。

二、原假设与备择假设

在t检验中,原假设(null hypothesis,通常表示为H0​)和备择假设(alternative hypothesis,通常表示为H1​)是根据研究目的和所要检验的具体问题来设定的。

通常,原假设是希望被推翻的假设,而备择假设是研究者认为如果原假设被推翻则可能成立的假设。在进行t检验时,如果计算出的p值小于事先设定的显著性水平(如0.05)则可以拒绝原假设,接受备择假设(在统计上显著的意义上)。

1、单样本t检验

原假设:样本所来自的总体的均值等于某个已知值μ0​,即μ=μ0​。

备择假设:样本所来自的总体的均值不等于μ0​(双侧检验);或者小于μ0​(左侧检验);或者大于μ0​(右侧检验)。

2、独立样本t检验

原假设:样本所来自的总体的均值等于某个已知值μ0​,即μ=μ0​。

备择假设:样本所来自的总体的均值不等于μ0​(双侧检验);或者小于μ0​(左侧检验);或者大于μ0​(右侧检验)

3、非独立样本t检验

原假设:配对样本差值所来自的总体的均值等于0,即μd​=0(其中μd​表示配对样本差值的均值)。

备择假设:配对样本差值所来自的总体的均值不等于0(双侧检验);或者小于0(左侧检验,即认为第一种情况或测量值平均上小于第二种);或者大于0(右侧检验,即认为第一种情况或测量值平均上大于第二种)

三、单样本t检验

单样本t检验(One-Sample t-Test)用于比较一个样本的均值与已知的总体均值是否存在显著差异。

sample_data <- c(5.1, 4.8, 5.0, 5.3, 4.9)# 总体均值
mu <- 5.0# 执行单样本t检验
t.test(sample_data, mu = mu)

四、独立样本t检验

独立样本t检验(Independent Two-Sample t-Test)用于比较两个独立样本的均值是否存在显著差异。

group1 <- c(5.1, 4.8, 5.0, 5.3, 4.9)
group2 <- c(5.5, 5.2, 5.4, 5.6, 5.0)# 执行独立样本t检验, 假设两组数据之间无差异
t.test(group1, group2, alternative="two.sided", var.equal = T)

 使用R语言中内置的数据集InsectSprays,检验杀虫剂C与杀虫剂D的效果差异。

head(InsectSprays)
library(dplyr)
# 筛选杀虫剂C和D
data1 <- dplyr::filter(InsectSprays, spray %in% c("C", "D"))
# 独立样本t检验:假设使用杀虫剂C后的昆虫数目比使用杀虫剂D后的少
t.test(count ~ spray, data1, alternative="less", var.equal=T)

从以下结果可知,使用杀虫剂C后的昆虫数目(2.08)显著(P=0.0027<0.05)小于使用杀虫剂D后的昆虫数目(4.917),即杀虫剂C的效果比D好。 

五、非独立样本t检验

非独立样本t检验又叫配对样本t检验。

下面用R语言内置的sleep数据集做演示:

现假设两种安眠药对增加睡眠时长的效果一样。

head(sleep)
t.test(extra ~ group, data = sleep, pair=TRUE, alternative="two.sided")

从下面结果可知,两种安眠药平均增加睡眠时长相差1.58小时,而且p=0.002833,远小于0.05,因此可以拒绝原假设,即这两种安眠药的效果不一样。 

需要注意的是,因为是进行的双侧t检验,所以可以假设“两种安眠药对增加睡眠时长的效果一样”。如果要假设“第1种安眠药的效果低于第2种安眠药的效果”,需要使用左侧检验,如下所示:

t.test(extra ~ group, data = sleep, pair=T, alternative="less")

原假设:两种安眠药的效果没有显著差异,或者第1种安眠药的效果不小于第2种。

备择假设:第1种安眠药的效果低于第2种安眠药的效果。

从下面结果可知,第1种安眠药效果确实显著低于第2种安眠药效果(p=0.0014)  

六、方差齐性检验

方差齐性检验(Equality of Variance Test)是指检验两个或两个以上样本的方差是否相等的统计检验方法。如果两个样本的方差相等,说明它们的变异情况相似,可以认为它们来自同一总体或具有可比性;反之,如果方差不相等,则可能说明它们的变异情况不相似,来自不同的总体或不具有可比性。

在t检验中,方差齐性是一个重要的前提假设,因为如果两个样本的方差差异显著,那么t检验的结果可能会受到严重影响,导致结论不准确。

在t检验中,进行方差齐性检验的目的是确保两个样本具有可比性。因为t检验是用于比较两个样本均值是否存在显著差异的统计方法,如果两个样本的方差差异显著,那么它们的均值比较可能就没有意义。因此,在进行t检验之前,通常需要先进行方差齐性检验

常用检验方法是Bartlett检验和Levene检验。

1、Bartlett检验

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/16701.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

TDengine 产品由哪些组件构成

目 录 背景产品生态taosdtaosctaosAdaptertaosKeepertaosExplorertaosXtaosX Agent应用程序或第三方工具 背景 了解一个产品&#xff0c;最好从了解产品包括哪些内容开始&#xff0c;我这里整理了一份儿 TDegnine 产品包括有哪些组件&#xff0c;每个组件作用是什么的说明&a…

实现限制同一个账号最多只能在3个客户端(有电脑、手机等)登录(附关键源码)

如上图&#xff0c;我的百度网盘已登录设备列表&#xff0c;有一个手机&#xff0c;2个windows客户端。手机设备有型号、最后登录时间、IP等。windows客户端信息有最后登录时间、操作系统类型、IP地址等。这些具体是如何实现的&#xff1f;下面分别给出android APP中采集手机信…

使用 Docker 安装 Open WebUI 并集成 Ollama 的 DeepSeek 模型

文章目录 使用 Docker 安装 Open WebUI 并集成 Ollama 的 DeepSeek 模型前提条件1. 安装ollama2. 拉取deepseek的模型3. Open-WebUI 说明4. 启动容器文档的方法如下优化命令&#xff08;可选&#xff09;1. 增加了健康检查机制&#xff08;--health-cmd&#xff09;2. 使 WebUI…

Untiy3d 铰链、弹簧,特殊的物理关节

&#xff08;一&#xff09;铰链组件 1.创建一个立方体和角色胶囊 2.给角色胶囊挂在控制脚本和刚体 using System.Collections; using System.Collections.Generic; using UnityEngine;public class plyer : MonoBehaviour {// Start is called once before the first execut…

【NLP 21、实践 ③ 全切分函数切分句子】

当无数个自己离去&#xff0c;我便日益坦然 —— 25.2.9 一、jieba分词器 Jieba 是一款优秀的 Python 中文分词库&#xff0c;它支持多种分词模式&#xff0c;其中全切分方式会将句子中所有可能的词语都扫描出来。 1.原理 全切分方式会找出句子中所有可能的词语组合。对于一…

团结引擎 OpenHarmony 平台全面支持 UAAL,实现引擎能力嵌入原生应用

团结引擎1.4版本已于近日正式发布&#xff01;在这一版本中&#xff0c;OpenHarmony 平台迎来了一个具有里程碑意义的更新&#xff1a;全面支持 Used as a Library&#xff08;UAAL&#xff09;。UAAL 这一技术方案&#xff0c;具有将引擎嵌入原生应用的独特能力&#xff0c;其…

自己部署DeepSeek 助力 Vue 开发:打造丝滑的标签页(Tabs)

前言&#xff1a;哈喽&#xff0c;大家好&#xff0c;今天给大家分享一篇文章&#xff01;并提供具体代码帮助大家深入理解&#xff0c;彻底掌握&#xff01;创作不易&#xff0c;如果能帮助到大家或者给大家一些灵感和启发&#xff0c;欢迎收藏关注哦 &#x1f495; 目录 自己…

DeepSeek渣机部署编程用的模型,边缘设备部署模型

DeepSeek渣机部署编程用的模型&#xff0c;边缘设备部署模型 文章目录 DeepSeek渣机部署编程用的模型&#xff0c;边缘设备部署模型前言一、python代码二、构建一个简单的前端来接入接口2.读入数据 总结 前言 也许大家伙都想完成一些部署DeepSeek的东西&#xff0c;不过部署并…

VS2019打开《喜缺全书算法册》附带代码的方法兼述单元测试

下载地址(大量的题目和测试用例) 下载:地址一&#xff0c;几乎实时更新 GitCode下载。 下载地址二&#xff0c;不定期更新csdn打包下载 如果这两个链接打不开&#xff0c;可能是这两个资源处于审核状态&#xff0c;快则几分钟&#xff0c;慢则2天。 可以加本文末的&#xff31…

急停信号的含义

前言&#xff1a; 大家好&#xff0c;我是上位机马工&#xff0c;硕士毕业4年年入40万&#xff0c;目前在一家自动化公司担任软件经理&#xff0c;从事C#上位机软件开发8年以上&#xff01;我们在开发C#的运动控制程序的时候&#xff0c;一个必要的步骤就是确认设备按钮的急停…

小白学网络安全难吗?需要具备哪些条件?

作为一名零基础小白&#xff0c;想要转行IT学习一门新技术&#xff0c;且上手难度低、就业前景好、薪资待遇高、入行门槛低&#xff0c;网络安全是最值得的选择&#xff0c;掌握它之后你可以获得一份收入不错的工作。那么零基础学网络安全好学吗?以下是具体内容介绍。 首先&am…

IEEE期刊Word导出PDF注意事项

在系统上提交论文时候一般要求PDF文档&#xff0c;但是word直接转PDF可能存在一些问题&#xff1a; 部分图片不清晰。字体未嵌入PDF。间距发生了变化。字体发生了变化。一张图片显示不完全。 下面介绍word转PDF最稳妥的技巧以及如何实现全部字体的嵌入。 1. 操作流程 ① 另…

CEF132 编译指南 MacOS 篇 - depot_tools 安装与配置 (四)

1. 引言 在 CEF132&#xff08;Chromium Embedded Framework&#xff09;的编译过程中&#xff0c;depot_tools 扮演着举足轻重的角色。这套由 Chromium 项目精心打造的脚本和工具集&#xff0c;专门用于获取、管理和更新 Chromium 及其相关项目&#xff08;包括 CEF&#xff…

【C++高并发服务器WebServer】-16:UDP简单实现

本文目录 一、UDP通信流程二、UDP API2.1 sendto()2.2 recvfrom() 一、UDP通信流程 UDP通信的流程比较简单&#xff0c;下面这张图可以总结。 二、UDP API 2.1 sendto() UDP相关头文件如下。 #include <sys/types.h> #include <sys/socket.h>ssize_t sendto(…

k8s管理工具之lens

什么是lens Lens 是当前市场上最强大的K8S IDE。它是一个独立的单机应用&#xff0c;可以同时运行在macOS、Windows和Linux上。 作为K8S IDE&#xff0c;该有的它基本都有了&#xff01; 集群管理 导入已有集群 首先&#xff0c;你需要在 Lens 中添加你的 Kubernetes 集群。点…

SynthDetoxM - 现代LLM是少样本的并行去毒化数据标注器

SynthDetoxM: Modern LLMs are Few-Shot Parallel Detoxification Data Annotators https://arxiv.org/html/2502.06394v1 1. 主要内容 这篇论文提出了一个 用于生成多语言平行去毒化数据的管道&#xff0c;并介绍了SynthDetoxM&#xff0c;一个包含16,000个高质量去毒化句子对…

云服务器流量包用尽(中病毒)

1. 情况 腾讯云提示我账号欠费了&#xff0c;服务器存在恶意文件。。。 一看&#xff0c;流量包用尽超额了&#xff0c;CPU直接爆了。 用iftop监测一下网络流量。可以看到向多个IP发送了大量的流量。 看来是中病毒了&#xff0c;被当成 “肉鸡”&#xff0c;纳入“僵尸网络”…

RK3588视觉控制器与AI 算法:开启工业视觉检测新境界

在实际应用中&#xff0c;工业相机拍摄产品的图像&#xff0c;RK3588 迅速接收并进行预处理。AI 算法随即对图像进行深入分析&#xff0c;提取特征并与预设的标准进行对比&#xff0c;从而准确判断是否存在缺陷。 例如&#xff0c;在电子元件生产线上&#xff0c;RK3588 和 AI…

android的ViewModel这个类就是业务逻辑层吗

android的ViewModel这个类就是业务逻辑层吗&#xff1f; 相似&#xff1a;业务逻辑代码应该放在ViewModel这个类吗&#xff1f; 嗯&#xff0c;我现在在学习Android架构组件&#xff0c;特别是ViewModel。用户问ViewModel是否就是业务逻辑层&#xff0c;我需要仔细思考这个问题…

Gui-Guider1.8.1 数字时钟控件找不到定义,无法编译

我们在Gui-Guider中使用的一些控件&#xff0c;生成后会发现在LVGL源码中找不到该控件的定义&#xff0c;这时因为Gui-Guider中的一些控件是其自己编写的而不是LVGL提供的&#xff0c;那么我们该如何应用呢&#xff1f;这里拿Digital Clock数字时钟控件举例&#xff1a; 这里我…