Wrapper for computing MSY from a "fisheries model"

msy_wrapper(model, newdata, hratevec, tsave)

Arguments

model

A fitGP_fish model.

newdata

Data frame for iterated prediction, such as made with makelags_iter.

hratevec

Vector of harvest rates

tsave

Number of timepoints to average over. Should not exceed nrow(newdata)