library(languageR)
library(ggplot2)
library(dplyr)
library(reshape2)
library(rms)
library(lattice)
library(grid)
library(lme4)
library(effects)
library(languageR)
library(lmerTest)

getwd()
setwd("")

RD2.all<- read.csv("", header=T, na.strings = "-999")
summary(RD2.all)

describe(RD2.all$TT)
describe(RD2.all$TD)
describe(RD2.all$TN)

RD2.all.correct <- filter(RD2.all, RT > 80)
RD2.all.correct1 <- filter(RD2.all.correct, RT < 700)

length(RD.all$LATENCY_CORRECT_SACCADE)
length(RD.all.correct1$LATENCY_CORRECT_SACCADE)

summary(RD.all.correct1)




RD2.all.lmer.1 <- lmer(LATENCY_CORRECT_SACCADE ~ 
                                
                                GROUP + (1|GROUP),
                      
                            
                          
                              
                              data = RD2.all.error, REML = FALSE)

summary(RD2.all.lmer.1)


RD.all.lmer.2 <- lmer(TD ~ 
                        
                        TT + TN + GROUP + (1|GROUP),
                      
                      
                      
                      
                      data = RD.all, REML = FALSE)

summary(RD.all.lmer.2)


RD2.all.lmer.3 <- lmer(TD_TN ~ 
                        
                        GROUP + (1|GROUP),
                      
                      
                      
                      
                      data = RD2.all, REML = FALSE)

summary(RD2.all.lmer.3)

RD.all.lmer.4 <- lmer(TN_TT ~ 
                        
                        TD_TN + GROUP + (1|GROUP),
                      
                      
                      
                      
                      data = RD.all, REML = FALSE)

summary(RD.all.lmer.4)