不完全区组设计和统计分析1

SSTATISTICSADVANCEDEXPERIMENTALDESIGNLECTURE01不完全区组设计和统计分析1高等试验设计课程·研究生/高年级本科BALANCEDINCOMPLETEBLOCKDESIGNv处理数b区组数k区组容量λ共现次数λ(v−1)=r(k−1)BBLOCKDESIGNCONTENTS课程内容概览不完全区组设计·六章结构导览01试验设计基础与不完全区组动机完全区组设计的局限性与不完全区组的必要性不完全区组设计的分类体系与应用场景02平衡不完全区组设计BIBDBIBD的五参数定义与基本关系式Fisher不等式与存在性必要条件03BIBD的构造方法对称BIBD、循环BIBD与差集法有限几何法与递归构造技术04其他不完全区组设计PBIBD、Youden设计与Lattice设计各类设计的适用条件与比较05统计分析与方差分析线性模型、调整处理均数与ANOVA表效率因子与多重比较方法06应用实例与总结农业、工业、医学领域真实案例课程核心要点回顾与展望CHAPTER01试验设计基础与不完全区组动机从完全区组到不完全区组的认知跃迁H不完全区组设计COMPARISON核心差异完全区组与不完全区组的核心差异完全区组设计假设每个区组能容纳全部处理，当处理数v超过区组容量k时该假设失效；不完全区组设计通过允许k<v来适应现实约束，但需以调整处理均数和更复杂的方差分析为代价换取可行性，其核心价值在于用结构化缺失替代不可行的完全观测。完全区组设计的理想与局限完全区组要求每个区组包含全部v个处理，当v较大或区组容量受物理限制时无法实现田间试验中地块面积、临床试验中患者数量等现实约束常使完全区组不可行强行扩大区组会导致区组内异质性增加，反而降低试验精度与估计效率不完全区组设计的应对策略允许每个区组仅包含k<v个处理，通过精心设计保证处理间比较的公平性以调整处理均数和更复杂的方差分析为代价，换取试验的现实可行性用结构化的缺失模式替代不可行的完全观测，是试验设计理论的重要进步ADVANTAGES不完全区组设计的四大核心优势不完全区组设计在区组容量受限的条件下，通过结构化安排同时实现经济性、平衡性、灵活性与计算严密性：经济性地解决试验成本问题，平衡性地保证处理比较公平，灵活性地适应不同区组规模，严密性地通过数学方法消除系统误差，四者缺一不可。经济性全部处理无需安排在同一区组内，对区组容量要求较低，经济地解决了试验成本与资源约束问题平衡性每个处理重复次数相同且每对处理相遇次数相同，保证了处理间比较的公平性与无偏性灵活性可根据实际区组容量k灵活调整设计方案，适应田间、临床、工业等不同场景的物理限制计算严密性有严格的数学方法有效消除系统误差，试验精度高，不因区组不完全而牺牲推断可靠性BBLOCKDESIGNCLASSIFICATIONTAXONOMY不完全区组设计的分类体系不完全区组设计按平衡程度分为BIBD、PBIBD与特殊设计三类：BIBD要求每对处理相遇次数λ恒定，平衡性最优但存在性条件严格；PBIBD放宽为多类关联，适应性更强；Youden与Lattice等特殊设计针对双向控制或特定处理数结构优化，各有专属应用场景。01平衡不完全区组设计BIBD每对处理恰好在λ个区组中相遇，平衡性最优，方差分析最简洁参数需满足vr=bk与λ(v-1)=r(k-1)，存在性条件严格广泛应用于农业品种比较、工业配方筛选等处理数适中的场景02部分平衡不完全区组设计PBIBD允许不同处理对相遇次数不同，通过关联方案维持部分结构当BIBD不存在或参数过大时提供可行的替代方案分析较BIBD复杂，需引入关联类与效率因子分层评价03Youden设计与Lattice设计Youden设计是不完全拉丁方，同时控制行与列两个区组因素Lattice设计要求v=k²或k(k+1)，利用格子结构简化构造与分析两者均为BIBD的特例或扩展，适用于特定约束条件下的试验CHAPTER02平衡不完全区组设计BIBD定义、参数关系与存在性条件DEFINITIONBIBD的三重定义条件BIBD由三个条件严格定义：每个区组含k<v个处理、每个处理重复r次、每对处理相遇λ次。三者共同保证设计的平衡性与可分析性，其中λ恒定是"平衡"的核心内涵，使所有处理对比较具有相同方差，是BIBD区别于其他不完全区组设计的本质特征。01区组容量约束处理重复均衡每个区组恰好包含k个处理且k<v，这是"不完全"的定义前提，若k=v则退化为完全区组设计每个处理在整个试验中恰好出现r次，保证所有处理的估计精度相同，避免某些处理被过度或不足采样02配对平衡性λ恒定原则任意两个不同处理恰好在λ个区组中同时出现，λ恒定是"平衡"的本质，使所有处理对比较方差相等平衡性意味着无论比较哪两个处理，其差值的方差均为2kσ²/(λv)，实现公平比较03参数符号体系五参数完整描述五个参数(v,b,r,k,λ)完整描述一个BIBD，后续推导均基于这三重定义条件的代数表达参数间存在两个必要等式：bk=vr与λ(v−1)=r(k−1)，确定任意三个即可推出其余PARAMETRICRELATIONS§2.1BIBD的基本参数关系式BIBD的两个基本参数关系vr=bk与λ(v-1)=r(k-1)是存在性的必要条件：前者源于总观测次数的双重计数，后者源于单个处理配对数的两种计算方式。满足这两个等式的参数组合仍可能无法构造实际设计，还需Fisher不等式等进一步约束。vr=bk总观测双重计数λ(v−1)=r(k−1)配对平衡代数表达关系式推导逻辑含义vr=bk处理视角总观测数=区组视角总观测数总试验单元数的双重计数恒等式λ(v−1)=r(k−1)单处理与其他v−1处理配对总数=该处理所在r个区组内配对数配对平衡性的代数表达r=λ(v−1)/(k−1)由第二式解出r重复数由v,k,λ唯一确定b=vr/k=λv(v−1)/[k(k−1)]结合两式解出b区组数由v,k,λ唯一确定THEORYFisher不等式与设计存在性Fisher不等式b≥v是BIBD存在的核心必要条件，表明区组数不得少于处理数。当b=v时设计对称且r=k，结构最优；该不等式限制了小区组设计的可行性，提醒试验者不能仅凭参数关系判断设计是否存在，还需验证此结构性约束。SECTION01Fisher不等式的表述与证明思路核心断言：Fisher不等式断言任何BIBD必须满足b≥v，即区组数不少于处理数，是存在性的关键必要条件证明思路：基于关联矩阵N的秩——NN'的特征值为rk（重数1）和r−λ（重数v−1），故rank(N)=v≤b对称情形：当b=v时设计称为对称BIBD，此时自动有r=k，且任意两个区组也恰好交于λ个处理SECTION02对试验设计选择的实践指导验证约束：不能仅凭vr=bk与λ(v−1)=r(k−1)判断设计可行，必须同时验证b≥v这一结构性约束参数调整：当计算出的b<v时，需增大r或调整k使b≥v成立，否则该参数组合无法实现平衡设计最优与存在：对称BIBD（b=v）在信息效率上最优，但并非所有参数都存在，需借助Bruck-Ryser-Chowla定理进一步判定BIBD·FISHERINEQUALITYMATRIXTHEORY关联矩阵的代数结构与性质BIBD的关联矩阵N是v×b的0-1矩阵，满足行和r、列和k、行间内积λ；NN'=(r-λ)I+λJ的简洁形式将组合设计转化为线性代数对象，其特征值rk与r-λ直接决定设计的秩、效率与信息量，是连接组合理论与统计分析的核心桥梁。关联矩阵N的核心代数性质性质数学表达统计/组合含义行和恒定∑ⱼnᵢⱼ=r(∀i)每个处理重复r次列和恒定∑ᵢnᵢⱼ=k(∀j)每个区组含k个处理行间内积恒定∑ₕnᵢₕnⱼₕ=λ(i≠j)每对处理共现λ次NN'的结构NN'=(r-λ)I+λJ特征值rk（1重）与r-λ（v-1重）秩与Fisher不等式rank(N)=v≤bb≥v的代数证明基础小结：关联矩阵将BIBD的组合定义完全编码为线性代数性质，是理论分析与统计推断的统一基础。B不完全区组设计SYMMETRICBIBDCOREPROPERTIES对称BIBD的对偶结构与存在性对称BIBD满足b=v且r=k，具有处理与区组的完美对偶性：每对区组也交于λ个点。其存在性除基本参数关系外还需满足Bruck-Ryser-Chowla定理的数论条件，有限射影平面是其λ=1的特例，体现了组合设计与代数几何的深刻联系。01对称性与对偶性质的体现定义特征：b=v且r=k是对称BIBD的核心标志，此时关联矩阵N为方阵且可逆，满足NN'=N'N=(r-λ)I+λJ，矩阵的对称性直接映射了设计的结构对称性对偶性质：不仅每对处理相遇λ次，每对区组也恰好交于λ个处理，处理与区组在组合结构中地位完全对称，对偶设计仍为同一参数的BIBD02Bruck-Ryser-Chowla定理与实例数论约束：v为偶数时n=k-λ必须是完全平方数；v为奇数时x²=n·y²+(-1)(v-1)/2λz²须有非零整数解，这是超越基本参数关系的深层存在性条件经典实例：(7,3,1)对应Fano平面、(13,4,1)对应射影平面阶3；而(22,7,2)被BRC定理证明不存在，展示了数论对组合设计的深刻约束H不完全区组设计SectionChapter03BIBD的构造方法差集、有限几何、递归与直接构造01差集法02有限几何法03递归构造04直接构造CONSTRUCTIONMETHODALGEBRAICDESIGN差集法构造对称BIBD差集法将对称BIBD构造转化为群论问题：(v,k,λ)-差集的全体平移构成对称(v,k,λ)-BIBD。该方法仅适用于对称设计，但提供了最简洁的代数构造途径，经典实例包括(7,3,1)、(13,4,1)、(21,5,1)等，其存在性与有限域上的二次剩余密切相关。差集的定义与核心定理差集定义：(v,k,λ)-差集是v阶群G中k元子集D，使{x−y|x,y∈D,x≠y}覆盖G\{0}中每个元素恰λ次核心定理：阿贝尔群G中的(v,k,λ)-差集D，其全体平移Dev(D)构成对称(v,k,λ)-BIBD必要条件：差集存在要求λ(v−1)=k(k−1)，与对称BIBD参数关系一致D={d1,d2,…,dk}⊂G⟹|{di−dj=g}|=λ∀g≠0经典实例与构造步骤(13,4,1)差集：Z₁₃中D={0,1,3,9}，12个非零差各出现1次，循环平移得13个区组更多实例：Z₇中D={0,1,3}为(7,3,1)-差集（Fano平面）；Z₂₁中存在非阿贝尔群差集构造步骤：选定群G→验证差集条件→生成全部平移→按行排列即为BIBD区组方案EXAMPLE·Z₁₃DIFFERENCESETD={0,1,3,9}→Dev(D)=13blocks{0,1,3,9}{1,2,4,10}{2,3,5,11}…{12,0,2,8}DIFFERENCESETMETHODBBIBDMETHODSFINITEGEOMETRYCONSTRUCTION有限几何法构造BIBD有限射影平面PG(2,q)与仿射平面AG(2,q)天然提供BIBD实例：前者给出对称(q²+q+1,q+1,1)-BIBD，后者给出可分解(q²,q,1)-BIBD。当q为素数幂时有限域保证存在性，几何直观（点=处理、线=区组、两点一线=λ=1）使抽象设计获得空间解释。PROJECTIVE射影平面与对称BIBDPG(2,q)有q²+q+1个点与等量条线，每条线含q+1个点，每对点恰在1条线上，构成对称(q²+q+1,q+1,1)-BIBDq=2得(7,3,1)即Fano平面，q=3得(13,4,1)，q=4得(21,5,1)，形成无穷族对称设计v=q²+q+1对称设计参数AFFINE仿射平面与可分解BIBDAG(2,q)有q²个点、q(q+1)条线，每条线含q个点，每对点恰在1条线上，构成(q²,q,1)-BIBD且可分解为q+1个平行类q=3得(9,3,1)，q=4得(16,4,1)；可分解性意味着区组可分组使每组覆盖全部处理，便于试验实施v=q²可分解设计参数几何直觉映射：点→处理线→区组两点一线→λ=1H学术讲义METHODOLOGYCONSTRUCTION循环BIBD的生成机制循环BIBD由初始区组模v循环平移生成全部区组，存储仅需记录初始区组而非完整设计。其存在性等价于循环差集的存在性，当v为素数时二次剩余常提供解。循环结构极大简化了试验实施的随机化与数据录入，是实际应用中最常用的BIBD类型之一。Z₇循环示意0123456B₀={0,1,3}初始区组（高亮节点）参数特征7个处理·7个区组λ=1每对恰相遇一次k=3区组容量TABLE(7,3,1)-循环BIBD的生成示例步骤操作生成区组初始区组选定基区组{0,1,3}{0,1,3}第1次平移{0,1,3}+1mod7{1,2,4}第2次平移{0,1,3}+2mod7{2,3,5}第3次平移{0,1,3}+3mod7{3,4,6}第4次平移{0,1,3}+4mod7{4,5,0}第5次平移{0,1,3}+5mod7{5,6,1}第6次平移{0,1,3}+6mod7{6,0,2}单个初始区组循环7次生成全部7个区组，每对处理恰相遇1次，验证λ=1METHODOLOGYRECURSIVECONSTRUCTION递归构造：从小设计到大设计递归构造通过直积、补集、嵌套等操作从已知小BIBD生成更大参数的设计，是建立无穷族BIBD存在性的核心策略。01直积与补集构造策略直积构造：将(v₁,k₁,λ₁)与(v₂,k₂,λ₂)-BIBD的区组做笛卡尔积，得到(v₁v₂,k₁k₂,λ₁λ₂)-BIBD补设计：(v,k,λ)-BIBD的补集是(v,v-k,b-2r+λ)-BIBD，可从已知设计快速获得互补参数两种操作均可保持平衡性，且能反复应用以生成参数越来越大的设计族02嵌套扩张与参数表编制嵌套构造在小设计中添加新点和新区组，逐步扩展参数，适用于无法直积的参数组合递归法是编制BIBD参数存在性表的理论基础，表中多数条目可通过少量种子设计递归生成实际应用中优先查表获取已知设计，仅在参数新颖时才需手动构造或计算机搜索BBIBD学术讲义METHODOLOGYComparisonBIBD构造方法对比与选择指南四种BIBD构造方法各有适用边界：差集法限于对称设计但最简洁，有限几何法提供无穷族但依赖素数幂阶，循环法实施便捷但存在性受限，递归法通用但构造复杂。实际选择应先查参数表，再据参数特征匹配方法，小参数优先查表，大参数依赖代数或递归。BIBD构造方法对比方法适用条件优点局限典型参数差集法对称BIBD，阿贝尔群构造最简洁，一区组生全设计仅限对称设计，差集存在性难判(7,3,1),(13,4,1),(21,5,1)有限几何法q为素数幂提供无穷族，几何直观强参数受限于q²+q+1或q²形式(7,3,1),(9,3,1),(13,4,1)循环法存在循环差集存储紧凑，随机化简便非所有BIBD都循环(7,3,1),(13,4,1),(19,9,4)递归法有种子设计可用通用性强，可生成无穷族构造过程复杂，需组合技巧(63,9,1),(21,5,1)等复合参数选择建议：实际工作优先查BIBD参数表，小参数（v≤20）查表即可，大参数依代数结构选方法。CHAPTER0419其他不完全区组设计PBIBD、Youden设计与Lattice设计部分平衡不完全区组设计Youden方设计Lattice设计不完全区组设计学术讲义DEFINITION部分平衡不完全区组设计PBIBDPBIBD放宽BIBD的λ恒定条件为多类关联：处理对被分为m个关联类，第i类处理对在λᵢ个区组中相遇。它通过关联方案维持部分结构，在BIBD不存在时提供可行替代，但分析复杂度随关联类数增加而上升，效率因子分层评价。01关联方案与分类体系处理对被分为m个关联类，第i类中任意一对处理恰在λᵢ个区组中相遇，λ₁,...,λₘ不必相等，打破了BIBD中所有处理对相遇次数恒定的限制常见类型包括一类关联（强正则图）、二类关联（群可分、三角、拉丁方型等），m越大分析越复杂，信息矩阵特征值层数增加，效率因子需分层评价02与BIBD的关系及应用场景BIBD是m=1的特例（所有λᵢ相等）；当BIBD参数不满足存在性条件时，PBIBD常是唯一结构化替代方案，保留了不完全区组设计的主要统计优势适用于处理数较多且无法构造BIBD的场景，如品种比较试验中v>20且k较小时常采用PBIBD，兼顾试验规模可控与统计分析可行DESIGNTHEORYYouden设计：不完全拉丁方Youden设计是b×k(b>k)矩形阵列，每处理每列出现一次、每行至多一次，行视角构成BIBD。它同时控制行与列两个区组因素，比单向BIBD更适合双向异质性场景，可通过删去拉丁方行或从对称BIBD构造，广泛用于实验室校准与工业过程优化。定义特征与结构性质阵列结构：b×k矩形阵列（b>k），每处理每列恰出现1次、每行至多1次，行视角构成(v,b,r,k,λ)-BIBD且v=b精度优势：同时消除行与列两个方向的系统变异，精度优于单向BIBD，是不完全拉丁方的标准形式构造途径与应用场景构造方法：可从完整拉丁方删去若干行获得，也可从对称BIBD直接构造（取关联矩阵的列为Youden方阵）应用场景：适用于实验室仪器校准、工业过程优化等需双向控制的场景，尤其当行列数不等时替代拉丁方EXPERIMENTALDESIGNLattice设计：格子结构的区组安排Lattice设计要求v=k²或k(k+1)，将处理排成k×k或k×(k+1)格子后按行、列、对角线方向生成区组。简单格子两重复、三重格子三重复、平衡格子多重复趋近BIBD。构造规则、实施简便、分析有简化公式，是品种比较试验中处理数为平方数时的首选设计。简单格子设计简单格子(v=k²)由两个重复组成：第一重复按行分组得k个区组，第二重复按列分组得k个区组，共2k个区组。每个区组包含k个处理，区组数恰好等于2k，适用于区组容量k受限而处理数v=k²的情形。三重格子设计三重格子在简单格子基础上增加对角线方向的第三重复，共3k个区组，提高了处理间比较的精度与信息量。三个方向的分组使处理间的有效重复数增加，方差分析精度优于简单格子设计。平衡格子与设计选择平衡格子通过足够多重复使每对处理相遇次数接近常数，趋近BIBD性质，分析可用BIBD近似公式。当v=k²且k≤10时优先选Lattice设计而非通用BIBD，因其构造规则、田间布局直观、分析简化。BBlockDesignCOMPARISON§Comparison四类不完全区组设计对比选择四类设计各有适用边界：BIBD平衡最优但存在性严，PBIBD灵活但分析复杂，Youden适合双向控制且行列不等，Lattice专攻v=k²情形且实施简便。选择时应据v与k的关系、区组因素个数、分析便利性及参数存在性综合权衡，无绝对优劣之分。不完全区组设计类型对比设计类型平衡程度适用条件分析复杂度典型应用BIBD完全平衡(λ恒定)vr=bk,λ(v-1)=r(k-1),b≥v中等，标准ANOVA品种比较、配方筛选PBIBD部分平衡(m类关联)BIBD不存在或参数过大时较高，需关联类分析大处理数品种试验Youden行BIBD+列完全需双向控制且行列数不等中等，双向调整实验室校准、工业优化Lattice趋近平衡(重复增多)v=k²或k(k+1)较低，专用简化公式品种比较(v为平方数)选择依据v与k关系、区组因素个数、分析便利性、参数存在性四维权衡INCOMPLETEBLOCKDESIGNSH不完全区组设计StatisticalAnalysisCHAPTER05统计分析与方差分析线性模型、调整均数与效率因子处理未在所有区组中出现，直接比较处理总和会产生偏倚——本章系统讲解BIBD的调整处理均数、ANOVA表构建与F检验方法。LINEARMODELBIBD的线性模型与假设体系BIBD线性模型Yᵢⱼ=μ+τᵢ+βⱼ+εᵢⱼ形式上与完全区组相同，但数据不完全导致最小二乘估计需解正规方程而非简单均值。核心假设为误差独立同正态分布，检验目标为H₀:τ₁=...=τᵥ，模型的可识别性依赖于Στᵢ=0或类似约束条件。01模型结构与参数含义Yᵢⱼ=μ+τᵢ+βⱼ+εᵢⱼ，仅当处理i出现在区组j时有观测，数据矩阵存在结构性缺失μ为总均值，τᵢ为第i个处理效应，βⱼ为第j个区组效应，εᵢⱼ~NID(0,σ²)为随机误差模型需加约束Στᵢ=0或Σβⱼ=0以保证可识别性，否则参数不唯一02假设检验目标与估计挑战核心检验H₀:τ₁=τ₂=...=τᵥ，备择假设为至少一对τᵢ≠τⱼ，采用F检验判断处理差异显著性因数据不完全，τᵢ的最小二乘估计不能直接用处理均值，需通过调整处理总量Qᵢ间接获得调整的本质是从处理总量中扣除区组效应的贡献，使不同处理的比较基于相同的区组基准FORMULATIONBALANCEDINCOMPLETEBLOCKDESIGN调整处理总量的计算与解释调整处理总量Qᵢ=Tᵢ−(1/k)∑Bᵒ从处理总和中扣除区组效应贡献，是消除区组偏差后的纯处理信息。E(Qᵢ)=(═v/k)τᵢ，故τ̂ᵢ=kQᵢ/(═v)为无偏估计。该公式的简洁性源于BIBD的平衡性，是不完全区组分析可行性的代数基础。调整处理总量Qᵢ的计算步骤步骤符号计算公式含义01处理总和TᵢTᵢ=∑Yᵢᵒ（对含处理i的区组求和）处理i的全部观测值之和，含区组效应02区组总和BᵒBᵒ=∑Yᵢᵒ（对区组j内所有处理求和）区组j的全部观测值之和03区组调整项(1/k)∑Bᵒ对含处理i的r个区组求和后除以k处理i所在区组效应的加权平均04调整处理总量QᵢQᵢ=Tᵢ−(1/k)∑Bᵒ扣除区组效应后的纯处理信息05处理效应估计τ̂ᵢτ̂ᵢ=kQᵢ/(═v)τᵢ的无偏最小二乘估计小结：Qᵢ是BIBD分析的核心中间量，其期望与τᵢ成正比，比例系数由设计参数(═,v,k)决定。公式的简洁性正是不完全区组设计平衡性的代数体现。ANOVABIBDANALYSISBIBD方差分析表的结构与计算BIBD的ANOVA将总变异分解为区组、调整处理与误差三部分，处理SS必须用Qi计算：SSTreat(adj)=(k/λv)ΣQi²。F=MSTreat(adj)/MSError检验处理差异，误差df=bk-b-v+1。直接使用未调整的Ti会导致处理SS有偏。BIBD方差分析表变异来源自由度(df)平方和(SS)均方(MS)F值区组（未调整）b−1SSBlock=(1/k)ΣBj²−CFMSBlock=SSBlock/(b−1)—处理（调整后）v−1SSTreat(adj)=(k/λv)ΣQi²MSTreat(adj)=SSTreat(adj)/(v−1)F=MSTreat(adj)/MSError误差bk−b−v+1SSError=SSTotal−SSBlock−SSTreat(adj)MSError=SSError/dfError总计bk−1SSTotal=ΣYij²−CFNOTECF=G²/(bk)为校正因子，G为总观测和；处理SS必须用调整后Qi，不可用原始Ti核心公式SSTreat(adj)=(k/λv)ΣQi²检验统计量F=MSTreat(adj)/MSError误差自由度dfError=bk−b−v+1TheoryEFFICIENCYFACTOR效率因子的定义与实践解读BIBD效率因子E=vλ/[r(k-1)]衡量相对于完全区组的信息保留比例，E=1为完全区组，E越低信息损失越大。它直接决定处理对比方差Var(τ̂ᵢ−τ̂ⱼ)=2σ²/(rE)，是设计选择与结果解释的关键指标。实践中E>0.8可接受，E<0.5需增加重复或换设计。01效率因子的理论定义与计算公式与含义：E=vλ/[r(k−1)]=v(k−1)/[k(v−1)]，表示BIBD处理对比方差与同等重复完全区组方差的比值倒数边界条件：E=1当且仅当k=v（完全区组）；k<v时E<1，且k越小E越低，信息损失越大矩阵视角：也可从信息矩阵C的非零特征值角度定义：E为规范效率因子的调和平均02实践解读与设计选择指导方差影响：Var(τ̂ᵢ−τ̂ⱼ)=2σ²/(rE)=2k·MS_Error/(λv)，E越低方差越大、置信区间越宽、检验功效越低经验阈值：实践中E>0.8通常可接受，0.5<E<0.8需谨慎解释，E<0.5建议增加重复数r或改用其他设计权衡判断：小区组虽降低E但提高区组内同质性（σ²减小），实际精度需综合E与σ²的变化权衡判断多重比较MULTIPLECOMPARISONSBIBD的多重比较方法与标准误BIBD所有处理对标准误相同：SE=√[2k·MSE/(λv)]，这是平衡性的