目前本人写的ggrcs包新的2.6版本已经在CRAN上线,增加了多元线性回归分析模块,现在支持逻辑回归(logistic回归)、cox回归和多元线性回归

需要的可以使用代码安装

install.packages("ggrcs")

如果原来安装了旧版本,可以通过Rstudio进行升级



这样就可以升级到最新版本了,新版本修正了一些前期的错误和增加了绘制多元线性回归分析RCS模块,以我的臭氧数据(公众号回复:臭氧数据,可以获得数据)为例子,演示一下,先导入R包和数据

library(rms)
library(ggplot2)
library(scales)
library(ggrcs)
library(foreign)
be <- read.spss("E:/r/test/ozone.sav",use.value.labels=F, to.data.frame=T)  #臭氧数据


数据中有七个变量,ozon每日臭氧水平为结局变量,Inversion base height(ibh)反转基准高度,Pressure gradient (mm Hg) 压力梯度(mm Hg),Visibility (miles) 能见度(英里),Temperature (degrees F) 温度(华氏度),Day of the year日期,vh我也不知道是什么,反正就是一参数,这里所有的变量都是连续的。
假设我们想了解温度和臭氧浓度的关系
整理数据

dd <- datadist(be)
options(datadist='dd')

建立模型

fit<-ols(ozon ~rcs(temp, 4)+dpg,data=be)

绘图

ggrcs(data=be,fit=fit,x="temp",histbinwidth=1)


假设我们想了解不同分组中,温度与臭氧浓度的关系,我们数据中没有分组变量,我自己生成一个

set.seed(1234)
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(temp, 4)+dpg+variables2,data=be)

绘图

ggrcs(data=be,fit=fit,x="temp",group="variables2",histbinwidth=1)


OK,演示完毕,其他操作同前两个模型一样,我就不演示了,有兴趣的可以看我的既往文章。值得一提的是,线性回归一定要注意数据的分布,最好是正态分布,不是正态分布的最好取个对数弄成近似正态分布,不然绘制出来的图形会很怪。新版本如有BUG欢迎可以发私信给我。

ggrcs包2.6版本已经上线,增加了绘制多元线性回归RCS(限制立方样条图)模块相关推荐

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

    ggrcs 包2.4版本已经发布一段时间了,大概几个月了吧,收到不少好评, 没听说太大的问题,最主要的问题有两个: 1.是说变量不是数字变量. 2.是说数据超过10万,无法处理 第一个问题非常好处理, ...

  2. 首发: ggrcs包,一个用于绘制直方图+限制立方样条+双坐标轴图的R包

    ggrcs包是我写的第一个R包,所有权和解释权归我本人.用于绘制直方图+限制立方样条+双坐标轴图的图形. ggrcs包可以支持cox回归和logistic回归.为什么是叫ggrcs,因为依托了ggpl ...

  3. 【高能】IM3.9.1版本正式上线!

    本次版本更新以音频的优化为主,主要体现在音频音质和音频通话的弱网表现. 对于音质,我们综合各种用户使用场景,归纳出了三种模式:基础通话模式.高清语音模式和高清音乐模式,三种模式分别有其不同的适应场景( ...

  4. Authing 低代码登录组件 Guard Web 3.0 版本全新上线

    Authing 低代码登录组件 Guard Web 3.0 版本全新上线 据相关数据显示,30% 的用户在遇到登录问题时会选择离开,超出 10% 的交易被放弃是由于缺乏安全信任. Authing 是以 ...

  5. 2299元 OPPO K3 8GB+256GB版本线上线下同步开售

    近日,OPPO K3的最高配8GB+256GB版本在线上线下正式开售,售价2299元. OPPO K3自开售以来,凭借硬核配置备受消费者喜爱.6.1-6.6日,K3在京东1500-1999价位段的累计 ...

  6. Pip/python-如何查看已安装的包有哪些版本?如何查看某个包存在哪些版本?pip list-pip freeze-pip show package

    目录 一.如何查看已安装的包有哪些版本? 1.在terminal终端模式或者cmd命令行模式下 1.1显示所有包pip list和pip freeze 1.2显示指定包findstr.show..__ ...

  7. 到maven仓库加载jar包(AS版本)

    一般开发app程序所需要的依赖,网上都会直接给了现成的代码,类似于: implementation 'org.apache.poi:poi:5.2.2' implementation 'org.apa ...

  8. 融云 CDN 播放器 2.0 版本正式上线

    近期,融云自研 CDN 播放插件 2.0 版本正式上线.在原有版本支持融云内置 CDN 直播流播放的基础上,开放了相关 API,满足开发者的多种业务需求.关注[融云全球互联网通信云]了解更多 融云 C ...

  9. android 调用搜狗地图api,搜狗地图API2.0版本正式上线:推一键搬家功能

    搜狗地图API2.0版本正式上线 7月20日,搜狗地图API2.0版本正式上线( http://map.sogou.com/api /),在此前API1.0的基础上,在易用性和灵活性方面全面优化,并将 ...

最新文章

  1. Handlebars模板引擎
  2. FlipView和自定义值转换器
  3. 【LeetCode从零单排】No.160 Intersection of Two Linked Lists
  4. Web Components 入门实例教程
  5. linux 条件变量
  6. 巴特沃斯滤波器python_如何用Scipy.signal.bu实现带通巴特沃斯滤波器
  7. 第十一届 蓝桥杯 单片机设计与开发项目 省赛 程序设计试题及源码
  8. java jqgrid treegrid_JqGrid中文文档之TreeGrid
  9. php订单金额分摊算法,PHP付钱找钱算法
  10. linux 11Gasm启动crs,11GR2重建CRS
  11. redis sentinel 哨兵
  12. 标准输入输出详解(C++)
  13. VC6.0(VC++6.0)下载地址和安装教程(图解)
  14. PCB叠层排布原则以及AD中如何设计叠层
  15. 深度解析网易严选和京东的会员体系,建议收藏
  16. 《说服力-让你的PPT会说话》9月上海公开课简章
  17. 北航研究生计算机网络实验报告,实验7_北航研究生计算机网络实验
  18. 第二周铁人战队学习总结
  19. PLC数据采集网关有哪些作用,PLC数据采集网关如何正确进行使用
  20. VL综述:视觉-语言智能:任务、表征学习、大模型

热门文章

  1. Cadence OrCAD Capture 属性显示的无用信息太长的解决方法
  2. 轨道运营管理专业自荐书_城市轨道运营管理专业的求职信(精选多篇)
  3. 定位相关属性(positon:relative absolution; z-index;right、left、bottom、top)
  4. 第二十四届CIO班开学典礼隆重举行
  5. Spring Boot 2.6.0 新特性默认禁止循环引用
  6. 匹配电阻帮助提高放大器性能-电子技术方案
  7. Apollo:决策模块
  8. LinkButton 报错
  9. 视觉/音质/佩戴体验大提升,Rokid Max产品解析
  10. 多商户java版小程序+公众号+h5+app秒杀拼团砍价分销商城源码