求助 系统分组试验的方差分析

Ynckmu

系统分组又称多层分组、嵌套设计、窝设计等等，最近在学习过程中发现有些不明白的地方，还请各位老师指点迷津。

数据：
FishMeal        Animal        Protein
A1        B11        82.50
A1        B11        82.40
A1        B12        87.10
A1        B12        86.50
A1        B13        84.00
A1        B13        83.90
A2        B21        86.60
A2        B21        85.80
A2        B22        86.20
A2        B22        85.70
A2        B23        87.00
A2        B23        87.60
A3        B31        82.00
A3        B31        81.50
A3        B32        80.00
A3        B32        80.50
A3        B33        79.50
A3        B33        80.30

此例是测定3种不同来源的鱼粉的蛋白质消化率的，由于同种鱼粉不能同时混于饲料中喂给同一种动物，不能交叉，所以采用系统分组试验的方差分析方法，下面是分析过程：

> datas <- read.table("clipboard", header=T)
> attach(datas)
> datas.aov <- aov(Protein~FishMeal+Animal+Animal%in%FishMeal)
> summary(datas.aov)
                Df   Sum Sq    Mean Sq   F value     Pr(>F)   
FishMeal     2   105.50     52.75       342.78     3.21e-09 ***
Animal       6    25.46       4.24         27.57       2.64e-05 ***
Residuals   9    1.38         0.15                     
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

书上说一级因素(FishMeal)的F值是MS_a/MS_b(a)，在此例中应是52.75/4.24=12.44104，二级因素的F值应是MS_b(a)/MS_e，此例中应是4.24/0.15=28.2667，这两个数字与书上的是一致的，但是R的结果前者是342.78，不知道这个是怎么算出来的，请了解这方面的老师指点一下。

Ynckmu

现在看来最接近的是R均将MS_e做为了F值的分母，这个有什么影响没？