四格表fisher检验

news2024/9/20 14:28:02

一、案例介绍

某医生用新旧两种药物治疗某病患者27人,治疗结果见下表,现在想知道两种两种药物的治疗效果有无差别?

二、问题分析

本案例的分析目的是探究两种治疗效果有无差异,总样本量为27<40,所以考虑使用四格表Fisher确切检验法进行分析,但是需要满足2个条件:

条件1:分组变量和观察变量均为二分类变量;本案例数据组别和治疗效果都是二分类变量,满足该条件。

条件2:观测值之间相互独立;本案例中27名患者之间都是独立的,治疗效果不存在互相干扰情况,满足该条件。

所以,可以使用四格表Fisher确切检验法进行分析。

三、软件操作及结果解读

(一)理论说明

四格表Fisher确切检验法的基本思想是:在周边合计数保持不变的条件下,计算表内4个实际频数变动时的各种组合之概率Pi;再按检验假设用单侧或双侧的累计概率P,依据所取的检验水准α做出推断。

Fisher确切检验法是由R.A.Fisher于1934年提出的,是基于超几何分布理论直接计算拒绝零假设的概率的,所以不会有卡方值。该方法依据超几何分布,利用排列的数学原理进行计算,计算量比较大,不适合手算,故使用软件计算直接得到p值进行判断。

(二)软件操作

在SPSSAU实验/医学研究模块,选择【Fisher卡方】。按照提示,在表格中填写案例数据,注意A1单元格一定要空着,且不需要填写汇总数据,操作如下图:

(三)结果解读

1、检验统计量表格

SPSSAU输出检验统计量结果如下:

从上表可知,Fisher确切检验对应的p值为0.3705,>0.05,所以不能拒绝原假设,即尚不能认为两种药物的治疗效果不同。

同时可以看到,除了Fisher卡方外,SPSSAU同时输出Pearson卡方值和连续校正卡方值以及分别对应的p值,从p值结果来看,3个值均呈现出一致性结论,即尚不能认为两种药物的治疗效果不同。

那么3种卡方值具体应该如何选择呢?需要看期望频数表格进行判断。

2、期望频数表格

SPSSAU输出案例期望频数如下:

通常情况下:

如果总样本量大于40且期望频数值全部均大于5时,一般使用Pearson卡方值;

如果总样本量大于40,但出现期望频数小于5的单元格时,可优先使用连续校正卡方,也或者使用fisher卡方值;

如果总样本量小于40,也或者出现期望频数小于1的单元格时,此时建议使用fisher卡方检验。

由此我们也明白了为什么在第二部分问题分析时提出,本案例总样本量为27小于40,所以使用Fisher确切检验法进行分析。

四、结论

本案例采用2*2Fisher确切检验法对新药和旧药的治疗效果进行比较分析,由于总样本量27<40,且满足Fisher确切检验法需要满足的条件,所以采用该方法进行分析。分析结果显示,Fisher确切检验法对应的p值为0,3705>0.05,所以不能拒绝原假设,两种药物的治疗效果差异无统计学意义,即尚不能认为新药和旧药两种药物的治疗效果不同。

五、知识小贴士

(1)三类卡方检验值应该看哪种?

通常情况下:如果总样本量>40且期望频数值全部均大于5时,一般使用pearshon卡方值;如果总样本量大于40,但出现期望频数小于5的单元格时,可优先使用连续校正卡方,也或者使用fisher卡方值;

如果总样本量小于40,也或者出现期望频数小于1的单元格时,此时建议使用fisher卡方检验。具体以文献为准即可,可能不同文献的标准不完全一致。

(2)期望频数如何计算?

我们称每个观测变量出现的次数为频数,期望频数是频数的期望值。对于r*c行列表,第r行第c列格子中观察频数O_{r c}对应的期望频数为E_{r c}=\frac{n_{r}n_{c}}{n},其中n_{r}为第r行的合计数,n_{c}为第c列的合计数,n为总样本量。

如本案例中,第1行1列旧药有效数为2,则其期望频数为5*16/27=2.963;其余计算类似。

参考文献:

[1]颜红,徐勇勇.医学统计学.第3版[M].人民卫生出版社,2015

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

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

相关文章

NB-IoT模块(BC系列—BC95)详解

NB-IoT模块&#xff08;BC系列—BC95&#xff09; 0. NB-IoT概述技术原理特点和优势应用领域 1. 常用的NB-IoT模块2. BC系列—BC95技术规格功能特点 3. STM32使用BC95方法BC95的AT指令示例代码 0. NB-IoT概述 NB-IoT&#xff08;Narrowband Internet of Things&#xff09;是一…

万字长文解析最常见的数据库恢复算法: ARIES

#万字长文解析最常见的数据库恢复算法: ARIES 首发地址&#xff1a; https://mp.weixin.qq.com/s/Kc13g8OHK1h_f7eMlnl4Aw Introduction 上图中为基于 WAL 的数据库的一种可能的架构情况。其中&#xff0c;In-Memory Data 为数据库数据在内存中的组织形式&#xff0c;可以是 B …

Element-ui 实现多个日期时间发范围查询

1、前端 <el-form-item label"生产时间"> <el-date-picker v-model"dateProduct" style"width: 240px" value-format"yyyy-MM-dd" type"daterange" range-separator"-" start-placeholder"生产开始…

(三)解析函数及其性质

本文主要内容如下&#xff1a; 1. 复变函数的导数与微分1.1. 复变函数可导、可微、解析与奇点的定义1.2. 复变函数可微的充要条件1.3. 关于复变函数可微性判定的其它形式1.4. 相关结论1.5. 解析函数的构造 2. 解析函数与调和函数2.1. 调和函数与共轭调和函数2.2. 解析函数与调和…

cglib bean复制报错:module java.base does not “opens java.lang“ to unnamed module

在使用cglib bean复制功能时&#xff0c;报下面的错误 Caused by: net.sf.cglib.core.CodeGenerationException: java.lang.reflect.InaccessibleObjectException-->Unable to make protected final java.lang.Class java.lang.ClassLoader.defineClass(java.lang.String,b…

牛客社区项目续

统一日志记录 我们的这个项目在很多地方都需要记录日志&#xff0c;比如帖子模块、评论模块、消息模块等&#xff0c;而以前我们记录日志都是在某一个功能点中使用日志工厂&#xff0c;像下面这样&#xff1a; 项目中很多地方都需要记录日志&#xff0c;像这样一个地方一个地方…

sumo的几种安装方法

sumo的几种安装方法 sumo有很多中安装方法&#xff0c;根据你需要的任务来自己选择&#xff1a; 采用官网的latest version来进行安装 sudo add-apt-repository ppa:sumo/stable sudo apt-get update sudo apt-get install sumo sumo-tools sumo-doc想要安装源码来进行自己b…

【Java可执行命令】(六)调试工具 jdb:深入解析应用程序调试工具jdb ~

Java可执行命令详解之jdb 1️⃣ 概念2️⃣ 优势和缺点3️⃣ 使用3.1 语法格式3.1.1 参数&#xff1a;-sourcepath < path>3.1.2 指令&#xff1a;run [class [args]]3.1.3 指令&#xff1a;print < expr>3.1.4 指令&#xff1a;stop at< class>:< line>…

如何利用Idea回滚代码以及Cherry-Pick部分代码

引言 大家在版本迭代过程中&#xff0c;是否遇到过开发好的需求&#xff0c;都已经合并到Master分支等待发布后&#xff0c;临时通知不需要上线了的情况。这个时候一般会要求只上一部分紧急功能或者别的新功能&#xff0c;那么这个时候就需要用到Git的Reset以及Cherry-Pick功能…

钉钉机器人用bitmap实现签到记录

现在是周五晚上&#xff0c;下面是一个二进制数字&#xff0c;其中&#xff0c;有16位&#xff0c;最后一位下标是15&#xff0c;今天晚上是14&#xff0c;我签到成功了

一定要收藏的30套可视化大屏制作模板!升职加薪不再是梦想!

前几天和朋友吃饭聊天&#xff0c;他吐槽说老板让他做可视化大屏&#xff0c;但他不会敲代码根本做不出来&#xff0c;老板动动嘴巴子根本不考虑技术难度只想看到成果&#xff0c;他焦虑得都睡不着觉。我给他分享了一套可视化大屏模板&#xff0c;10分钟就制作完成了老板要求的…

了解浏览器缓存

什么是HTTP缓存&#xff0c;如何工作的&#xff1f; 当我们打开一个页面时&#xff0c;会向服务端发起很多次请求&#xff0c;如下图打开百毒首页&#xff0c;发起了HTML、各种图片、JS、CSS等资源共101次请求。这里面很多资源并不会频繁变化&#xff0c;每次打开页面都重新请…

利用Python构建宁德时代、比亚迪、隆基绿能股票时间序列预测模型

存货 import tushare as ts # 导包 import numpy as np import matplotlib.pyplot as plt from scipy.signal import find_peaks from scipy.stats import norm import datetime import pandas as pd import seaborn as sns # pip install seaborn import matplotlib.patches …

Oracle报错:“Error in invoking target ‘agent nmhs’ of makefile...”

前言&#xff1a;Oracle在安装过程中的报错一定要重视&#xff0c;这决定你后续是否能完成安装以及是否能使用。我这边会陆续汇总一些报错现象以及解决方案共享。 ##Install Product 86%报错信息 &#xff1a;“Error in invoking target agent nmhs of makefile...”解决方案…

意大利语翻译哪个公司比较专业?

据了解&#xff0c;意大利语除了通行于意大利之外&#xff0c;还通行于美国、加拿大、阿根廷和巴西等29个国家的&#xff0c;其应用极为广泛。随着意大利语翻译需求量的日益增加&#xff0c;也促进了意大利语与其它语言间的交流、转化和发展。那么&#xff0c;意大利语翻译难吗…

抖音短视频矩阵号系统开源部署搭建分享(二)

开发背景&#xff1a;抖音seo源码&#xff0c;抖音矩阵系统源码。抖音获客系统源码&#xff0c;短视频矩阵系统源码开源搭建&#xff08;MySQL数据库&#xff0c;NGINX&#xff0c;PHP7.4&#xff0c;MySQL5.7&#xff0c;redis&#xff09; 技术要点&#xff1a; 服务器配置服…

vue-router.esm.js:2248 Error: Cannot find module ‘@/views/dylife/ 报错解决

具体是展示 一直加载 控制台报找不到模块 webpack版本问题&#xff0c;webpack4 不支持变量方式的动态 import &#xff0c;新版本需要使用 require() 来解决此问题。 return () > import(/views/${view}) 改写成 return (resolve) > require([/views/${view}], reso…

python接口自动化(四)--接口测试工具介绍(详解)

简介 “工欲善其事必先利其器”&#xff0c;通过前边几篇文章的介绍&#xff0c;大家大致对接口有了进一步的认识。那么接下来让我们看看接口测试的工具有哪些。 目前&#xff0c;市场上有很多支持接口测试的工具。利用工具进行接口测试&#xff0c;能够提供测试效率。例如&…

【C++】浅析C++中的虚函数

关于虚函数 Q1&#xff1a;观察一个类引入虚函数后&#xff0c;类会发生什么变化&#xff1f; 首先&#xff0c;创建一个空类A&#xff0c;并实例化出A的一个对象a&#xff0c;计算一下这个对象占用多少字节&#xff1a; #include<iostream> using namespace std;clas…