# inspect the structure of the dataframe
head(data_low_vision, 10)
#------
# restrict dataset to one MNREAD test per subject (regular polarity only)
data_regular <- data_low_vision %>%
filter (polarity == "regular")
# run the NLME model for data grouped by subject
nlme_model <- nlmeModel(data_regular, ps, vd, rt, err, subject)
#------
# extract the critical print size required
# to achieve 40 words/min (ie. spot reading) according to the NLME fit
nlmePredict_PS(nlme_model, 40)
#------
# extract the critical print size required
# to achieve 80 words/min (ie. fluent reading) according to the NLME fit
nlmePredict_PS(nlme_model, 80)
Run the code above in your browser using DataLab