> MangroveToads <- read.csv(file.choose(), header=TRUE)
> MangroveToads.mod=aov(BP~Diet+factor(Location), data=MangroveToads)
> anova(MangroveToads.mod)
Analysis of Variance Table

Response: BP
                 Df  Sum Sq Mean Sq F value    Pr(>F)
Diet              2 1049.20  524.60 24.1011 0.0004105 ***
factor(Location)  4  523.07  130.77  6.0077 0.0155690 *
Residuals         8  174.13   21.77
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> TukeyHSD(MangroveToads.mod)
  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = BP ~ Diet + factor(Location), data = MangroveToads)

$Diet
              diff       lwr      upr     p adj
Diet_B-Diet_A  4.0 -4.431478 12.43148 0.4065532
Diet_C-Diet_A 19.4 10.968522 27.83148 0.0004488
Diet_C-Diet_B 15.4  6.968522 23.83148 0.0020439

$`factor(Location)`
           diff        lwr       upr     p adj
2-1  -4.0000000 -17.160335  9.160335 0.8260367
3-1 -15.0000000 -28.160335 -1.839665 0.0261261
4-1 -11.0000000 -24.160335  2.160335 0.1085312
5-1  -0.6666667 -13.827002 12.493669 0.9997400
3-2 -11.0000000 -24.160335  2.160335 0.1085312
4-2  -7.0000000 -20.160335  6.160335 0.4163798
5-2   3.3333333  -9.827002 16.493669 0.8983862
4-3   4.0000000  -9.160335 17.160335 0.8260367
5-3  14.3333333   1.172998 27.493669 0.0329869
5-4  10.3333333  -2.827002 23.493669 0.1376592