2. 单因素实验设计基础

在单因素实验中，“因素”是被主动改变的那个变量（如温度、催化剂种类、饲料配方），它的不同取值称为“水平”。实验设计要解决的核心问题是：如何把有限次试验安排到实验单元上，使得水平效应能被干净地估计出来，同时把已知或未知的干扰因素尽量隔离。完全随机设计、随机区组设计和拉丁方设计，正是按“对干扰源的控制程度”由弱到强排列的三种常见方案。

一、完全随机设计：最基础、自由度最大

完全随机设计（Completely Randomized Design, CRD）只做一件事：把全部实验单元随机分配给各处理水平，不设任何区组。例如考察三种饲料配方 A、B、C 对仔猪日增重的影响，共 18 头仔猪，则每种配方随机分配 6 头，谁用哪种配方完全由随机数决定。

它的适用前提是实验单元彼此同质：品种、日龄、初始体重、栏舍条件尽量一致。此时处理间的差异可以较有把握地归因于水平差异。CRD 的优点是安排简单、缺失数据不致命、误差自由度最大（本例误差 df = 18 − 3 = 15），检验灵敏度高；缺点是当单元本身存在明显梯度（如体重差异大、土壤肥力不均、设备老化程度不同）时，这些差异会全部混入误差项，导致误差方差被抬高、检验功效下降。

分析上，CRD 对应单因素方差分析，模型可写为：

y_ij = μ + τ_i + ε_ij
i = 1..a（水平数），j = 1..n_i（重复数）
ε_ij ~ N(0, σ²) 独立
检验 H0: τ_1 = τ_2 = ... = τ_a = 0

实操中要特别注意随机化的落实：不能“前 6 头给 A、中间 6 头给 B”，那等于把顺序或位置变成了隐含区组。应使用随机数表或随机数生成器产生分配序列，并保留随机化记录以便复核。

二、随机区组设计：把已知干扰“扣”出来

当实验单元存在一个明确的、可事先测量的系统梯度时，就应使用随机区组设计（Randomized Complete Block Design, RCBD）。做法是：先按该梯度把单元分成若干个“区组”，使区组内尽量同质、区组间差异明显；然后在每个区组内，把各水平各自随机分配一次。这样每个区组都完整包含所有水平，故称“完全区组”。

仍以饲料配方为例：若 18 头仔猪初始体重从 25 kg 到 45 kg 不等，可按体重分为 6 个区组，每区组 3 头体重接近的仔猪，区组内随机分配 A、B、C。此时模型变为：

y_ij = μ + τ_i + β_j + ε_ij
τ_i：处理（配方）效应
β_j：区组（体重层次）效应
ε_ij ~ N(0, σ²)

关键变化在于：原本混在误差里的体重差异被 β_j 吸收，误差平方和随之减小，处理效应的检验更有力。代价是区组占用了自由度（本例误差 df = (3−1)(6−1) = 10），因此区组不宜过多，一般以 4～8 个为宜，且区组划分必须依据实验前可测的变量，不能事后按结果分组。

RCBD 有一个硬性约束：每个区组必须能容纳全部水平。若水平数很多而区组容量有限（如一块地只能种 5 个品种却有 8 个品种要比较），完全区组就做不成，此时需转向不完全区组或拉丁方。

三、拉丁方设计：同时控制两个方向的干扰

拉丁方设计（Latin Square Design, LSD）在行、列两个方向上同时设置区组，要求水平数 = 行数 = 列数，每个水平在每一行、每一列中各出现且仅出现一次。典型场景是田间试验：土壤肥力沿东西方向有梯度（行），沿南北方向也有梯度（列），用 4×4 拉丁方安排 4 个品种，就能把两个方向的肥力差异同时扣除。

其模型为：

y_ijk = μ + τ_i + ρ_j + κ_k + ε_ijk
τ_i：处理效应；ρ_j：行区组；κ_k：列区组
误差 df = (t−1)(t−2)，t 为水平数

以 t = 4 为例，总试验次数 16，误差 df = 3×2 = 6。可以看到，拉丁方用两个方向的区组换来了更强的干扰控制，但误差自由度下降很快：t = 3 时误差 df 仅为 2，检验几乎无力；t = 5 时误差 df = 12 才比较可用。因此拉丁方通常要求 t ≥ 4，且行、列与处理之间不存在交互作用——这是它最容易被忽视的假设。如果行因素（如肥力梯度）与品种存在交互，拉丁方的估计就会有偏。

四、如何选择：看干扰源的数量与形态

三种设计的选择逻辑可以归结为一句判断：数一数有几个方向上的系统差异需要排除。没有可识别的梯度，用 CRD；有一个明确梯度（体重、批次、时间先后、操作员），用 RCBD；有两个相互独立且都不可忽略的梯度（田间两个方向的肥力、培养箱内位置与批次），且水平数不少于 4，才考虑拉丁方。

下表给出一个便于对照的决策参考（示例参数，用于说明量级关系）：

设计类型	控制方向	示例规模	误差自由度	典型适用场景
完全随机 CRD	无	a=3，n=6，N=18	15	同批同质动物、均匀试验台
随机区组 RCBD	1 个	a=3，b=6，N=18	10	初始体重分层、设备批次分层
拉丁方 LSD	2 个	t=4，N=16	6	双向肥力梯度、双维度位置效应