在前面的章节中,我们已经探讨了2种SCI单因素表的制作方法,今天我们来将第三种表,其实这三种表已经涵盖了绝大部分的SCI的单因素表,只要您有心,绝对可以做出来。今天我们来看看这第三种,这种就像是第二种的加强版,先把产妇分为未生育的和已经生育的,再在里面比较两个人群前置胎盘患者的关系。

我们还是使用既往的乳腺癌的数据
age表示年龄,pathsize表示病理肿瘤大小(厘米),lnpos表示腋窝淋巴结阳性,histgrad表示病理组织学等级,er表示雌激素受体状态,pr表示孕激素受体状态,status结局事件是否死亡,pathscat表示病理肿瘤大小类别(分组变量),ln_yesno表示是否有淋巴结肿大,time是生存时间,后面的agec是我们自己设定的,不用管它。
library(foreign)
library(survival)
bc <- read.spss(“E:/r/Breast cancer survival agec.sav”,
use.value.labels=F, to.data.frame=T)

删除缺失值,查看数据
bc <- na.omit(bc)
head(bc)
这是数据格式

我们先把是否有淋巴结转移,分成两个数据表,再分别进行分析
ln_yesno0<-subset(bc,bcKaTeX parse error: Expected 'EOF', got '#' at position 12: ln_yesno<1)#̲我们告诉subset 函数在l…ln_yesno>=1)
这样我们就得到了ln_yesno0,ln_yesno1这两个分类数据集,我们就可以用它来分别进行统计了。
我们把年龄分为3个类别,小于40岁,40-60岁,大于60岁的,比较在这三个年龄区间孕激素状态对对无淋巴结转移人群对死亡结局的影响
先来做一个小于40岁人群中,孕激素状态(pr)在无淋巴结转移人群对死亡结局的影响
建立模型
f.age1 <- glm(status ~ pr, family = binomial(link = “logit”), data = ln_yesno0,subset = (age<=40))
对模型解析
summary(f.age1)

已经得出P值和标准误,已经可以做表了,如果你要求更高一点,还可以继续求95%CI和OR(参看我原来的帖子)
再来来做一个小于40岁人群中,孕激素状态(pr)在有淋巴结转移人群对死亡结局的影响(注意这里就要换一个数据集)
f.age2<- glm(status ~ pr, family = binomial(link = “logit”), data = ln_yesno1,subset = (age<=40))

再来一个做一个小于40-60岁人群中,孕激素状态(pr)在无淋巴结转移人群对死亡结局的影响
f.age3 <- glm(status ~ pr, family = binomial(link = “logit”), data = ln_yesno0,subset = (age>40 & age<=60))
再来一个做一个小于40-60岁人群中,孕激素状态(pr)在有淋巴结转移人群对死亡结局的影响(换个数据集而已)
f.age4 <- glm(status ~ pr, family = binomial(link = “logit”), data = ln_yesno0,subset = (age>40 & age<=60))

这样就可以依次做出表中所有的数据了,其实代码和原来很多都是重复,也不难做,但是你的思路一定要清晰,明白它是怎么做出来的。我们已经把大部分的单因素表格做出来了,下次将继续介绍怎么做SCI的多因素表。

手把手教你使用R语言做出SCI论文中的表二(单因素分析表)(3)相关推荐

  1. 手把手教你使用R语言做出SCI论文中的表二(单因素分析表)(2)

    在前面的章节,我们已经用R语言分析了一种SCI单因素表格的制作方法,今天我们来继续分析,这种分类的,其实方法都是差不多的 我们先要对表进行分析,要看它表达的是什么意思,其实就是通过回归分析对子集的分析 ...

  2. 手把手教你使用R语言做出SCI论文中的表二(单因素分析表)(1)

    在SCI论文中,我们经常可以看见一些这样的表格,大多数命名表格2,主要用来表示原因和结果的单因素分析的关系或者是分组变量的关系,如下图 这样论文中的表格数不胜数,今天我们通过一个实例数据演示告诉大家怎 ...

  3. R语言绘制SCI论文中常见的箱线散点图,并自动进行方差分析计算显著性水平

    显著性标记箱线散点图 本篇笔记的内容是在R语言中利用ggplot2,ggsignif,ggsci,ggpubr等包制作箱线散点图,并计算指定变量之间的显著性水平,对不同分组进行特异性标记,最终效果如下 ...

  4. 手把手教你入门R语言--数据分析篇(1)——数据导入与基础绘图

    R语言入门系列课程(6) 前面已将R语言的基础全部交给大家了,这两节将教大家用R语言进行数据分析的基础.数据读入(数据源包括txt文件.Excel文件.MySLQ数据库.MongoDB数据库).数据写 ...

  5. 手把手教你使用R语言建立COX回归并画出列线图(Nomogram)

    列线图,又称诺莫图(Nomogram),它是建立在回归分析的基础上,使用多个临床指标或者生物属性,然后采用带有分数高低的线段,从而达到设置的目的:基于多个变量的值预测一定的临床结局或者某类事件发生的概 ...

  6. 独家 | 手把手教你学习R语言(附资源链接)

    作者:NSS 翻译:杨金鸿 术语校对:韩海畴 全文校对:林亦霖 本文约3000字,建议阅读7分钟. 本文为带大家了解R语言以及分段式的步骤教程! 人们学习R语言时普遍存在缺乏系统学习方法的问题.学习者 ...

  7. 超级干货 :手把手教你学习R语言(附资源链接)

    作者:NSS:翻译:杨金鸿:校对:韩海畴,林亦霖: 本文约3000字,建议阅读7分钟. 本文为带大家了解R语言以及分段式的步骤教程! 人们学习R语言时普遍存在缺乏系统学习方法的问题.学习者不知道从哪开 ...

  8. 手把手教你学习R语言

    本文为带大家了解R语言以及分段式的步骤教程! 人们学习R语言时普遍存在缺乏系统学习方法的问题.学习者不知道从哪开始,如何进行,选择什么学习资源.虽然网络上有许多不错的免费学习资源,然而它们多过了头,反 ...

  9. 手把手教你使用R语言制作临床决策曲线

    DCA(Decision Curve Analysis)临床决策曲线是一种用于评价诊断模型诊断准确性的方法,在2006年由AndrewVickers博士创建,我们通常判断一个疾病喜欢使用ROC曲线的A ...

最新文章

  1. IAR编译提示could not open file lnkstm8s003f3.icf
  2. 【HDU - 3466 】Proud Merchants(dp,背包问题,巧妙排序)
  3. 谷爱凌惊“险”一跳,最少价值10个亿!
  4. WSPSRV.EXE错误导致ISA中断
  5. 006 认识BeanNameAware
  6. 系统设计和机器学习算法
  7. keras调用自己训练的模型,并去掉全连接层
  8. 前端实现成绩管理系统
  9. 计算机里怎样打字,电脑知识:如何电脑快速打字
  10. python类和对象的应用:烤地瓜
  11. 怎么知道跟交换机互联的交换机_电话交换机日常怎么维护和保养?
  12. Python数据分析入门笔记10——简单案例练习(学生信息分析)
  13. 小白算法积累——单链表6#带头结点单链表+递增有序
  14. 51单片机+DS18B20+数码管显示+串口通讯+Proteus仿真
  15. 图片链接打不开检测工具-免费链接失败检测软件
  16. 详解跨境电商企业信息化建设的关键点:业财一体化
  17. 卡迪夫大数据专业排名_2019QS排名出炉啦!看看卡迪夫大学那些世界前百强的学科...
  18. git学习之Linux安装Git
  19. 深入Spring底层透析Bean创建过程之拨云见日篇
  20. WinCC 7.0 安装及使用过程部分注意事项

热门文章

  1. bde连接oracle失败,BDE联接出错,求助
  2. JVM|记一次生产环境 CPU 占用飙高问题解决
  3. Android Hook 实战--替换第三方sdk的类
  4. 计算机研究生学什么课本,计算机研究生到底该怎么读?
  5. linux系统fasta程序,Linux生信练习2--fastq/fasta
  6. Photoshop设计精讲精练笔记
  7. 四位电脑故障诊断卡使用说明
  8. python webpy启动服务
  9. git clone加速(实测推荐)
  10. AI独角兽泡沫破裂:高管离职、数据打架、上市中止、营收崩溃