Open sean-rohan-NOAA opened 3 months ago
@MargaretSiple-NOAA Is there a GAM for estimating spread? If not, how are missing spread values filled in the GOA/AI?
From our DPR:
Is there some code/query for retrieving or calculating vessel speed and total catch? Aiming to ensure the queries and data sets used in the models match.
Developed some missing spread/height fill code with GOA and AI models on dev
Still need to verify that EBS/NBS outputs match the original code since there were quite a few changes. Also need to check that it's OK for GOA/AI models to be fit to data from a single vessel.
Issue
GOA/AI models:
with.spread.model.gam <- gam(net.height ~ factor(vessel) + s(net.spread) + s(depth) + s(speed) + s(scope.ratio)
s(total.catch), data = with.spread.model.data)
without.spread.model.gam <- gam(net.height ~ factor(vessel) + s(depth) + s(speed) + s(scope.ratio) + s(total.catch), data = without.spread.model.data)
without.height.with.spread.data.gam <- gam(net.height ~ factor(vessel) + s(net.spread) + s(depth) + s(speed) + s(scope.ratio)
without.height.without.spread.data.gam <- gam(net.height ~ factor(vessel) + s(depth) + s(speed) + s(scope.ratio)