#install packages
install.packages("data.table")
install.packages("BayesFactor")


#load packages
library(data.table)
library(BayesFactor)


########## data preparation (if not already done with "SNP-B5_R-Script_MainText.R") ##########


#read data
dat <- as.data.frame(fread("https://madata.bib.uni-mannheim.de/364/2/SNP-B5_Data.csv", header = T, sep = ","))


#apply data-exclusion criteria (Footnote 7)
dat <- subset(dat, subset = bigfive_incomplete == 0)
dat <- subset(dat, subset = sample < 5 | c(sample == 5 & external_predictors_incomplete == 0) | c(sample == 6 & external_predictors_incomplete == 0))
dat <- subset(dat, subset = double_participation == 0)
dat <- subset(dat, subset = instruction_check_item_part2 != 2)
dat <- subset(dat, subset = not_followed_instructions == 0)
dat <- subset(dat, subset = c(chinese_characters_known != 1 & recall_correct_total > 0) | c(chinese_characters_known == 1 & recall_correct_values > 0))
dat <- subset(dat, subset = participated_seriously != 2)
dat <- subset(dat, subset = not_use_data != 1)


#compute variables

##compute BFI/BFI-2 domains
dat_BFI <- subset(dat, subset = sample == 1 | sample == 2)
dat_BFI2 <- subset(dat, subset = sample == 3 | sample == 4 | sample == 5 | sample == 6)
dat_BFI$agr <- rowMeans(dat_BFI[c("agr01r","agr02","agr03r","agr04","agr05","agr06r","agr07","agr08r","agr09")], na.rm=TRUE)
dat_BFI$cns <- rowMeans(dat_BFI[c("cns01","cns02r","cns03","cns04r","cns05r","cns06","cns07","cns08","cns09r")], na.rm=TRUE)
dat_BFI$opn <- rowMeans(dat_BFI[c("opn01","opn02","opn03","opn04","opn05","opn06","opn07r","opn08","opn09r","opn10")], na.rm=TRUE)
dat_BFI$ext <- rowMeans(dat_BFI[c("ext01","ext02r","ext03","ext04","ext05r","ext06","ext07r","ext08")], na.rm=TRUE)
dat_BFI$neu <- rowMeans(dat_BFI[c("neu01","neu02r","neu03","neu04","neu05r","neu06","neu07r","neu08")], na.rm=TRUE)
dat_BFI2$agr <- rowMeans(dat_BFI2[c("agr01","agr02","agr03r","agr04r","agr05r","agr06","agr07","agr08r","agr09r","agr10r","agr11","agr12")], na.rm=TRUE)
dat_BFI2$cns <- rowMeans(dat_BFI2[c("cns01r","cns02r","cns03","cns04","cns05r","cns06r","cns07","cns08","cns09","cns10r","cns11","cns12r")], na.rm=TRUE)
dat_BFI2$opn <- rowMeans(dat_BFI2[c("opn01r","opn02","opn03","opn04","opn05r","opn06r","opn07","opn08","opn09r","opn10r","opn11r","opn12")], na.rm=TRUE)
dat_BFI2$ext <- rowMeans(dat_BFI2[c("ext01","ext02","ext03r","ext04r","ext05","ext06r","ext07r","ext08r","ext09","ext10","ext11r","ext12")], na.rm=TRUE)
dat_BFI2$neu <- rowMeans(dat_BFI2[c("neu01r","neu02r","neu03","neu04","neu05r","neu06r","neu07","neu08","neu09r","neu10r","neu11","neu12")], na.rm=TRUE)

##compute BFI-2 facets
dat_BFI2$compassion <- rowMeans(dat_BFI2[c("agr01","agr04r","agr07","agr10r")], na.rm=TRUE)
dat_BFI2$respectfulness <- rowMeans(dat_BFI2[c("agr02","agr05r","agr08r","agr11")], na.rm=TRUE)
dat_BFI2$trust <- rowMeans(dat_BFI2[c("agr03r","agr06","agr09r","agr12")], na.rm=TRUE)
dat_BFI2$organization <- rowMeans(dat_BFI2[c("cns01r","cns04","cns07","cns10r")], na.rm=TRUE)
dat_BFI2$productiveness <- rowMeans(dat_BFI2[c("cns02r","cns05r","cns08","cns11")], na.rm=TRUE)
dat_BFI2$responsibility <- rowMeans(dat_BFI2[c("cns03","cns06r","cns09","cns12r")], na.rm=TRUE)
dat_BFI2$aesthetic <- rowMeans(dat_BFI2[c("opn01r","opn04","opn07","opn10r")], na.rm=TRUE)
dat_BFI2$creative <- rowMeans(dat_BFI2[c("opn03","opn06r","opn09r","opn12")], na.rm=TRUE)
dat_BFI2$intellectual <- rowMeans(dat_BFI2[c("opn02","opn05r","opn08","opn11r")], na.rm=TRUE)
dat_BFI2$assertiveness <- rowMeans(dat_BFI2[c("ext02","ext05","ext08r","ext11r")], na.rm=TRUE)
dat_BFI2$energy <- rowMeans(dat_BFI2[c("ext03r","ext06r","ext09","ext12")], na.rm=TRUE)
dat_BFI2$sociability <- rowMeans(dat_BFI2[c("ext01","ext04r","ext07r","ext10")], na.rm=TRUE)
dat_BFI2$anxiety <- rowMeans(dat_BFI2[c("neu01r","neu04","neu07","neu10r")], na.rm=TRUE)
dat_BFI2$depression <- rowMeans(dat_BFI2[c("neu02r","neu05r","neu08","neu11")], na.rm=TRUE)
dat_BFI2$emotional <- rowMeans(dat_BFI2[c("neu03","neu06r","neu09r","neu12")], na.rm=TRUE)

##merge datasets
dat_BFI[ ,c("compassion","respectfulness","trust","organization","productiveness","responsibility","aesthetic","creative","intellectual",
            "assertiveness","energy","sociability","anxiety","depression","emotional")] <- NA
dat <- rbind(dat_BFI,dat_BFI2)

##compute external predictors in Experiment 3
dat3 <- subset(dat, subset = sample == 5 | sample == 6)
dat3$rational <- rowMeans(dat3[c("premed01","premed02","premed03","premed04","premed05","premed06","premed07","premed08","premed09","premed10","premed11")], na.rm=TRUE)
dat3$nfc <- rowMeans(dat3[c("nfc01","nfc02","nfc03r","nfc04r","nfc05r","nfc06","nfc07r","nfc08r","nfc09r","nfc10","nfc11","nfc12r","nfc13","nfc14","nfc15","nfc16r","nfc17r","nfc18")], na.rm=TRUE)
dat3$attention <- rowMeans(dat3[c("sa01","sa02","sa03","sa04","sa05","sa06","sa07","sa08","sa09","sa10","sa11","sa12r","sa13r")], na.rm=TRUE)

##merge datasets
dat3 <- subset(dat3, select = c(id,pair_index_within_sample,rational,nfc,attention))
dat <- merge(dat, dat3, by = c("id","pair_index_within_sample"), all = TRUE)
rm(dat_BFI,dat_BFI2,dat3)


#create datasets for experiments

##Experiment 1
dat1 <- subset(dat, subset = sample == 1 | sample == 2)
##Experiment 2
dat2 <- subset(dat, subset = sample == 3 | sample == 4)
##Experiment 3
dat3 <- subset(dat, subset = sample == 5 | sample == 6)


#grand-mean centering of BFI/BFI-2 domains (level-2 predictors)

##Experiment 1
dat1$agr_c <- dat1$agr - mean(dat1$agr)
dat1$cns_c <- dat1$cns - mean(dat1$cns)
dat1$opn_c <- dat1$opn - mean(dat1$opn)
dat1$ext_c <- dat1$ext - mean(dat1$ext)
dat1$neu_c <- dat1$neu - mean(dat1$neu)
##Experiment 2
dat2$agr_c <- dat2$agr - mean(dat2$agr)
dat2$cns_c <- dat2$cns - mean(dat2$cns)
dat2$opn_c <- dat2$opn - mean(dat2$opn)
dat2$ext_c <- dat2$ext - mean(dat2$ext)
dat2$neu_c <- dat2$neu - mean(dat2$neu)
##Experiment 3
dat3$agr_c <- dat3$agr - mean(dat3$agr)
dat3$cns_c <- dat3$cns - mean(dat3$cns)
dat3$opn_c <- dat3$opn - mean(dat3$opn)
dat3$ext_c <- dat3$ext - mean(dat3$ext)
dat3$neu_c <- dat3$neu - mean(dat3$neu)


##grand-mean centering of BFI-2 facets (level-2 predictors)

##Experiment 2
dat2$compassion_c <- dat2$compassion - mean(dat2$compassion)
dat2$respectfulness_c <- dat2$respectfulness - mean(dat2$respectfulness)
dat2$trust_c <- dat2$trust - mean(dat2$trust)
dat2$organization_c <- dat2$organization - mean(dat2$organization)
dat2$productiveness_c <- dat2$productiveness - mean(dat2$productiveness)
dat2$responsibility_c <- dat2$responsibility - mean(dat2$responsibility)
dat2$aesthetic_c <- dat2$aesthetic - mean(dat2$aesthetic)
dat2$creative_c <- dat2$creative - mean(dat2$creative)
dat2$intellectual_c <- dat2$intellectual - mean(dat2$intellectual)
dat2$assertiveness_c <- dat2$assertiveness - mean(dat2$assertiveness)
dat2$energy_c <- dat2$energy - mean(dat2$energy)
dat2$sociability_c <- dat2$sociability - mean(dat2$sociability)
dat2$anxiety_c <- dat2$anxiety - mean(dat2$anxiety)
dat2$depression_c <- dat2$depression - mean(dat2$depression)
dat2$emotional_c <- dat2$emotional - mean(dat2$emotional)
##Experiment 3
dat3$compassion_c <- dat3$compassion - mean(dat3$compassion)
dat3$respectfulness_c <- dat3$respectfulness - mean(dat3$respectfulness)
dat3$trust_c <- dat3$trust - mean(dat3$trust)
dat3$organization_c <- dat3$organization - mean(dat3$organization)
dat3$productiveness_c <- dat3$productiveness - mean(dat3$productiveness)
dat3$responsibility_c <- dat3$responsibility - mean(dat3$responsibility)
dat3$aesthetic_c <- dat3$aesthetic - mean(dat3$aesthetic)
dat3$creative_c <- dat3$creative - mean(dat3$creative)
dat3$intellectual_c <- dat3$intellectual - mean(dat3$intellectual)
dat3$assertiveness_c <- dat3$assertiveness - mean(dat3$assertiveness)
dat3$energy_c <- dat3$energy - mean(dat3$energy)
dat3$sociability_c <- dat3$sociability - mean(dat3$sociability)
dat3$anxiety_c <- dat3$anxiety - mean(dat3$anxiety)
dat3$depression_c <- dat3$depression - mean(dat3$depression)
dat3$emotional_c <- dat3$emotional - mean(dat3$emotional)


#grand-mean centering of external predictors in Experiment 3 (level-2 predictors)
dat3$rational_c <- dat3$rational - mean(dat3$rational)
dat3$nfc_c <- dat3$nfc - mean(dat3$nfc)
dat3$attention_c <- dat3$attention - mean(dat3$attention)


#exclude pairs of Chinese characters and social values for which the sociocultural norm was not recalled correctly

##Experiment 1
dat1 <- subset(dat1, subset = recall_correct == 1)
##Experiment 2
dat2 <- subset(dat2, subset = recall_correct == 1)
##Experiment 3
dat3 <- subset(dat3, subset = recall_correct == 1)


#exclude pairs of Chinese characters for participants who knew the meaning of at least one of the characters 

##Experiment 1
dat1 <- subset(dat1, subset = chinese_characters_known != 1 | c(chinese_characters_known == 1 & pair_index_across_samples > 18))
##Experiment 2 
dat2 <- subset(dat2, subset = chinese_characters_known != 1 | c(chinese_characters_known == 1 & pair_index_across_samples > 18))
##Experiment 3
dat3 <- subset(dat3, subset = chinese_characters_known != 1 | c(chinese_characters_known == 1 & pair_index_across_samples > 18))


##############################################################################


#conduct mixed-effects models with the BayesFactor package (note: running time of a model can be up to several days)

##convert id and sociocultural_norm to factor

###Experiment 1
dat1$id <-as.factor(dat1$id)
dat1$sociocultural_norm <-as.factor(dat1$sociocultural_norm)
###Experiment 2
dat2$id <-as.factor(dat2$id)
dat2$sociocultural_norm <-as.factor(dat2$sociocultural_norm)
###Experiment 3
dat3$id <-as.factor(dat3$id)
dat3$sociocultural_norm <-as.factor(dat3$sociocultural_norm)

##set the seed for reproducibility
set.seed(123)

##conduct domains-as-predictors models 

###Experiment 1
exp1_main_effects_model = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c + opn_c
                               + ext_c + neu_c + id + sociocultural_norm:id, 
                               data = dat1, whichRandom = c('id','sociocultural_norm:id'))

exp1_model_no_agr = lmBF(personal_preference ~ sociocultural_norm + cns_c + opn_c
                         + ext_c + neu_c + id + sociocultural_norm:id, 
                         data = dat1, whichRandom = c('id','sociocultural_norm:id'))

exp1_model_no_cns = lmBF(personal_preference ~ sociocultural_norm + agr_c + opn_c
                         + ext_c + neu_c + id + sociocultural_norm:id, 
                         data = dat1, whichRandom = c('id','sociocultural_norm:id'))

exp1_model_no_opn = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c
                         + ext_c + neu_c + id + sociocultural_norm:id, 
                         data = dat1, whichRandom = c('id','sociocultural_norm:id'))

exp1_model_no_ext = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c + opn_c
                         + neu_c + id + sociocultural_norm:id, 
                         data = dat1, whichRandom = c('id','sociocultural_norm:id'))

exp1_model_no_neu = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c + opn_c
                         + ext_c + id + sociocultural_norm:id, 
                         data = dat1, whichRandom = c('id','sociocultural_norm:id'))

exp1_model_no_norms = lmBF(personal_preference ~ agr_c + cns_c + opn_c
                           + ext_c + neu_c + id, 
                           data = dat1, whichRandom = c('id'))

exp1_full_model = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                       + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                       data = dat1, whichRandom = c('id','sociocultural_norm:id'))

exp1_model_no_IEagr = lmBF(personal_preference ~ sociocultural_norm + agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                           + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat1, whichRandom = c('id','sociocultural_norm:id'))

exp1_model_no_IEcns = lmBF(personal_preference ~ sociocultural_norm * agr_c + cns_c + sociocultural_norm * opn_c
                           + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat1, whichRandom = c('id','sociocultural_norm:id'))

exp1_model_no_IEopn = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + opn_c
                           + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat1, whichRandom = c('id','sociocultural_norm:id'))

exp1_model_no_IEext = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                           + ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat1, whichRandom = c('id','sociocultural_norm:id'))

exp1_model_no_IEneu = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                           + sociocultural_norm * ext_c + neu_c + id + sociocultural_norm:id, 
                           data = dat1, whichRandom = c('id','sociocultural_norm:id'))

###Experiment 2
exp2_main_effects_model = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c + opn_c
                               + ext_c + neu_c + id + sociocultural_norm:id, 
                               data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_agr = lmBF(personal_preference ~ sociocultural_norm + cns_c + opn_c
                         + ext_c + neu_c + id + sociocultural_norm:id, 
                         data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_cns = lmBF(personal_preference ~ sociocultural_norm + agr_c + opn_c
                         + ext_c + neu_c + id + sociocultural_norm:id, 
                         data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_opn = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c
                         + ext_c + neu_c + id + sociocultural_norm:id, 
                         data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_ext = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c + opn_c
                         + neu_c + id + sociocultural_norm:id, 
                         data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_neu = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c + opn_c
                         + ext_c + id + sociocultural_norm:id, 
                         data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_norms = lmBF(personal_preference ~ agr_c + cns_c + opn_c
                           + ext_c + neu_c + id, 
                           data = dat2, whichRandom = c('id'))

exp2_full_model = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                       + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                       data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEagr = lmBF(personal_preference ~ sociocultural_norm + agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                           + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEcns = lmBF(personal_preference ~ sociocultural_norm * agr_c + cns_c + sociocultural_norm * opn_c
                           + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEopn = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + opn_c
                           + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEext = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                           + ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEneu = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                           + sociocultural_norm * ext_c + neu_c + id + sociocultural_norm:id, 
                           data = dat2, whichRandom = c('id','sociocultural_norm:id'))

###Experiment 3
exp3_main_effects_model = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c + opn_c
                               + ext_c + neu_c + id + sociocultural_norm:id, 
                               data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_agr = lmBF(personal_preference ~ sociocultural_norm + cns_c + opn_c
                         + ext_c + neu_c + id + sociocultural_norm:id, 
                         data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_cns = lmBF(personal_preference ~ sociocultural_norm + agr_c + opn_c
                         + ext_c + neu_c + id + sociocultural_norm:id, 
                         data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_opn = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c
                         + ext_c + neu_c + id + sociocultural_norm:id, 
                         data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_ext = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c + opn_c
                         + neu_c + id + sociocultural_norm:id, 
                         data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_neu = lmBF(personal_preference ~ sociocultural_norm + agr_c + cns_c + opn_c
                         + ext_c + id + sociocultural_norm:id, 
                         data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_norms = lmBF(personal_preference ~ agr_c + cns_c + opn_c
                           + ext_c + neu_c + id, 
                           data = dat3, whichRandom = c('id'))

exp3_full_model = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                       + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                       data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEagr = lmBF(personal_preference ~ sociocultural_norm + agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                           + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEcns = lmBF(personal_preference ~ sociocultural_norm * agr_c + cns_c + sociocultural_norm * opn_c
                           + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEopn = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + opn_c
                           + sociocultural_norm * ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEext = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                           + ext_c + sociocultural_norm * neu_c + id + sociocultural_norm:id, 
                           data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEneu = lmBF(personal_preference ~ sociocultural_norm * agr_c + sociocultural_norm * cns_c + sociocultural_norm * opn_c
                           + sociocultural_norm * ext_c + neu_c + id + sociocultural_norm:id, 
                           data = dat3, whichRandom = c('id','sociocultural_norm:id'))


##conduct facets-as-predictors models

###Experiment 2
exp2_main_effects_model_facets = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                                      + organization_c + productiveness_c + responsibility_c 
                                      + aesthetic_c + creative_c + intellectual_c 
                                      + assertiveness_c + energy_c + sociability_c 
                                      + anxiety_c + depression_c + emotional_c
                                      + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_comp = lmBF(personal_preference ~ sociocultural_norm + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_respect = lmBF(personal_preference ~ sociocultural_norm + compassion_c + trust_c 
                             + organization_c + productiveness_c + responsibility_c 
                             + aesthetic_c + creative_c + intellectual_c 
                             + assertiveness_c + energy_c + sociability_c 
                             + anxiety_c + depression_c + emotional_c 
                             + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_trust = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c 
                           + organization_c + productiveness_c + responsibility_c 
                           + aesthetic_c + creative_c + intellectual_c 
                           + assertiveness_c + energy_c + sociability_c 
                           + anxiety_c + depression_c + emotional_c 
                           + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_orga = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_prod = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_respo = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                            + organization_c + productiveness_c 
                            + aesthetic_c + creative_c + intellectual_c 
                            + assertiveness_c + energy_c + sociability_c 
                            + anxiety_c + depression_c + emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_aest = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + creative_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_crea = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_intel = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                           + organization_c + productiveness_c + responsibility_c 
                           + aesthetic_c + creative_c 
                           + assertiveness_c + energy_c + sociability_c 
                           + anxiety_c + depression_c + emotional_c 
                           + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_asse = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_ener = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c
                          + assertiveness_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_soci = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + assertiveness_c + energy_c 
                          + anxiety_c + depression_c + emotional_c
                          + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_anxi = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c
                          + assertiveness_c + energy_c + sociability_c 
                          + depression_c + emotional_c
                          + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_depr = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c
                          + anxiety_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_emot = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c 
                          + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_norms_facets = lmBF(personal_preference ~ compassion_c + respectfulness_c + trust_c 
                                  + organization_c + productiveness_c + responsibility_c 
                                  + aesthetic_c + creative_c + intellectual_c 
                                  + assertiveness_c + energy_c + sociability_c 
                                  + anxiety_c + depression_c + emotional_c 
                                  + id, data = dat2, whichRandom = c('id'))

exp2_full_model_facets = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                              + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                              + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                              + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                              + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                              + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEcomp = lmBF(personal_preference ~ sociocultural_norm + compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IErespect = lmBF(personal_preference ~ sociocultural_norm * compassion_c + respectfulness_c + sociocultural_norm * trust_c 
                               + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                               + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                               + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                               + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                               + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEtrust = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + trust_c 
                             + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                             + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                             + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                             + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                             + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEorga = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEprod = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IErespo = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                              + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + responsibility_c 
                              + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                              + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                              + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                              + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEaest = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEcrea = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEintel = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                             + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                             + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + intellectual_c 
                             + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                             + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                             + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEasse = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEener = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEsoci = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEanxi = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEdepr = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

exp2_model_no_IEemot = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + emotional_c 
                            + id + sociocultural_norm:id, data = dat2, whichRandom = c('id','sociocultural_norm:id'))

###Experiment 3
exp3_main_effects_model_facets = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                                      + organization_c + productiveness_c + responsibility_c 
                                      + aesthetic_c + creative_c + intellectual_c 
                                      + assertiveness_c + energy_c + sociability_c 
                                      + anxiety_c + depression_c + emotional_c
                                      + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_comp = lmBF(personal_preference ~ sociocultural_norm + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_respect = lmBF(personal_preference ~ sociocultural_norm + compassion_c + trust_c 
                             + organization_c + productiveness_c + responsibility_c 
                             + aesthetic_c + creative_c + intellectual_c 
                             + assertiveness_c + energy_c + sociability_c 
                             + anxiety_c + depression_c + emotional_c 
                             + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_trust = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c 
                           + organization_c + productiveness_c + responsibility_c 
                           + aesthetic_c + creative_c + intellectual_c 
                           + assertiveness_c + energy_c + sociability_c 
                           + anxiety_c + depression_c + emotional_c 
                           + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_orga = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_prod = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_respo = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                            + organization_c + productiveness_c 
                            + aesthetic_c + creative_c + intellectual_c 
                            + assertiveness_c + energy_c + sociability_c 
                            + anxiety_c + depression_c + emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_aest = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + creative_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_crea = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_intel = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                           + organization_c + productiveness_c + responsibility_c 
                           + aesthetic_c + creative_c 
                           + assertiveness_c + energy_c + sociability_c 
                           + anxiety_c + depression_c + emotional_c 
                           + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_asse = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + energy_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_ener = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c
                          + assertiveness_c + sociability_c 
                          + anxiety_c + depression_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_soci = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + assertiveness_c + energy_c 
                          + anxiety_c + depression_c + emotional_c
                          + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_anxi = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c
                          + assertiveness_c + energy_c + sociability_c 
                          + depression_c + emotional_c
                          + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_depr = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c 
                          + assertiveness_c + energy_c + sociability_c
                          + anxiety_c + emotional_c 
                          + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_emot = lmBF(personal_preference ~ sociocultural_norm + compassion_c + respectfulness_c + trust_c 
                          + organization_c + productiveness_c + responsibility_c 
                          + aesthetic_c + creative_c + intellectual_c
                          + assertiveness_c + energy_c + sociability_c 
                          + anxiety_c + depression_c 
                          + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_norms_facets = lmBF(personal_preference ~ compassion_c + respectfulness_c + trust_c 
                                  + organization_c + productiveness_c + responsibility_c 
                                  + aesthetic_c + creative_c + intellectual_c 
                                  + assertiveness_c + energy_c + sociability_c 
                                  + anxiety_c + depression_c + emotional_c 
                                  + id, data = dat3, whichRandom = c('id'))

exp3_full_model_facets = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                              + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                              + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                              + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                              + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                              + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEcomp = lmBF(personal_preference ~ sociocultural_norm + compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IErespect = lmBF(personal_preference ~ sociocultural_norm * compassion_c + respectfulness_c + sociocultural_norm * trust_c 
                               + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                               + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                               + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                               + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                               + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEtrust = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + trust_c 
                             + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                             + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                             + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                             + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                             + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEorga = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEprod = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IErespo = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                              + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + responsibility_c 
                              + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                              + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                              + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                              + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEaest = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEcrea = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEintel = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                             + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                             + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + intellectual_c 
                             + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                             + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                             + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEasse = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEener = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEsoci = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEanxi = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + anxiety_c + sociocultural_norm * depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEdepr = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + depression_c + sociocultural_norm * emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEemot = lmBF(personal_preference ~ sociocultural_norm * compassion_c + sociocultural_norm * respectfulness_c + sociocultural_norm * trust_c 
                            + sociocultural_norm * organization_c + sociocultural_norm * productiveness_c + sociocultural_norm * responsibility_c 
                            + sociocultural_norm * aesthetic_c + sociocultural_norm * creative_c + sociocultural_norm * intellectual_c 
                            + sociocultural_norm * assertiveness_c + sociocultural_norm * energy_c + sociocultural_norm * sociability_c 
                            + sociocultural_norm * anxiety_c + sociocultural_norm * depression_c + emotional_c 
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))


##conduct externals-as-predictors model in Experiment 3
exp3_main_effects_model_externals = lmBF(personal_preference ~ sociocultural_norm + rational_c + nfc_c + attention_c
                         + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_rat = lmBF(personal_preference ~ sociocultural_norm + nfc_c + attention_c
                         + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_nfc = lmBF(personal_preference ~ sociocultural_norm + rational_c + attention_c
                         + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_att = lmBF(personal_preference ~ sociocultural_norm + rational_c + nfc_c
                         + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_norms_externals = lmBF(personal_preference ~ rational_c + nfc_c + attention_c
                         + id, data = dat3, whichRandom = c('id'))

exp3_full_model_externals = lmBF(personal_preference ~ sociocultural_norm * rational_c + sociocultural_norm * nfc_c + sociocultural_norm * attention_c
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IErat = lmBF(personal_preference ~ sociocultural_norm + rational_c + sociocultural_norm * nfc_c + sociocultural_norm * attention_c
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEnfc = lmBF(personal_preference ~ sociocultural_norm * rational_c + nfc_c + sociocultural_norm * attention_c
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))

exp3_model_no_IEatt = lmBF(personal_preference ~ sociocultural_norm * rational_c + sociocultural_norm * nfc_c + attention_c
                            + id + sociocultural_norm:id, data = dat3, whichRandom = c('id','sociocultural_norm:id'))


#compute Bayes factors for the Big Five domains (Table 1)

##Experiment 1
exp1_agr_BF <- extractBF(exp1_main_effects_model/exp1_model_no_agr, onlybf = TRUE)
exp1_cns_BF <- extractBF(exp1_main_effects_model/exp1_model_no_cns, onlybf = TRUE)
exp1_opn_BF <- extractBF(exp1_main_effects_model/exp1_model_no_opn, onlybf = TRUE)
exp1_ext_BF <- extractBF(exp1_main_effects_model/exp1_model_no_ext, onlybf = TRUE)
exp1_neu_BF <- extractBF(exp1_main_effects_model/exp1_model_no_neu, onlybf = TRUE)
exp1_norms_BF <- "2e+845" #value from (exp1_main_effects_model/exp1_model_no_norms)
exp1_IEagr_BF <- extractBF(exp1_full_model/exp1_model_no_IEagr, onlybf = TRUE)
exp1_IEcns_BF <- extractBF(exp1_full_model/exp1_model_no_IEcns, onlybf = TRUE)
exp1_IEopn_BF <- extractBF(exp1_full_model/exp1_model_no_IEopn, onlybf = TRUE)
exp1_IEext_BF <- extractBF(exp1_full_model/exp1_model_no_IEext, onlybf = TRUE)
exp1_IEneu_BF <- extractBF(exp1_full_model/exp1_model_no_IEneu, onlybf = TRUE)
##Experiment 2
exp2_agr_BF <- extractBF(exp2_main_effects_model/exp2_model_no_agr, onlybf = TRUE)
exp2_cns_BF <- extractBF(exp2_main_effects_model/exp2_model_no_cns, onlybf = TRUE)
exp2_opn_BF <- extractBF(exp2_main_effects_model/exp2_model_no_opn, onlybf = TRUE)
exp2_ext_BF <- extractBF(exp2_main_effects_model/exp2_model_no_ext, onlybf = TRUE)
exp2_neu_BF <- extractBF(exp2_main_effects_model/exp2_model_no_neu, onlybf = TRUE)
exp2_norms_BF <- "2e+845" #value from (exp2_main_effects_model/exp2_model_no_norms)
exp2_IEagr_BF <- extractBF(exp2_full_model/exp2_model_no_IEagr, onlybf = TRUE)
exp2_IEcns_BF <- extractBF(exp2_full_model/exp2_model_no_IEcns, onlybf = TRUE)
exp2_IEopn_BF <- extractBF(exp2_full_model/exp2_model_no_IEopn, onlybf = TRUE)
exp2_IEext_BF <- extractBF(exp2_full_model/exp2_model_no_IEext, onlybf = TRUE)
exp2_IEneu_BF <- extractBF(exp2_full_model/exp2_model_no_IEneu, onlybf = TRUE)
##Experiment 3
exp3_agr_BF <- extractBF(exp3_main_effects_model/exp3_model_no_agr, onlybf = TRUE)
exp3_cns_BF <- extractBF(exp3_main_effects_model/exp3_model_no_cns, onlybf = TRUE)
exp3_opn_BF <- extractBF(exp3_main_effects_model/exp3_model_no_opn, onlybf = TRUE)
exp3_ext_BF <- extractBF(exp3_main_effects_model/exp3_model_no_ext, onlybf = TRUE)
exp3_neu_BF <- extractBF(exp3_main_effects_model/exp3_model_no_neu, onlybf = TRUE)
exp3_norms_BF <- "2e+845" #value from (exp3_main_effects_model/exp3_model_no_norms)
exp3_IEagr_BF <- extractBF(exp3_full_model/exp3_model_no_IEagr, onlybf = TRUE)
exp3_IEcns_BF <- extractBF(exp3_full_model/exp3_model_no_IEcns, onlybf = TRUE)
exp3_IEopn_BF <- extractBF(exp3_full_model/exp3_model_no_IEopn, onlybf = TRUE)
exp3_IEext_BF <- extractBF(exp3_full_model/exp3_model_no_IEext, onlybf = TRUE)
exp3_IEneu_BF <- extractBF(exp3_full_model/exp3_model_no_IEneu, onlybf = TRUE)
rm(exp1_main_effects_model,exp1_model_no_agr,exp1_model_no_cns,exp1_model_no_opn,exp1_model_no_ext,exp1_model_no_neu,exp1_model_no_norms,
   exp1_full_model,exp1_model_no_IEagr,exp1_model_no_IEcns,exp1_model_no_IEopn,exp1_model_no_IEext,exp1_model_no_IEneu,
   exp2_main_effects_model,exp2_model_no_agr,exp2_model_no_cns,exp2_model_no_opn,exp2_model_no_ext,exp2_model_no_neu,exp2_model_no_norms,
   exp2_full_model,exp2_model_no_IEagr,exp2_model_no_IEcns,exp2_model_no_IEopn,exp2_model_no_IEext,exp2_model_no_IEneu,
   exp3_main_effects_model,exp3_model_no_agr,exp3_model_no_cns,exp3_model_no_opn,exp3_model_no_ext,exp3_model_no_neu,exp3_model_no_norms,
   exp3_full_model,exp3_model_no_IEagr,exp3_model_no_IEcns,exp3_model_no_IEopn,exp3_model_no_IEext,exp3_model_no_IEneu)


#compute Bayes factors for the Big Five facets (Table 2)

##Experiment 2
exp2_comp_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_comp, onlybf = TRUE)
exp2_respect_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_respect, onlybf = TRUE)
exp2_trust_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_trust, onlybf = TRUE)
exp2_orga_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_orga, onlybf = TRUE)
exp2_prod_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_prod, onlybf = TRUE)
exp2_respo_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_respo, onlybf = TRUE)
exp2_aest_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_aest, onlybf = TRUE)
exp2_crea_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_crea, onlybf = TRUE)
exp2_intel_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_intel, onlybf = TRUE)
exp2_asse_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_asse, onlybf = TRUE)
exp2_ener_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_ener, onlybf = TRUE)
exp2_soci_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_soci, onlybf = TRUE)
exp2_anxi_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_anxi, onlybf = TRUE)
exp2_depr_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_depr, onlybf = TRUE)
exp2_emot_BF <- extractBF(exp2_main_effects_model_facets/exp2_model_no_emot, onlybf = TRUE)
exp2_norms_facets_BF <- "3e+1104" #value from (exp2_main_effects_model_facets/exp2_model_no_norms_facets)
exp2_IEcomp_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEcomp, onlybf = TRUE)
exp2_IErespect_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IErespect, onlybf = TRUE)
exp2_IEtrust_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEtrust, onlybf = TRUE)
exp2_IEorga_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEorga, onlybf = TRUE)
exp2_IEprod_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEprod, onlybf = TRUE)
exp2_IErespo_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IErespo, onlybf = TRUE)
exp2_IEaest_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEaest, onlybf = TRUE)
exp2_IEcrea_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEcrea, onlybf = TRUE)
exp2_IEintel_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEintel, onlybf = TRUE)
exp2_IEasse_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEasse, onlybf = TRUE)
exp2_IEener_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEener, onlybf = TRUE)
exp2_IEsoci_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEsoci, onlybf = TRUE)
exp2_IEanxi_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEanxi, onlybf = TRUE)
exp2_IEdepr_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEdepr, onlybf = TRUE)
exp2_IEemot_BF <- extractBF(exp2_full_model_facets/exp2_model_no_IEemot, onlybf = TRUE)
##Experiment 3
exp3_comp_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_comp, onlybf = TRUE)
exp3_respect_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_respect, onlybf = TRUE)
exp3_trust_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_trust, onlybf = TRUE)
exp3_orga_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_orga, onlybf = TRUE)
exp3_prod_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_prod, onlybf = TRUE)
exp3_respo_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_respo, onlybf = TRUE)
exp3_aest_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_aest, onlybf = TRUE)
exp3_crea_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_crea, onlybf = TRUE)
exp3_intel_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_intel, onlybf = TRUE)
exp3_asse_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_asse, onlybf = TRUE)
exp3_ener_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_ener, onlybf = TRUE)
exp3_soci_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_soci, onlybf = TRUE)
exp3_anxi_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_anxi, onlybf = TRUE)
exp3_depr_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_depr, onlybf = TRUE)
exp3_emot_BF <- extractBF(exp3_main_effects_model_facets/exp3_model_no_emot, onlybf = TRUE)
exp3_norms_facets_BF <- "3e+1104" #value from (exp3_main_effects_model_facets/exp3_model_no_norms_facets)
exp3_IEcomp_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEcomp, onlybf = TRUE)
exp3_IErespect_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IErespect, onlybf = TRUE)
exp3_IEtrust_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEtrust, onlybf = TRUE)
exp3_IEorga_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEorga, onlybf = TRUE)
exp3_IEprod_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEprod, onlybf = TRUE)
exp3_IErespo_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IErespo, onlybf = TRUE)
exp3_IEaest_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEaest, onlybf = TRUE)
exp3_IEcrea_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEcrea, onlybf = TRUE)
exp3_IEintel_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEintel, onlybf = TRUE)
exp3_IEasse_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEasse, onlybf = TRUE)
exp3_IEener_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEener, onlybf = TRUE)
exp3_IEsoci_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEsoci, onlybf = TRUE)
exp3_IEanxi_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEanxi, onlybf = TRUE)
exp3_IEdepr_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEdepr, onlybf = TRUE)
exp3_IEemot_BF <- extractBF(exp3_full_model_facets/exp3_model_no_IEemot, onlybf = TRUE)
rm(exp2_main_effects_model_facets,exp2_model_no_comp,exp2_model_no_respect,exp2_model_no_trust,exp2_model_no_orga,exp2_model_no_prod,exp2_model_no_respo,exp2_model_no_aest,exp2_model_no_crea,exp2_model_no_intel,exp2_model_no_asse,
   exp2_model_no_ener,exp2_model_no_soci,exp2_model_no_anxi,exp2_model_no_depr,exp2_model_no_emot,exp2_model_no_norms_facets,exp2_full_model_facets,exp2_model_no_IEcomp,exp2_model_no_IErespect,exp2_model_no_IEtrust,exp2_model_no_IEorga,
   exp2_model_no_IEprod,exp2_model_no_IErespo,exp2_model_no_IEaest,exp2_model_no_IEcrea,exp2_model_no_IEintel,exp2_model_no_IEasse,exp2_model_no_IEener,exp2_model_no_IEsoci,exp2_model_no_IEanxi,exp2_model_no_IEdepr,exp2_model_no_IEemot,
   exp3_main_effects_model_facets,exp3_model_no_comp,exp3_model_no_respect,exp3_model_no_trust,exp3_model_no_orga,exp3_model_no_prod,exp3_model_no_respo,exp3_model_no_aest,exp3_model_no_crea,exp3_model_no_intel,exp3_model_no_asse,
   exp3_model_no_ener,exp3_model_no_soci,exp3_model_no_anxi,exp3_model_no_depr,exp3_model_no_emot,exp3_model_no_norms_facets,exp3_full_model_facets,exp3_model_no_IEcomp,exp3_model_no_IErespect,exp3_model_no_IEtrust,exp3_model_no_IEorga,
   exp3_model_no_IEprod,exp3_model_no_IErespo,exp3_model_no_IEaest,exp3_model_no_IEcrea,exp3_model_no_IEintel,exp3_model_no_IEasse,exp3_model_no_IEener,exp3_model_no_IEsoci,exp3_model_no_IEanxi,exp3_model_no_IEdepr,exp3_model_no_IEemot)


#compute Bayes factors for the external predictors in Experiment 3
exp3_rat_BF <- extractBF(exp3_main_effects_model_externals/exp3_model_no_rat, onlybf = TRUE)
exp3_nfc_BF <- extractBF(exp3_main_effects_model_externals/exp3_model_no_nfc, onlybf = TRUE)
exp3_att_BF <- extractBF(exp3_main_effects_model_externals/exp3_model_no_att, onlybf = TRUE)
exp3_norms_externals_BF <- "1e+1172" #value from (exp3_main_effects_model_externals/exp3_model_no_norms_externals)
exp3_IErat_BF <- extractBF(exp3_full_model_externals/exp3_model_no_IErat, onlybf = TRUE)
exp3_IEnfc_BF <- extractBF(exp3_full_model_externals/exp3_model_no_IEnfc, onlybf = TRUE)
exp3_IEatt_BF <- extractBF(exp3_full_model_externals/exp3_model_no_IEatt, onlybf = TRUE)
rm(exp3_main_effects_model_externals,exp3_model_no_rat,exp3_model_no_nfc,exp3_model_no_att,exp3_model_no_norms_externals,
   exp3_full_model_externals,exp3_model_no_IErat,exp3_model_no_IEnfc,exp3_model_no_IEatt)