月下醉光遛鸟

R语言 数据列出现X. 检查列名有效性 check.names参数问题 改为FRUE去除检查

R语言reshape2包中的melt和dcast函数进行长宽转换，dcast长数据转换为宽数据中如何选择多列

R语言关于散点图的边框和颜色填充问题

ggplot 图例如何进行分段

Error in `vec_as_location()`:! `...` must be empty.x Problematic argument:* call = call重新安装dplyr包

ggplot2()绘制geom_ribbon()和geom_line()的画法，以及注意点

ggplot2绘制折线图：factor和group

ggplot2如何在填充颜色，主要使用的是geom_ribbon()函数；如何给折线添加置信区间和拟合线条。

ggplot2 fact_grid()按行和列分泌按后如何进行背景、字体、间距的设置

R语言绘制次级坐标轴、双坐标轴

R Error in match.names(clabs, names(xi)) : 名字同原来已有的名字不相对

par、ts.plot、R、

par、画图、Rstudio、R

git remote出现错误

博主在使用SPSS统计软件做pearson相关性分析的时候，发现用热图可视化可能效果会更好，于是，博主就去度娘了，然后用R代码绘制，但是发现做显著性的时候（P值检验），p值计算出来竟然不一样，惊呆了！! >_< ...

person相关系数、相关性分析、显著性检验

相关系数如何绘制呐？

相关系数如何绘制呐？

相关系数如何绘制呐？

查询死因数据的编码方法：使用网址查询：>>点击这儿<<疾病名称+ICD-10/11查询,如下所示：

数据篇corrplot绘制相关热图，csv数据里面含有NA值 or 数值源不是数值类型的，使用 read.csv() 读取的时候可以读取出来，但是运行完毕的时候，会出现?，如下图：解决办法就是将NA值改为0，对应数据修改成数值类型...

点击source，不出现图，解决办法：在画图代码前面添加– dev.new()，使用此函数，创建新的图形窗。销毁叼用下面函数或者点击叉叉– dev.off()，使用此函数，销毁创建的图形窗。补充：– dev.list() 计算有多少图形窗– dev.curt() 当前是哪个图形窗– dev.next() 下一个窗口...

RR估计值得时候会看显著性的问题，把今天学到的知识点记录一下：p值小于0.05,还有RR所有的值大于1就是正显著，负显著就是p值小于0.05,RR所有的结果小于1，其他情况就是无显著解释： 这个是跑出来的结果图，其中黑色的就是不显著，记为0，即所有结果（效应值和置信区间）骑在1.00这条横线上。红色的是正显著，记为1，即所有结果在1.00线以上；蓝色的是负显著，记为-1，即所有结果在1.00以下。...

关于cox比例风险模型的一些讲解比较清楚的网站的备注，方便自己查找

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)...

1. summary()函数可以获取描述性统计量可以提供最小值、最大值、四分位数和数值型变量的均值，以及因子向量和逻辑型向量的频数统计2. misc包中的describe()函数可返回变量和观测的数量、缺失值和唯一值的数目、平均值、分位数，以及五个最大的值和五个最小的值3.psych包中的describe()函数psych包也拥有一个名为describe()的函数，它可以计算非缺失值的数量、平均数、标准差、中位数、截尾均值、绝对中位差、最小值、最大值、值域、偏度、峰度和平均值的标准误4.past

分面写的很详细的文章：关于网格分面，封装分面，标度控制，分组与分面，并列与分面，连续型变量的分面。请移步

R语言中关于as.Date,时间函数POSIXct与POSIXlt,以及日期相关的提取，计算，请看这篇这个博客文章内容，写的挺好的。R相关时间处理，请移步！

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==

R语言针对调整列的位置方式：

具体的信息请移步：AIC和BIC讲解

经常用数据分析，有时不同的文件的小数位数不一样，但是我们可以让它们的位数保持一致的，下面的介绍就是设置小数位数。#使用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

水一下：提取nc内容最重要的是要确认好变量，时间变量这些格式有几个，有几层。然后经纬度循环，时间循环一次读写出来。注意：conc格式、xx.lat_lon格式都是nc格式

有一个数据框名叫df，R语言获取数据框行数的方法dim(df)[1]nrow(df)length(df[,1])

提取特定列或者行，使用c()时候提示错误：longer object length is not a multiple of shorter object length改成使用 %in%%in%的使用：

主要是站点监测值如何分组提取单一的站点和对应的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 =

主要是通过日期将站点监测值和PM2.5模拟值对应的日期匹配连接起来：主要是通过merge()函数来实现具体的实现步骤：（1）首先将日期在excel里面修改一下格式，修改成20130101这样的，不然的画2013/01/01这样直接匹配的话，匹配的模拟值都是NA注： 具体如何实现的可以百度。在这里简单说一下，就是选中日期列，右击选择设置单元格格式，选择自定义，手动修改成yyyymmdd,确定就好了。日期就是这样修改的，当然也可以把20130101修改成2013/01/01，一样的操作，还有一种方法就

重新跑了一遍，但是貌似还是缺了一个区县，gan。步骤：一、先把给的数据展示到gis中，然后把这个数据和不缺的标准区县的地图进行空间连接（这次没有做FID_county，主要是这个都是唯一的，没有重复的，其实也不需要用这个来，用这个只是简单一点）二、把数据导出来删除不要的列，然后留做合并数据，通过代码跑出来不缺的区县的数据tq <- function(input,output){ c <- dir(input) cd <-

就是说站点的值要和对应的数据集的格网的值要对应起来。或者说从区县角度来看，就是所有的值都用区县对应的PAC来匹配站点监测值，空间上匹配上具体的步骤：第一步，首先要将站点展示到GIS中（13-17年），保存成shp文件，然后将shp文件和标准区县文件空间连接（站点做目标文件，区县做连接文件），保存成1.csv.第二步，使用代码将1.csv和数据集的对应年份的污染值使用merge函数提取连接合并。fun <- function(input,station_TAP_spation,output){

拿到手的站点监测数据是小时的文件，并且有很多不需要的信息，主要是清理成自己想要的数据格式，主要用到的是aggregate（）函数，具体的做法自己百度。清理前：代码：#主要是清理并合并tq <- function(input,output){ bb <- dir(input) cd <- paste(input,bb,sep = "") m <- NULL for(i in 1:length(cd)){ aa <- read.c