# \donttest{
# Load the dataset
data(vulture_diet)
head(vulture_diet)
summary(vulture_diet)
# Examine isotope signatures by species
boxplot(d13C ~ species, data = vulture_diet,
main = "Carbon Isotope Signatures by Vulture Species",
xlab = "Species", ylab = "d13C (per mill)")
# Create isotope biplot
plot(vulture_diet$d13C, vulture_diet$d15N,
col = as.numeric(vulture_diet$species),
pch = 16, cex = 1.2,
xlab = "d13C (per mill)", ylab = "d15N (per mill)",
main = "Vulture Feather Isotope Signatures")
legend("topright", legend = levels(vulture_diet$species),
col = 1:nlevels(vulture_diet$species), pch = 16)
# Summary statistics by species
aggregate(. ~ species, data = vulture_diet[,1:4], FUN = mean)
# }
Run the code above in your browser using DataLab