. anova testm2 c.testm1 i.hptreat c.cortm1 c.cortm2 i.female if inelig == 0 & an > youtv1 == 0

Number of obs =39    R-squared =0.7048
Root MSE= 16.0144    Adj R-squared =0.6601

 Source | Partial SS         df         MS        F    Prob>F
---------------------------------------------------------------
  Model |  20209.281          5   4041.8563     15.76  0.0000
 testm1 |  3516.6527          1   3516.6527     13.71  0.0008
 hptreat|  1183.5007          1   1183.5007      4.61  0.0391
 cortm1 |  8.5753841          1   8.5753841      0.03  0.8560
 cortm2 |  2810.9353          1   2810.9353     10.96  0.0023
 female |  2557.3444          1   2557.3444      9.97  0.0034
Residual|  8463.2532         33   256.46222
----------------------------------------------------------------
  Total |  28672.535         38   754.54038

FosseTest<-aov(testm2~testm1+hptreat+cortm1+cortm2+female,data=X2data)
summary(FosseTest)

            Df Sum Sq Mean Sq F value   Pr(>F)
testm1       1  15121   15121  58.962 7.68e-09 ***
hptreat      1    524     524   2.043  0.16228
cortm1       1     23      23   0.089  0.76715
cortm2       1   1984    1984   7.735  0.00888 **
female       1   2557    2557   9.972  0.00339 **
Residuals   33   8463     256
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1