But I don't want the original mean y ,but the "corrected mean y" of analysis of 
covariant.
Also,the difference of original mean y is not 136.83,which is the difference of 
the "corrected meany".




 


ÔÚ 2013-03-09 20:11:11£¬"R. Michael Weylandt" <michael.weyla...@gmail.com> 
дµÀ£º



On Sat, Mar 9, 2013 at 10:35 AM, meng <laomen...@163.com> wrote:
Hi all:
My data is in the attachment.
I want to analysis the mean difference of y between 2 sex.

My code:
result_lm<-lm(y~factor(sex) + x1 + x2)
summary(result_lm)

The result of "factor(sex)m" 136.83, is the mean difference of y between 2 
sex,and the corresponding p value is 0.07618.

My question is: how to get the mean y of sex(m) and sex(f) respectively via lm 
function?



lm() would be overkill no? I'd just use tapply(), by(), ave(), or the like: 
e.g., 


tapply(dat$y, dat$sex, mean)


Cheers, 


MW 
        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to