Please recreate Figure 9.7 from your text. Make the figures look tidy, and include a vector field and nullclines. Make up your own parameter values but also look for some guidance within the book if you don’t know where to start.
If you want to create a multi-panel plot, then use the function
par()
and set the argument mfrow = c(n, m)
with n
being the number of rows and m
being
the number of columns (e.g., par(mfrow = c(2, 3))
creates a
plot with 2 panels for rows and 3 columns).
There is often a lot of white space as a default around
plot()
. You can control that with
mar = c(a, b, c, d)
) as an argument in par()
.
a
—d
respectfully adjust the whitespace at the
bottom, left, top, and right of each panel.
It might be useful to make the graphs that go in the first column first (they use linear interaction terms) then the second column (they have the nonlinear interaction terms).