> guppies<- read.csv(file.choose(), header=TRUE)
> guppies.mod=aov(Offspring~Feedings+factor(Temperature), data=guppies)
> anova(guppies.mod)
Analysis of Variance Table

Response: Offspring
                    Df  Sum Sq Mean Sq F value    Pr(>F)
Feedings             1 1333.33 1333.33  55.843 3.573e-09 ***
factor(Temperature)  2 1059.73  529.87  22.192 2.943e-07 ***
Residuals           41  978.93   23.88
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> TukeyHSD(guppies.mod)
  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = Offspring ~ Feedings + factor(Temperature), data = guppies)

$`factor(Temperature)`
             diff       lwr       upr     p adj
75F-70F  5.333333 0.9946773  9.671989 0.0128112
80F-70F 11.866667 7.5280106 16.205323 0.0000002
80F-75F  6.533333 2.1946773 10.871989 0.0020105