我已经按照惯例回答了我的部分问题。这是为了在MASS软件包中为厌食症数据找到可加线性模型。
对于第二部分，我有一个问题：
写出设计矩阵的一部分，该矩阵包括每个治疗级别的三名患者，并使用角点约束来赋予其权重。
我理解这涉及控制处理功能，但是我很少有示例，不足以说明这种情况在这种情况下如何工作。我将如何为我的数据集执行此操作？library(MASS)
data(anorexia)
names(anorexia)
anorexia$Treat<-factor(anorexia$Treat) ##make factores
anorex.1 <- lm(Postwt ~ Prewt + Treat + offset(Prewt), data = anorexia)
plot(anorex.1)
summary(anorex.1)
anova(anorex.1)
max(cooks.distance(anorex.1))
boxplot(residuals(anorex.1))