emulating pax (catch at age calculation) in R