# Supplementary Material ESM2 for:

# Behavioural traits modulate the use of heterospecific social information for nest site selection: experimental evidence from a wild bird population
# Jennifer Morinay*, Jukka T. Forsman, Marion Germain and Blandine Doligez

# Proceedings of the Royal Society B: Biological Sciences, 2020, doi: 10.1098/rspb.2020.0265

# * Correspondance: jennifer.morinay@gmail.com

#---------------------#
#### PRELIMINARIES ####
#---------------------#
library(lme4)


#------------------------#
#### DATA PREPARATION ####
#------------------------#
data <- read.table("TableS5_dataset_Morinay_procB.csv",header=T,sep=";",na.string="NA",dec=".")
data$copying.binary <- as.factor(data$copying.binary)

#-------------------------#
#### Models for MALES  ####
#-------------------------#

### Aggressiveness M ####
#-----------------------#

dat.m.a <- data[!is.na(data$aggressiveness.m),]

dat.m.a$tit.density <- scale(dat.m.a$tit.density)
dat.m.a$tit.clutch.brood.size <- scale(dat.m.a$tit.clutch.brood.size)
dat.m.a$aggressiveness.m <- scale(dat.m.a$aggressiveness.m)
dat.m.a$dev.symbol <- scale(dat.m.a$dev.symbol)

m.a.full <- glmer(copying.binary ~ 1 + age.m + aggressiveness.m + tit.clutch.brood.size  + tit.density + dev.symbol +
                         age.m:aggressiveness.m + 
                         tit.clutch.brood.size:aggressiveness.m + tit.density:aggressiveness.m + 
                         tit.clutch.brood.size:age.m + tit.density:age.m +
                         (1|forest.patch)+(1|year),
                       family = "binomial",data = dat.m.a)

### Boldness M ####
#-----------------#

dat.m.b <- data[!is.na(data$boldness.m),]

dat.m.b$tit.density <- scale(dat.m.b$tit.density)
dat.m.b$tit.clutch.brood.size <- scale(dat.m.b$tit.clutch.brood.size)
dat.m.b$boldness.m <- scale(dat.m.b$boldness.m)
dat.m.b$dev.symbol <- scale(dat.m.b$dev.symbol)

m.b.full <- glmer(copying.binary ~ 1 + age.m + boldness.m + tit.clutch.brood.size  + tit.density + dev.symbol +
                         age.m:boldness.m + 
                         tit.clutch.brood.size:boldness.m + tit.density:boldness.m + 
                         tit.clutch.brood.size:age.m + tit.density:age.m +
                         (1|forest.patch)+(1|year),
                       family = "binomial",data = dat.m.b)

### Neophobia M ####
#------------------#

dat.m.n <- data[!is.na(data$neophobia.m),]

dat.m.n$tit.density <- scale(dat.m.n$tit.density)
dat.m.n$tit.clutch.brood.size <- scale(dat.m.n$tit.clutch.brood.size)
dat.m.n$neophobia.m <- scale(dat.m.n$neophobia.m)
dat.m.n$dev.symbol <- scale(dat.m.n$dev.symbol)

m.n.full <- glmer(copying.binary ~ 1 + age.m + neophobia.m + tit.clutch.brood.size  + tit.density + dev.symbol +
                         age.m:neophobia.m + 
                         tit.clutch.brood.size:neophobia.m + tit.density:neophobia.m + 
                         tit.clutch.brood.size:age.m + tit.density:age.m +
                         (1|forest.patch)+(1|year),
                       family = "binomial",data = dat.m.n)


#---------------------------#
#### Models for FEMALES  ####
#---------------------------#


### Aggressiveness F ####
#-----------------------#

dat.f.a <- data[!is.na(data$aggressiveness.f) & !is.na(data$age.f),]

dat.f.a$tit.density <- scale(dat.f.a$tit.density)
dat.f.a$tit.clutch.brood.size <- scale(dat.f.a$tit.clutch.brood.size)
dat.f.a$aggressiveness.f <- scale(dat.f.a$aggressiveness.f)
dat.f.a$dev.symbol <- scale(dat.f.a$dev.symbol)

f.a.full <- glmer(copying.binary ~ 1 + age.f + aggressiveness.f + tit.clutch.brood.size  + tit.density + dev.symbol +
                         age.f:aggressiveness.f + 
                         tit.clutch.brood.size:aggressiveness.f + tit.density:aggressiveness.f + 
                         tit.clutch.brood.size:age.f + tit.density:age.f +
                         (1|forest.patch)+(1|year),
                       family = "binomial",data = dat.f.a)

### Boldness F ####
#-----------------#

dat.f.b <- data[!is.na(data$age.f) & !is.na(data$boldness.f),]

dat.f.b$tit.density <- scale(dat.f.b$tit.density)
dat.f.b$tit.clutch.brood.size <- scale(dat.f.b$tit.clutch.brood.size)
dat.f.b$boldness.f <- scale(dat.f.b$boldness.f)
dat.f.b$dev.symbol <- scale(dat.f.b$dev.symbol)

f.b.full <- glmer(copying.binary ~ 1 + age.f + boldness.f + tit.clutch.brood.size  + tit.density + dev.symbol +
                         age.f:boldness.f + 
                         tit.clutch.brood.size:boldness.f + tit.density:boldness.f + 
                         tit.clutch.brood.size:age.f + tit.density:age.f +
                         (1|forest.patch)+(1|year),
                       family = "binomial",data = dat.f.b)


### Neophobia F ####
#------------------#

dat.f.n <- data[!is.na(data$age.f) & !is.na(data$neophobia.f),]

dat.f.n$tit.density <- scale(dat.f.n$tit.density)
dat.f.n$tit.clutch.brood.size <- scale(dat.f.n$tit.clutch.brood.size)
dat.f.n$neophobia.f <- scale(dat.f.n$neophobia.f)
dat.f.n$dev.symbol <- scale(dat.f.n$dev.symbol)

f.n.full <- glmer(copying.binary ~ 1 + age.f + neophobia.f + tit.clutch.brood.size  + tit.density + dev.symbol +
                         age.f:neophobia.f + 
                         tit.clutch.brood.size:neophobia.f + tit.density:neophobia.f + 
                         tit.clutch.brood.size:age.f + tit.density:age.f +
                         (1|forest.patch)+(1|year),
                       family = "binomial",data = dat.f.n)


#----------------------------------#
#### Models for MALES + FEMALES ####
#----------------------------------#

## Aggressiveness
##----------------##

dat.fm.a <- data[!is.na(data$age.f) & !is.na(data$aggressiveness.f) &
                     !is.na(data$age.m) & !is.na(data$aggressiveness.m),]

dat.fm.a$tit.density <- scale(dat.fm.a$tit.density)
dat.fm.a$tit.clutch.brood.size <- scale(dat.fm.a$tit.clutch.brood.size)
dat.fm.a$aggressiveness.f <- scale(dat.fm.a$aggressiveness.f)
dat.fm.a$aggressiveness.m <- scale(dat.fm.a$aggressiveness.m)
dat.fm.a$dev.symbol <- scale(dat.fm.a$dev.symbol)

fm.a.full <- glmer(copying.binary ~ 1 + age.f + aggressiveness.f + 
                         age.m + aggressiveness.m +
                         aggressiveness.m:aggressiveness.f+
                         tit.clutch.brood.size  + tit.density + dev.symbol +
                         age.m:aggressiveness.m +
                         (1|forest.patch) + (1|year),
                       family = "binomial",data = dat.fm.a)

## Boldness
##---------##

dat.fm.b <- data[!is.na(data$age.f) & !is.na(data$boldness.f) &
                !is.na(data$age.m) & !is.na(data$boldness.m),]

dat.fm.b$tit.density <- scale(dat.fm.b$tit.density)
dat.fm.b$tit.clutch.brood.size <- scale(dat.fm.b$tit.clutch.brood.size)
dat.fm.b$boldness.f <- scale(dat.fm.b$boldness.f)
dat.fm.b$boldness.m <- scale(dat.fm.b$boldness.m)
dat.fm.b$dev.symbol <- scale(dat.fm.b$dev.symbol)

fm.b.full <- glmer(copying.binary ~ 1 + age.f + boldness.f + 
                         age.m + boldness.m +
                         boldness.m:boldness.f+
                         tit.clutch.brood.size  + tit.density + dev.symbol +
                         age.f:boldness.f +
                         age.m:boldness.m +
                         (1|forest.patch) + (1|year),
                       family = "binomial",data = dat.fm.b)

## Neophobia
##-----------##

dat.fm.n <- data[!is.na(data$age.f) & !is.na(data$neophobia.f) &
                !is.na(data$age.m) & !is.na(data$neophobia.m),]

dat.fm.n$tit.density <- scale(dat.fm.n$tit.density)
dat.fm.n$tit.clutch.brood.size <- scale(dat.fm.n$tit.clutch.brood.size)
dat.fm.n$neophobia.f <- scale(dat.fm.n$neophobia.f)
dat.fm.n$neophobia.m <- scale(dat.fm.n$neophobia.m)
dat.fm.n$dev.symbol <- scale(dat.fm.n$dev.symbol)

fm.n.full <- glmer(copying.binary ~ 1 + age.f + neophobia.f + 
                         age.m + neophobia.m +
                         neophobia.m:neophobia.f+
                         tit.clutch.brood.size  + tit.density + dev.symbol +
                         tit.clutch.brood.size:neophobia.f +
                         (1|forest.patch) + (1|year),
                       family = "binomial",data = dat.fm.n)