美文网首页
亚组分析交互作用的P值计算及其可视化

亚组分析交互作用的P值计算及其可视化

作者: 灵活胖子的进步之路 | 来源:发表于2022-03-19 22:51 被阅读0次

首先从一篇JCO 的图表说起

举例文章

文章的目的是为比较HAIC与TACE对于大肝癌的治疗效果,研究设计为前瞻性随机对照三期临床试验,最后获得了一个阳性结果,既HAIC疗效优于TACE。

在得到阳性结果后,为了进一步观察结论在不同亚组中是否和总人群一样,既结果的稳健性。作者进行了亚组分析并绘制了森林图,如下图

举例文章的亚组分析森立图

可以从图中看到,在年龄的亚组分析中,虽然小于50岁和大于50岁患者都可以从HAIC治疗获益(HR值小于1且95%CI不过1),但是其交互作用p for interaction为0.01。HAIC治疗治疗在不同年龄分组之间的差异是有统计学意义,换句话说,就是说HAIC治疗对于年青组获益更多。

亚组分析的交互作用P值可以确定干预在不同亚组的区分是否有统计学意义,这种分析可以从一个侧面反应研究结果的可靠性。另外,其可以用于探索某个干预的浅在获益人群。

以下利用R语言的QualInt计算这种交互作用的P值

R包官方说明 形式参数 形式参数定义

以下进行代码举例

以下绘制连续变量为结果变量

##以下绘制连续变量为结果变量
ynorm <- rnorm(300)
trtment <- sample(c(0, 1), 300, prob = c(0.4, 0.6),replace = TRUE)
subgrp <- sample(c(0, 1, 2), 300, prob = c(1/3, 1/3, 1/3),replace = TRUE)
test1 <- qualint(ynorm, trtment, subgrp)
plot(test1)
print(test1)
连续变量绘图结果
连续变量统计学结果

以下进行二分类变量的计算

##以下进行二分类变量的计算
ybin <- sample(c(0, 1), 300, prob = c(0.3, 0.7),
               replace = TRUE)
test5 <- qualint(ybin, trtment, subgrp,type = "binary",
                 scale = "OR", test = "IBGA")

plot(test5)
print(test5)
二分类变量结果

以下进行生存变量为结果变量计算

##以下进行生存变量为结果变量计算
time <- rpois(300, 200)
censor <- sample(c(0, 1), 300, prob = c(0.7, 0.3),
                 replace = TRUE)
test7 <- qualint(Surv(time, censor), trtment, subgrp,
                 type = "survival", test = "IBGA")
plot(test7)
print(test7)
结果图

相关文章

  • 亚组分析交互作用的P值计算及其可视化

    首先从一篇JCO 的图表说起 文章的目的是为比较HAIC与TACE对于大肝癌的治疗效果,研究设计为前瞻性随机对照三...

  • Qiime2 (六)-Alpha和beta多样性分析

    计算核心多样性 Alpha多样性组间显著性分析和可视化

  • P值转化为Padj

    把计算出来的储存成一个数组,然后调用函数计算, 一定是所有p值(要计算的p值集合,不能单独只计算一个p值),因为p...

  • Python 资源大全(七)

    目录: 用来进行科学计算和数据分析的库进行数据可视化的库计算机视觉库机器学习库MapReduce 框架和库使用 P...

  • 11统计基础- 计算P值

    计算P值 有两种p值,单侧和双侧。双侧p值是最常见的,重点是计算它们。相反,单侧p值是很少被使用。 p-value...

  • T-test——T检验

    前面几节内容,我们了解了在回归分析中,如何判断变量之间的相关性——计算R2,如何判断相关的真实性——计算F值和P值...

  • Oracle开窗函数笔记及应用场景

    介绍Oracle的开窗函数之前先介绍一下分析函数,因为开窗函数也属于分析函数 分析函数用于计算基于组的某种聚合值,...

  • over (Partition by...) of oracle

    over()分析函数用于计算基于组的某种聚合值,它和聚合函数的不同之处是:对于每个组返回多行,而聚合函数对于每个组...

  • 卡方检验P值 两组间百分比差异性

    最近有项目需要计算p值,比较两组百分比的p值。 比较 分别在药物情况A、药物情况B时,BD I 与 BD II 之...

  • ggplot2自定义添加P值

    数据分析中通常都需要对数据进行统计检验,仿佛加个P值就能显得论文数据更加专业一点,本节不介绍如何计算P值,我们来介...

网友评论

      本文标题:亚组分析交互作用的P值计算及其可视化

      本文链接:https://www.haomeiwen.com/subject/ubbndrtx.html