bombus_subset <- bombus[bombus$colony %in% c("17", "104", "20", "24"), ]
results <- bumbl(bombus_subset, colonyID = colony, t = week,
formula = mass ~ week)
library(ggplot2)
autoplot(results)
Run the code above in your browser using DataLab