推断统计(独立样本t检验)

 

这里我们是采用假设检验中的独立样本t 检验来比较两个独立正态总体均值之间是否存在显著性差异,以比较城市与农村孩子的心理素质是否有显著差异为例 。

 

 

 这里我们首先是假设城市孩子与农村孩子心理素质无显著差异,但是此时方差是否齐性是未知的,所以首先对方差是否齐性做一个判定,然后根据方差是否齐性的不同结果做出不同的决策。

'''step1 调用包'''
from  scipy.stats import ttest_ind, levene
import numpy as np'''step2 读(导)入数据'''
x = [4.75,6.40,2.62,3.44,6.50,5.20,5.60,3.80,4.30,5.78,3.76,4.15]
y = [2.38,2.60,2.10,1.80,1.90,3.65,2.30,3.80,4.60,4.85,5.80,4.25,4.22,3.84]
'''Step3 前提性检验:方差齐性检验	'''	
a =  levene(x,y)
print('方差齐性检验结果:')
print('\n T统计量的值为:',np.round(a.statistic,3))
print('\n P值为:',np.round(a.pvalue,3))
print('----------------')'''Step4 独立样本T检验,默认方差齐性'''
res = ttest_ind(x, y)
print('方差齐性情况下,均值比较结果(P值):\n',np.round(res.pvalue,3))
print('----------------')
'''Step4' 如果方差不齐性,则equal_var=False'''
res2 =ttest_ind(x,y,equal_var=False)
print('方差不具有齐性情况下,均值比较结果(P值):\n',np.round(res2.pvalue,3))

此时我们首先需要将所要用到的包从 scipy 导入进来,将要用于检验的数据读取进来,

`a = levene(x, y)`是进行莱文氏检验(Levene's test)的代码。莱文氏检验是检验两个或多个独立样本的方差是否相等的统计学方法。其中:- x和y是两个独立样本的数据,可以是数组或dataframe。- a是检验的统计量和p值,通常会以一个tuple的形式返回。- 如果p值小于显著性水平(如0.05),则拒绝原假设,认为两个样本的方差不相等。- 如果p值大于显著性水平,则不能拒绝原假设,认为两个样本的方差相等。这个检验经常在进行t检验或方差分析前,用来检查方差齐性的假设是否成立。所以莱文氏检验是一个重要的假设检验方法,用来检验方差齐性这一前提条件。

在独立样本t检验中,T统计量的作用是:

1. 用来检验两个独立样本的均值是否有显著差异。

2. T统计量的计算公式为:T = (样本1均值 - 样本2均值) / (两样本均值的标准误)其中,标准误反映了样本均值的误差范围。

3. 检验时,将计算所得的T统计量与理论T分布进行比较,以得到p值。

4. 如果p值小于显著性水平α(如0.05),则认为两样本均值之间存在显著差异,即拒绝原假设。

5. 如果p值大于显著性水平α,则不能拒绝原假设,认为两样本均值之间无显著差异。

6. 所以T统计量反映了两样本差异的程度,p值反映了这种差异是否显著。

7. 通过T统计量和p值,可以对两样本均值间的差异进行推断,这就是T统计量在独立样本t检验中的作用。

此时我们可以看到根据方差齐性检验结果,p值为 0.906,接受原假设方差齐性,此时根据方差是否齐性检验的结果判断此时方差齐性,在方差齐性的条件下比较均值是否有显著性差异,此时得到的p值为0.017小于0.05,所以此时拒绝原假设,认为有显著差异,即城市和农村孩子的心理素质有显著差异。

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

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

相关文章

题目:2566.替换一个数字后的最大差值

​​题目来源: leetcode题目,网址:2566. 替换一个数字后的最大差值 - 力扣(LeetCode) 解题思路: 将从左到右第一个非 9 数字全部修改为 9 以得到最大值。将从左到右第一个数字全部修改为 0 以得到最小值&a…

c#实现命令模式

下面是一个使用C#实现命令模式的示例代码: using System; using System.Collections.Generic;// 命令接口 public interface ICommand {void Execute();void Undo(); }// 具体命令:打开文件 public class OpenFileCommand : ICommand {private FileMana…

【MySQL】MySQL不走索引的情况分析

未建立索引 当数据表没有设计相关索引时,查询会扫描全表。 create table test_temp (test_id int auto_incrementprimary key,field_1 varchar(20) null,field_2 varchar(20) null,field_3 bigint null,create_date date null );expl…

ffmpeg命令行是如何打开vf_scale滤镜的

前言 在ffmpeg命令行中,ffmpeg -i test -pix_fmt rgb24 test.rgb,会自动打开ff_vf_scale滤镜,本章主要追踪这个流程。 通过gdb可以发现其基本调用栈如下: 可以看到,query_formats()中创建的v…

Unity框架学习--2

接上文 IOC 容器是一个很方便的模块管理工具。 除了可以用来注册和获取模块,IOC 容器一般还会有一个隐藏的功能,即: 注册接口模块 抽象-实现 这种形式注册和获取对象的方式是符合依赖倒置原则的。 依赖倒置原则(Dependence I…

maven install

maven install maven 的 install 命令,当我们的一个 maven 模块想要依赖其他目录下的模块时,直接添加会找不到对应的模块,只需要找到需要引入的模块,执行 install 命令,就会将该模块放入本地仓库,就可以进…

Linux tar包安装 Prometheus 和 Grafana(知识点:systemd Unit/重定向)

0. 介绍 用tar包的方式安装 Prometheus 和 Grafana Prometheus:开源的监控方案Grafana:将Prometheus的数据可视化平台 Prometheus已经有了查询功能为什么还需要grafana呢?Prometheus基于promQL这一SQL方言,有一定门槛!Grafana基于浏览器的操作与可视化图表大大降低了理解难…

Vue3 setup tsx 子组件向父组件传值 emit

需求:Vue3 setup 父组件向子组件传值,子组件接收父组件传入的值;子组件向父组件传值,父组件接收的子组件传递的值。 父组件:parent.tsx: import { defineComponent, ref, reactive } from vue; import To…

【Android】okhttp爆java.lang.IllegalStateException: closed的解决方法

解决 java.lang.IllegalStateException: closed异常通常是由于OkHttp中的Response对象在调用response.body().string()后被关闭而导致的。 在代码中,在onResponse()方法中如果两次调用了response.body().string(),每次调用都会消耗掉响应体并关闭Respo…

如何优化PHP Smarty模板的性能?

Smarty模板是一种非常强大的模板引擎,但是如果不正确地使用,可能会导致你的网站慢得像一只树懒! 那么,如何优化Smarty模板的性能呢? 减少Smarty对象的创建 你可能会在代码中多次创建Smarty对象。但是,每次…

Server - 文字转语音 (Text to Speech) 的在线服务 TTSMaker

欢迎关注我的CSDN:https://spike.blog.csdn.net/ 本文地址:https://spike.blog.csdn.net/article/details/132287193 TTSMaker 是一款免费的文本转语音工具,提供语音合成服务,支持多种语言,包括英语、法语、德语、西班…

什么是冒烟测试?

冒烟测试,刚进公司就接触到了。只是刚开始一直没有体会到冒烟的含义和精髓,一直以为是冒烟测试就是把待测产品的主要功能测试一下就行了。后面回想一下,不是那么回事的。 冒烟测试源自硬件行业,对一个硬件或者硬件组件改动后&…

Exams/ece241 2013 q4

蓄水池问题 S3 S2 S1 例如:000 代表 无水 ,需要使FR3, FR2, FR1 都打开(111) S3 S2 S1 FR3 FR2 FR1 000 111 001 011 011 001 111 000 fr代表水变深为…

快手商品详情数据API 抓取快手商品价格、销量、库存、sku信息

快手商品详情数据API是用来获取快手商品详情页数据的接口,请求参数为商品ID,这是每个商品唯一性的标识。返回参数有商品标题、商品标题、商品简介、价格、掌柜昵称、库存、宝贝链接、宝贝图片、商品SKU等。 接口名称:item_get 公共参数 名…

【PostgreSQL的CLOG解析】

同样还是这张图,之前发过shared_buffer和os cache、wal buffer和work mem的文章,今天的主题是图中的clog,即 commit log,PostgreSQL10之前放在数据库目录的pg_clog下面。PostgreSQL10之后修更名为xact,数据目录变更为pg_xact下面&…

WPF 本地化的最佳做法

WPF 本地化的最佳做法 资源文件英文资源文件 en-US.xaml中文资源文件 zh-CN.xaml 资源使用App.xaml主界面布局cs代码 App.config辅助类语言切换操作类资源 binding 解析类 实现效果 应用程序本地化有很多种方式,选择合适的才是最好的。这里只讨论一种方式&#xff0…

pytorch单机多卡后台运行

nohup sh ./train_chat.sh > train_chat20230814.log 2>1&参考资料 Pytorch单机多卡后台运行的解决办法

kafka-2.12使用记录

kafka-2.12使用记录 安装kafka 2.12版本 下载安装包 根据你的系统下载rpm /deb /zip包等等, 这里我使用的是rpm包 安装命令 rpm -ivh kafka-2.12-1.nfs.x86_64.rpm启动内置Zookeeper 以下命令要写在同一行上 /opt/kafka-2.12/bin/zookeeper-server-start.sh /opt/kafka-2…

实验二十八、三角波发生电路参数的确认

一、题目 利用 Multisim 确定图1所示电路中各元件的参数,使输出电压的频率为 500 Hz 500\,\textrm{Hz} 500Hz、幅值为 6 V 6\,\textrm{V} 6V 的三角波。 图 1 三角波发生电路 图1\,\,三角波发生电路 图1三角波发生电路 2、仿真电路 A 1 \textrm A_1 A1​ 采用…

深入解析 MyBatis 中的 lt;foreachgt; 标签:优雅处理批量操作与动态 SQL

在当今的Java应用程序开发中&#xff0c;数据库操作是一个不可或缺的部分。MyBatis作为一款颇受欢迎的持久层框架&#xff0c;为我们提供了一种优雅而高效的方式来管理数据库操作。在MyBatis的众多特性中&#xff0c;<foreach>标签无疑是一个强大的工具&#xff0c;它使得…