bivariate_response(), store and restore graphical
parameters only when plotting area is modified.response_curves(), omit graphical parameter storing
and resetting to allow plotting in multipanel figures.projection_changes(), create data.frame
to set levels representing changes outside of the loop. This ensures it
works correctly even when by_gcm is set to
FALSE.