概念

随机游程检验是检验一个序列中的数据出现是否与顺序无关的一种方法。值得注意的是,随机游程检验的零假设是数据出现顺序随机,备择假设是数据出现顺序不随机。

实例 & 代码

在试验设计中,经常要关心试验误差 ( experiment error ) 是否与序号无关。假设有 A, B, C 三个葡萄品种,用完全试验设计需要重复测量4次,安排在12个试验田中栽种,共得到12组数据,每个试验田试验结果收成 ( 单位:kg ) 如下表所示。试问按试验田的序号,检查误差分布是否按序号随机?显著性水平为0.05。

(1) B (2) C (3) B (4) B (5) C (6) A (7) A (8) C (9) A (10) B (11) C (12) A
23 24 18 23 19 11 6 22 14 22 27 15

解答:
首先计算每个品种的均值 mean( A )=11.5 ,mean( B )=21.5 ,mean( C )=23。各试验田实际收成与各自误差成分之间出现顺序为正和负的记录为:

(1) B (2) C (3) B (4) B (5) C (6) A (7) A (8) C (9) A (10) B (11) C (12) A
+ + - + - - - - + + + +
import numpy as np
import pandas as pd
import scipy.stats as stats
from scipy.special import comb
def run_test(list_2e):arr=pd.factorize(list_2e)[0]if (2 in arr) or (1 not in arr):return '输入的列表不是二分类数据'count={}for i in arr:count[i]=count.get(i,0)+1n,n0,n1=len(arr),count[0],count[1]r=1for i in range(len(arr)-1):if arr[i+1]!=arr[i]:r+=1if r%2==0:p=2*comb(n1-1,r/2-1)*comb(n0-1,r/2-1)/comb(n,n1)else:p=(comb(n1-1,(r-1)/2-1)*comb(n0-1,(r-1)/2)+comb(n1-1,(r-1)/2)*comb(n0-1,(r-1)/2-1))/comb(n,n1)return {'游程数':r,'p值':p}
run_test(['+','+','-','+','-','-','-','-','+','+','+','+'])

结果如下:
{‘游程数’: 5,‘p值’: 0.121}
由于p值大于0.05,因此认为试验误差分布随机。

非参数统计的Python实现——随机游程检验相关推荐

  1. 【非参数统计02】单一样本的位置推断问题:符号检验、分位数推断、Cox-Staut趋势存在性检验、随机游程检验、Wilcoxon符号秩检验、位置参数区间估计、正态记分、分布一致性检验、稳健性评价

    目录导引 2 单一样本的位置推断 2.1 符号检验 2.1.2 大样本场合 2.1.3 配对样本 2.1.4 分位数检验 2.2 Cox-Staut 趋势存在性检验 2.2.1 最优权重 2.2.2 ...

  2. 大数据时代的“小数据”系列2 随机游程检验

    什么是游程检验 游程检验:是根据样本标志表现排列所形成的游程的多少进行判断的检验方法,主要用于检验一件事件的发生是否是随机的.游程检验广泛应用于工业产品的质量检验和生物制药检验等领域.随着大数据技术的 ...

  3. Python生成随机五位数——模仿手机验证码

    使用Python生成随机的五位手机验证码. # -*- coding:utf-8 -*-#生成五位随机数,模仿手机验证码#导入random库,可以生成随机数 import randomdef ran( ...

  4. python 生成随机图片

    python numpy生成随机图 # -*- coding: utf-8 -*-import numpy as np import cv2grayImage=np.random.randint(0, ...

  5. gbdt 算法比随机森林容易_用Python实现随机森林算法

    CDA数据分析师 出品 拥有高方差使得决策树(secision tress)在处理特定训练数据集时其结果显得相对脆弱.bagging(bootstrap aggregating 的缩写)算法从训练数据 ...

  6. Python 创建随机名字的文件夹/文件

    Python 创建随机名字的文件夹/文件 导入库 创建文件名 创建文件 导入库 import random import string import os 创建文件名 dir_name = ''.jo ...

  7. 用Python生成随机的中文验证码图片

    大家都在 Python论坛里讨论.交流,你还在等什么!赶快加入! 在登录很多网站的时候,他们已经不在使用简单的英文和数字的验证码,为了防止恶心注册和群发软件的侵袭,现在都开始使用中文的验证码了. 今天 ...

  8. python随机生成100内的10个整数_用python随机生成数字教程_如何用Python编程随机产生10个随机整数,并输出这10个整数的和�9�3...

    如何用Python编程随机产生10个随机整数,并输出这10个整数的和 用Python编程产生10个随机整数,并整数的和,方法如下 准备材料python.电脑 1.需要加载的模块是:random. 2. ...

  9. Python模拟随机漫步

    本文用Python模拟随机漫步行为. 1 使用内建的的random模块 import randomposition = 0 walk = [position] steps = 1000 for i i ...

  10. python大学随机点名名_python用tkinter实现一个简易能进行随机点名的界面

    今天学到Python 有一个名为 random 的内置模块,可用于生成随机数, 这个好玩~可以用来做为上课随机点名的程序了哈哈.. def main_cor(): import random i = ...

最新文章

  1. 5 种将死的编程语言!
  2. Linux基础——linux和gcc的介绍
  3. 轨迹系列8——记某真实项目中轨迹展示查询效率优化方案一(初步设计)
  4. [转]MyBatis的foreach语句详解
  5. python替代php,Python架构的PHP替代方案
  6. 南昌网络赛E Interesting Series
  7. charles代理手机调试_H5开发 移动端 调试之 Charles 抓包 和 Map Remote
  8. 三角剖分求多边形面积的交 HDU3060
  9. 减治法解决八枚硬币问题/假币问题(JAVA)----二分,三分,不知轻重的情况
  10. 前端怎么画三角形_WEB前端-CSS三角形绘制方法
  11. 分区工具parted的详解及常用分区使用方法
  12. linux下gdb调试查看源代码
  13. Groovy 教程-- Groovy 入门
  14. Spark学习:Win10中编译Spark源码
  15. geoserver发布shp格式的图层 实现步骤(含图).doc
  16. 选手的名次 C语言实现
  17. 【web性能】获取web各个阶段响应时间:服务器响应时间、首页白屏时间、dom渲染完成时间等
  18. 沈航C语言上机实验题答案,2017年沈阳航空航天大学航空航天工程学部823C语言程序设计考研题库...
  19. 对创业团队的一点想法
  20. 百度网盘不限速版正式推出,35M/S,附iOS/安卓双端下载

热门文章

  1. Linux常用命令(超详细)
  2. jmeter录制 过滤_Jmeter录制pc脚本
  3. 常规英文字体 电商_电商设计中常见的中英字体的应用
  4. android View的层级结构
  5. java 无限循环语句_java学习之循环语句
  6. 2021-07-02 swift大礼包
  7. unity三维地形制作
  8. 趣谈网络协议——UDP协议
  9. java 16进制与汉字_Java汉字与16进制数相互转换
  10. C语言从入门到精通需要多久?你真的知道吗!