![](https://img-blog.csdnimg.cn/20201014180756738.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
R语言相关
存放和R语言相关的学习经验、技巧
初于青丝mc终于白发
学习如逆水行舟,不进则退
展开
-
R语言在读取数据的时候自动加上X. 检查列名的有效性
R语言 数据列出现X. 检查列名有效性 check.names参数问题 改为FRUE去除检查原创 2022-11-15 23:44:53 · 1552 阅读 · 0 评论 -
【R语言】实现长宽数据转换——dcast()长转宽
R语言reshape2包中的melt和dcast函数进行长宽转换,dcast长数据转换为宽数据中如何选择多列原创 2022-10-15 12:30:09 · 2085 阅读 · 0 评论 -
【R语言】ggplot2---散点图的边框和颜色填充问题
R语言关于散点图的边框和颜色填充问题原创 2022-09-18 15:22:58 · 8398 阅读 · 0 评论 -
【R语言】ggplot绘制图形的图例如何进行分段标记
ggplot 图例如何进行分段原创 2022-09-14 15:52:38 · 551 阅读 · 0 评论 -
【Rstudio】出现的错误:Error in `vec_as_location()`: `...` must be empty.
Error in `vec_as_location()`:! `...` must be empty.x Problematic argument:* call = call重新安装dplyr包原创 2022-09-11 22:10:47 · 2364 阅读 · 0 评论 -
ggplot2:geom_ribbon()去除置信区间的框线
ggplot2()绘制geom_ribbon()和geom_line()的画法,以及注意点原创 2022-07-29 20:20:44 · 3762 阅读 · 0 评论 -
x轴年份只有完整的时候,使用factor的问题
ggplot2绘制折线图:factor和group原创 2022-07-28 18:10:24 · 156 阅读 · 0 评论 -
R语言:ggplot2如何在两条折线之间添加颜色
ggplot2如何在填充颜色,主要使用的是geom_ribbon()函数;如何给折线添加置信区间和拟合线条。原创 2022-07-28 17:34:38 · 1807 阅读 · 0 评论 -
R语言:ggplot2的facet_grid按行和列分面后的设置
ggplot2 fact_grid()按行和列分泌按后如何进行背景、字体、间距的设置原创 2022-07-24 17:40:18 · 1775 阅读 · 0 评论 -
R语言 ggplot2绘制双坐标轴(次坐标轴)
R语言绘制次级坐标轴、双坐标轴原创 2022-07-15 18:35:13 · 2313 阅读 · 0 评论 -
Error in match.names(clabs, names(xi)) : 名字同原来已有的名字不相对
R Error in match.names(clabs, names(xi)) : 名字同原来已有的名字不相对原创 2022-06-19 20:32:08 · 513 阅读 · 0 评论 -
ts.plot()绘图后,坐标轴修改失败,使用par函数进行修改
par、ts.plot、R、原创 2022-06-17 15:08:26 · 956 阅读 · 0 评论 -
Error in plot.new() : figure margins too large
par、画图、Rstudio、R原创 2022-06-17 14:30:02 · 564 阅读 · 0 评论 -
Git推送remote的时候报错原因
git remote出现错误原创 2022-06-12 10:19:27 · 5862 阅读 · 0 评论 -
R代码和SPSS绘制相关性热图/矩阵的差别说明
博主在使用SPSS统计软件做pearson相关性分析的时候,发现用热图可视化可能效果会更好,于是,博主就去度娘了,然后用R代码绘制,但是发现做显著性的时候(P值检验),p值计算出来竟然不一样,惊呆了!! >_< ...原创 2022-06-06 14:37:42 · 4320 阅读 · 0 评论 -
person相关性分析-显著性检验
person相关系数、相关性分析、显著性检验原创 2022-06-05 16:27:05 · 1954 阅读 · 0 评论 -
绘制相关性热图/矩阵
相关系数如何绘制呐?原创 2022-06-02 11:56:34 · 217 阅读 · 0 评论 -
绘制相关性热图/矩阵
相关系数如何绘制呐?原创 2022-06-02 11:56:04 · 397 阅读 · 0 评论 -
绘制相关性热图/矩阵
相关系数如何绘制呐?原创 2022-06-01 23:28:48 · 262 阅读 · 0 评论 -
ICD-10/11 查询疾病编码方法
查询死因数据的编码方法:使用网址查询:>>点击这儿<<疾病名称+ICD-10/11查询,如下所示:原创 2022-05-30 15:31:40 · 8933 阅读 · 0 评论 -
R语言运行出现的小问题
数据篇corrplot绘制相关热图,csv数据里面含有NA值 or 数值源不是数值类型的,使用 read.csv() 读取的时候可以读取出来,但是运行完毕的时候,会出现?,如下图:解决办法就是将NA值改为0,对应数据修改成数值类型...原创 2022-05-29 10:54:43 · 477 阅读 · 0 评论 -
关于使用Rstudio不能画图问题(source点击不出现图)ggplot2
点击source,不出现图,解决办法:在画图代码前面添加– dev.new(),使用此函数,创建新的图形窗。销毁叼用下面函数或者点击叉叉– dev.off(),使用此函数,销毁创建的图形窗。补充:– dev.list() 计算有多少图形窗– dev.curt() 当前是哪个图形窗– dev.next() 下一个窗口...原创 2022-05-22 22:19:47 · 2242 阅读 · 0 评论 -
RR正显著-不显著 -负显著
RR估计值得时候会看显著性的问题,把今天学到的知识点记录一下:p值小于0.05,还有RR所有的值大于1就是正显著,负显著就是p值小于0.05,RR所有的结果小于1,其他情况就是无显著解释: 这个是跑出来的结果图,其中黑色的就是不显著,记为0,即所有结果(效应值和置信区间)骑在1.00这条横线上。红色的是正显著,记为1,即所有结果在1.00线以上;蓝色的是负显著,记为-1,即所有结果在1.00以下。...原创 2022-05-13 17:20:42 · 2655 阅读 · 0 评论 -
R语言—cox比例风险模型
关于cox比例风险模型的一些讲解比较清楚的网站的备注,方便自己查找原创 2022-04-14 16:13:47 · 570 阅读 · 0 评论 -
R语言 formula()
R 语言回归函数里面的公式函数, as.formula(). 其作用就是将字符串转换成公式。> aa = "ReadCount~Age+BMI+Sex+HAMD+1+(1|Sex)"> aa[1] "ReadCount~Age+BMI+Sex+HAMD+1+(1|Sex)"> as.formula(aa)ReadCount ~ Age + BMI + Sex + HAMD + 1 + (1 | Sex)...原创 2022-04-08 10:00:40 · 4046 阅读 · 0 评论 -
R语言中描述统计量的多种方法
1. summary()函数可以获取描述性统计量可以提供最小值、最大值、四分位数和数值型变量的均值,以及因子向量和逻辑型向量的频数统计2. misc包中的describe()函数可返回变量和观测的数量、缺失值和唯一值的数目、平均值、分位数,以及五个最大的值和五个最小的值3.psych包中的describe()函数psych包也拥有一个名为describe()的函数,它可以计算非缺失值的数量、平均数、标准差、中位数、截尾均值、绝对中位差、最小值、最大值、值域、偏度、峰度和平均值的标准误4.past原创 2022-04-07 21:01:18 · 3088 阅读 · 0 评论 -
R 数据可视化 —— ggplot2 分面
分面写的很详细的文章:关于网格分面,封装分面,标度控制,分组与分面,并列与分面,连续型变量的分面。请移步原创 2022-03-29 17:39:42 · 427 阅读 · 0 评论 -
R语言中的时间处理
R语言中关于as.Date,时间函数POSIXct与POSIXlt,以及日期相关的提取,计算,请看这篇这个博客文章内容,写的挺好的。R相关时间处理,请移步!原创 2022-03-27 14:14:45 · 308 阅读 · 0 评论 -
R语言 subset()函数用法
subset() 函数:subset(dataset , subset , select )dataset 是要进行操作的数据集subset 是对数据的某些字段进行操作select 选取要显示的字段```{r spdep}hightechaa = subset(hightech,year==2001&province=="北京",select = c("year","province","PAT"))bb = subset(hightech,year==2001|province==原创 2022-03-07 10:11:07 · 6472 阅读 · 0 评论 -
R语言如何进行列的调整
R语言针对调整列的位置方式:原创 2022-03-02 14:00:39 · 3730 阅读 · 0 评论 -
AIC(赤池信息准则)和BIC(贝叶斯信息准则)
具体的信息请移步:AIC和BIC讲解原创 2022-02-08 16:34:46 · 3529 阅读 · 0 评论 -
R语言中对小数点的位数的设置
经常用数据分析,有时不同的文件的小数位数不一样,但是我们可以让它们的位数保持一致的,下面的介绍就是设置小数位数。#使用options函数options(digits) #默认为7位> a=0.234333323#9位> a[1] 0.2343333 #默认为7位#下面开始设置下> options(digits=3) #设置为3位> a=0.34434434#8位> a[1] 0.344 #3位#最大位数是22位> a = 0.11原创 2022-02-07 17:34:39 · 2848 阅读 · 0 评论 -
关于nc文档提取的那些事
水一下:提取nc内容最重要的是要确认好变量,时间变量这些格式有几个,有几层。然后经纬度循环,时间循环一次读写出来。注意:conc格式、xx.lat_lon格式都是nc格式原创 2022-01-22 17:06:20 · 993 阅读 · 0 评论 -
R语言获取数据框行数
有一个数据框名叫df,R语言获取数据框行数的方法dim(df)[1]nrow(df)length(df[,1])原创 2022-01-04 15:29:35 · 5524 阅读 · 0 评论 -
R语言 使用 c() 遇到的问题的解决---%in%
提取特定列或者行,使用c()时候提示错误:longer object length is not a multiple of shorter object length改成使用 %in%%in%的使用:原创 2021-11-21 23:39:00 · 1508 阅读 · 0 评论 -
站点监测值如何分组合并
主要是站点监测值如何分组提取单一的站点和对应的PM2.5值的问题:如下辣么多主要是用函数aggregate()来实现:tq <- function(input,output){ bb <- dir(input) cd <- paste(input,bb,sep = "") m <- NULL for(i in 1:length(cd)){ aa <- read.csv(cd[i],header = T,fileEncoding =原创 2021-12-13 18:29:56 · 286 阅读 · 0 评论 -
如何实现站点监测值和模拟值通过日期匹配
主要是通过日期将站点监测值和PM2.5模拟值对应的日期匹配连接起来:主要是通过merge()函数来实现具体的实现步骤:(1)首先将日期在excel里面修改一下格式,修改成20130101这样的,不然的画2013/01/01这样直接匹配的话,匹配的模拟值都是NA注: 具体如何实现的可以百度。在这里简单说一下,就是选中日期列,右击选择设置单元格格式,选择自定义,手动修改成yyyymmdd,确定就好了。日期就是这样修改的,当然也可以把20130101修改成2013/01/01,一样的操作,还有一种方法就原创 2021-12-13 18:40:23 · 284 阅读 · 0 评论 -
TAP_restart
重新跑了一遍,但是貌似还是缺了一个区县,gan。步骤:一、先把给的数据展示到gis中,然后把这个数据和不缺的标准区县的地图进行空间连接(这次没有做FID_county,主要是这个都是唯一的,没有重复的,其实也不需要用这个来,用这个只是简单一点)二、把数据导出来删除不要的列,然后留做合并数据,通过代码跑出来不缺的区县的数据tq <- function(input,output){ c <- dir(input) cd <-原创 2021-12-14 13:10:16 · 242 阅读 · 0 评论 -
如何画R2(拟合散点图)
就是说站点的值要和对应的数据集的格网的值要对应起来。或者说从区县角度来看,就是所有的值都用区县对应的PAC来匹配站点监测值,空间上匹配上具体的步骤:第一步,首先要将站点展示到GIS中(13-17年),保存成shp文件,然后将shp文件和标准区县文件空间连接(站点做目标文件,区县做连接文件),保存成1.csv.第二步,使用代码将1.csv和数据集的对应年份的污染值使用merge函数提取连接合并。fun <- function(input,station_TAP_spation,output){原创 2021-12-20 15:59:52 · 1774 阅读 · 0 评论 -
如何清理站点监测数据
拿到手的站点监测数据是小时的文件,并且有很多不需要的信息,主要是清理成自己想要的数据格式,主要用到的是aggregate()函数,具体的做法自己百度。清理前:代码:#主要是清理并合并tq <- function(input,output){ bb <- dir(input) cd <- paste(input,bb,sep = "") m <- NULL for(i in 1:length(cd)){ aa <- read.c原创 2021-12-20 16:06:44 · 264 阅读 · 0 评论