手把手教你使用R语言做出一篇20分SCI文章的交互效应表(p for Interaction)
高分SCI中几乎必备有交互效应分析表,通常为文章的表三或者表五,上一篇文章我们已经讲过怎么使用SPSS做出交互效应分析,今天我们R手把手实战做出下图这篇20分的文章的交互效应表。
我们先要看他是怎么分析的,他把血小板进行4分组和2分组后,在进行血小板和叶酸的交互作用分析,从各个层面分析了血小板对叶酸的影响,所以别人能发高分不是没有原因的。在模型中,他还使用了调整模型和未调整模型进行比较,支持自己的证据力度。我们还是沿用上次的数据来进行分析
把数据导入并查看数据
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)
假如我们研究的是乳腺癌各个pr(孕激素受体状态)和死亡的关系,我们想知道pr(孕激素受体状态)和年龄是否产生交互作用。
先使用CUT函数进行分割
age2<-cut(bc$age,breaks=c(0,40,60,100),include.lowest=T,
labels = c(1,2,3))#把age划分为0-20,20-60,60到100这样3个区间
dd<-cbind(bc,age2)#把变量加入表格
刚才是3分组的,再来一个2分组的
Age3<-cut(bc$age,breaks=c(0,40, 100),include.lowest=T,labels = c(1,2))
df<-cbind(bc,age2)
然后通过函数subset和if函数设定子集,每个分类变量都要设定一个子集
我们先来做3分类变量的age2,它的3个标签是,1,2,3,简单设置一下
目前age2是因子变量,我们把它转换成数字的
n.age2<-as.numeric(as.character(df$age2))
df<-cbind(df,n.age2)
然后再对它取其中是1的子集
agec1<-subset(df,df$n.age2<=1) #就是age2中取1的子集
得出这个子集后就对他进行分析
先来个广义线性方程得个回归系数,
f1<-glm(status ~pr,family = binomial,)
summary(f1)
计算OR和95%CI
exp(confint(f1))
exp(coef(f1))
由此已经算出年龄0-40岁的亚组的OR,95%CI和P值,同理可以算出其他的组
现在我们来做2分组的交互效应
交互效应就不能拿子集了,要用原来的数据
f2<-glm(status ~pr+age3,family = binomial,data =df)
summary(f2)
求交互效应其实很简单,加入一个a*b就可以了,箭头所示就是交互的P值了
f3<-glm(status ~pr+age3+pr*age3,family = binomial,data =df)
调整模型就是在模型中加入其他变量就可以了,我这里就调整了er,其实做出来用到R的代码并不难,关键是你要思路清晰,明白它是怎么做出来的,
动动小手关注一下吧,更多精彩文章尽在零基础说科研
手把手教你使用R语言做出一篇20分SCI文章的交互效应表(p for Interaction)相关推荐
- 手把手教你使用R语言做出SCI论文中的表二(单因素分析表)(3)
在前面的章节中,我们已经探讨了2种SCI单因素表的制作方法,今天我们来将第三种表,其实这三种表已经涵盖了绝大部分的SCI的单因素表,只要您有心,绝对可以做出来.今天我们来看看这第三种,这种就像是第二种 ...
- 手把手教你使用R语言做出SCI论文中的表二(单因素分析表)(2)
在前面的章节,我们已经用R语言分析了一种SCI单因素表格的制作方法,今天我们来继续分析,这种分类的,其实方法都是差不多的 我们先要对表进行分析,要看它表达的是什么意思,其实就是通过回归分析对子集的分析 ...
- 手把手教你使用R语言做出SCI论文中的表二(单因素分析表)(1)
在SCI论文中,我们经常可以看见一些这样的表格,大多数命名表格2,主要用来表示原因和结果的单因素分析的关系或者是分组变量的关系,如下图 这样论文中的表格数不胜数,今天我们通过一个实例数据演示告诉大家怎 ...
- 手把手教你入门R语言--数据分析篇(1)——数据导入与基础绘图
R语言入门系列课程(6) 前面已将R语言的基础全部交给大家了,这两节将教大家用R语言进行数据分析的基础.数据读入(数据源包括txt文件.Excel文件.MySLQ数据库.MongoDB数据库).数据写 ...
- 手把手教你使用R语言建立COX回归并画出列线图(Nomogram)
列线图,又称诺莫图(Nomogram),它是建立在回归分析的基础上,使用多个临床指标或者生物属性,然后采用带有分数高低的线段,从而达到设置的目的:基于多个变量的值预测一定的临床结局或者某类事件发生的概 ...
- 独家 | 手把手教你学习R语言(附资源链接)
作者:NSS 翻译:杨金鸿 术语校对:韩海畴 全文校对:林亦霖 本文约3000字,建议阅读7分钟. 本文为带大家了解R语言以及分段式的步骤教程! 人们学习R语言时普遍存在缺乏系统学习方法的问题.学习者 ...
- 超级干货 :手把手教你学习R语言(附资源链接)
作者:NSS:翻译:杨金鸿:校对:韩海畴,林亦霖: 本文约3000字,建议阅读7分钟. 本文为带大家了解R语言以及分段式的步骤教程! 人们学习R语言时普遍存在缺乏系统学习方法的问题.学习者不知道从哪开 ...
- 手把手教你学习R语言
本文为带大家了解R语言以及分段式的步骤教程! 人们学习R语言时普遍存在缺乏系统学习方法的问题.学习者不知道从哪开始,如何进行,选择什么学习资源.虽然网络上有许多不错的免费学习资源,然而它们多过了头,反 ...
- 手把手教你使用R语言制作临床决策曲线
DCA(Decision Curve Analysis)临床决策曲线是一种用于评价诊断模型诊断准确性的方法,在2006年由AndrewVickers博士创建,我们通常判断一个疾病喜欢使用ROC曲线的A ...
- 独家 | 手把手教你用R语言做回归后的残差分析(附代码)
作者:Abhijit Telang 翻译:张睿毅 校对:丁楠雅 本文约2600字,建议阅读10分钟. 本文介绍了做残差分析的方法及其重要性,以及利用R语言实现残差分析. 在这篇文章中,我们通过探索残差 ...
最新文章
- arduino项目_quot;Arduino智能车项目quot;实战班开班!
- 关于使用fastjson统一序列化响应格式。
- RPA如何助力电信业数字化转型?
- html中图片的属性优化,Html标签元素在SEO中的优化方式(二)
- 获取系统信息1——linux系统中的时间
- 构思新巧的173dmba飞鸽
- CVPR 2019 | 微软亚研院提出用于语义分割的结构化知识蒸馏
- c语言添加收支情况,C语言编写一个计算个人所得税的程序,要求输入收入金额,能够输...
- 序列化_SpreadJS序列化与反序列化表格Serialization+Deserialization
- python如何将天数转换为日期字符串
- jQuery实现倒计时重新发送短信验证码功能示例
- 个人知识管理能解决什么问题?
- (PTA)java基础编程题汇总
- Android HttpURLConnection 获取下载文件大小
- ios 判断打开相机权限_iOS中判断照片和相机权限
- 洛谷P3537 [POI2012]SZA-Cloakroom(背包)
- Win10的Linux子系统Ubuntu安装图形界面
- ZYNQ学习之旅--PS_AXI_VDMA(利用VDMA实现将PS端的数据显示在PL端的HDMI上)
- win7 dll怎么在xp运行_微信DLL劫持反弹shell复现
- Java基础知识之静态