Dear R user:
   
  I have a question about doing ANCOVA in S-plus or R. 
   
  I know that many users use lm to do the regression and check the ANCOVA. But is there a way to get the traditional Table form of the ANCOVA test through S-plus (like what we would get from SPSS or SAS)?
   
  The problem I’m interested in is whether or not there is a treatment effect on some medical measurement. I will have the endpoint (E) and the baseline (B) of the measurement,   together with a dummy variable (T) of what treatment the patients get (control, or new). I’m thinking of using lm( (E-B)  ~ B + T, data = ..)  
   
  Question 1: if I get the regression results (E-B) = b0 + b1*B + b2*T and in order to check whether or not there is a treatment effect, all I need to do is to check whether the p-value associated with b2 is less than the significant level I set or there’s more I need to do?
   
  Question 2: If I put  anova( lm( (E-B)  ~ B + T, data = ..) ), is the table I get here the one I want for an ANCOVA? My concern is that ANCOVA is supposed to work only on residual compared to ANOVA, but we are using the same command here. If this gives me the ANCOVA table, does the order of the variables in regression matter? Will 
  anova( lm( (E-B)  ~  T + B, data = ..) ) gives the same result?
   
  Question 3: if I have another dummy variable to indicate the location of hospital each patient goes to (L)  and want to include it to my model, should I do
   lm( (E-B)  ~ B + T + L, data = ..)   or  I need to consider the interactions?
   
  Thanks a lot!
   
  Cindy

 __________________________________________________



	[[alternative HTML version deleted]]

