my_data <- read.csv("C:\\Sire half-sibs.csv", header=TRUE) # Add in the appropriate path to this CSV file for your computer

results <- subset(my_data, select=c(2,3))

# Compute the analysis of variance
res.aov <- aov(Weight ~ Sire, data = results)
# Summary of the analysis
summary(res.aov)