#################################################################################################### #################################################################################################### ################################ Code Experiment 3 ################################################# #################################################################################################### #################################################################################################### # R version 3.5.1 # lme4 package version 1.1-17 # lmerTest package version 3.0-1 # set working directory setwd("") data = read.table("Nahon_et_al_Truth feels easy_SUPPL_Data_Experiment_3.csv", header = TRUE, sep = ";", as.is = TRUE, na.strings = c("-77", "-99", "-66")) newdata <- subset(data, andere_Sprache == 2, select = c(lfdn, group, Noise1_DE, Noise2_DE, Noise3_DE, Noise4_DE, Noise5_DE, Noise6_DE, Noise7_DE, Noise8_DE, Noise1_SVK, Noise2_SVK, Noise3_SVK, Noise4_SVK, Noise5_SVK, Noise6_SVK, Noise7_SVK, Noise8_SVK, gr1_f_unknown1, gr1_f_unknown2, gr1_f_unknown3, gr1_f_unknown4, gr1_f_unknown5, gr1_f_unknown6, gr1_f_unknown7, gr1_f_unknown8, gr1_t_unknown1, gr1_t_unknown2, gr1_t_unknown3, gr1_t_unknown4, gr1_t_unknown5, gr1_t_unknown6, gr1_t_unknown7, gr1_t_unknown8, gr1_f_known1, gr1_f_known2, gr1_f_known3, gr1_f_known4, gr1_f_known5, gr1_f_known6, gr1_f_known7, gr1_f_known8, gr1_t_known1, gr1_t_known2, gr1_t_known3, gr1_t_known4, gr1_t_known5, gr1_t_known6, gr1_t_known7, gr1_t_known8, gr2_f_unknown1, gr2_f_unknown2, gr2_f_unknown3, gr2_f_unknown4, gr2_f_unknown5, gr2_f_unknown6, gr2_f_unknown7, gr2_f_unknown8, gr2_t_unknown1, gr2_t_unknown2, gr2_t_unknown3, gr2_t_unknown4, gr2_t_unknown5, gr2_t_unknown6, gr2_t_unknown7, gr2_t_unknown8, gr2_f_known1, gr2_f_known2, gr2_f_known3, gr2_f_known4, gr2_f_known5, gr2_f_known6, gr2_f_known7, gr2_f_known8, gr2_t_known1, gr2_t_known2, gr2_t_known3, gr2_t_known4, gr2_t_known5, gr2_t_known6, gr2_t_known7, gr2_t_known8, audio_gr1_1, audio_gr1_2, audio_gr1_3, audio_gr1_4, audio_gr1_5, audio_gr1_6, audio_gr1_7, audio_gr1_8, audio_gr1_9, audio_gr1_10, audio_gr1_11, audio_gr1_12, audio_gr1_13, audio_gr1_14, audio_gr1_15, audio_gr1_16, audio_gr1_17, audio_gr1_18, audio_gr1_19, audio_gr1_20, audio_gr1_21, audio_gr1_22, audio_gr1_23, audio_gr1_24, audio_gr1_25, audio_gr1_26, audio_gr1_27, audio_gr1_28, audio_gr1_29, audio_gr1_30, audio_gr1_31, audio_gr1_32, audio_gr2_1, audio_gr2_2, audio_gr2_3, audio_gr2_4, audio_gr2_5, audio_gr2_6, audio_gr2_7, audio_gr2_8, audio_gr2_9, audio_gr2_10, audio_gr2_11, audio_gr2_12, audio_gr2_13, audio_gr2_14, audio_gr2_15, audio_gr2_16, audio_gr2_17, audio_gr2_18, audio_gr2_19, audio_gr2_20, audio_gr2_21, audio_gr2_22, audio_gr2_23, audio_gr2_24, audio_gr2_25, audio_gr2_26, audio_gr2_27, audio_gr2_28, audio_gr2_29, audio_gr2_30, audio_gr2_31, audio_gr2_32)) library(reshape2) newdata_long <- melt(newdata, id.vars = c("lfdn", "group", "Noise1_DE", "Noise2_DE", "Noise3_DE", "Noise4_DE", "Noise5_DE", "Noise6_DE", "Noise7_DE", "Noise8_DE", "Noise1_SVK", "Noise2_SVK", "Noise3_SVK", "Noise4_SVK", "Noise5_SVK", "Noise6_SVK", "Noise7_SVK", "Noise8_SVK", "gr1_f_unknown1", "gr1_f_unknown2", "gr1_f_unknown3", "gr1_f_unknown4", "gr1_f_unknown5", "gr1_f_unknown6", "gr1_f_unknown7", "gr1_f_unknown8", "gr1_t_unknown1", "gr1_t_unknown2", "gr1_t_unknown3", "gr1_t_unknown4", "gr1_t_unknown5", "gr1_t_unknown6", "gr1_t_unknown7", "gr1_t_unknown8", "gr1_f_known1", "gr1_f_known2", "gr1_f_known3", "gr1_f_known4", "gr1_f_known5", "gr1_f_known6", "gr1_f_known7", "gr1_f_known8", "gr1_t_known1", "gr1_t_known2", "gr1_t_known3", "gr1_t_known4", "gr1_t_known5", "gr1_t_known6", "gr1_t_known7", "gr1_t_known8", "gr2_f_unknown1", "gr2_f_unknown2", "gr2_f_unknown3", "gr2_f_unknown4", "gr2_f_unknown5", "gr2_f_unknown6", "gr2_f_unknown7", "gr2_f_unknown8", "gr2_t_unknown1", "gr2_t_unknown2", "gr2_t_unknown3", "gr2_t_unknown4", "gr2_t_unknown5", "gr2_t_unknown6", "gr2_t_unknown7", "gr2_t_unknown8", "gr2_f_known1", "gr2_f_known2", "gr2_f_known3", "gr2_f_known4", "gr2_f_known5", "gr2_f_known6", "gr2_f_known7", "gr2_f_known8", "gr2_t_known1", "gr2_t_known2", "gr2_t_known3", "gr2_t_known4", "gr2_t_known5", "gr2_t_known6", "gr2_t_known7", "gr2_t_known8")) newdata_long <- subset(newdata_long, value > 0) newdata_long$truth <- ifelse(newdata_long$variable == "audio_gr1_1" | newdata_long$variable == "audio_gr1_2" |newdata_long$variable == "audio_gr1_3" | newdata_long$variable == "audio_gr1_4" | newdata_long$variable == "audio_gr1_5" | newdata_long$variable == "audio_gr1_6" | newdata_long$variable == "audio_gr1_7" | newdata_long$variable == "audio_gr1_8" | newdata_long$variable == "audio_gr1_17" | newdata_long$variable == "audio_gr1_18" |newdata_long$variable == "audio_gr1_19" | newdata_long$variable == "audio_gr1_20" | newdata_long$variable == "audio_gr1_21" | newdata_long$variable == "audio_gr1_22" | newdata_long$variable == "audio_gr1_23" | newdata_long$variable == "audio_gr1_24" | newdata_long$variable == "audio_gr2_1" | newdata_long$variable == "audio_gr2_2" |newdata_long$variable == "audio_gr2_3" | newdata_long$variable == "audio_gr2_4" | newdata_long$variable == "audio_gr2_5" | newdata_long$variable == "audio_gr2_6" | newdata_long$variable == "audio_gr2_7" | newdata_long$variable == "audio_gr2_8" | newdata_long$variable == "audio_gr2_17" | newdata_long$variable == "audio_gr2_18" |newdata_long$variable == "audio_gr2_19" | newdata_long$variable == "audio_gr2_20" | newdata_long$variable == "audio_gr2_21" | newdata_long$variable == "audio_gr2_22" | newdata_long$variable == "audio_gr2_23" | newdata_long$variable == "audio_gr2_24", -1, 1) newdata_long$known <- ifelse(newdata_long$variable == "audio_gr1_1" | newdata_long$variable == "audio_gr1_2" |newdata_long$variable == "audio_gr1_3" | newdata_long$variable == "audio_gr1_4" | newdata_long$variable == "audio_gr1_5" | newdata_long$variable == "audio_gr1_6" | newdata_long$variable == "audio_gr1_7" | newdata_long$variable == "audio_gr1_8" | newdata_long$variable == "audio_gr1_9" | newdata_long$variable == "audio_gr1_10" |newdata_long$variable == "audio_gr1_11" | newdata_long$variable == "audio_gr1_12" | newdata_long$variable == "audio_gr1_13" | newdata_long$variable == "audio_gr1_14" | newdata_long$variable == "audio_gr1_15" | newdata_long$variable == "audio_gr1_16" | newdata_long$variable == "audio_gr2_1" | newdata_long$variable == "audio_gr2_2" |newdata_long$variable == "audio_gr2_3" | newdata_long$variable == "audio_gr2_4" | newdata_long$variable == "audio_gr2_5" | newdata_long$variable == "audio_gr2_6" | newdata_long$variable == "audio_gr2_7" | newdata_long$variable == "audio_gr2_8" | newdata_long$variable == "audio_gr2_9" | newdata_long$variable == "audio_gr2_10" |newdata_long$variable == "audio_gr2_11" | newdata_long$variable == "audio_gr2_12" | newdata_long$variable == "audio_gr2_13" | newdata_long$variable == "audio_gr2_14" | newdata_long$variable == "audio_gr2_15" | newdata_long$variable == "audio_gr2_16", -1, 1) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_17" | newdata_long$variable == "audio_gr2_17" | newdata_long$variable == "audio_gr1_25" | newdata_long$variable == "audio_gr2_25", newdata_long$Noise1_DE, 0) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_18" | newdata_long$variable == "audio_gr2_18" | newdata_long$variable == "audio_gr1_26" | newdata_long$variable == "audio_gr2_26", newdata_long$Noise2_DE, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_19" | newdata_long$variable == "audio_gr2_19" | newdata_long$variable == "audio_gr1_27" | newdata_long$variable == "audio_gr2_27", newdata_long$Noise3_DE, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_20" | newdata_long$variable == "audio_gr2_20" | newdata_long$variable == "audio_gr1_28" | newdata_long$variable == "audio_gr2_28", newdata_long$Noise4_DE, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_21" | newdata_long$variable == "audio_gr2_21" | newdata_long$variable == "audio_gr1_29" | newdata_long$variable == "audio_gr2_29", newdata_long$Noise5_DE, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_22" | newdata_long$variable == "audio_gr2_22" | newdata_long$variable == "audio_gr1_30" | newdata_long$variable == "audio_gr2_30", newdata_long$Noise6_DE, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_23" | newdata_long$variable == "audio_gr2_23" | newdata_long$variable == "audio_gr1_31" | newdata_long$variable == "audio_gr2_31", newdata_long$Noise7_DE, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_24" | newdata_long$variable == "audio_gr2_24" | newdata_long$variable == "audio_gr1_32" | newdata_long$variable == "audio_gr2_32", newdata_long$Noise8_DE, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_1" | newdata_long$variable == "audio_gr2_1" | newdata_long$variable == "audio_gr1_9" | newdata_long$variable == "audio_gr2_9", newdata_long$Noise1_SVK, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_2" | newdata_long$variable == "audio_gr2_2" | newdata_long$variable == "audio_gr1_10" | newdata_long$variable == "audio_gr2_10", newdata_long$Noise2_SVK, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_3" | newdata_long$variable == "audio_gr2_3" | newdata_long$variable == "audio_gr1_11" | newdata_long$variable == "audio_gr2_11", newdata_long$Noise3_SVK, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_4" | newdata_long$variable == "audio_gr2_4" | newdata_long$variable == "audio_gr1_12" | newdata_long$variable == "audio_gr2_12", newdata_long$Noise4_SVK, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_5" | newdata_long$variable == "audio_gr2_5" | newdata_long$variable == "audio_gr1_13" | newdata_long$variable == "audio_gr2_13", newdata_long$Noise5_SVK, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_6" | newdata_long$variable == "audio_gr2_6" | newdata_long$variable == "audio_gr1_14" | newdata_long$variable == "audio_gr2_14", newdata_long$Noise6_SVK, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_7" | newdata_long$variable == "audio_gr2_7" | newdata_long$variable == "audio_gr1_15" | newdata_long$variable == "audio_gr2_15", newdata_long$Noise7_SVK, newdata_long$Noise) newdata_long$Noise <- ifelse(newdata_long$variable == "audio_gr1_8" | newdata_long$variable == "audio_gr2_8" | newdata_long$variable == "audio_gr1_16" | newdata_long$variable == "audio_gr2_16", newdata_long$Noise8_SVK, newdata_long$Noise) library(plyr) newdata_long$Noise <- revalue(newdata_long$Noise, c(" " = "1", " "= "2", " " = "3", " "= "4", " " = "5", " "= "6", " " = "7", " "= "8", " " = "1", " "= "2", " " = "3", " "= "4", " " = "5", " "= "6", " " = "7", " "= "8")) newdata_long$Noise <- as.factor(newdata_long$Noise) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_1", newdata_long$gr1_f_unknown1, 0) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_2", newdata_long$gr1_f_unknown2, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_3", newdata_long$gr1_f_unknown3, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_4", newdata_long$gr1_f_unknown4, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_5", newdata_long$gr1_f_unknown5, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_6", newdata_long$gr1_f_unknown6, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_7", newdata_long$gr1_f_unknown7, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_8", newdata_long$gr1_f_unknown8, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_9", newdata_long$gr1_t_unknown1, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_10", newdata_long$gr1_t_unknown2, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_11", newdata_long$gr1_t_unknown3, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_12", newdata_long$gr1_t_unknown4, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_13", newdata_long$gr1_t_unknown5, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_14", newdata_long$gr1_t_unknown6, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_15", newdata_long$gr1_t_unknown7, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_16", newdata_long$gr1_t_unknown8, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_17", newdata_long$gr1_f_known1, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_18", newdata_long$gr1_f_known2, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_19", newdata_long$gr1_f_known3, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_20", newdata_long$gr1_f_known4, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_21", newdata_long$gr1_f_known5, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_22", newdata_long$gr1_f_known6, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_23", newdata_long$gr1_f_known7, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_24", newdata_long$gr1_f_known8, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_25", newdata_long$gr1_t_known1, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_26", newdata_long$gr1_t_known2, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_27", newdata_long$gr1_t_known3, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_28", newdata_long$gr1_t_known4, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_29", newdata_long$gr1_t_known5, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_30", newdata_long$gr1_t_known6, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_31", newdata_long$gr1_t_known7, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr1_32", newdata_long$gr1_t_known8, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_1", newdata_long$gr2_f_unknown1, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_2", newdata_long$gr2_f_unknown2, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_3", newdata_long$gr2_f_unknown3, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_4", newdata_long$gr2_f_unknown4, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_5", newdata_long$gr2_f_unknown5, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_6", newdata_long$gr2_f_unknown6, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_7", newdata_long$gr2_f_unknown7, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_8", newdata_long$gr2_f_unknown8, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_9", newdata_long$gr2_t_unknown1, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_10", newdata_long$gr2_t_unknown2, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_11", newdata_long$gr2_t_unknown3, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_12", newdata_long$gr2_t_unknown4, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_13", newdata_long$gr2_t_unknown5, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_14", newdata_long$gr2_t_unknown6, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_15", newdata_long$gr2_t_unknown7, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_16", newdata_long$gr2_t_unknown8, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_17", newdata_long$gr2_f_known1, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_18", newdata_long$gr2_f_known2, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_19", newdata_long$gr2_f_known3, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_20", newdata_long$gr2_f_known4, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_21", newdata_long$gr2_f_known5, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_22", newdata_long$gr2_f_known6, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_23", newdata_long$gr2_f_known7, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_24", newdata_long$gr2_f_known8, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_25", newdata_long$gr2_t_known1, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_26", newdata_long$gr2_t_known2, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_27", newdata_long$gr2_t_known3, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_28", newdata_long$gr2_t_known4, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_29", newdata_long$gr2_t_known5, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_30", newdata_long$gr2_t_known6, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_31", newdata_long$gr2_t_known7, newdata_long$Stimuli) newdata_long$Stimuli <- ifelse(newdata_long$variable == "audio_gr2_32", newdata_long$gr2_t_known8, newdata_long$Stimuli) newdata_long$Stimuli <- as.factor(newdata_long$Stimuli) newdata_long$lfdn <- as.factor(newdata_long$lfdn) newdata_long$truth_fact <- as.factor(newdata_long$truth) newdata_long$known_fact <- as.factor(newdata_long$known) library(lme4) library(lmerTest) model <- lmer(value ~ 1 + truth_fact + known_fact + truth_fact * known_fact + (1 + truth_fact + known_fact + truth_fact * known_fact | lfdn) + (1 + truth_fact + known_fact + truth_fact * known_fact | Noise) + (1 | Stimuli), control=lmerControl(optimizer="bobyqa"), newdata_long) anova(model) library(phia) testInteractions(model, fixed = "known_fact", across = "truth_fact", adjustment = "none")