活鱼现杀-CSDN博客

原创 R语言计算时间序列（CST）交集、并集

针对union()、intersect()并集、交集等基础函数在处理时间序列（如CST北京时）时，时间序列发生变化情况的若干改进与处理####生成时间数据序列####series_time_a<-as.POSIXct("2022-05-02 00:00:00")+1:5series_time_aseries_time_b<-as.POSIXct("2022-05-02 00:00:03")+1:5series_time_b####计算a,b时间序列并集####series_t.

2022-05-02 11:09:36 597

原创 R语言判断是否是闰年-自定义函数

R语言判断是否是闰年-自定义函数定义含参函数判断年份是否为闰年，代码如下：leap<-function(x){ if(x%%4==0){ if(((x%%100==0)&(x%%400==0))|(x%%100!=0)) print(TRUE) else print(FALSE) } else

2020-11-02 09:52:06 2423 1

原创 R语言计算带有绝对值符号的函数偏导并输出偏导结果

R语言对带有绝对值符号的函数求偏导并输出导数结果求偏导后的导数表达式网上已有众多教程，这里主要对导数结果的输出及含有绝对值符号的的函数偏导运算进行方法阐述对函数f=x*|y|*(z-x)求对x、y、z的偏导，输出在点（1,2,3）和（1,-2,3）的导数值，代码如下：#当y>0时，函数f如下：fpo<-expression(x*y*(z-x))#当y<0时，函数f如下：fne<-expression(x*(-y)*(z-x))#构造求偏导函数fun = functio

2020-07-13 16:07:54 2808

原创 R语言求和上三角矩阵

R语言求和上三角矩阵举例如下w<-matrix(20:31,3,3)w_up<-w[which(upper.tri(w))]w_upsum<-sum(w_up)ww_upw_upsum结果如下w[,1] [,2] [,3][1,] 20 23 26[2,] 21 24 27[3,] 22 25 28w_up[1] 23 26 27w_upsum[1] 76...

2020-06-11 11:18:40 2681

原创利用R语言进行MK检验以检测序列变化趋势

关于MK检验，目前因为搜索不到相关的R语言程序包，所以本人按照魏凤英老师在其著作《现代气候统计诊断与预测技术》中的MK检验方法，利用R语言写下如下代码，希望能够帮到大家，相互切磋，相互交流。本文所用数据也为魏老师书中1900-1990年上海年平均气温序列注：本文只是为大家提供思路，文中所画图形与魏凤英老师书中对应图像有所差异，但我觉得的这个程序没有错，本人才疏学浅，如果看到这篇文章的您觉得我程...

2019-10-12 17:32:16 14564 53

原创 R语言二进制grd文件读取&合并&导出

数据格式：1* 31* 10，共17个文件，存在D:\cyh\zpmi\路径下setwd("D:/cyh/zpmi/")path <- "D:/cyh/zpmi"fileNames <- dir(path)filePath <- sapply(fileNames, function(x){ paste(path,x,sep='/')})wind<-array(...

2019-04-19 22:40:55 2194 3

原创 R语言气象相关分析方法

R语言气象相关分析方法（包括二进制文件R语言读取）######读取热带印度洋海温场资料#################read.filename <- file("D:\\99\\NCEP_IOSST_30y_Wt.dat", "rb")bindata <- readBin(read.filename, numeric(),16*7*30,size = 4)IOSST<...

2019-02-03 12:06:32 2923

weixin_39634034的博客