attach(S3_Recruit_1997.2008_R_ready)
str(S3_Recruitment_1997.2008_R_ready)

library(lme4)
library(MuMIn)

fm1<-lmer(logB_std~1+(1|Bay))
fm2<-lmer(logB_std~1+(1|Bay)+(1|Site))
fm3<-lmer(logB_std~1+(1|Year))
fm4<-lm(logB_std~Size)
fm5<-lmer(logB_std~1+(1|Bay)+(1|Year))
fm6<-lmer(logB_std~1+(1|Bay)+(1|Year)+(1|Bay:Year))
fm7<-lmer(logB_std~Size+(1|Bay))
fm8<-lmer(logB_std~Size+(1|Bay)+(1|Bay:Size))
fm9<-lmer(logB_std~Size+(1|Year))
fm10<-lmer(logB_std~Size+(1|Year)+(1|Year:Size))
fm11<-lmer(logB_std~Size+(1|Bay)+(1|Year))
fm12<-lmer(logB_std~Size+(1|Bay)+(1|Year)+(1|Bay:Year))
fm13<-lmer(logB_std~Size+(1|Bay)+(1|Year)+(1|Bay:Size))
fm14<-lmer(logB_std~Size+(1|Bay)+(1|Year)+(1|Year:Size))
fm15<-lmer(logB_std~Size+(1|Bay)+(1|Year)+(1|Bay:Year)+(1|Bay:Size))
fm16<-lmer(logB_std~Size+(1|Bay)+(1|Year)+(1|Bay:Year)+(1|Year:Size))
fm17<-lmer(logB_std~Size+(1|Bay)+(1|Year)+(1|Bay:Size)+(1|Year:Size))
fm18<-lmer(logB_std~1+(1|Bay)+(1|Site)+(1|Year))
fm19<-lmer(logB_std~1+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year))
fm20<-lmer(logB_std~1+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Year:Site))
fm21<-lmer(logB_std~Size+(1|Bay)+(1|Site))
fm22<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Bay:Size))
fm23<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Bay:Size)+(1|Site:Size))
fm24<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year))
fm25<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year))
fm26<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Size))
fm27<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Site:Size))
fm28<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Bay:Size))
fm29<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Bay:Size)+(1|Year:Site))
fm30<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Bay:Size)+(1|Site:Size))
fm31<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Year:Size))
fm32<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Year:Size)+(1|Year:Site))
fm33<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Year:Size)+(1|Site:Size))
fm34<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Size)+(1|Year:Size))
fm35<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Size)+(1|Year:Size)+(1|Site:Size))
fm36<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Bay:Size)+(1|Year:Size))
fm37<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Bay:Size)+(1|Year:Size)+(1|Year:Bay:Size))
fm38<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Bay:Size)+(1|Year:Size)+(1|Year:Site))
fm39<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Bay:Size)+(1|Year:Size)+(1|Site:Size))
fm40<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Bay:Size)+(1|Year:Size)+(1|Year:Bay:Size)+(1|Year:Site))
fm41<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Bay:Size)+(1|Year:Size)+(1|Year:Bay:Size)+(1|Site:Size))
fm42<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Bay:Size)+(1|Year:Size)+(1|Year:Site)+(1|Site:Size))
fm43<-lmer(logB_std~Size+(1|Bay)+(1|Site)+(1|Year)+(1|Bay:Year)+(1|Bay:Size)+(1|Year:Size)+(1|Year:Bay:Size)+(1|Year:Site)+(1|Site:Size))

modelsel <- model.sel(fm1,  fm2,  fm3,  fm4,	fm5,	fm6,	fm7,	fm8,	fm9,	fm10,	fm11,	fm12,
                     fm13,	fm14,	fm15,	fm16,	fm17,	fm18,	fm19,	fm20,	fm21,	fm22,	fm23,	fm24,
                     fm25,	fm26,	fm27,	fm28,	fm29,	fm30,	fm31,	fm32,	fm33,	fm34,	fm35,	fm36,
                     fm37,	fm38,	fm39,	fm40,	fm41,	fm42,	fm43)
modelsel