gdyang
一个数据集中有身高(h)、体重(w)、性别(s),按性别不同分别做身高相对体重的回归分析,这个怎么写呢,看了半天书也没有找到。
是lm(h~w|s)吗,还是用tapply?
yihui
lm的formula没有|这种写法。你可以用subset参数。
比如先建立一个lm对象:lm1=lm(w~h, data = somedata, subset = s=='M') #对Male
然后update这个模型,把subset改掉:update(lm1, subset = s=='F') #对Female
gdyang
谢谢老大,假如分类变量不是两个,而是很多,那不是很麻烦?有别的方法吗?
abel
?by
这个命令就可以了。
TTG
by(df,S,function(x) lm(W~H,data=x))->lm12
summary(lm12[1])
gdyang
谢谢楼上各位!
rtist
But this will not pool the error variances.