ggrcs包2.6版本已经上线,增加了绘制多元线性回归RCS(限制立方样条图)模块
目前本人写的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(限制立方样条图)模块相关推荐
- ggrcs 包2.4绘制RCS(限制立方样条图)实际操作演示(1)
ggrcs 包2.4版本已经发布一段时间了,大概几个月了吧,收到不少好评, 没听说太大的问题,最主要的问题有两个: 1.是说变量不是数字变量. 2.是说数据超过10万,无法处理 第一个问题非常好处理, ...
- 首发: ggrcs包,一个用于绘制直方图+限制立方样条+双坐标轴图的R包
ggrcs包是我写的第一个R包,所有权和解释权归我本人.用于绘制直方图+限制立方样条+双坐标轴图的图形. ggrcs包可以支持cox回归和logistic回归.为什么是叫ggrcs,因为依托了ggpl ...
- 【高能】IM3.9.1版本正式上线!
本次版本更新以音频的优化为主,主要体现在音频音质和音频通话的弱网表现. 对于音质,我们综合各种用户使用场景,归纳出了三种模式:基础通话模式.高清语音模式和高清音乐模式,三种模式分别有其不同的适应场景( ...
- Authing 低代码登录组件 Guard Web 3.0 版本全新上线
Authing 低代码登录组件 Guard Web 3.0 版本全新上线 据相关数据显示,30% 的用户在遇到登录问题时会选择离开,超出 10% 的交易被放弃是由于缺乏安全信任. Authing 是以 ...
- 2299元 OPPO K3 8GB+256GB版本线上线下同步开售
近日,OPPO K3的最高配8GB+256GB版本在线上线下正式开售,售价2299元. OPPO K3自开售以来,凭借硬核配置备受消费者喜爱.6.1-6.6日,K3在京东1500-1999价位段的累计 ...
- Pip/python-如何查看已安装的包有哪些版本?如何查看某个包存在哪些版本?pip list-pip freeze-pip show package
目录 一.如何查看已安装的包有哪些版本? 1.在terminal终端模式或者cmd命令行模式下 1.1显示所有包pip list和pip freeze 1.2显示指定包findstr.show..__ ...
- 到maven仓库加载jar包(AS版本)
一般开发app程序所需要的依赖,网上都会直接给了现成的代码,类似于: implementation 'org.apache.poi:poi:5.2.2' implementation 'org.apa ...
- 融云 CDN 播放器 2.0 版本正式上线
近期,融云自研 CDN 播放插件 2.0 版本正式上线.在原有版本支持融云内置 CDN 直播流播放的基础上,开放了相关 API,满足开发者的多种业务需求.关注[融云全球互联网通信云]了解更多 融云 C ...
- android 调用搜狗地图api,搜狗地图API2.0版本正式上线:推一键搬家功能
搜狗地图API2.0版本正式上线 7月20日,搜狗地图API2.0版本正式上线( http://map.sogou.com/api /),在此前API1.0的基础上,在易用性和灵活性方面全面优化,并将 ...
最新文章
- Handlebars模板引擎
- FlipView和自定义值转换器
- 【LeetCode从零单排】No.160 Intersection of Two Linked Lists
- Web Components 入门实例教程
- linux 条件变量
- 巴特沃斯滤波器python_如何用Scipy.signal.bu实现带通巴特沃斯滤波器
- 第十一届 蓝桥杯 单片机设计与开发项目 省赛 程序设计试题及源码
- java jqgrid treegrid_JqGrid中文文档之TreeGrid
- php订单金额分摊算法,PHP付钱找钱算法
- linux 11Gasm启动crs,11GR2重建CRS
- redis sentinel 哨兵
- 标准输入输出详解(C++)
- VC6.0(VC++6.0)下载地址和安装教程(图解)
- PCB叠层排布原则以及AD中如何设计叠层
- 深度解析网易严选和京东的会员体系,建议收藏
- 《说服力-让你的PPT会说话》9月上海公开课简章
- 北航研究生计算机网络实验报告,实验7_北航研究生计算机网络实验
- 第二周铁人战队学习总结
- PLC数据采集网关有哪些作用,PLC数据采集网关如何正确进行使用
- VL综述:视觉-语言智能:任务、表征学习、大模型
热门文章
- Cadence OrCAD Capture 属性显示的无用信息太长的解决方法
- 轨道运营管理专业自荐书_城市轨道运营管理专业的求职信(精选多篇)
- 定位相关属性(positon:relative absolution; z-index;right、left、bottom、top)
- 第二十四届CIO班开学典礼隆重举行
- Spring Boot 2.6.0 新特性默认禁止循环引用
- 匹配电阻帮助提高放大器性能-电子技术方案
- Apollo:决策模块
- LinkButton 报错
- 视觉/音质/佩戴体验大提升,Rokid Max产品解析
- 多商户java版小程序+公众号+h5+app秒杀拼团砍价分销商城源码