Please enable JavaScript to view this site.

包含定量因子的方差分析

当其中一个因子为定量因子时,有时会使用双因素方差分析,例如在比较时间进程或剂量反应曲线时。在这些情况下,其中一个因子是剂量或时间,您将其设定为若干(或许多)离散值之一。

方差分析完全不考虑时间点(或剂量)的顺序。请仔细思考这一点。您的实验目的可能正是为了观察某种趋势或剂量-反应关系。  但方差分析的计算完全忽略了时间点或剂量的顺序。如果您随机打乱时间点或剂量,双因素方差分析得出的结果将完全相同。方差分析对待不同的时间点或不同的剂量,与对待不同的药物、不同的基因型或不同的国家完全一样。

既然当其中一个因素是定量因素时,方差分析会忽略实验的全部意义,建议考虑采用替代分析方法(如回归分析)。在某些情况下,若数据不足或理论基础不够充分以拟合曲线,方差分析仍可能是一种合理的分析方法。

含定量因子的 P 值解读

假设您在六个时间点上比较两种处理方案。  

双因素方差分析将报告三个P值:

其中一个 P 值用于检验“时间对结果无影响”的零假设。在许多情况下,您知道结果会随时间变化,这正是您进行时间序列分析的原因。由于预期时间效应的 P 值较小,因此该结果并无太大参考价值。

另一个 P 值用于检验“平均而言,两种处理方式无差异”这一零假设。在某些情况下,检验该假设或许有意义。但在多数情况下,您预期早期时间点无差异,仅关注后期时间点的差异。在此类情境下,检验平均处理效应可能帮助不大。

第三个P值用于检验交互作用。其零假设是:各时间点上不同处理组之间的差异完全相同。但若在时间零点或早期时间点收集数据,您并不期望此时发现任何差异。 您的实验设计真正关注的是后期时间点。在此情境下,您预期存在交互作用,因此交互作用的P值较小并不能帮助您理解数据。如果处理组之间的差异在某些时间点变大,而在后期时间点又变小,这种P值就更无用了。

含定量因子的多重比较检验的解读

那么,关于多重比较检验又该如何处理?

有些科学家希望确定在哪个最低剂量(或时间点)下,响应的变化具有统计学显著性。 多重比较检验可以给出答案,但这个答案依赖于样本量。增加受试者数量,或为每条曲线增加更多剂量或时间点,答案就会改变。如果样本量足够大(在每个剂量或时间点),您会发现极小的剂量或非常早的时间点上存在统计学显著(但在生物学上微不足道)的效果。 若每个剂量或时间点的重复次数较少,则需等到更高剂量或更晚的时间点才会观察到统计学显著性。由于追问能产生“显著”效应的最小剂量并非针对该系统提出根本性问题,因此结果可能无实际帮助。

若想了解最小有效剂量,不妨考虑寻找能产生大于某阈值效应的最小剂量 - 该阈值应基于生理学(或其他科学背景)设定。例如,寻找能使脉搏率每分钟增加超过10次跳动的最小剂量。这种方法能得出有用的答案,而寻找导致“显著”效应的最小剂量则无法做到。

如果您考察所有多重比较检验(而不仅仅是询问哪个剂量或时间点能产生“显著”效应),可能会得到毫无意义的结果。 您可能会发现,在第3、5、6和9个时间点差异具有统计学显著性,但在第1、2、4、7、8和10个时间点则不显著。该如何解读这种现象?仅知晓治疗在哪些剂量或时间点具有统计学显著性,通常既无法帮助您理解系统的生物学机制,也难以指导新实验的设计。

双因素方差分析的替代方法

双因素方差分析的替代方案是什么?

若采用重复测量设计,可考虑使用这种重复测量方差分析的替代方案,Will G Hopkins将其称为“受试者内建模”

首先,以某种具有生物学意义的方式量化每个受试者的数据。这可能是曲线下面积,也可能是峰值水平,或是达到峰值所需的时间。或许您可以通过非线性回归拟合曲线,从而确定速率常数或斜率。

接下来,将这些数值(面积或速率常数等)用于受试者组间的比较,若仅含两组处理,可采用t检验;若含三组或更多,则采用单因素方差分析(ANOVA)。与双因素方差分析不同,此类分析遵循实验的科学逻辑,因此得出的结果易于理解,并能引导您迈向下一步(设计更优的实验)。

若未采用重复测量设计,仍可针对每种处理方案拟合曲线。随后作为线性和非线性回归的一部分,比较斜率、EC50值或滞后时间。

请深入思考您的科学目标,并设法让统计检验与科学目标相匹配。在许多情况下,您会发现比使用双因素方差分析更合适的方法。

倾向测试

在单因素方差分析之后进行多重比较检验时,线性趋势检验是一种可选方案。当然,该检验确实会考虑处理顺序。其他软件(但不包括Prism)提供多项式后检验,该方法同样会考虑处理顺序。