install.packages("foreign")
install.packages("lme4")
install.packages("MuMIn")

library(foreign)
library(lme4)
library(MuMIn)


library(lmerTest)
library(multilevel)


#### Open the data set. Please add file path (e.g., of the spss-data set) and file name
setwd("C:/file")

diagnose <- read.spss(file = "file_name.sav", to.data.frame = T, use.value.labels = F)
diagnose$correct.n <- as.numeric(diagnose$correct)



########## Correctness of the diagnosis

kglm <- glmer(correct.n ~ 1 + format +  vis    + schoolexamgrade +   order + clinical_year +(1 |id),data= diagnose, family = "binomial")

summary(kglm)

r.squaredGLMM(kglm)

########## Speed of diagnosis - without eight outliers with time > 3xinterquartile range

diagnose2 <- read.spss(file = "file_name.sav", to.data.frame = T, use.value.labels = F)

diagnose2$time_exclude <- as.numeric(diagnose2$time_exclude)

k2glm <- lmer(scale(time_exclude) ~ 1 + scale(format)+ scale(vis) +scale(context_Cushing)+ scale(context_FHH) +scale(context_thyroid)   +scale(schoolexamgrade) +scale(order)    +(1 | id),data= diagnose2)
summary(k2glm)