3. 实验单元与处理结构
在生物实验设计中,统计推断的合法性并不取决于用了多复杂的模型,而取决于一个更基础的问题:独立重复到底发生在哪一层。如果这一点判断错误,即使方差分析、混合模型或置换检验的代码完全正确,得到的 p 值也可能严重失真。本节围绕实验单元与观测单元的区分,梳理处理、对照、因子与水平的定义,并重点分析伪重复、嵌套与交叉结构对推断的影响。
一、实验单元与观测单元:独立性的边界在哪里
实验单元(experimental unit)是接受一次独立处理分配的最小实体。它决定了处理效应在统计上能获得多少个独立重复。若一个实验单元被随机分配了某个处理,那么该单元内所有测量值都共享同一处理条件,彼此之间不再独立。
观测单元(observational unit)是实际产生一个测量值的最小实体。它可能等于实验单元,也可能小于实验单元。例如,把 12 只小鼠随机分为 3 组,每组 4 只,分别接受药物 A、药物 B、生理盐水;每只小鼠取 3 个组织切片测量荧光强度。这里实验单元是小鼠(n=12),观测单元是组织切片(n=36)。处理分配发生在小鼠层面,因此处理效应的独立重复数是 12,而不是 36。若把 36 个切片当作独立样本做 t 检验,就犯了典型的伪重复错误。
判断实验单元的关键不是“测量了几次”,而是“随机化发生在哪一层”。随机化单位、处理施加单位和独立重复单位通常应保持一致;当它们不一致时,统计模型必须显式引入随机效应来刻画层级结构。
二、处理、对照、因子与水平
处理(treatment)是施加给实验单元的具体条件组合。它可以是单一因素的一个水平,也可以是多个因素水平的组合。因子(factor)是研究者主动操纵的自变量,例如温度、药物种类、基因型;水平(level)是因子的具体取值,例如温度因子的 20℃、25℃、30℃。
对照(control)是用于估计处理效应基准的条件。常见类型包括:空白对照(不施加处理)、溶剂对照(施加载体但不含活性成分)、阳性对照(施加已知有效处理)以及假手术对照。对照的选择直接决定效应量的解释:若缺少溶剂对照,药物效应可能被载体效应混淆。
以一个双因子实验为例:因子 A 为光照(水平:强光、弱光),因子 B 为氮素(水平:高氮、低氮),共 4 个处理组合。每个组合分配 5 个实验单元,则总实验单元数为 20。若每个单元再测 2 个观测值,总观测值为 40,但处理效应的独立重复仍为 5。这一区分是后续所有推断的基础。
三、伪重复:最常见的推断陷阱
伪重复(pseudoreplication)指把非独立的观测值当作独立重复使用,从而人为放大自由度、缩小标准误、抬高假阳性率。它有两种典型形式:
- 简单伪重复:同一实验单元内多次测量被当作多个独立样本。例如同一培养瓶取 5 次吸光度,仍只是 1 个实验单元。
- 牺牲性伪重复:处理施加在较大单元上,但测量在子单元上完成,且子单元未独立接受处理。例如整块田施一种肥料,田内取 10 株测量,独立重复仍是田块数。
识别伪重复的实用方法是画一张“随机化—处理—测量”三层图:随机化发生在哪一层,处理施加在哪一层,测量发生在哪一层。若测量层比随机化层更细,且细层内观测不独立,就必须使用嵌套随机效应或先聚合到实验单元层再分析。
四、嵌套与交叉:两种基本结构
嵌套结构(nested structure)指一个因子的水平只出现在另一个因子的特定水平之内。例如“班级”嵌套于“学校”:A 校的 1 班与 B 校的 1 班不是同一个班。嵌套结构中,子因子的效应无法与父因子完全分离,统计模型需写成 y ~ 处理 + (1|学校/班级) 或等价形式。
交叉结构(crossed structure)指两个因子的所有水平组合都出现。例如“温度”与“pH”两个因子,每个温度都测试了每个 pH,则二者交叉。交叉结构可以估计交互作用,而嵌套结构通常不能。
判断嵌套还是交叉,可以问:因子 B 的某个水平是否在因子 A 的每个水平下都出现?若都出现,是交叉;若只出现在部分水平下,是嵌套。这一判断直接决定模型公式和自由度分配。
五、结构误判对推断的实际影响
把嵌套当交叉,会把本应属于上层单元的变异错误地归入误差项,导致分母自由度虚增、标准误偏小、p 值偏小。把交叉当嵌套,则会丢失可估计的交互作用,降低检验功效。伪重复的后果更直接:它把“同一处理下的多个测量”当作“多个独立处理重复”,使误差方差被严重低估。
以一个具体示例说明:某实验比较 3 种饲料对仔猪增重的影响,每种饲料分配 4 个猪栏,每栏 6 头猪,共 72 头。若以猪为实验单元做单因素方差分析,误差 df = 72 − 3 = 69;但随机化发生在猪栏层面,正确的误差 df 应为 12 − 3 = 9。前者会把标准误缩小约 √(69/9) ≈ 2.77 倍,p 值被显著夸大。正确做法是使用 y ~ 饲料 + (1|猪栏) 的混合模型,或先按猪栏求均值再分析。
六、设计检查清单
| 检查项 | 判断依据 | 处理方式 |
|---|---|---|
| 实验单元是否明确 | 随机化发生在哪一层 | 以随机化单位作为独立重复数 |
| 是否存在伪重复 | 测量层比随机化层更细 | 引入嵌套随机效应或先聚合 |
| 因子是嵌套还是交叉 | B 的水平是否在 A 的每个水平下都出现 | 嵌套用 /,交叉用 * |
| 对照是否充分 | 是否包含溶剂对照或阳性对照 | 补充缺失对照,避免效应混淆 |
| 自由度是否匹配 | 误差 df 是否按实验单元数计算 | 重新核算 df,必要时改用混合模型 |
实验单元与观测单元的区分、处理结构的正确识别,是后续方差分析、混合模型和置换检验能够成立的前提。设计阶段多花十分钟画清随机化层级,往往能避免分析阶段数小时的返工和结论层面的系统性偏差。