library(readxl) library(dplyr) library(tidyr) library(ordinal)

df <- read_excel(“~/Desktop/DigitalvsPhysical.xlsx”)

df <- df %>% mutate( Score = factor(Score, ordered = TRUE), Environment = factor(Environment, levels = c(“digital”, “physical”)), Panelist = factor(Panelist), Item = factor(Item) )

str(df)

clmm_model <- clmm( Score ~ Environment + (1 | Panelist) + (1 | Item), data = df, Hess = TRUE)

summary(clmm_model)

confint(clmm_model) citation() getRversion()