TO-DO list September 2025 Check the code all-round to see where we can find bottlenecks in computation and make things quicker. Should we make ggplot the default for graphs (instead of base)?