本文通过加载数据并使用线性模型，对比了不同对比编码方法（如Treatment、Helmert及Sum对比）对性别和年龄因素在体重变化研究中的影响。
Ancovacontrasts <- read.table("c:\\temp\\Ancovacontrasts.txt",header=T)
attach(Ancovacontrasts)
attach(Ancovacontrasts)
names(Ancovacontrasts)
lm(weight[sex=="male"]~age[sex=="male"])
lm(weight~age,subset=(sex=="female"))
lm(weight~age)
options(contrasts=c("contr.treatment", "contr.poly"))
model1 <- lm(weight~age*sex)
summary(model1)
options(contrasts=c("contr.helmert", "contr.poly"))
model2 <- lm(weight~age*sex)
summary(model2)
options(contrasts=c("contr.sum", "contr.poly"))
model3 <- lm(weight~age*sex)
summary(model3)
扫一扫
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
