#Models with 2 dimensional interactions M2D_LST_DTDB <- gam(PreyWeight ~ te(Length, SolarTime,k=c(3,5), bs=c('tp','cc'))+s(TrawlDepth,k=3)+s(BottomDepth, k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_LST_DT <- gam(PreyWeight ~ te(Length, SolarTime,k=c(3,5), bs=c('tp','cc'))+s(BottomDepth , k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_LST_DB <- gam(PreyWeight ~ te(Length, SolarTime,k=c(3,5), bs=c('tp','cc'))+s(TrawlDepth,k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_LST <- gam(PreyWeight ~ te(Length, SolarTime,k=c(3,5), bs=c('tp','cc')), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_LDT_STDB <- gam(PreyWeight ~ te(Length, TrawlDepth,k=c(3,3), bs=c('tp','tp'))+s(SolarTime, bs ='cc', k=5)+s(BottomDepth , k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_LDT_ST <- gam(PreyWeight ~ te(Length, TrawlDepth,k=c(3,3), bs=c('tp','tp'))+s(SolarTime, bs ='cc', k=5), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_LDT_DB <- gam(PreyWeight ~ te(Length, TrawlDepth,k=c(3,3), bs=c('tp','tp'))+s(BottomDepth , k=3), data=dfGAMData, family = tw()) M2D_LDT <- gam(PreyWeight ~ te(Length, TrawlDepth,k=c(3,3), bs=c('tp','tp')), data=dfGAMData, family = tw()) M2D_STDT_LDB <- gam(PreyWeight ~ te(TrawlDepth, SolarTime,k=c(3,5), bs=c('tp','cc'))+s(Length,k=3)+s(BottomDepth , k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_STDT_L <- gam(PreyWeight ~ te(TrawlDepth, SolarTime,k=c(3,5), bs=c('tp','cc'))+s(Length,k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_STDT_DB <- gam(PreyWeight ~ te(TrawlDepth, SolarTime,k=c(3,5), bs=c('tp','cc'))+s(BottomDepth , k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_STDT <- gam(PreyWeight ~ te(TrawlDepth, SolarTime,k=c(3,5), bs=c('tp','cc')), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_STDB_LDT <- gam(PreyWeight ~ te(BottomDepth, SolarTime,k=c(3,5), bs=c('tp','cc'))+s(Length,k=3)+s(TrawlDepth, k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_STDB_L <- gam(PreyWeight ~ te(BottomDepth, SolarTime,k=c(3,5), bs=c('tp','cc'))+s(Length,k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_STDB_DT <- gam(PreyWeight ~ te(BottomDepth, SolarTime,k=c(3,5), bs=c('tp','cc'))+s(TrawlDepth, k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M2D_STDB <- gam(PreyWeight ~ te(BottomDepth, SolarTime,k=c(3,5), bs=c('tp','cc')), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) #Without interactions M_LSTDTDB <- gam(PreyWeight ~ s(Length, k=3) + s(SolarTime, bs ='cc', k=5) + s(TrawlDepth , k=3)+s(BottomDepth , k=3) , data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M_LSTDT <- gam(PreyWeight ~ s(Length, k=3) + s(SolarTime, bs ='cc', k=5) + s(TrawlDepth , k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M_LSTDB <- gam(PreyWeight ~ s(Length, k=3) + s(SolarTime, bs ='cc', k=5) +s(BottomDepth , k=3) , data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M_LDTDB <- gam(PreyWeight ~ s(Length, k=3) + s(TrawlDepth , k=3)+s(BottomDepth , k=3) , data=dfGAMData, family = tw()) M_STDTDB <- gam(PreyWeight ~ s(SolarTime, bs ='cc', k=5) + s(TrawlDepth , k=3)+s(BottomDepth , k=3) , data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M_LST <- gam(PreyWeight ~ s(Length, k=3) + s(SolarTime, bs ='cc', k=5) , data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M_LDT <- gam(PreyWeight ~ s(Length, k=3) + s(TrawlDepth , k=3), data=dfGAMData, family = tw()) M_LDB <- gam(PreyWeight ~ s(Length, k=3) +s(BottomDepth , k=3) , data=dfGAMData, family = tw()) M_STDT <- gam(PreyWeight ~ s(SolarTime, bs ='cc', k=5) + s(TrawlDepth , k=3), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M_STDB <- gam(PreyWeight ~ s(SolarTime, bs ='cc', k=5) +s(BottomDepth , k=3) , data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M_L <- gam(PreyWeight ~ s(Length, k=3), data=dfGAMData, family = tw()) M_ST <- gam(PreyWeight ~ s(SolarTime, bs ='cc', k=5), data=dfGAMData, family = tw(), knots=c(list(SolarTime=seq(0,24, length=5)))) M_DT <- gam(PreyWeight ~ s(TrawlDepth , k=3), data=dfGAMData, family = tw()) M_DB <- gam(PreyWeight ~ s(BottomDepth , k=3) , data=dfGAMData, family = tw())