library(sjstats) Area.Lp<-readRDS("models/AreaLp.rds") Area.Tv<-readRDS("models/AreaTv.rds") Prey.Lp<-readRDS("models/PreyLp.rds") Prey.Tv<-readRDS("models/PreyTv.rds") Prey_Ephem.Lp<-readRDS("models/PreyEphemLp.rds") Prey_Ephem.Tv<-readRDS("models/PreyEphemTv.rds") BCR.Lp<-readRDS("models/BCRLp.rds") BCR.Tv<-readRDS("models/BCRTv.rds") Trans.Lp<-readRDS("models/TransectLp.rds") Trans.Tv<-readRDS("models/TransectTv.rds") Table<-rbind( data.frame(tidy_stan(Area.Lp,prob=c(.9,.8)),Model="Capture area",species="L. patagiatus"), data.frame(tidy_stan(Area.Tv,prob=c(.9,.8)),Model="Capture area",species="T. versicolor"), data.frame(tidy_stan(Prey.Lp,prob=c(.9,.8)),Model="Prey (Diptera)",species="L. patagiatus"), data.frame(tidy_stan(Prey.Tv,prob=c(.9,.8)),Model="Prey (Diptera)",species="T. versicolor"), data.frame(tidy_stan(Prey_Ephem.Lp,prob=c(.9,.8)),Model="Prey (Ephemeroptera)",species="L. patagiatus"), data.frame(tidy_stan(Prey_Ephem.Tv,prob=c(.9,.8)),Model="Prey (Ephemeroptera)",species="T. versicolor"), data.frame(tidy_stan(BCR.Lp,prob=c(.9,.8)),Model="Body condition",species="L. patagiatus"), data.frame(tidy_stan(BCR.Tv,prob=c(.9,.8)),Model="Body condition",species="T. versicolor"), data.frame(tidy_stan(Trans.Lp,prob=c(.9,.8)),Model="Abundance",species="L. patagiatus"), data.frame(tidy_stan(Trans.Tv,prob=c(.9,.8)),Model="Abundance",species="T. versicolor") ) Table<-Table[,-c(1,3)] names(Table)<-c("Parameter","Estimate","SE", "80% CI","90% CI","PD","ESS", "Rhat","MCSE","Response","species") Table<-Table[,c(10,11,1:9)] unique(Table$Parameter) Table$Parameter[Table$Parameter=="dB.s"]<-"Sound pressure level" Table$Parameter[Table$Parameter=="SOUNDRIVER"]<-"River playback" Table$Parameter[Table$Parameter=="SOUNDSHIFTED"]<-"Shifted River playback" Table$Parameter[Table$Parameter=="YearYear2"]<-"Year: 2018" Table$Parameter[Table$Parameter=="YearYear3"]<-"Year: 2019" Table$Parameter[Table$Parameter=="TASS.s"]<-"Time after sunset" Table$Parameter[Table$Parameter=="OrdinalDay.s"]<-"Ordinal day" Table$Parameter[Table$Parameter=="Temp.s"]<-"Temperature" Table$Parameter[Table$Parameter=="WS.s"]<-"Wind speed" Table$Parameter[Table$Parameter=="Ele.s"]<-"Elevation" Table$Parameter[Table$Parameter=="reciprocal_dispersion"]<-"Reciprocal dispersion" Table$Parameter[Table$Parameter=="TreatRIVER"]<-"River playback" Table$Parameter[Table$Parameter=="TreatSHIFTED"]<-"Shifted River playback" Table$Parameter[Table$Parameter=="Stream.s"]<-"Stream width" Table$Parameter[Table$Parameter=="ele.s"]<-"Elevation" Table$Parameter[Table$Parameter=="YEARYEAR2"]<-"Year: 2018" Table$Parameter[Table$Parameter=="YEARYEAR3"]<-"Year: 2019" head(Table) write.csv(Table,"Model_Output_Table.csv")