ggrcs 包2.4绘制RCS(限制立方样条图)实际操作演示(1)

news2025/1/17 3:52:21

ggrcs 包2.4版本已经发布一段时间了,大概几个月了吧,收到不少好评,
在这里插入图片描述
没听说太大的问题,最主要的问题有两个:
1.是说变量不是数字变量。
在这里插入图片描述
2.是说数据超过10万,无法处理
在这里插入图片描述
第一个问题非常好处理,这个问题是RMS包的报错,主要是数据格式必须是data.frame,不是这个格式容易报错,我们导入数据的时候最好以csv格式读入,这个格式比较稳定,不容易出错,如果你是其他格式导入的,可以使用as.data.frame()把它进行强制转换就可以了,在后面的版本中我也会对它进行强制转换。
第二个问题是数据量太大超过10万个,rms包的函数处理不了,就会报错,我们可以改一下公式给它生成一个模拟数据,画出来的RCS图和原始数据差不多,但是这样直方图就会和原始数据的不一样,这个暂时没有很好的办法解决,别人的R包我们改不了。
第三个其实不算问题,就是粉丝画出来的图不怎么好看,下图

library(ggrcs)
fit <- coxph(Surv(y1t,y1==1)~x1+age+gender,data = mydata)
ggrcs(mydata,fit,"x1")

在这里插入图片描述
造成这个原因是数据造成的,这个不是正态分布,而是泊松分布,分布过于集中所致,我们只要对它的x变量取个对数,转换一下,让它变成正态分布,从新做图,变成如下

mydata$x99<-log(mydata$x1)
fit2 <- cph(Surv(y1t,y1==1) ~ rcs(x99,4)+age+gender, x=TRUE, y=TRUE,data=mydata)
ggrcs(mydata,fit2,"x99",histbinwidth=0.5)

在这里插入图片描述
另外一个例子也是一模一样的,做的时候是下图,也不是正态分布
在这里插入图片描述
对x变量取对数以后,做图如下(代码就不贴了)
在这里插入图片描述
取对数后两个图的趋势都很明显,图也很漂亮,发个文章应该是没问题的。因此可以总结出,绘制直方图和RCS图,最好数据x是正态分布,如果偏态偏得太厉害,我们最好给它取个对数或者其他方法转换一下,让它变成近似正态分布。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/44691.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

【一文秒懂——Profile配置】

目录 1. Profile配置 2. 实例 1. Profile配置 Spring框架允许使用Profile配置&#xff0c;即某些“个性化配置文件”&#xff0c;这些配置文件默认并不会被应用&#xff0c;需要“激活”后才生效&#xff01; 在Spring Boot项目中&#xff0c;简化了Profile配置的使用&…

智能型电动机保护器(马保)可应用在水泥厂,保证生产过程的连续不间断

安科瑞 李可欣 具体可咨询&#xff1a;Acrel_lkx 一、水泥厂的要求 水泥生产是一项较庞大的工程&#xff0c;生产环境比较恶劣&#xff0c;粉尘污染严重&#xff0c;受生产工艺的限制使得其工作场所的通风条件差且环境温度较高&#xff0c;而且需要设定大量的工作站点&…

NetCore路由的Endpoint模式

IdentityServer里有各种Endpoint,如TokenEndpoint,UserInfoEndpoint,Authorize Endpoint,Discovery Endpoint等等。Endpoint从字面意思来看是“终端节点"或者“终节点”的意思。无独有偶NetCore的路由也有Endpoint的概念。那么我们提出一个问题来&#xff0c;究竟什么是En…

R语言MCMC:Metropolis-Hastings采样用于回归的贝叶斯估计

全文链接&#xff1a;http://tecdat.cn/?p19664 MCMC是从复杂概率模型中采样的通用技术。蒙特卡洛马尔可夫链Metropolis-Hastings算法&#xff08;点击文末“阅读原文”获取完整代码数据&#xff09;。问题如果需要计算有复杂后验pdf p&#xff08;θ| y&#xff09;的随机变量…

中远通创业板IPO过会:上半年营收7.3亿 拟募资2.3亿

雷递网 雷建平 11月28日深圳市核达中远通电源技术股份有限公司&#xff08;简称&#xff1a;“中远通”&#xff09;日前IPO过会&#xff0c;准备在深交所创业板上市。中远通计划募资2.3亿元。其中&#xff0c;1.3亿元用于研发中心改造提升项目&#xff0c;3248万元用于企业信息…

规则引擎Drools在贷后催收业务中的应用

作者&#xff1a;vivo 互联网服务器团队- Feng Xiang 在日常业务开发工作中我们经常会遇到一些根据业务规则做决策的场景。为了让开发人员从大量的规则代码的开发维护中释放出来&#xff0c;把规则的维护和生成交由业务人员&#xff0c;为了达到这种目的通常我们会使用规则引擎…

(十八)Spring6集成MyBatis3.5

文章目录实现步骤具体实现第一步&#xff1a;准备数据库表第二步&#xff1a;IDEA中创建一个模块&#xff0c;并引入依赖第三步&#xff1a;基于三层架构实现&#xff0c;所以提前创建好所有的包第四步&#xff1a;编写pojo第五步&#xff1a;编写mapper接口第六步&#xff1a;…

怎么建网站?【模版建站】

关于怎么建网站&#xff0c;除了公司企业&#xff0c;甚至有些个人用户都想了解。大家印象中的建站网站都是要会编程&#xff0c;不然就是找外包公司解决。其实现在建网站也是比较简单的&#xff0c;模版建站一般都能解决基本的建站需求。下面我们一起来看看怎么建网站吧。 一…

考阿里云ACE需要准备什么?考试内容难不难?

最近几年云计算技术发展得越来越好&#xff0c;市场上大多数企业已经选择转型&#xff0c;使用云计算技术来发展自己的新业务&#xff0c;这样一来就需要大量的人才来维持市场的运行。另一方面&#xff0c;为了在现在内卷的社会中的脱颖而出&#xff0c;获得一份稳定、高薪的工…

目标级联分析法( Analytical Target Cascading , ATC )理论matlab程序

目标级联分析法&#xff08; Analytical Target Cascading &#xff0c; ATC &#xff09;理论matlab程序 目标级联分析法&#xff08;Analytical Target Cascading&#xff0c;ATC&#xff09;是一种采用并行思想解决复杂系统的设计方法&#xff0c;最初由密执安大学研究人员…

微服务框架 SpringCloud微服务架构 4 Ribbon 4.1 负载均衡原理

微服务框架 【SpringCloudRabbitMQDockerRedis搜索分布式&#xff0c;系统详解springcloud微服务技术栈课程|黑马程序员Java微服务】 SpringCloud微服务架构 文章目录微服务框架SpringCloud微服务架构4 Ribbon4.1 负载均衡原理4.1.1 负载均衡流程4 Ribbon 4.1 负载均衡原理 …

基于STM32单片机电子相册设计全套资料

资料编号&#xff1a;188 功能介绍&#xff1a; 本系统采用STM32f103单片机通过SPI接口读取sd卡模块中的图片数据。并在单片机为sd卡模块生成fat文件系统。方便读取sd卡中的文件信息。将Bmp格式的图片存放到sd卡的picture文件夹中&#xff0c;然后单片机进行Bmp解码&#xff0…

GEE批量下载 Python本地快速下载GEE数据(比网页版保存到网盘再下载快几十倍,尤其是在下载几十年的长时间系列数据时,速度提升更加明显)

前言 可根据研究区直接裁剪数据以及进行一些计算处理后再下载&#xff0c;GEE成为了大家下载数据的重要途径&#xff0c;然而直接通过官网网页将数据先保存到网盘再下载的下载方法速度太慢&#xff0c;新号速度还好&#xff0c;越用速度越来越慢&#xff0c;本文提供了一种直接…

DolphinScheduler 机器学习工作流预测今年 FIFA 世界杯冠军大概率是荷兰!

点击蓝字&#xff0c;关注我们作者 | DolphinScheduler Committer 周捷光2022 FIFA 世界杯火热进行中&#xff01;这段时间&#xff0c;这场盛宴吸引了全球球迷的目光。除了让人心跳加快的赛况和被大家调侃像馄饨皮的吉祥物之外&#xff0c;预测和投注哪支队伍将会夺冠绝对是球…

Antd中Table列表行默认包含修改及删除功能的封装

一、前言 ant-design是非常不错、方便的一款前端组件库&#xff0c;而这次用到的ProComponents则是在 Ant Design 上进行了自己的封装&#xff0c;更加易用&#xff0c;与 Ant Design 设计体系一脉相承&#xff0c;无缝对接 antd 项目&#xff0c;样式风格与 antd 一脉相承&am…

Java:多线程基础(二)-线程生命周期

目录 线程生命周期 Thread类的常用方法 构造方法 静态方法 常用实例方法 线程生命周期 线程有其创建、就绪、运行、阻塞、死亡的过程&#xff0c;将其称之为“线程的生命周期”。如下图所示&#xff0c; 对应以上5个状态&#xff0c;jdk-Thread类的源码中定义了枚举类Stat…

计算机网络第五章知识点回顾(自顶向下)

1. 网络层控制面 1.1 网络层功能 1.2选路问题 选路问题的描述&#xff1a; 给定一组路由器和连接路由器的链路&#xff0c;寻找一条从源路由器到目的路由器的最佳路径。 1.3 什么是最佳路径&#xff1f; 1.4 图抽象 1.5 选路算法分类 1.6 链路状态&#xff08;LS&#xff0…

“生成式技术”正在颠覆人类创作!

整理 | 王启隆在过去的半年里&#xff0c;AI 写小说、绘画和剪视频等热点新闻火爆全球&#xff0c;现在只需要在键盘上敲几个关键词&#xff0c;AI 就能在烧着我们显卡的同时画出一幅幅优美的图画&#xff0c;一个全新的应用世界向未来的初创公司敞开了大门。人类现在拥有着一大…

碳中和专利创新专题:各省市县专利面板(原始文件)、低碳专利授权数等多指标数据

一、各省市县专利面板含原始文件 1、数据来源&#xff1a;国家知识产权局 2、时间跨度&#xff1a;1985-2019年 3、区域范围&#xff1a;全国 4、指标说明&#xff1a; 来源 序号 标题 合享价值度 链接到incoPat 公开&#xff08;公告&#xff09;号 公开&#xff08…

数据结构之线性表中的栈和队列【详解】

文章目录引言&#xff1a;栈和队列的讲解&#xff08;一、&#xff09;什么是栈1.栈的概念、结构和图解&#xff1a;&#xff08;1.&#xff09;顺序表和链表的对比&#xff08;严格来说这两个结构是相辅相成的&#xff09;&#xff08;2.&#xff09;栈的概念和结构&#xff0…