如何在r中的时间序列数据中执行mann kendall和sen的斜率

时间:2017-01-01 05:34:41

标签: r plot time-series

https://drive.google.com/file/d/0B3fQ9_46L-O0ZXVYb3lzZDBZaHM/view?usp=sharing https://drive.google.com/file/d/0B3fQ9_46L-O0TjJwYmF6UThNSGs/view?usp=sharing

这里,我有一个CSV格式的同一个地方的两个站的降水的时间序列数据。我试图找到这两个站的降水趋势以及趋势的大小。为此,我使用了Mann-Kendall和Sen的斜坡。我试着这种特殊的方式,我想知道我做得对不对?在下面的代码中,我单独计算趋势和幅度,但我的目标是一起计算两个站的趋势和幅度。我如何能够同时对两个站点进行预白化,Mann-Kendall和Sen的斜率测试以及它们的情节?

library(xts)
library(Kendall)
library(zyp)
library(trend)
setwd("~/Desktop/pawan thapa 2/pawan thapa/rautahat")
ram = read.csv('ramoli_preci.csv',header=TRUE, stringsAsFactors=FALSE)
gaur2= read.csv('preci_Gaur.csv',header=TRUE, stringsAsFactors=F)
# Get data properties
dim(gaur2)
dim(ram)
head(ram)
date <- as.Date(as.character(gaur2$Year), "%Y")
date_r <- as.Date(as.character(ram$Year), "%Y")
gaur.ts<-as.ts(gaur2)
ram.ts<-as.ts(ram)
head(ram.ts)
class(ram.ts)
head(ram.ts)
hist(ram.ts)
acf(ram.ts)

mk.test(ram.ts)
sens.slope(ram.ts)

head(date)
tail(date)
class(date)
class(gaur2)
str(gaur2)

class(gaur.ts)
head(gaur.ts)
mk.test(gaur.ts)
sens.slope(gaur.ts)
hist(gaur.ts)

0 个答案:

没有答案