恳请帮助：我现在有E，T，R，B四个因素，表型为y，模型中包含E因素，T因素，T与E的互作，R因素（R嵌套在E下），B因素（B嵌套在R下），请问如何用III型平方和计算各因素平方和，XE,XT,XTE,XR,XB分别为对应的设计矩阵
X0 <- cbind(1,XE,XT,XTE,XR,XB)
X1 <- cbind(1,XE,XT,XTE,XR)
X2 <- cbind(1,XE,XT,XR,XB)
X3 <- cbind(1,XE,XT,XTE)
X4 <- cbind(1,XE,XR,XB)
X5 <- cbind(1,XR,XB)
SSall <- sum((y-mean(y))^2)
##lse（y，a）是用来计算y=a模型的残差平方和的
SSe0 <- lse(y,X0)
SSe1 <- lse(y,X1)
SSe2 <- lse(y,X2)
SSe3 <- lse(y,X3)
SSe4 <- lse(y,X4)
SSe5 <- lse(y,X5)
SSb <- SSe1-SSe0
SSte <- SSe2-SSe0
SSr <- SSe3-SSe0-SSb
SSt <- SSe4-SSe0-SSte
SSee <- SSe5-SSe4
我算的T与E的互作，R因素（R嵌套在E下），B因素（B嵌套在R下）的平方和（分别为SSte，SSr，SSb）都与SAS计算的III型平方和一致，但是就是T因素(SSt)，E因素(SSee)老是算不对，请教一下，我该如何修改以上的使得SSt和SSee计算的对呢？
