library(tidyverse)
library(here)
library(lmerTest)

Load data

trial_metadata <- read.csv(here("data","metadata","trial_metadata.csv"))
trial_summary_data <- read.csv(here("data","processed_data","trial_summary_data.csv"))

usable_trials <- trial_summary_data |>
  filter(exclude_participant_insufficient_data == 0 & trial_exclusion == 0 & exclude_participant == 0) 

# Merging with similarity information and mean-centering main effects
trials_with_effect_vars <- usable_trials |>
  left_join(trial_metadata) |>
  mutate(text_similarity_centered = (text_similarity - mean(text_similarity, na.rm = TRUE)),
         # Converting scale of age to months from days before centering
         test_age_centered = (SubjectInfo.testAge / 30 - mean(SubjectInfo.testAge / 30 , na.rm = TRUE))
  )

Joining with `by = join_by(Trials.trialID, Trials.targetImage,
Trials.distractorImage, Trials.imagePair)`

Run mixed-effects model

main_effect <- lmer(corrected_target_looking ~ (text_similarity_centered * test_age_centered) 
                    + (text_similarity_centered | SubjectInfo.subjID) 
                    + (1|Trials.targetImage) 
                    + (1|Trials.imagePair), 
                    data = trials_with_effect_vars)