1、第一步:安装读取excel文件包"openxlsx"install.packages("openxlsx")library(openxlsx)
2、第二步:载入自己的数据data <- read.xlsx("F:\\你自己数据的路径.xlsx")#查看数据前几行 head(data)
3、第三步:安装趋势分析包“trend”install.packages("trend")library(trend)
4、第四步:进行MK计算使用方法:mk.test(x, alternative = c("two.sided", "greater", "less"), continuity = TRUE)示例:mk.test(data$径流深, continuity = TRUE)
5、第五步:计算结果解读。若z值小于0,表明该时间序列呈下降趋势,若z值大于0,表明该时间序列拂痒缰萼呈上升趋势。若p>0.01,接受零假设,表明数据趋势变化不显著。p<0.01,拒绝零假设,表明数据趋势变化显著。