欢迎光临散文网 会员登陆 & 注册

R语言:autoReg包,自动输出整洁的回归模型结果

2023-04-13 15:32 作者:郑老师妙趣横生统计学  | 我要投稿

2023年以来浙中医大学郑老师开设了一系列医学科研统计课程,零基础入门医学统计包括R语言、meta分析、临床预测模型、真实世界临床研究、问卷与量表分析、医学统计与SPSS、临床试验数据分析、重复测量资料分析、结构方程模型、孟德尔随机化等10门课,如果您有需求,不妨点击下方跳转查看:

2023年郑老师多门科研统计课程:多次直播,含孟德尔随机化方法 

回归建模是论文中常见的统计方法,今天来学习一个R包——autoReg,使用这个R包可以快速输出回归模型的统计结果。

1. 安装和加载R包

可以从CRAN上安装R包。

install.packages("autoReg")
library(autoReg)

2.加载数据

使用survival包的colon数据集进行演示。

data(colon, package="survival")
colon

这个数据集收集了B/C 期结肠癌患者辅助化疗后的生存时间数据。

数据集中的变量解释:

id # 患者编号
study # 所有患者都是1
rx # 表示治疗方式,有三种:观察、Levamisole、Levamisole + 5-FU
sex # 性别,男性为1,女性为0
age # 年龄
obstruct # 肿瘤是否阻塞结肠,1 为有,0 为无
perfor # 结肠是否穿孔,1 为有,0 为无
adhere # 肿瘤是否粘附邻近器官,1 为有,0 为无
nodes # 检出淋巴结的数目
status # 生存状态,1 为发生感兴趣终点事件,0 为删失
differ # 肿瘤的分化程度(1=well, 2=moderate, 3=poor)
extent # 局部转移程度(1=submucosa, 2=muscle, 3=serosa, 4=contiguous structures)
surg # 从手术到登记注册的时间(0=short, 1=long)
node4 # 超过4 个阳性淋巴结
time # 直至发生感兴趣终点事件或删失的时间
etype # 事件类型: 1= 复发,2= 死亡

查看数据类型。

str(colon)

有些数字型变量其实是分类变量,需要先转换数据类型。

转换后构建新的数据集。

mycolon <- colon %>% # 创建新数据集新变量
  transmute(time,
            status,
            Age = age,
            Sex = factor(sex, levels = c(0, 1),
                         labels = c("Female", "Male")),
            Obstruct = factor(colon$obstruct),
            Differ = factor(colon$differ),
            Extent = factor(colon$extent))
str(mycolon) # 查看数据集结构

3.构建逻辑回归模型

使用glm()函数构建逻辑回归模型,可以使用summary()函数输出模型的摘要信息。

fit <- glm(status ~ Age + Sex + Obstruct + Differ + Extent,
           data=mycolon,
           family="binomial")
summary(fit)

详情请点击下方:

https://mp.weixin.qq.com/s?__biz=MzAwOTYyMDY3OQ==&mid=2650405111&idx=4&sn=d5a4f5991d7835811522d50cf6954e98&chksm=8351855fb4260c49bcb7130ec5f767d75b157322fc410bf5cdc044926139f201cb6b8032d5fc&token=482467525&lang=zh_CN#rd

vx关注“医学论文与统计分析”,获取更多精彩内容!

2023年统计服务

2023年,我们将开展从科研设计、数据分析、统计学报告等医学科研研究方法咨询与服务多项服务,若您有课题经费可以支持,欢迎您提前和我们联系,2022底前采用预付方式与我们开展合作。

2023年统计服务开启!欢迎提前洽谈数据分析、科研合作服务 

R语言:autoReg包,自动输出整洁的回归模型结果的评论 (共 条)

分享到微博请遵守国家法律