医学案例|单因素重复测量方差

news2025/2/23 12:08:45

一、案例介绍

为研究某诱导方法对患者收缩压的影响,将手术要求基本相同的15名患者,在手术过程中采用A麻醉诱导方法,并在T0(诱导前)、T1、T2、T3、T4五个时刻测量患者的收缩压,试着进行方差分析。

二、问题分析

可以从案例得到,样本数据包括诱导前T0、T1、T2、T3、T4五个时刻,所以考虑使用重复测量方差进行单因素重复测量方差,对于重复测量实验数据的方差分析一般需要考虑组内的影响。

三、软件操作及结果解读

(一) 数据导入

1.数据格式

重复测量方差数据的特殊之处在于需要有ID号,比如病例号等,以及时间点数据,同一个ID(病例)会有多个不同时间点的数据,比如10个样本分别有5个相同测量时间点,一个样本有5个,那么就有10*5=50行数据。数据格式如下:

2.导入数据

将整理好的数据上传至SPSSAU系统内,如下:

上传结果如下:

(二) 适用条件判断

1.异常值

和大多数分析一样,一般在分析前都需要进行一些检验,重复测量方差需要检查数据是否存在异常值,以及因变量是否满足正态分布。

异常值分析数据中不正常的值,也称离群值。异常值的检验方法有很多,比如可以采用描述统计判断法(一般认为介于3倍标准差外的数据),或者箱线图等。这里使用箱线图进行判断。利用SPSSAU描述分析,结果如下:

从箱线图可以粗略的判断没有异常值。接下来查看因变量是否满足正态分布。

2. 正态性检验

判断数据是否符合正态分布有很多种方法,比如正态性检验(最为严格),描述统计法(查看偏度和峰度),以及图示法等等,这里使用直方图进行判断。

发现“收缩压”呈现倒扣的钟形,所以满足正态性分布标准,数据满足方差分析的前提条件,可进行重复测量方差。

(三)单因素重复测量方差分析

  1. 软件操作
    重复测量方差的分析路径为点击【实验医学/研究】→【重复测量方差】然后进行分析:

  1. 结果解读

数据满足分析条件,接下来进行重复测量方差,首先查看球形度检验的结果:

球形度检验可以看到p值大于0.05,所以通过球形度检验不需要校正p值。【特别提示:进行球形度检验p值大于0.05则通过检验不需要校正p值,如果p值小于0.05,则没有通过球形度检验需要校正p值,接下来在看球形度W值,若w值小于0.75则使用GG校正,反之使用HF校正】。

组内分析:

根据球形度检验发现需要校正p值,从上表可以看出主效应的p值小于0.05,呈现0.05水平显著,使用A麻醉诱导方法不同时间点的收缩压有显著性差异分析,最终结果如下:

  1. 均方(效应)
    平方和SS(效应)/df(效应),例如:1110.28/4=277.57;
  2. 均方(误差)
    平方和SS(误差)/df(误差),例如:480.92 /4=13.359;
  3. F
    F值=均方(效应)/均方(误差),例如:277.57/13.359=20.778;
  4. P
    P值是由F值得到的。
  5. ges(generalizedeta-squared)
    平方和SS(效应)/【平方和SS(效应)+平方和SS(误差组内)+平方和SS(误差组间)】,例如:1110.28 /(1110.28+1136.48+480.92)=0.407;
  6. 偏Eta方

偏Eta方=平方和SS(效应)/[平方和SS(效应)+平方和SS(误差)],例如:745176.32/(745176.32+1136.48)=0.998;

如果主效应显著想要进一步两两之间比较,可以使用事后多重比较进行分析。SPSSAU提供多种事后多重比较方法,对于他们之间的差异可以查看帮助手册,这里选择Tukey方法(各组别的样本数量相同时使用),假如想要了解不同时间点之间抑郁程度是否有差别。结果如下:

通过事后多重比较发现T0和T1、T0和T3、T0和T4等(标黄区域)的数据p值小于0.05,即存在差异。说明T0和T1等两个时间测量的收缩压有明显差异。

四、结论

对案例数据进行重复测量方差分析,首先整理成正确的数据格式,然后分析前检查数据是否满足相关检验,接着对数据进行分析,分析结果得到数据不满足球形度检验发现需要校正p值,并且发现使用A麻醉诱导方法不同时间点的收缩压有显著性差异分析,利用事后多重比较进行具体查看,发现T0和T1、T0和T3、T0和T4等p值小于0.05,即存在差异。

五、知识小贴士

1.SPSSAU对事后多重比较,简单效应的表格输出说明如下:

第一、事后多重比较规则为:针对组内项和组间项(如果有)均输出事后多重比较;

第二、简单效应输出规则为:仅针对第1个组内项,分别与组间项(如果有)进行简单效应表格输出,如果组内项超过1个,可通过互换组内项的位置输出其余组内项与组间项的简单效应表格。

2.事后多重比较与‘单独进行事后多重比较’结果不一致?

单独进行事后多重比较(进阶方法->事后多重比较法)时,模型实质上为单因素方差,仅考虑1个X的情况,标准误差的计算并不一致,因此结果会不一致,但通常情况下结论会保持一致;以及此处事后多重比较使用的是边际估计均值(偏最小二乘均值)与一般意义上的平均值有所区别,类似于SPSS软件的EMMEANS功能。SPSSAU默认针对组内项和组间项均进行事后多重比较。

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

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

相关文章

小程序相较于APP,广告变现有哪些优势?

对于开发者而言,微信小程序开发门槛相对较低,难度不及APP,能够满足简单的基础应用,适合生活服务类线下商铺以及非刚需低频应用的转换。 和早期相比,今天小程序在产品功能、UI设计、交互体验等方面,越来越精…

python_day16_装饰器

装饰器基础写法 def outer(func):def inner():print("睡眠")func()print("起床")return innerdef sleep():import timeprint("睡眠中、、、")time.sleep(5)f1 outer(sleep) f1()装饰器语法糖 def outer(func):def inner():print("休息一…

Inpaint Anything: 自动化抹除视频元素

自动化抹除视频元素 不用逐帧抠图,直接SAM Tracking Video Inpainting就能实现自动化抹除奔跑吧idol。 https://github.com/geekyutao/Inpaint-Anything 目录 网站演示参考文献 网站 https://huggingface.co/spaces/InpaintAI/Inpaint-Anything 演示 原理就是&a…

Bert经典变体学习

ALBert ALBERT就是为了解决模型参数量大以及训练时间过长的问题。ALBERT最小的参数只有十几M, 效果要比BERT低1-2个点,最大的xxlarge也就200多M。可以看到在模型参数量上减少的还是非常明显的,但是在速度上似乎没有那么明显。最大的问题就是这种方式其实…

力扣算法递归类—(链表)两数相加

目录 2. 两数相加 题解: 知识点: 链表定义: 链表赋值: 链表题小技巧: 代码: 结果: 给你两个 非空 的链表,表示两个非负的整数。它们每位数字都是按照 逆序 的方式存储的&…

idea 里 controller service impl mapper xml 切换跳转快捷键

首先在controller层,对着接口点方法的方法上按着ctrl和鼠标左键,你会进入service层。 对着方法ctrlaltb不按鼠标,你会进入impl层。service层的方法上按ctrl和鼠标左键会回到controller,ctrlaltb不按鼠标也会进入到impl层,impl上的…

0基础学习VR全景平台篇 第73篇:VR直播-如何自定义邀请二维码(直播邀请)

自定义直播邀请二维码是自定义直播间邀请卡上显示的二维码,若上传,那么便会替换掉邀请卡上原有的二维码,原二维码为本场直播活动的二维码。 建议上传的尺寸为300px*300px,可选择开启二维码的弹出效果,开启后&#xff0…

第三章 数据链路层

第三章 数据链路层 3.1 数据链路层的几个基本概念 数据发送模型 数据链路层主要的两种信号类型 点对点信号:这种信道使用一对一的点对点通信方式;广播信道:这种信道使用一对多的广播方式,因此过程比较复杂。广播信道上连接的主机…

管理类联考——写作——真题篇——论说文——企业管理=2011+2013+2014+2015+2016+2017+2020真题论说文

2020年探究的就是管理者做决策的时候,应该多听听专家的意见; 2017年探究的是企业在研发新产品和扩大生产之间的抉择问题; 2016年也可以看成企业的问题,就是企业组织的多样化和一致性。 2015年则说明企业在积累财富的过程中&#x…

又双叒叕被调查?竟是Elsevier旗下中科院TOP?如何避雷看这几个指标!

近期小编发现处于“On Hold”状态的期刊愈发的多了,总在不明所以的情况下,科睿唯安官网就将期刊拉入“审查”名单,包括前几天发现的同为Elsevier旗下的中科院TOP期刊Chemosphere (👉参考:慎投!…

科技云报道:是时候全员FinOps了吗?

科技云报道原创。 在论坛上,国外某企业的真实案例引发了热议。一开始该企业只顾技术创新,积极上云,不顾成本。 直到有一天,高层介入喊停:“这个云不能再上了,成本已经远大于收益了”。该企业因为成本失控…

数据库优化器设计穿越探索之旅

作者:王晨(道客) 一、前言 引用来自百度百科的话术:在数据库技术发展历史上,1970 年是发生伟大转折的一年,因为这一年的6月,IBM的圣约瑟研究实验室的高级研究员Edgar Frank Codd在Communications of ACM 上发表了《A…

【JavaEE初阶】Servlet (二) Servlet中常用的API

文章目录 HttpServlet核心方法 HttpServletRequest核心方法 HttpServletResponse核心方法 Servlet中常用的API有以下三个: HttpServletHttpServletRequestHttpServletResponse HttpServlet 我们写 Servlet 代码的时候, 首先第一步就是先创建类, 继承自 HttpServlet, 并重写其…

springboot拿到yml文件中的配置,决绝将配置写在代码里

springboot拿到yml文件中的配置 写你对应的文件 package com.gym.vueea.config;import org.springframework.boot.context.properties.ConfigurationProperties; import org.springframework.stereotype.Component;/*** author ymgu* Package com.gym.vueea.config* date 2023…

230717期班5G网络优化就业班正式开班!命运的齿轮开始转动~

当前,无论是个人生活还是商业领域,5G网络已经成为了不可或缺的一部分,所以各位同学看好这个行业,在这个炎炎夏日的七月,经历了8天的基础班之后,也愿意投身到5G网络优化这个行业中,所以继续参与2…

logback中文一直是乱码,logback中文问号

logback一直是乱码 方案一加上UTF-8 方案二我这边方案一不行 在启动参数加上 -Dfile.encodingutf-8 这个竟然就可以了

IDEA如何快捷创建serialVersionUID【详细图解】

在Java Bean,快速创建serialVersionUID,详细操作流程如下: 1.File->Settings->Editor->Inspections 在搜索框里搜索 Uid,选择下图中勾选的选择 2.如何使用 双击选中需要序列化Uid的类名,使用Alt enter,快捷…

软件测试基本知识

安全测试 安全防护策略?(漏洞扫描、入侵检查、安全日志、隔离防护) 安全日志:用于记录非法用户的登录名称、操作时间及内容等信息,以便发现问题并提出解决措施;安全日志仅记录相关信息,不对非…

openssl3.1.1关于国标支持的验证笔记

openssl3.1.1关于国标支持的验证笔记 openssl的版本差异日志 openssl虽然有3个大分支,我们就以3.1大分支查看关于国密的差异日志。 Changes between 1.1.0i and 1.1.1 [11 Sep 2018] Changes between 1.1.1b and 1.1.1c [28 May 2019] Changes between 1.1.1k and…

kotlin 编写一个简单的天气预报app(四)

编写界面来显示返回的数据 用户友好性&#xff1a;通过界面设计和用户体验优化&#xff0c;可以使天气信息更易读、易理解和易操作。有效的界面设计可以提高用户满意度并提供更好的交互体验。 增加城市名字的TextView <TextViewandroid:id"id/textViewCityName"…