# \donttest{
# Load the dataset
data(badger_energy)
# Basic exploration
head(badger_energy)
summary(badger_energy)
# Compare energy expenditure by disease status
boxplot(daily_energy ~ disease, data = badger_energy,
main = "Daily Energy Expenditure by TB Status",
ylab = "DEE (kJ/day)", xlab = "TB Status")
# Examine relationship between body mass and energy expenditure
plot(daily_energy ~ body_mass, data = badger_energy,
col = as.numeric(disease), pch = 16,
main = "Energy Expenditure vs Body Mass",
xlab = "Body Mass (kg)", ylab = "DEE (kJ/day)")
legend("topright", levels(badger_energy$disease),
col = 1:3, pch = 16)
# }
Run the code above in your browser using DataLab