setwd("/home1/p309444/models/AIC")

library(mgcv)
library(chron)
library(lubridate)
library(lme4)
library(tidyverse)

physio<-read.csv("physio_28_5.csv",stringsAsFactors = F)

physio$Observer<-as.factor(physio$Observer)
physio$BirdID<-as.factor(physio$BirdID)
physio$SexEstimate<-as.factor(physio$SexEstimate)
physio$newstat<-as.factor(physio$newstat)


wing_g<-gam(WingLength_z~s(age_year,k=15)+s(newlifespan)+RightTarsus_z+newstat+
                SexEstimate+avg_invert+
                s(birthyear, bs="re")+s(Observer, bs="re")+s(BirdID, bs="re"),
              data=physio, method = "REML", family = gaussian)

wing_gs<-gam(WingLength_z~s(age_year,k=15)+s(newlifespan)+RightTarsus_z+newstat+
                   s(age_year, SexEstimate, k=15, bs='fs', m=2)+avg_invert+
                   s(birthyear, bs="re")+s(Observer, bs="re")+s(BirdID, bs="re"),
                 data=physio, method = "REML", family = gaussian)

wing_s<-gam(WingLength_z~s(age_year, SexEstimate, k=15, bs='fs')+s(newlifespan)+RightTarsus_z+newstat+
                   avg_invert+s(birthyear, bs="re")+s(Observer, bs="re")+s(BirdID, bs="re"),
                 data=physio, method = "REML", family = gaussian)


saveRDS(wing_g, "wing_g.rds")
saveRDS(wing_gs, "wing_gs.rds")
saveRDS(wing_s, 'wing_s.rds')



#######