当前位置：首页 > news >正文

Time to event ：Kaplan-Meier曲线、Log Rank检验与Shiny R

news 来源：原创 2025/7/28 7:50:57

代码：

# 创建数据框
data_a <- data.frame(
usubjid = c('1- 1', '1- 2', '1- 3', '1- 4', '1- 5',
'1- 6', '1- 7', '1- 8', '1- 9', '1-10',
'2- 1', '2- 2', '2- 3', '2- 4', '2- 5',
'2- 6', '2- 7', '2- 8', '2- 9', '2-10'),
cnsr = c(0,1,0,1,0,1,0,0,0,1,
1,0,1,0,1,0,0,1,0,1),
time = c(5,7,9,1,10,2,9,11,8,4,
11,11,9,11,10,12,9,9,9,13),
arm = c(rep('A', 10), rep('B', 10))
)

# 查看数据结构
str(data_a)

# 显示前6行
head(data_a)

# 检查观测数
nrow(data_a) # 应返回20

# 检查分组比例
table(data_a$arm) # 应显示 A:10, B:10

# 加载必要的包
library(survival)
library(survminer) # 用于可视化（可选）

# 创建生存对象：事件定义为cnsr=0（需反转cnsr）
Surv_obj <- Surv(time = data_a$time, event = 1 - data_a$cnsr)

# 分层分析（按arm分组）并进行Log-Rank检验
logrank_test <- survdiff(Surv_obj ~ arm, data = data_a)

# 提取原始p值并应用Bonferroni校正（仅一次比较，校正后p值不变）
p_value <- 1 - pchisq(logrank_test$chisq, df = 1)
adjusted_p_value <- p.adjust(p_value, method = "bonferroni")

# 输出Log-Rank检验结果
cat("Log-Rank Test Results:\n")
cat("----------------------\n")
cat("Chi-Square Statistic:", logrank_test$chisq, "\n")
cat("Degrees of Freedom:", logrank_test$df, "\n")
cat("Raw p-value:", p_value, "\n")
cat("Bonferroni-adjusted p-value:", adjusted_p_value, "\n\n")

# 拟合Kaplan-Meier曲线（对照组为B组）
km_fit <- survfit(Surv_obj ~ arm, data = data_a)

# 输出生存曲线摘要（含中位生存时间及置信区间）
cat("Kaplan-Meier Survival Summary:\n")
cat("------------------------------\n")
print(summary(km_fit))

# 加载包
library(survival)
library(survminer) # 确保已加载！

# 绘制生存曲线
ggsurvplot(km_fit,
data = data_a,
pval = TRUE, # 显示p值
conf.int = TRUE, # 显示置信区间

相关文章：

Time to event ：Kaplan-Meier曲线、Log Rank检验与Shiny R

代码： # 创建数据框 data_a <- data.frame( usubjid = c(1- 1, 1- 2, 1- 3, 1- 4, 1- 5, 1- 6, 1- 7, 1- 8, 1- 9, 1-10, 2- 1, 2- 2, 2- 3, 2- 4, 2- 5, 2- 6, 2- 7, 2- 8, 2- 9, 2-10), cnsr = c(0,1,0,1,0,1,0,0,0,1,…...

编程日记 2025/7/28 7:50:57

线上救急-AWS限频

线上救急-AWS限频问题在一个天气炎热的下午，我正喝着可口可乐，悠闲地看着Cursor生成代码，忽然各大群聊中出现了加急➕全体的消息，当时就心里一咯噔，点开一看，果然，线上服务出问题&#xff0…...

编程日记 2025/7/27 5:30:59

JavaWeb学习打卡-Day1-分层解耦、Spring IOC、DI

三层架构 Controller（控制层）：接收前端发送的请求，对请求进行处理，并响应数据。Service（业务逻辑层）：处理具体的业务逻辑。DAO（数据访问层/持久层）&#xff…...

编程日记 2025/7/27 5:30:57

【LeetCode】1.两数之和

目录 📚 题目概要🧰 前置知识🚧 问题难点🔑 关键思路步骤拆解 💻 代码实现代码注释 📊 复杂度分析❗ 易错点与测试案例易错点测试案例 🔗 总结与扩展模式归纳核心思维 📚 题目概要在…...

编程日记 2025/7/27 5:30:55

mongodb 存储数据的具体实现方式

MongoDB 存储数据的具体实现方式涉及数据模型、存储引擎、分片机制等多个核心模块，以下是其实现原理的详细分析： 一、数据模型 1.1 文档型数据模型‌ MongoDB 使用 BSON格式存储数据，支持键值对、嵌套文档和数组等复杂结构。 1.2 无模式设…...

编程日记 2025/7/28 7:50:02

【手机】vivo手机应用声音分离方案

文章目录前言方案前言尝试分离vivo手机音乐与其他应用的声音方案最佳方案：网易云音乐设置内关闭音量均衡上传不同的白噪音，成功 goodlock，主要适用于三星手机，vivo不一定适用 app volume control ，可行...

编程日记 2025/7/27 5:30:53

多级缓存架构，让系统更快的跑起来！

大家好，今天，咱们来聊聊一个超级实用的话题——多级缓存架构。别一听“架构”俩字就头大，我保证，这篇文章既有趣又易懂，让你秒变缓存小达人！一、多级缓存，为啥这么火？在互联网的汪洋大海里，数据就是咱们的宝藏。但每次从数据库里捞数据，都跟挖宝藏似的，慢得很！…...

编程日记 2025/7/27 5:30:51

Vibracostic EDI 需求分析

Vibracostic 是德国Freudenberg集团旗下全球领先的减振与噪音控制技术公司，专注于为汽车及工业领域提供高效振动管理和隔音解决方案，客户涵盖宝马、奔驰、特斯拉等主流车企。 Vibracostic EDI 需求分析供应商接收Vibracostic发来的DELFOR交付预测报文…...

编程日记 2025/7/28 7:50:57

基于超启发鲸鱼优化算法的混合神经网络多输入单输出回归预测模型 HHWOA-CNN-LSTM-Attention

基于超启发鲸鱼优化算法的混合神经网络多输入单输出回归预测模型 HHWOA-CNN-LSTM-Attention 随着人工智能技术的飞速发展，回归预测任务在很多领域得到了广泛的应用。尤其在金融、气象、医疗等领域，精确的回归预测模型能够为决策者提供宝贵的参考信息。为…...

编程日记 2025/7/28 7:48:22

Linux卸载删除gitlab

1、停止 gitlab服务 gitlab-ctl stop 2、卸载 gitlab（社区版） rpm -e gitlab-ce 或者 yum remove gitlab-ce 3、查看 gitlab 进程 ps aux | grep gitlab 4、杀掉gitlab service进程，该进程与runsvdir相关（带有好多..........…...

编程日记 2025/7/27 5:30:46

高品质性价比之王-特伦斯便携钢琴V10

在电子钢琴选购过程中，预算与品质的平衡常常让消费者感到纠结。但特伦斯 V10 88 键可折叠电子钢琴的出现，为广大音乐爱好者带来了惊喜，亲民的价格实现了高品质的音乐体验。先看便携性，同价位的电子钢琴大多体型庞大&#xff0c…...

编程日记 2025/7/27 7:10:18

解决方案评测｜告别复杂配置！基于阿里云云原生应用开发平台CAP快速部署Bolt.diy

写在前面的话突然看到上线了关于Bolt.new开源版本的解决方案测评，其实心里还是挺高兴的，我最早接触到Bolt.new的时候应该是在去年的11月份，当时是撰写了一篇名为一种基于通义千问prompt辅助Qwen2.5-coder-32bBolt.newv0Cursor的无代码对话…...

编程日记 2025/7/27 7:10:16

python测试框架之pytest

Pytest pytest 基础使用pytest安装pytest的测试case收集规则pytest - fixture的使用skip and xfailpytest - 属性标记测试函数pytest - 参数化测试pytest - mock/monkeypatch的使用pytest - 运行方式pytest - 运行方式/命令pytest - 处理测试失败的case pytest - 测试输出捕获 …...

编程日记 2025/7/27 7:10:13

uni-app 开发企业级小程序课程

课程大小：7.7G 课程下载：https://download.csdn.net/download/m0_66047725/90616393 更多资源下载：关注我备注：缺少两个视频5-14 tabs组件进行基本的数据展示和搜索历史处理searchData的删除操作 1-1导学.mp4 2-10小程序内…...

编程日记 2025/7/27 7:10:11

深度图可视化

import cv2# 1.读取一张深度图 depth_img cv2.imread("Dataset_depth/images/train/1112_0-rgb.png", cv2.IMREAD_UNCHANGED) print(depth_img.shape) cv2.imshow("depth", depth_img) # (960, 1280) print(depth_img)# 读取一张rgb的图片做对比 input_p…...

编程日记 2025/7/27 7:10:09

Java实现希尔排序算法

1. 希尔排序原理图解希尔排序是插入排序的一种高效改进版本，通过比较和交换间隔较远的元素来减少数据的移动次数。以下是希尔排序的步骤： 1. 选择初始间隔：通常选择数组长度的一半作为初始间隔。 2. 分组和插入排序：将数组分成若…...

编程日记 2025/7/27 7:10:07

Apache RocketMQ 荣获 2024 开源创新榜单“年度开源项目

持续深耕 RocketMQ 领域，更多实战避坑指南、源码解析、最佳实践，已整理至「RocketMQ 中文社区」 🔍 近日，以“新纪天工、开物焕彩——致敬开源的力量”为活动主题的“重大科技成就发布会（首场）”在国家科技…...

编程日记 2025/7/27 7:10:05

2025年智慧交通与能源技术国际会议（ITNET 2025）

The 2nd International Conference on Intelligent Transportation and Energy Technology 一、大会信息会议简称：ITNET 2025 大会时间：(以官网为准） 大会地点：中国南京收录检索：提交Ei Compendex,CPCI,CNKI,Google…...

编程日记 2025/7/27 7:10:03

Spring Boot 集成 Redis 实战总结

Spring Boot 集成 Redis 实战总结一、基础集成与配置依赖引入 <dependency><groupId>org.springframework.boot</groupId><artifactId>spring-boot-starter-data-redis</artifactId> </dependency>默认使用 Lettuce 作为连接池&#xf…...

编程日记 2025/7/26 13:45:29

ubuntu18.04启动不了修复

参考: 虚拟机里的Ubuntu18.4启动时进入到grub rescue救援模式（无法正常进入到系统），ls查看后只有一个硬盘和分区，且无法找到/boot/grub文件【已解决】_ubuntu grub rescue-CSDN博客本人fdisk错误使用,导致了grub启动不了第一步…...

编程日记 2025/7/27 7:09:59

人类行为的原动力是自我保存-来自ChatGPT

自我保存（Self-Preservation）确实可以说是人类行为最原始、最底层的驱动力。简单来说： 无论我们做什么，表面看动机五花八门，实际上归根到底都绕不开活下去、保护自己。 💡 从不同层面理解这个观点&#…...

编程日记 2025/7/27 7:09:57

退役淘汰的硬盘数据安全处置不可忽视-硬盘数据抹除清零

在信息化时代，硬盘作为数据存储的核心载体，其退役处理直接关系到信息安全，淘汰硬盘中往往留存大量敏感信息，若处置不当，极易引发数据泄露风险，造成难以挽回的损失。因此，建立规范的硬盘数据销毁…...

编程日记 2025/7/27 7:09:55

创建 Node.js Playwright 项目：从零开始搭建自动化测试环境

一、环境准备在开始创建 Playwright 项目之前，确保你的电脑上已经安装了以下工具： Node.js：Playwright 依赖于 Node.js 环境，确保你已经安装了最新版本的 Node.js。可以通过以下命令检查是否安装成功： node -v npm -…...

编程日记 2025/7/27 7:09:53

1+X： Python程序开发职业技能等级要求（初级）练习资料分享

以下将根据表1中Python程序开发职业技能等级要求（初级）的不同工作领域和任务，分别给出练习代码或操作步骤。 Python基础编程开发环境搭建操作步骤：以Windows系统为例，首先从Python官方网站（https://www…...

编程日记 2025/7/27 7:09:50

【漏洞复现】CVE-2024-38856（ApacheOfbiz RCE）

【漏洞复现】CVE-2024-38856（ApacheOfbiz RCE） 1. 漏洞描述 Apache OFBiz 是一个开源的企业资源规划（ERP）系统。它提供了一套企业应用程序，用于集成和自动化企业的许多业务流程。这个漏洞是由于对 CVE-2023-51467 的…...

编程日记 2025/7/27 7:09:48

【React】获取元素距离页面顶部的距离

文章目录代码实现代码实现 import { useEffect, useRef, useState } from react;const DynamicPositionTracker () > {const [distance, setDistance] useState(0);const divRef useRef(null);useEffect(() > {const targetDiv divRef.current;if (!targetDiv) re…...

编程日记 2025/7/27 7:09:46

【Markdown】【HTML】在Markdown中实现康奈尔笔记模式（右侧留白）

文章目录 HTML表格——康奈尔笔记法标签解释<table> 相关标签文本内容标签注释标签最低限度应了解的HTML标签 (Essential HTML Tags)结构标签 (Structure Tags)文本标签 (Text Tags)列表标签 (List Tags)链接和图片 (Links and Images)格式化 (Formatting)表格 (Tables)…...

编程日记 2025/7/27 7:09:44

Go语言中的Context

目录 Go语言中的Context 1. Context的基本概念 1.1 Context的核心作用 2. Context的基本用法 2.1 创建Context 背景Context 可取消的Context 带有超时的Context 2.2 在Goroutine间传递Context 2.3 获取Context的值为Context添加自定义数据访问Context中的值 3. C…...

编程日记 2025/7/27 7:09:42

基于Python的图片/签名转CAD小工具开发方案

基于Python的图片/签名转CAD工具开发方案一、项目背景传统设计流程中，设计师常常需要将手写签名或扫描图纸转换为CAD格式。本文介绍如何利用Python快速开发图像矢量化工具，实现： 📷 图像自动预处理✏️ 轮廓精确提取⚙️ 参数…...

编程日记 2025/7/27 7:09:40

uniapp微信小程序：WIFI设备配网之TCP/UDP开发AP配网

一、AP配网技术原理 1.1 配网模式选择 AP配网（SoftAP模式）是IoT设备配网成功率最高的方案之一 1、其核心原理： 设备端：启动AP模式（如SSID格式YC3000_XXXX，默认IP192.168.4.1）手…...

编程日记 2025/7/27 7:09:38

docker部署seata

seata版本：2.0.0 使用eureka作为注册中心，mysql作为存储拉取镜像 docker pull seataio/seata-server:2.0.0mysql建表语句可从 github下载sql脚本 -- seata_server.branch_table definition CREATE TABLE branch_table (branch_id bigint NOT NULL,…...

编程日记 2025/7/27 7:09:36

JVM 内存分布详解

/*** 示例类，展示各种成员在JVM中的存储位置*/ public class MemoryDistributionExample {// 1. 实例变量（普通成员变量） - 存储在堆内存的对象实例中private int instanceVar 10;private String instanceRef "Instance Reference&quo…...

编程日记 2025/7/27 7:09:34

AR行业应用案例与NXP架构的结合

1. 工业巡检AR头盔场景示例：宁德核电基地使用AR智能头盔进行设备巡检，通过实时数据叠加和远程指导，将工作效率提升35%。头盔需处理传感器数据、图像渲染和低延迟通信1。 NXP架构支持： 协处理器角色：NXP i.MX RT系列M…...

编程日记 2025/7/28 7:20:37

OpenCV 图形API（49）颜色空间转换-----将 NV12 格式的图像数据转换为 BGR 颜色空间函数NV12toBGR()

操作系统：ubuntu22.04 OpenCV版本：OpenCV4.9 IDE:Visual Studio Code 编程语言：C11 算法描述将图像从NV12（YUV420p）颜色空间转换为BGR。该函数将输入图像从NV12颜色空间转换为RGB。Y、U和V通道值的常规范围是0到25…...

编程日记 2025/7/27 7:09:29

Transformer中Post-Norm和Pre-Norm如何选择？

pre-norm: x n 1 x n f ( n o r m ( x n ) ) x_{n1}x_nf(norm(x_n)) xn1xnf(norm(xn)) 其中第二项的方差由于有 norm 是不随层数变化的，于是 x 的方差会在主干上随层数积累。到了深层以后，单层对主干的影响可以视为小量，而不同层的 …...

编程日记 2025/7/27 7:09:27

jvm-描述符与特征签名的区别

在Java虚拟机（JVM）中，存储的是方法签名，而不是仅仅方法描述符。方法签名包含了方法的参数类型和返回值类型的信息，而方法描述符通常指的是仅包含参数类型的那部分信息。为了更清晰地理解这两者的区别以及它们如何在JVM…...

编程日记 2025/7/27 7:09:25

JVM有什么调优参数？

-Xms：堆最小值 -Xmx：堆最大值。-Xms与-Xmx 的单位默认字节都是以k、m做单位的。通常这两个配置参数相等，避免每次空间不足，动态扩容带来的影响。 -Xmn：新生代大小 -Xss：每个线程池的栈大小。在jdk5以上…...

编程日记 2025/7/27 7:09:23

QT：Qt5 串口模块 (QSerialPort) 在 VS2015 中正确关闭串口避免被占用

以下是使用 Qt5 串口模块 (QSerialPort) 在 VS2015 中正确关闭串口避免被占用的完整示例代码： #include <QSerialPort> #include <QDebug>// 创建全局或类成员变量（推荐使用智能指针） QSerialPort *serialPort nullptr; // 打开…...

编程日记 2025/7/27 7:09:21

14.QT-多元素控件|QListWidget|QTableWidget|QTreeWidget(C++)

Qt中提供的多元素控件有: QListWidgetQListViewQTableWidgetQTableViewQTreeWidgetQTreeView 以QTableWidget和QTableView为例. QTableView是基于MVC设计的控件.QTableView⾃⾝不持有数据.使⽤QTableView的时候需要⽤⼾创建⼀个Model对象(⽐如 QStandardModel ),并且把Model…...

编程日记 2025/7/27 7:09:19

单页面应用的特点,什么是路由,VueRouter的下载,安装和使用,路由的封装抽离,声明式导航的介绍和使用

文章目录一.什么是单页面应用?二.什么是路由?生活中的路由和Vue中的路由三.VueRouter(重点)0.引出1.介绍2.下载与使用(5个基本步骤2个核心步骤)2.1 五个基本步骤2.2 两个核心步骤一.什么是单页面应用? 定义和特点:所有功能都在一个HTML中实现示例:网易云音乐多页面应…...

编程日记 2025/7/27 7:09:17

Java并发编程-线程通讯

Java并发编程-线程通讯线程通讯常见场景实现方法等待通知信号量实现栅栏实现锁机制实现 CountDownLatchCountDownLatch 常用方法说明模拟报表统计场景其他样例 CyclicBarrier 线程通讯线程通讯指的是多个线程之间通过共享内存或消息传递等方式来协调和同步它们的执行。在多线…...

编程日记 2025/7/27 7:09:15

ubuntu学习day4

4 Linux 文件操作 4.1 基于文件指针的文件操作 Linux 中对目录和设备的操作都是文件操作，文件分为普通文件，目录文件，链接文件和设备文件。普通文件：也称磁盘文件，并且能够进行随机的数据存储(能够自由 seek 定位到某…...

编程日记 2025/7/26 7:19:26

百度暑期实习岗位超3000个，AI相关岗位占比87%，近屿智能携AIGC课程加速人才输出

今年3月，百度重磅发布3000暑期实习岗位，聚焦大模型、机器学习、自动驾驶等AI方向的岗位比例高达87%。此次实习岗位涉及技术研发、产品策划、专业服务、管理支持、政企解决方案等四大类别，覆盖超300个岗位细分方向。值得一提的是，百…...

编程日记 2025/7/27 7:09:11

vin码识别技术-车辆vin识别代码-Java接口集成

快速识别车辆vin码有助于提升汽车保险行业、4s店车辆出入库、二手车交易、车辆租赁公司等众多业务流程的管理。vin码识别接口，为汽车后市场实现车辆vin信息的快速提取提供了高效且可靠的解决方案。基于先进的ocr技术，能够迅速定位图片中的VIN码区域&…...

编程日记 2025/7/27 7:09:08

随机深林算法是分类还是回归？

随机森林算法既可以用于分类也可以用于回归。一、随机森林用于分类的情况原理随机森林是一种集成学习算法，它通过构建多棵决策树来提高预测性能。在分类任务中，每棵决策树都是基于随机选择的特征子集和随机抽样的数据子集（通常是通过自助…...

编程日记 2025/7/27 4:46:07

【现代深度学习技术】循环神经网络05：循环神经网络的从零开始实现

【作者主页】Francek Chen 【专栏介绍】 ⌈ ⌈ ⌈PyTorch深度学习 ⌋ ⌋ ⌋ 深度学习 (DL, Deep Learning) 特指基于深层神经网络模型和方法的机器学习。它是在统计机器学习、人工神经网络等算法模型基础上，结合当代大数据和大算力的发展而发展出来的。深度学习最重…...

编程日记 2025/7/26 10:39:27

北京大学｜第8讲：《生成未必理解：基于扩散模型能否实现视觉世界模型》

导读INTRODUCTION 随着人工智能技术的飞速发展，AI在图像和视频生成领域取得了显著成就，但生成和理解之间的统一仍然是一个挑战。北京大学深圳研究生院助理教授袁粒的团队正在探索如何实现这一目标，他们的研究可能会改变我们对AI生成和理解能…...

编程日记 2025/7/27 7:09:02

C++如何处理多线程环境下的异常？如何确保资源在异常情况下也能正确释放

多线程编程的基本概念与挑战多线程编程的核心思想是将程序的执行划分为多个并行运行的线程，每个线程可以独立处理任务，从而充分利用多核处理器的性能优势。在C中，开发者可以通过std::thread创建线程，并使用同步原语如std::mutex、…...

编程日记 2025/7/27 5:30:43

Gmssl实战

最近项目要引入gm算法，实现aes和sm4数据加密，研究了一下两个系统的安装和测试用例 linux，从github下载GmSSL库（ https://github.com/guanzhi/GmSSL.git ），以下为安装方法。 ## GmSSL库安装编译,需要确认…...

编程日记 2025/7/27 10:09:48

C语言函数(上)

目录 1.函数的概念 2.库函数 3.自定义函数 4.形参和实参 5.return语句 1.函数的概念数学中我们其实就见过函数的概念，比如：一次函数y kx b，k和b都是常数，给一个任意的x， 就得到一个y值。其实在C语言也引入…...

编程日记 2025/7/27 7:09:00