非参数统计的Python实现—— Durbin 不完全区组分析法

概念

当处理组非常大，而区组中可允许样本量有限时，在一个区组中很难包含所有处理。于是出现了不完全的数据设计结构，其中较为常见的是均衡不完全区组 BIB 设计。Durbin 检验能用于不完全区组设计中。
均衡不完全区组需满足以下条件：

1. 每个区组包含 t 个水平——区组大小相等；
2. 每个处理在 r 个区组中出现——处理重复数相等；
3. 每对处理在 λ 个区组中相遇——相遇数相等。

实例 & 代码

设需要对4种饲料（处理）的养猪效果进行试验，用以比较饲料的质量。选4胎母猪所生的小猪进行试验，每头小猪体重相当，选择3头进行实验。3个月后测量所有小猪增加的体重如下表所示，试比较4种饲料品质有无差别。显著性水平为0.05。

饲料	胎别1	胎别2	胎别3	胎别4
A	73	74		71
B		75	67	72
C	74	75	68
D	75		72	75

解答：
t = 3, k = 4, r = 3, b = 4, λ = 2

import pandas as pd
import numpy as np
import scipy.stats as stats
A=[73,74,np.nan,71]
B&#