mfc绘制bezier b样条三种曲线_生存曲线(二):SPSS和Origin绘图教程及相关问题

上期生存分析推送后,有粉丝在后台问:有发生率的情况,为什么要做生存曲线分析?

举个例子:临床试验中,共招募30位胃腺癌患者,均分为3组,分别使用了A/B/C三套不同的治疗方案,治疗期为2年。2年后统计3个组每组均死亡3人,但A组0.5年时死亡3人,B组是1.5年时死亡3人,C组是接近2年时死亡3人。此时,如果单纯分析死亡率,3个组的死亡率完全相同。但是,如果用生存分析,就能看出3个组的患者生存情况明显不同。


进入今日正题......

上一期,我们提到了如何使用GraphPad Prism绘制生存曲线。

实际上,很多软件都可以绘制生存曲线并加以统计分析有这种功能,包括SPSS、Origin、Stata以及R语言等。

综合看来,GraphPad Prism更容易上手,简单有效。Stata和R语言涉及一丢丢编程语言,可能相对不太容易上手。

由于,有相当一部分人喜欢使用SPSS和Origin,因此今天就拿这两个软件说一说如何绘制生存曲线。

4a6790751a056dec5a2f902ff9935788.png

cdc233d9844b697ffdbe745f7cbc3279.png

还是使用上一期的数据作为示例。

eg:一共有36只裸鼠,施加某因素用于制作模型后全部成模,并随机均分为A、B、C、D组,每组9只。随后A组被设定为control组,不做治疗。B、C、D组动物分别给予药物1、药物2、药物3治疗。生存时间小于等于50天作为终点事件(即50天内死亡),以此比较几组之间是否存在差异?

那么常见的原始数据如下图,数值为动物存活天数。

557adb187edf6ba2f8037ee06ff8e3ae.png

SPSS篇

1. 将上表中的数据重新编码为以下格式。state中的1代表死亡(裸鼠存活≤50天),而0代表存活(裸鼠存活>50天)。treatment中的1代表control,2/3/4分别代表3种药物。

f58e3c5f71d5b2c11dc5bd1958080314.png

2.仅复制数据(不要复制标题)到SPSS数据表中,然后再定义列名称和值标签。

dfb6a7cef1ead412cf728b4c475c3cff.png

3. 数据放好后,选择分析 → 生存分析 → Kaplan-Meier。

注:Kaplan-Meier是生存函数,允许有一个分组变量进行生存率的组间比较,还可容许一个分层变量,是生存分析的常用统计方法。

aba12288e61c946acc22967690ca6c2d.png

4. 在弹出的对话框中,按照下图将数据放入相应框中。

00e627df7be82fbe8525983c0339c208.png

5. 注意了。我们要看不同组之间生存率的差异,关键事件为动物死亡,编码数据时就将死亡编码为1,而存活则编码为0。所以,点击“定义事件”,单值填1,1代表了动物死亡这个事件已发生,再点击继续。

033053589ce031aa06f2356b7d540d07.png

有人会问,如果我选择单值填0呢?大多数动物生存时间小于50天,仅少量动物生存大于50天。因此“≤50天”内的生存差异成为关注的焦点。如果选择单值填0,此时软件认为你关注的焦点在大于50天之后,即50天后存活的动物存活期差别,生存曲线也会彻底改变。大家可以试试,比较一下。

6. 点击“选项”,选择“生存分析函数”,其它选项默认,点击继续后,再点击确定。

83430ed02fb8f330560d36458b806ce2.png

7. 在输出文档中,可以看到生存曲线。样子是有了,就是不太美观。

a1e4ec5d9e388ad8d99586facc174953.png

Origin篇

1.Origin的数据录入格式与SPSS一样,只不过需要将标题一起复制粘贴到Origin数据表之中。

7aba98ea8f5336b4f5c5052f614e653f.png

2. 粘贴到F(x)这一栏下方,并修改列标题名称,如下图。

bf05f6424cf791d2305542d1ddebc1e6.png

3. 继续按照下图选择Kaplan-Meier分析

1e84f3568847c865e3f8fe7d439c9882.png

4. 在弹出的对话框中,按照下图设置。点击小三角符号,将标题分别填入。

642977cc5abf9ad4e82f91c78b67f68c.png

5. 注意了。Censoring Value(删失值)处设置为0。解释一下,在生存分析中,只要没有发生关键事件(如意外、病例脱落等)的情况,这种值都要定义为删失值,此值需定义为0,因为它没有发生死亡这个关键事件。当然了,删失值不能太多,否则会影响生存分析结果的可信度。

7fa5c931723ddc4a6f548af227953c82.png

6.将 Test Equality处的三种差异分析方法都选上。所有其它选项默认,点击OK。

f30789f8aaf7397aa7880e600dff5432.png

7.输出的生存曲线效果如下。当然了,Origin支持对图形就行修改和美化,此处暂且不表。

317eecf786d48063fd6a4110228e5f8a.png

Ending

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

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

相关文章

8位可控加减法电路设计_高级工程师:相同的地线GND,却有不同的电路设计含义...

问一个简单而又很难回答的电路问题:电路中的地线GND,它的本质是什么?工程师,在研发设计一个电路项目时,一般会经历三个阶段:电路项目PCBA板1,项目方案论证项目方案论证,是在项目前期…

word刷子刷格式_Excel技巧—开始菜单之格式刷六大功能

点赞再看,养成习惯千里之行,始于足下今天我们主要学习的是“第一区块-剪贴板”中的“格式刷”。“格式刷”主要将一个对象的颜色、字体样式、字体大小、边框样式等所有格式复制到目标对象上,我们可以把它理解为格式的复制粘贴。说到“格式刷”…

html/css题库,DIV+CSS题库

CD2、在CSS语言中下列哪一项是"左边框"的语法?(C )A、border-left-width:B、border-top-width:C、border-left:D、border-top-width:3. 选出你认为最合理的定义标题的方法( C )A. 文章标题B.文章标题C.文章标题D. 文章标题4、在CSS语言中下列哪一项的适用…

紫金计算机网络,南京理工大学紫金学院《计算机网络技术》考试复习试题试题(含答案解析)2...

一、选择题:(每题2分,共40分)1.在TCP/IP参考模型中,TCP协议工作在()。A、应用层B、传输层C、互联网层D、网络接口层2.下列()是合法的E-mail地址。A、mnetworkhttp://www.doczj.com/doc/fe7e6e475ebfc77da26925c52cc58b…

go语言和python结合_Go+Python双语言混合开发-第1章 【阶段一:Go语言基础】

1-1 GoPython双语言混合开发-课程导学 (11:35)1-2 课程资源获取方式 - 非常重要!!! (06:07)1-3 课程中需要用到的开发环境介绍 (03:47)1-4 windows中安装centos7虚拟机 (19:16)1-5 xshell的安装和使用 (04:14)1-6 git的安装和配置 (06:22)1-7…

cwinthread*线程指针怎么销毁结束_最新版Web服务器项目详解 01 线程同步机制封装类...

点 击 关 注 上 方"两猿社"设 为"置 顶 或 星 标",干 货 第 一 时 间 送 达。互 联 网 猿 | 两 猿 社基础知识RAIIRAII全称是“Resource Acquisition is Initialization”,直译过来是“资源获取即初始化”.…

您的计算机无法访问dota2,解决win7系统DOTA2无法连接至Steam网络的设置方法

随着电脑的使用率越来越高,我们有时候可能会遇到win7系统DOTA2无法连接至Steam网络问题,如果我们遇到了win7系统DOTA2无法连接至Steam网络的问题,要怎么处理win7系统DOTA2无法连接至Steam网络呢?接下来给大家带来win7系统DOTA2无法…

根据经纬度计算范围_高中地理必修一二三思维导图+计算公式全汇总!能用3年...

必修一必修二必修三计算公式1.经纬度计算:经度差与地方时差算经度——地方时每相差1小时,经度相差1;纬差法与正午太阳高度算纬度——正午太阳相差多小,纬度相差多少;北极星的仰角即地平高度等于当地地理纬度;经纬线上长度算经纬度——1经线长…

linux 系统 可视化工具实时监控_Linux上的实时监控平台-你需要这样做

应朋友们要求,介绍上Linux系统下的实时监控平台,在上次提到了glances,它提供了较多的监控指标,那如果我们要看历史数据呢?某一时间段的回放呢?显然glances是做不到的。因此,实时监控平台就应运而…

职称类计算机论文,计算机类职称论文发表(2)

计算机类职称论文发表篇二基于计算机通信网络的采矿传输技术【摘 要】煤矿开采是人类合理开发自然资源的主要途径,对国民产业经营收益水平也有很大的影响。信息科技改革指导下,分析了基于计算机通信网络的采矿信息传输方案,指导了煤炭开采新模…

linux 普通用户touch权限不够_Linux 开启指令

简单指令1.1、 ls 指令语法1:# ls [路径]表示列出指定路径下的文件夹和文件 的名字,如果路径没有指定则列出当前路径下的(list)列如:在root用户的家目录中输入ls命令,则会看到以下的效果:要求列…

5、this调用语句必须是构造函数中的第一个可执行语句_ES6中的Promise和Generator详解...

简介ES6中除了上篇文章讲过的语法新特性和一些新的API之外,还有两个非常重要的新特性就是Promise和Generator,今天我们将会详细讲解一下这两个新特性。Promise什么是PromisePromise 是异步编程的一种解决方案,比传统的解决方案“回调函数和事件”更合理和…

华中师范大学本科生计算机课教师,华中师范大学第八届大学生计算机设计大赛的通知...

为了激发大学生学习计算机知识和技能的兴趣和潜能,培养其创新能力和团队合作意识,以及运用信息技术解决实际问题的综合实践能力,切实提高计算机教学质量,教育部高等学校计算机类专业教学指导委员会、软件工程专业教学指导委员会、…

labelimg如何调整框的颜色_如何制作摄影集(下)

如果说上篇讲的是有关排版设计,那么下篇就是有关具体实施了。如果你正对作品集的制作跃跃欲试,那么下篇的内容仍旧非常重要,而且非常干货。 目录:软件推荐(Id)字体推荐纸张选择打印方式(激光、喷…

highcharts 显示平均值数值_拼多多评价多久能显示,有什么出评价技巧吗?

很多商家都知道,拼多多的评价,有时不一定会在买家评价后就立马出现,一般得等一段时间才会显示,,这个时间一般不会很久,24小时之内,评论一般在审核通过后就会展示出来。如果长时间不展示&#xf…

辽宁省大学生计算机系统与程序设计竞赛,2019CCF大学生计算机系统与程序设计竞赛(华东分赛区)在我校顺利举办...

5月18日,2019CCF大学生计算机系统与程序设计竞赛(CCSP)华东分赛区比赛及颁奖会在我校举办。今年是CCSP大赛首次采用区域分赛区的比赛,共分为7个赛区,包括东北区(哈尔滨工业大学承办),华北区(北京邮电大学承办)、华东区(南京航空航…

arm b bl 地址无关码_ARM汇编语言入门(六)

Part 6:条件状态和分支在探讨CPSR时我们已经接触了条件状态。我们通过跳转(分支)或者一些只有满足特定条件才执行的指令来控制程序在运行时的执行流。通过CPSR寄存器中的特定bit位来表示条件状态。这些位根据指令每次执行的结果而不断变化。例…

bcp out 带列名导出_从零开始学习 MySQL 系列索引、视图、导入和导出

阅读本文大概需要 8 分钟前言上篇文章我们学习了数据库和数据表操作语句,今天我们学习下数据库索引,视图,导入和导出的知识。作为基础篇,不会涉及到关于索引和视图的高级应用和核心概念,但是基本操作大家会了解&#x…

四年级计算机课程,信息技术(四年级)全部课程PPT课件.ppt

信息技术(四年级)全部课程PPT课件调整窗口的大小 最小化 最大化 关闭 水平的双向箭头可以调整窗口的宽度,让窗口变宽或变窄。 调整窗口的大小 最小化 最大化 关闭 倾斜的双向箭头可同时调整窗口的高和宽。 做一做 打开我的文档,对窗口进行最小化、最大化…

python等腰梯形_简单空实心图形打印|Python练习系列[3]

def shape_print(n): #实心等腰三角形 for y in range(n): for x in range(n-y-1):#先循环打印空格 形成一个倒直角三角形 range()中的值是形成的规律 print( ,end) for z in range(y*21):#再循环打印X 形成一个等腰三角形 range()中的值是形成的规律 print(X,end) print() pri…