ggrcs包4.0版本发布—重新对密度图宽度进行了设计

目前本人写的ggrcs包新的4.0版本已经在CRAN上线,目前支持逻辑回归(logistic回归)、cox回归和多元线性回归。
在这里插入图片描述
需要的可以使用代码安装

install.packages("ggrcs")

如果原来安装了旧版本,重新在安装一次就可以升级到新版本。本次主要是修复了之前的两个问题,下面我来演示一下。
第一个是有柱子在负轴的问题,既往偶有柱子在负轴的问题,我用一个粉丝的数据来演示一下,代码不解释了,可以看我既往的文章

library(rms)
library(ggplot2)
library(scales)
library(cowplot)

setwd("E:/公众号文章2024年/ggrcs报错")
df <- read.csv('rscdata.csv',sep = ',',header = T)
str(df)
covars = c('age','Gender','Race','edu','marital',
           'PIR','alcohol.user','smoke')

df[,covars] <- lapply(df[,covars],as.factor)

dd <- datadist(df) 
options(datadist='dd') 
fit<-ols(PAA ~rcs(X2.OHNAP,4)+age+Gender+Race+marital+edu+PIR+bmi+
           smoke+alcohol.user,data=df)
summary(fit)
ggrcs(data=df,fit=fit,x="X2.OHNAP",histbinwidth=0.01,px=1.1)

在这里插入图片描述
这是因为线性回归种可能会产生负值,可能会导致这种情况,新版本解决了这个问题

library(ggrcs)
ggrcs(data=df,fit=fit,x="X2.OHNAP",px=1.1)

在这里插入图片描述
第二个就是线性回归中因为宽度问题导致绘图不美观的问题,使用我的臭氧数据来演示这个情况(公众号回复:臭氧数据,可以获得这个数据)

library(ggrcs)
library(rms)
library(ggplot2)
library(scales)
library(cowplot)
library(foreign)
be <- read.spss("E:/r/test/ozone.sav",
                use.value.labels=F, to.data.frame=T)  #????
names(be)
be$variables2<-sample(0:1,size=330,replace=TRUE)
be$variables2<-as.factor(be$variables2)
dd <- datadist(be) 
options(datadist='dd') 
fit<-ols(ozon ~rcs(vh, 4)+temp+variables2,data=be)
ggrcs(data=be,fit=fit,x="vh")

在这里插入图片描述
绘制出的这个图不怎么美观是因为宽度没有控制好,新版本对宽度重新进行自动计算后,还是原来的代码

ggrcs(data=be,fit=fit,x="vh")

在这里插入图片描述
这样图片就比原来美观多啦,多谢各位粉丝提出的意见。

  • 1
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

天桥下的卖艺者

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值