1.提出科学问题
2.收集临床资料
3.选择合适统计学方法
4.选择合适的统计软件
5.评估结果,结合专业知识回答科学问题
这个框架真的很好用的,我在间接性残废地学习统计一年多以来,觉着这样的框架才是真正的功夫。将自己的统计工作,时刻比对这这五点进行下去,才不至于工作做了一半发现少了数据,然后重头返工,要么就是做到一半不知道怎么结束这项工作。
坑,是开始时候自己就给自己挖好了。我要做的就是选择一条好走的路,没有坑的路。
这个框架里,最难的是第一条,提出一个科学问题,往往别人忽略。
第三第四是结果的展示阶段,出了一张图,p值<0.05然后就可以万事大吉,所以也是各种帖子热衷展示的。当然,就像自己的内在不能总是拿出来展示,只能通过拍照的方式进行外表展示。所以,方法类的学习帖子一般作为自己学习笔记真的非常棒,但是为了技术而技术,为了搞出一个模型去收集数据,有点本末倒置了。我们是为了想解决一个实实在在临床问题,去收集数据,用统计学方法得出结论,从而提示医生如何对患者采取治疗措施,从而减少患者住院时间,减轻患者的病情,延长患者总生存时间。(这是我的初心吧)
闲话聊完,开始水这一篇的主角临床预测模型。我是看到公众号上《简单易懂》,觉得写的比较好,就忍不住想自己也写一点东西,毕竟自己摸鱼一年多了,一篇像样的文章都没出来,总得要做点什么来缓解自己无处释放的才华,所以打算对照自己学习路径以及课题需要,进行一系列临床模型构建的笔记输出。
《简单易懂》上对于临床预测模型的解释和举例,真的很通俗易懂了。临床模型是什么嘞:
分类和回归的区别在于输出变量的类型。
定量输出称为回归,或者说是连续变量预测; 定性输出称为分类,或者说是离散变量预测。
举个例子: 预测明天的气温是多少度,这是一个回归任务; 预测明天是阴、晴还是雨,就是一个分类任务。
(详见参考资料3)
临床模型,把它看成y=ax+b,通过身高预测孩子的升高,就像我们在儿科考试中遇到,给你有几个骨头预测小孩多大一样。而这个公式怎么来的,以及这个公式靠不靠谱就是模型构建的大部分内容了。
学术一点话就是模型构建和模型评价,在医学类临床资料统计中,我们常用的有线性回归,逻辑回归,cox回归等。
在我们拿到收集好的数据时候,我们要开始构建预测模型了,我们首先面对的一个问题就是该选择什么的因素纳入分析?
周支瑞老师对这种情况给出了三条原则:
1.当自己样本量足够大(每一个研究因素至少有20个患者,比如我想研究20个因素,那么我有400例以上的患者),那就SPSS软件中在Logistic回归和Cox回归中给出了7种变量筛选的方法,任君挑选;
2.如果不满足样本足够大,那么这个时候比较通用的做法是单因素cox回归(详见参考资料4,5)。用spss点点点,把所有因素都放进去,然后点cox回归,看p值,将p值<0.1的因素都挑选出来,进行多因素cox回归;
3.这个时候会出现,统计学不显著但是临床意义明确,例如,前列腺癌的gleason评分,单因素分析结果与已知临床专业知识决定纳入回归方程的变量。
这张周支瑞老师总结的“三种预测模型的研究思路”,真的大而全。

步骤1. 首先筛选影响低出独立影响因素,构建Logistic回归模型;
步骤2.绘制Nomogram;
步骤3.计算模型的区分度 C-Statistics。有三种方法可以计算其C-Statistics。
方法1. 利用 {rms} 包中的 lrm 函数构建Logistic回归模型,直接读取模型Rank Discrim.参数 C,即为C-Statistics。
方法2. 构建Logistic回归模型,predict函数计算模型预测概率,然后利用ROCR包根据此预测概率画ROC曲线,并计算曲线下面积AUC,此即为C-Statistics。注:此方法与SPSS中的计算方法一致。
方法3. 构建Logistic回归模型,predict函数计算模型预测概率,利用Hmisc包中somers2函数直接计算ROC曲线下面积AUC。注:此方法与SPSS中的计算方法一致。
不管那种方法都未给出标准误,所以可信区间的计算就很麻烦,如果一定要报告C-Statistics可信区间,可以考虑使用SPSS软件进行ROC分析,软件可以直接给出AUC的可信区间。
步骤4.重抽样的方法进行模型验证,并绘制Calibration曲线
步骤 1. 我们首先使用Cox回归基于构建预测模型并筛选独立预后因素(用于建模的数据集一般称为训练集或者内部数据集)。
步骤 2. 我们就以这三个独立的预后因素绘制Nomogram,建模完成。
步骤 3. 对上述两步所构建的预测模型的区分能力 (Discrimination) 进行评价,并计算C-index。
Cox回归的C-index计算方法:
方法1:直接从survival包的函数coxph结果中输出,需要R的版本高于2.15.需要提前安装survival包可以看出这种方法输出了C-index (对应模型参数C),也输出了标准误,95%可信区间就可以通过C加减1.96*se得到。并且这种方法也适用于很多指标联合。
方法2:利用rms包中的cph函数和validate函数,可提供un-adjusted和bias adjusted C指数两种。
步骤 4. 对模型进行验证,可通过外部数据集进行验证,如果无法获得外部数据集,笔者推荐采用Bootstrap冲抽样法基于训练集验证模型并绘制校正曲线(Calibration plot)。
Cox回归模型中C-Index计算,Nomogram绘制过程,Bootstrap法验证模型及绘制标准曲线
净重新分类指数(NRI)这个指标最初用于评价诊断试验中新的诊断指标较旧诊断指标把研究对象进行正确分类在数量上的变化。
综合判别改善指数(Integrated Discrimination Improvement, IDI)用于判断预测模型改善情况
决策曲线(Decision Curve Analysis):寻找一个预测净受益最大的模型
外部数据验证
内部数据重抽样验证
这个是真功夫,就像单细胞测序,有的人4例结肠癌单细胞测序可以发在NC,有了人只能一直砸在手里,总是分析不出来。所以,功夫在日常积累中。
我是看着周支瑞老师的书长大的,从《统计学》三本,到《临床预测模型》都是从买的正版书,如果大家想要专门深入学习,他的书真的可以买来反复演练,一定会有收获。还有他的公众号《 临床研究与医学统计》也是干活满满,安利一波
参考资料:
我有一个模型:classItem项目有一个属性“商店”基于存储的值,我希望Item对象对特定方法具有不同的行为。Rails中是否有针对此的通用设计模式?如果方法中没有大的if-else语句,这是如何干净利落地完成的? 最佳答案 通常通过Single-TableInheritance. 关于ruby-on-rails-Rails-子类化模型的设计模式是什么?,我们在StackOverflow上找到一个类似的问题: https://stackoverflow.co
我需要从一个View访问多个模型。以前,我的links_controller仅用于提供以不同方式排序的链接资源。现在我想包括一个部分(我假设)显示按分数排序的顶级用户(@users=User.all.sort_by(&:score))我知道我可以将此代码插入每个链接操作并从View访问它,但这似乎不是“ruby方式”,我将需要在不久的将来访问更多模型。这可能会变得很脏,是否有针对这种情况的任何技术?注意事项:我认为我的应用程序正朝着单一格式和动态页面内容的方向发展,本质上是一个典型的网络应用程序。我知道before_filter但考虑到我希望应用程序进入的方向,这似乎很麻烦。最终从任何
我有一个包含模块的模型。我想在模块中覆盖模型的访问器方法。例如:classBlah这显然行不通。有什么想法可以实现吗? 最佳答案 您的代码看起来是正确的。我们正在毫无困难地使用这个确切的模式。如果我没记错的话,Rails使用#method_missing作为属性setter,因此您的模块将优先,阻止ActiveRecord的setter。如果您正在使用ActiveSupport::Concern(参见thisblogpost),那么您的实例方法需要进入一个特殊的模块:classBlah
我有一个表单,其中有很多字段取自数组(而不是模型或对象)。我如何验证这些字段的存在?solve_problem_pathdo|f|%>... 最佳答案 创建一个简单的类来包装请求参数并使用ActiveModel::Validations。#definedsomewhere,atthesimplest:require'ostruct'classSolvetrue#youcouldevencheckthesolutionwithavalidatorvalidatedoerrors.add(:base,"WRONG!!!")unlesss
我想向我的Controller传递一个参数,它是一个简单的复选框,但我不知道如何在模型的form_for中引入它,这是我的观点:{:id=>'go_finance'}do|f|%>Transferirde:para:Entrada:"input",:placeholder=>"Quantofoiganho?"%>Saída:"output",:placeholder=>"Quantofoigasto?"%>Nota:我想做一个额外的复选框,但我该怎么做,模型中没有一个对象,而是一个要检查的对象,以便在Controller中创建一个ifelse,如果没有检查,请帮助我,非常感谢,谢谢
我有一些非常大的模型,我必须将它们迁移到最新版本的Rails。这些模型有相当多的验证(User有大约50个验证)。是否可以将所有这些验证移动到另一个文件中?说app/models/validations/user_validations.rb。如果可以,有人可以提供示例吗? 最佳答案 您可以为此使用关注点:#app/models/validations/user_validations.rbrequire'active_support/concern'moduleUserValidationsextendActiveSupport:
对于Rails模型,是否可以/建议让一个类的成员不持久保存到数据库中?我想将用户最后选择的类型存储在session变量中。由于我无法从我的模型中设置session变量,我想将值存储在一个“虚拟”类成员中,该成员只是将值传递回Controller。你能有这样的类(class)成员吗? 最佳答案 将非持久属性添加到Rails模型就像任何其他Ruby类一样:classUser扩展解释:在Ruby中,所有实例变量都是私有(private)的,不需要在赋值前定义。attr_accessor创建一个setter和getter方法:classUs
我有一个正在构建的应用程序,我需要一个模型来创建另一个模型的实例。我希望每辆车都有4个轮胎。汽车模型classCar轮胎模型classTire但是,在make_tires内部有一个错误,如果我为Tire尝试它,则没有用于创建或新建的activerecord方法。当我检查轮胎时,它没有这些方法。我该如何补救?错误是这样的:未定义的方法'create'forActiveRecord::AttributeMethods::Serialization::Tire::Module我测试了两个环境:测试和开发,它们都因相同的错误而失败。 最佳答案
ruby如何管理内存。例如:如果我们在执行过程中采用C程序,则以下是内存模型。类似于这个ruby如何处理内存。C:__________________|||stack|||------------------||||------------------|||||Heap|||||__________________|||data|__________________|text|__________________Ruby:? 最佳答案 Ruby中没有“内存”这样的东西。Class#allocate分配一个对象并返回该对象。这就是程序
我正在使用Rails3.1并在一个论坛上工作。我有一个名为Topic的模型,每个模型都有许多Post。当用户创建新主题时,他们也应该创建第一个Post。但是,我不确定如何以相同的形式执行此操作。这是我的代码:classTopic:destroyaccepts_nested_attributes_for:postsvalidates_presence_of:titleendclassPost...但这似乎不起作用。有什么想法吗?谢谢! 最佳答案 @Pablo的回答似乎有你需要的一切。但更具体地说...首先改变你View中的这一行对此#