登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x

本帖最后由 sfhstcn2 于 2017-10-19 07:29 编辑

mutation_test这是一包代码,是基于《现代气候统计诊断与预测技术》(魏凤英)书中的气候突变检测一章的7种检测方法改写+编写的用于Matlab计算的代码。里面一共有8个function,都可直接调用,正好对应了书中的8种检验方法:

1 Movet  滑动t检验

2 Cramert  克拉默法

3 Yamamoto  山本法

4 MannKendall  M-K法

5 Pettitt  佩蒂特法

6 LePage  勒帕热法

7 BGSA  B-G分割算法

8 wvlt_trans 小波变换法

每个函数的介绍和用法都有精心撰写的详细注释,help 函数名 即可看到,不过由于想让程序更国(zhuang)际(bi)一点,因而注释都写成了英文。导入原序列后,每个程序都能算出用于检验的序列,还能直接画出检验图(当然如果不想画图,只要算出来的序列,可以把程序里面figure前面注释掉的return注释取消)。这里附一张B-G分割算法画的图。用法:下载->解压->放在某个路径下->set path->可以直接调用了。这里补了一个我测试用的序列,比如对于BGSA,导入数据后输入:[ FLAGS ] = BGSA( year,zscore(series),.95,20 );

当然,突变检验远不止这些方法,如果大家还有什么其它好的检验方法可以提在下面,附上链接什么的,有程序当然更好,我争取整合进来,壮大这一包代码。

mutation_test_V2.0修改了之前的一些错误;增进了函数调用的灵活性(可以不必对所有输入的变量赋值,比如对于BGSA,直接[ FLAGS ] = BGSA( series );都行,程序会给这些参数赋上缺省值default。当然,具体设置视情况而定)接下来再来个高级一点的,mutation_test_menu这包增加了一个可以做突变检验的下拉菜单(如图2所示)。设置好路径后,命令窗口输入mt_menu; 就出来了。直接点击操作就可以,很简单。弹出的对话框用于设置各种参数。(注意:在计算画图之前一定要将序列变量名赋成y,对应的时间变量名赋成x!)

pettitt检验r语言_[原创]Matlab气候突变检测程序合集(更新了)相关推荐

  1. pettitt检验r语言_基于Kriging方法和Pettitt检验的数据冲突判别

    基于 Kriging 方法和 Pettitt 检验的数据冲突判别 李建勋 ; 唐子豪 ; 张永进 ; 王小江 [期刊名称] <计算机工程与应用> [年 ( 卷 ), 期] 2020(056 ...

  2. R语言构建xgboost模型并评估模型(测试集、训练集每一轮):误分类率指标(misclassification rate)、logloss

    R语言构建xgboost模型并评估模型(测试集.训练集每一轮):误分类率指标(misclassification rate).logloss 目录

  3. mk突变点检测_气候突变检测

    1.气候突变是什么意思 科学家们所说的气候变化,通常指的是"气候渐变".也就是我们经常说的,"地球在逐渐变暖",或者是"气候在慢慢变暖".但 ...

  4. MATLAB有趣或有用小程序合集

    最近编辑时间:2021-5-11 这是一些比较有趣,但是感觉不值得单独发文章的小程序合集: 不定时更新,大家可以闲的没事来淘个宝 这里写目录标题 1.字幕画 2.鱼形曲线 3.鼠标沙盘 4.地图寻路演 ...

  5. stata代码命令全集:跟检验do文件、stata常用命令、代码大合集、面板熵值法

     一.常用模型代码整理 1.数据来源: 2.时间跨度:无 3.区域范围:无 4.指标说明: 包含如下模型代码: l OLS模型 l Heckman两阶段模型 l PSM+DID模型 l 固定效应模型( ...

  6. 二元置信椭圆r语言_医学统计与R语言:多分类logistic回归HosmerLemeshow拟合优度检验...

    微信公众号:医学统计与R语言如果你觉得对你有帮助,欢迎转发 输入1:multinominal logistic regression "nnet") 结果1: test (mult ...

  7. 检验杜宾 瓦森检验法R语言_应用统计学与R语言实现学习笔记(九)——线性回归...

    Chapter 9 Linear Regression 本篇是第九章,内容是回归分析(主要以线性回归为主).回归分析是数理统计.数理分析中最基础(也可以说是最重要)的一个分析,所以这一章内容相对来说也 ...

  8. 泰尔指数r语言_【数学建模】泰尔指数及分解的计算方法与Matlab实现

    补充更方便使用的R语言版本:张敬信:[R语言]泰尔指数及其分解​zhuanlan.zhihu.com 前言 最近查一点泰尔指数的资料,发现无论是公式还是软件实现,都说的特别乱,看不出所以然. 特整理了 ...

  9. pycharm中配置r语言_【R语言】R语言中的循环

    编程中减少代码重复的两个工具,一是循环,一是函数. 循环,用来处理对多个同类输入做相同事情(即迭代),如对不同列做相同操作.对不同数据集做相同操作. R语言有三种方式实现循环: (1)for循环.wh ...

最新文章

  1. pytorch torch.item()(返回此张量的值作为标准Python数字。 这仅适用于具有一个元素的张量。)
  2. mysql中一个表怎么查询多以上的信息,MySQL怎么样实现多个表的或查询?
  3. 深度学习之表面缺陷识别
  4. linux集群流程运行,linux 怎么配置集群
  5. python自动答题软件_广东开放大学(广开)线上作业自动答题python-selenium
  6. Git标签tag及tag远程同步
  7. 机器学习-朴素贝叶斯分类
  8. 35 万行代码,旷视重磅开源天元深度学习框架 ,四大特性实现简单开发
  9. android中搭建NDK环境及使用JNI技术
  10. fatal: empty string is not a valid pathspec, please use . instead if you meant to match all paths
  11. MT2503处理器性能介绍,MT2503/MT2503A/MT2503D芯片资料下载
  12. crazybox路由器解决授权码问题
  13. 原神ios android,原神ios和安卓数据互通吗 原神ios和安卓能一起玩吗
  14. 微信小程序请求本地服务器测试成功
  15. AI之路最近的一些思考
  16. import * as 用法导入模块变量
  17. 页面使用html生成一个n行n列表格,HTML静态网页:表格、表单
  18. 【命名规则】驼峰命名法
  19. 各种语言的特点和介绍
  20. 矢量切片工具:tippecanoe

热门文章

  1. 解决冲突(避免冲突)
  2. 统计学习导论(三)习题
  3. #Geek Talk# 庄明浩:直播到底是什么鬼?
  4. 10--File类字节流异常
  5. vb.net 获取打印机工作状态,是否离线,打印机名称,端口号等打印机属性信息,windows 7下亲测可用
  6. 利用python编写祝福_我用Scratch和Python编程祝福送给您Merry Christmas!
  7. 如何理解LTV(单个用户整个生命周期价值)的5种计算方法?
  8. 常用温度测量技术及其接口电路
  9. 极客日报:腾讯、字节展开拉锯战;谷歌支付 380 万美元和解此前被指不公平对待女性和亚裔;罗永浩称做锤子科技时太业余...
  10. 除了鲁大师,还有什么软件可以验机?