![](https://img-blog.csdnimg.cn/20201014180756916.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
R
文章平均质量分 56
schneesnow
这个作者很懒,什么都没留下…
展开
-
R语言rle()函数
计算向量x中等值游程的长度和值。rle显示一个向量x中每一个连续重复元素及该元素对应的值。x可以是不同的格式,如数值,字符串,逻辑变量等。适合用于游程理论的一些编程。例子1,x为数值向量x <- rev(rep(6:10, 1:5))## x<-c(10,10,10,10,10,9,9,9,9,8,8,8,7,7,6)rle(x)## lengths [1:5] 5 4 3 2 1## values [1:5] 10 9 8 7 6例子1,x为逻辑向量z <原创 2021-08-13 22:55:08 · 2122 阅读 · 0 评论 -
ggplot2实现色阶散点图+曲线拟合图
library(ggplot2)library(scales)# humidity以百分比的形式出现T=sort(runif(500,min=24,max=38))O3=8*10^(-11)*T^7.99O3=jitter(O3,amount=30)O3=O3-min(O3)+10# 保证O3没有负值Humidity=sort(runif(500,0,1),decreasing = TRUE)df=data.frame(T=T,O3=O3,Humidity=Humidity)#你有自己的数据原创 2021-04-06 08:52:33 · 1136 阅读 · 0 评论 -
Diebold-Mariano检验
给定两个预测的预测结果,我们希望比较他们的预测结果,以用于预测模型精度的比较。Diebold-Mariano检验本质是一个t检验,用于检验替代预测的两个损失序列的平均值是否相等。即,它是一系列损失差的零均值的t检验。在存在自相关的情况下,它使用损失差分时间序列标准差的自相关一致性估计。注意,在使用DM检验式时,假设损失序列是平稳的。另外,DM检验在小样本数据时往往会拒绝零假设。对于小样本数据,推荐Harvey, Leybourne and Newbold (HLN)检验【1】;参考文献【1】htt原创 2021-01-20 16:07:39 · 5240 阅读 · 4 评论 -
R语言无法从GitLAB/GITHUB远程安装
问题描述我以ggeconodist包安装为例,一般地,安装的命令可如下形式。install.packages("ggeconodist", repos = "https://cinc.rud.is")# orremotes::install_git("https://git.rud.is/hrbrmstr/ggeconodist.git")# orremotes::install_git("https://git.sr.ht/~hrbrmstr/ggeconodist")# orremote原创 2021-01-08 16:40:07 · 1175 阅读 · 1 评论 -
Circular-Circular 和 Circular-Linear Regression
Circular-Circular的例子在一个海边的小岛上有着某种稀缺鱼类,科学家关注这种鱼的产卵时间(TS)是否取决于低潮时间(TLT)。在产卵期,科学家甩了几个小硕士在岛上记录了一波实验数据[TSi,TLTi]。为了研究TS对TLT的依赖关系,科学家将这波数据甩给了某个秃头的延期博士。该博士将TS的0 ~ 24小时的周期转换为[0,24];将TLT转换为相位区间[0,2π]。因此配对的TS和TLT是二元循环数据。如下图,它们被绘制成圆圈。根据文献【1】,可分析TS对TLT的统计相干性,若存在显原创 2020-08-20 15:49:22 · 548 阅读 · 1 评论 -
利用R.matlab包将MATLAB数据导入到R语言
目录1 安装R.matlab包2.在MATLAB中保存.mat文件3. 在R环境中读取.mat文件1 安装R.matlab包安装完R.matlab后,加载该包的过程中可能会提示加载不成功。按照提示继续装提示缺乏的包环境,加载R.methodsS3、R.oo、R.utils。2.在MATLAB中保存.mat文件下列示例代码将ER矩阵保存为ER.mat。注意‘’-v7’’为可用的MAT文件版本,该版本与R.matlab包相匹配。save('ER.mat','ER','-v7');3. 在R环境中读原创 2020-08-17 15:20:13 · 2718 阅读 · 0 评论 -
ggplot2实现多组箱线图上添加均值点(线)
正文library(ggplot2)n <- 1000df <- data.frame(var1=sample(LETTERS[1:3],n, T), var2=sample(LETTERS[1:2],n,T), var3=sample(LETTERS[1:3],n,T), y=rnorm(n))p1 <- ggplot(df, aes(x=var1, y=y)) +原创 2020-06-27 21:04:29 · 10434 阅读 · 2 评论 -
circular包绘图笔记
官方画图流程:1.initialize layout -> 2.create track -> 3.add graphics -> 4.create track -> 5.add graphics - … ->6. clear流程的土味讲解,1.初始化布局–定义factor,也就是定义我一共有几个扇区,扇区的名字和属性 。2.创建轨道–定义circular一共有几圈。无需太在意每个扇区图的完整程度,因为每个圈的每个扇区(即cell)后面都可以单独修改。3.添加图的信原创 2020-05-26 22:31:02 · 974 阅读 · 0 评论 -
ggplot2点图+线性趋势+公式+$R^2$+p值
先看效果。R语言代码如下,data("faithful")library(ggplot2)p <- ggplot(faithful,aes(x=eruptions,y=waiting)) + geom_point() + stat_smooth(method='lm',formula = y~x,colour='red')lm(formula = waiting ~ erupt...原创 2020-04-06 00:28:58 · 10173 阅读 · 6 评论 -
R语言下载气象数据
近来气象数据共享网的账户到期,失去日序列数据的下载权限。总结了以下R语言气象数据的下载包。package一览下列所有包的来源均属于 rOpenSci 项目.1. GSODTools2. GSOD3. rnoaa4. bomrang5. riem6. weathercan7. CliFlo8. weatherData上述包,简介如下:The GSODTools by Flo...原创 2020-03-28 17:02:37 · 2354 阅读 · 0 评论 -
Cramér‐von Mises检验
文章目录0写在前面1.单样本CVM检验2.双样本CVM检验参考文献0写在前面Cramér‐von Mises检验用于检验一列数据是否与另一组数据有着同样的分布一种非参数检验方法,是一种比KS检验更加严谨的检验方法。单样本CVM检验检验未知分布是否来自某参数分布,如正态分布等。若要比较两组未知分布的数据则需要双样本CVM检验。1.单样本CVM检验2.双样本CVM检验公式以及详细介绍...原创 2019-11-15 22:51:50 · 5426 阅读 · 0 评论 -
ggplot2默认字体修改
文章目录1.查看当前的ggplot2默认字体2.查看Windows系统下的字体3.修改成Time New Roman4.参考文献1.查看当前的ggplot2默认字体theme_get()$text显示一个list,family为默认空白,默认是sans(Arial字体)。List of 11$ family : chr “”$ face : chr “pl...原创 2019-11-13 10:16:59 · 34942 阅读 · 2 评论 -
pheatmap绘制热力图
pheatmap绘制热力图(1)坐标轴顺序pheatmap包在绘制热力图时,行名不与data的行名一致,如何保留原始data的行名,控制参数,cluster_rows=F, cluster_cols=F.(2)调色板函数①colorRamp {grDevices}函数创建颜色梯度,1000表示离散成一千格。要显示的数据就可以插值对应1000种颜色中的某一种。colorRampPal...原创 2019-11-08 21:36:59 · 1507 阅读 · 0 评论 -
ggplot坐标轴形式转换
**0.写在前边**在画图的时候,数据可能集中在图的某一局部范围内或者我们论文描述的重点是图的某一特定范围。因此在图中清晰的表达出重点关注的范围需要进行坐标轴转换,将图片在坐标轴上适当的压缩拉伸。下面简单的介绍了ggplot2两类坐标轴转换形式。1.将x轴设置成平方根转换,coord_trans(x = “sqrt”)*任意选用一数据包,观察数据在[0,1000]范围内集中,可通过...原创 2019-04-19 16:42:04 · 2646 阅读 · 0 评论 -
R语言零碎整理
R语言零碎整理1.ggplot2画柱形图时,柱子顺序的调整默认情况下是按照英文字母递增顺序排序,即b,d,g,k,r,w。需要对date.frame调整,修改的命令参考如下。ma$classify <- factor(ma$classify, levels=c( "b", "r","g","w","d","k"), ordered=TRUE)2.ggplot2调整调色板颜色(1).........原创 2019-10-30 16:54:04 · 5838 阅读 · 0 评论