Installing and Checking Packages

library("lavaan")
library("effects")
library("psych")
library("lme4")
library("lmerTest")
library("sjstats")
library("sjmisc")
library("lattice")

Importing and Attaching Data

library("sjmisc")
setwd("C:/Users/pickrese/Dropbox/Vanderbilt/Cutting Lab/Reading Error Measure (RED)/PCRC/Poster Data/Wide Data Y2")
REDmasterPassageCombined <- read.csv("C:/Users/pickrese/Dropbox/Vanderbilt/Cutting Lab/Reading Error Measure (RED)/PCRC/Poster Data/Wide Data Y2/REDmasterPassageCombined.csv")
View(REDmasterPassageCombined)
attach(REDmasterPassageCombined)

Basic Reading Models

NullE <- lm(rc3_wj_br_ss_age ~ CWPM_E,  data = REDmasterPassageCombined)
summary(NullE)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_E, data = REDmasterPassageCombined)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -22.6877  -5.2226  -0.1916   5.6572  17.0930 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 81.27224    1.74766   46.50   <2e-16 ***
## CWPM_E       0.28467    0.01766   16.12   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.76 on 150 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.6341, Adjusted R-squared:  0.6317 
## F-statistic:   260 on 1 and 150 DF,  p-value: < 2.2e-16
M1 <- lm(rc3_wj_br_ss_age~ CWPM_E + DeltaDSyM_E_Sum ,  data = REDmasterPassageCombined)
summary(M1)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_E + DeltaDSyM_E_Sum, data = REDmasterPassageCombined)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -21.0403  -5.2204  -0.5376   5.2748  16.7166 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     89.36673    2.18579  40.885  < 2e-16 ***
## CWPM_E           0.22596    0.01945  11.619  < 2e-16 ***
## DeltaDSyM_E_Sum -0.28984    0.05326  -5.442 2.12e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.111 on 149 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.6948, Adjusted R-squared:  0.6907 
## F-statistic: 169.6 on 2 and 149 DF,  p-value: < 2.2e-16
NullN <- lm(rc3_wj_br_ss_age ~ CWPM_N,  data = REDmasterPassageCombined)
summary(NullN)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N, data = REDmasterPassageCombined)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -21.7427  -6.0862   0.3508   5.6966  21.6311 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 82.44366    1.85521   44.44   <2e-16 ***
## CWPM_N       0.22430    0.01545   14.52   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 8.271 on 150 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.5843, Adjusted R-squared:  0.5815 
## F-statistic: 210.8 on 1 and 150 DF,  p-value: < 2.2e-16
M2 <- lm(rc3_wj_br_ss_age~ CWPM_N +  DeltaDSyM_N_Sum ,  data = REDmasterPassageCombined)
summary(M2)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N + DeltaDSyM_N_Sum, data = REDmasterPassageCombined)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -20.8016  -6.3350  -0.1194   4.7234  30.5875 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     89.66602    2.36256  37.953  < 2e-16 ***
## CWPM_N           0.18276    0.01718  10.641  < 2e-16 ***
## DeltaDSyM_N_Sum -0.26286    0.05797  -4.534 1.18e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.78 on 149 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.6347, Adjusted R-squared:  0.6298 
## F-statistic: 129.5 on 2 and 149 DF,  p-value: < 2.2e-16
M3 <- lm(rc3_wj_br_ss_age~ CWPM_E*DeltaDSyM_E_Sum ,  data = REDmasterPassageCombined)
summary(M3)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_E * DeltaDSyM_E_Sum, data = REDmasterPassageCombined)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -21.0294  -5.1819  -0.5504   5.0560  15.3662 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            89.793397   2.175652  41.272   <2e-16 ***
## CWPM_E                  0.234371   0.019727  11.881   <2e-16 ***
## DeltaDSyM_E_Sum        -0.144267   0.090708  -1.590   0.1139    
## CWPM_E:DeltaDSyM_E_Sum -0.004284   0.002172  -1.973   0.0504 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.043 on 148 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.7026, Adjusted R-squared:  0.6966 
## F-statistic: 116.5 on 3 and 148 DF,  p-value: < 2.2e-16
plot(allEffects(M3), multiline=TRUE)

M4 <- lm(rc3_wj_br_ss_age~ CWPM_N*DeltaDSyM_N_Sum ,  data = REDmasterPassageCombined)
summary(M4)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N * DeltaDSyM_N_Sum, data = REDmasterPassageCombined)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -22.817  -5.221  -0.723   5.474  27.380 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            89.161153   2.332015  38.234   <2e-16 ***
## CWPM_N                  0.202187   0.018627  10.854   <2e-16 ***
## DeltaDSyM_N_Sum        -0.052294   0.102495  -0.510   0.6107    
## CWPM_N:DeltaDSyM_N_Sum -0.004679   0.001893  -2.472   0.0146 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.65 on 148 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.6492, Adjusted R-squared:  0.6421 
## F-statistic:  91.3 on 3 and 148 DF,  p-value: < 2.2e-16
plot(allEffects(M4), multiline=TRUE)

Model fit

anova(NullE,M1)
anova(NullN,M2)

Regions of significance

library("reghelper")
## Warning: package 'reghelper' was built under R version 3.5.3
## 
## Attaching package: 'reghelper'
## The following object is masked from 'package:psych':
## 
##     ICC
## The following object is masked from 'package:base':
## 
##     beta
simple_slopes(M4)

Analysis of poor readers <25%ile

detach(REDmasterPassageCombined)
REDmasterPassage25percentileE <- read.csv("C:/Users/pickrese/Dropbox/Vanderbilt/Cutting Lab/Reading Error Measure (RED)/PCRC/Poster Data/Wide Data Y2/REDmasterPassage25percentileE.csv")
View(REDmasterPassage25percentileE)

attach(REDmasterPassage25percentileE)

P25 <- lm(rc3_wj_br_ss_age ~ CWPM_E,  data = REDmasterPassage25percentileE)
summary(P25)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_E, data = REDmasterPassage25percentileE)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -16.9693  -5.0653   0.7464   5.4387  17.1698 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 76.41592    3.75615   20.34  < 2e-16 ***
## CWPM_E       0.35383    0.06911    5.12 5.36e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.324 on 48 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.3532, Adjusted R-squared:  0.3397 
## F-statistic: 26.21 on 1 and 48 DF,  p-value: 5.36e-06
M25 <- lm(rc3_wj_br_ss_age~ CWPM_E + DeltaDSyM_E_Sum ,  data = REDmasterPassage25percentileE)
summary(M25)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_E + DeltaDSyM_E_Sum, data = REDmasterPassage25percentileE)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -9.691 -5.113 -1.248  5.228 16.279 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     89.57388    4.36392  20.526  < 2e-16 ***
## CWPM_E           0.18568    0.06996   2.654   0.0108 *  
## DeltaDSyM_E_Sum -0.24267    0.05491  -4.419 5.81e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 6.221 on 47 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.5431, Adjusted R-squared:  0.5236 
## F-statistic: 27.93 on 2 and 47 DF,  p-value: 1.016e-08
M25interaction <- lm(rc3_wj_br_ss_age~ CWPM_E*DeltaDSyM_E_Sum ,  data = REDmasterPassage25percentileE)
summary(M25interaction)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_E * DeltaDSyM_E_Sum, data = REDmasterPassage25percentileE)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -9.627 -4.999 -1.624  5.380 15.889 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            88.833341   4.993587  17.789   <2e-16 ***
## CWPM_E                  0.203802   0.091081   2.238   0.0301 *  
## DeltaDSyM_E_Sum        -0.198151   0.151737  -1.306   0.1981    
## CWPM_E:DeltaDSyM_E_Sum -0.001283   0.004072  -0.315   0.7540    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 6.281 on 46 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.544,  Adjusted R-squared:  0.5143 
## F-statistic:  18.3 on 3 and 46 DF,  p-value: 5.899e-08
plot(allEffects(M25))

below 25 for narrative

detach(REDmasterPassage25percentileE)
REDmasterPassageNbelow25 <- read.csv("C:/Users/pickrese/Dropbox/Vanderbilt/Cutting Lab/Reading Error Measure (RED)/PCRC/Poster Data/Wide Data Y2/REDmasterPassageNbelow25.csv")
View(REDmasterPassageNbelow25)
attach(REDmasterPassageNbelow25)

PN25 <- lm(rc3_wj_br_ss_age ~ CWPM_N,  data = REDmasterPassageNbelow25)
summary(P25)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_E, data = REDmasterPassage25percentileE)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -16.9693  -5.0653   0.7464   5.4387  17.1698 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 76.41592    3.75615   20.34  < 2e-16 ***
## CWPM_E       0.35383    0.06911    5.12 5.36e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.324 on 48 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.3532, Adjusted R-squared:  0.3397 
## F-statistic: 26.21 on 1 and 48 DF,  p-value: 5.36e-06
MN25 <- lm(rc3_wj_br_ss_age~ CWPM_N + DeltaDSyM_N_Sum ,  data = REDmasterPassageNbelow25)
summary(MN25)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N + DeltaDSyM_N_Sum, data = REDmasterPassageNbelow25)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -11.383  -4.760  -0.696   2.540  25.851 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     80.56754    8.83136   9.123 7.01e-10 ***
## CWPM_N           0.30735    0.14218   2.162   0.0393 *  
## DeltaDSyM_N_Sum -0.14747    0.08257  -1.786   0.0849 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.925 on 28 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.4101, Adjusted R-squared:  0.3679 
## F-statistic: 9.732 on 2 and 28 DF,  p-value: 0.0006182
anova(PN25, MN25)
plot(allEffects(MN25))

between 25-50 percentile

detach(REDmasterPassageNbelow25)
REDmasterPassageNbelow50 <- read.csv("C:/Users/pickrese/Dropbox/Vanderbilt/Cutting Lab/Reading Error Measure (RED)/PCRC/Poster Data/Wide Data Y2/REDmasterPassageNbelow50.csv")
View(REDmasterPassageNbelow50)
attach(REDmasterPassageNbelow50)

PN50 <- lm(rc3_wj_br_ss_age ~ CWPM_N,  data = REDmasterPassageNbelow50)
summary(PN50)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N, data = REDmasterPassageNbelow50)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -12.5978  -4.5407  -0.1293   5.4022  15.8136 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  53.3350    13.9913   3.812 0.000726 ***
## CWPM_N        0.5886     0.1666   3.533 0.001502 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 6.939 on 27 degrees of freedom
## Multiple R-squared:  0.3161, Adjusted R-squared:  0.2908 
## F-statistic: 12.48 on 1 and 27 DF,  p-value: 0.001502
MN50 <- lm(rc3_wj_br_ss_age~ CWPM_N + DeltaDSyM_N_Sum ,  data = REDmasterPassageNbelow50)
summary(MN50)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N + DeltaDSyM_N_Sum, data = REDmasterPassageNbelow50)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -12.6175  -3.1372  -0.7115   4.2071  12.9171 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)      75.1447    13.4299   5.595 7.04e-06 ***
## CWPM_N            0.3880     0.1525   2.544  0.01724 *  
## DeltaDSyM_N_Sum  -0.5512     0.1607  -3.430  0.00202 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 5.867 on 26 degrees of freedom
## Multiple R-squared:  0.5292, Adjusted R-squared:  0.493 
## F-statistic: 14.61 on 2 and 26 DF,  p-value: 5.583e-05
anova(PN50, MN50)

between 50 and 75 percentile

detach(REDmasterPassageNbelow50)
REDmasterPassageNbelow75 <- read.csv("C:/Users/pickrese/Dropbox/Vanderbilt/Cutting Lab/Reading Error Measure (RED)/PCRC/Poster Data/Wide Data Y2/REDmasterPassageNbelow75.csv")
View(REDmasterPassageNbelow75)
attach(REDmasterPassageNbelow75)

PN75 <- lm(rc3_wj_br_ss_age ~ CWPM_N,  data = REDmasterPassageNbelow75)
summary(PN75)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N, data = REDmasterPassageNbelow75)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -16.0485  -5.2016  -0.2306   3.9687  13.2246 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)   
## (Intercept)  50.0223    18.6684   2.680  0.01221 * 
## CWPM_N        0.5228     0.1688   3.098  0.00441 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.392 on 28 degrees of freedom
## Multiple R-squared:  0.2552, Adjusted R-squared:  0.2286 
## F-statistic: 9.596 on 1 and 28 DF,  p-value: 0.004405
MN75 <- lm(rc3_wj_br_ss_age~ CWPM_N + DeltaDSyM_N_Sum ,  data = REDmasterPassageNbelow75)
summary(MN75)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N + DeltaDSyM_N_Sum, data = REDmasterPassageNbelow75)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -12.1595  -4.3262   0.2361   4.8705  11.7125 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)   
## (Intercept)      62.2253    17.7043   3.515  0.00157 **
## CWPM_N            0.4528     0.1565   2.893  0.00746 **
## DeltaDSyM_N_Sum  -0.7340     0.2865  -2.562  0.01631 * 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 6.752 on 27 degrees of freedom
## Multiple R-squared:  0.4009, Adjusted R-squared:  0.3565 
## F-statistic: 9.032 on 2 and 27 DF,  p-value: 0.0009923
anova(PN75, MN75)

above 75%ile narrative

detach(REDmasterPassageNbelow75)
REDmasterPassageNabove75 <- read.csv("C:/Users/pickrese/Dropbox/Vanderbilt/Cutting Lab/Reading Error Measure (RED)/PCRC/Poster Data/Wide Data Y2/REDmasterPassageNabove75.csv")
View(REDmasterPassageNabove75)
attach(REDmasterPassageNabove75)

PNa75 <- lm(rc3_wj_br_ss_age ~ CWPM_N,  data = REDmasterPassageNabove75)
summary(PNa75)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N, data = REDmasterPassageNabove75)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -19.6198  -6.7506  -0.4346   7.2981  19.4543 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 96.12441    8.85411  10.856 8.69e-16 ***
## CWPM_N       0.13426    0.05633   2.384   0.0203 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 8.67 on 60 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0865, Adjusted R-squared:  0.07127 
## F-statistic: 5.681 on 1 and 60 DF,  p-value: 0.02033
MNa75 <- lm(rc3_wj_br_ss_age~ CWPM_N + DeltaDSyM_N_Sum ,  data = REDmasterPassageNabove75)
summary(MNa75)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N + DeltaDSyM_N_Sum, data = REDmasterPassageNabove75)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -20.3453  -6.2957  -0.9518   7.3035  19.5774 
## 
## Coefficients:
##                  Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     104.22885    9.54550  10.919 8.61e-16 ***
## CWPM_N            0.10122    0.05741   1.763   0.0830 .  
## DeltaDSyM_N_Sum  -0.59790    0.29919  -1.998   0.0503 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 8.462 on 59 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.1444, Adjusted R-squared:  0.1154 
## F-statistic: 4.979 on 2 and 59 DF,  p-value: 0.01004
anova(PNa75, MNa75)

between 75-90 narrative

detach(REDmasterPassageNabove75)
REDmasterPassageNbelow90 <- read.csv("C:/Users/pickrese/Dropbox/Vanderbilt/Cutting Lab/Reading Error Measure (RED)/PCRC/Poster Data/Wide Data Y2/REDmasterPassageNbelow90.csv")
View(REDmasterPassageNbelow90)
attach(REDmasterPassageNbelow90)

PN90 <- lm(rc3_wj_br_ss_age ~ CWPM_N,  data = REDmasterPassageNbelow90)
summary(PN90)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N, data = REDmasterPassageNbelow90)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -15.5748  -5.5240  -0.9251   5.0523  16.3744 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)  
## (Intercept)  66.8857    36.1657   1.849   0.0763 .
## CWPM_N        0.3503     0.2649   1.323   0.1980  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 8.428 on 25 degrees of freedom
## Multiple R-squared:  0.06539,    Adjusted R-squared:  0.028 
## F-statistic: 1.749 on 1 and 25 DF,  p-value: 0.198
MN90 <- lm(rc3_wj_br_ss_age~ CWPM_N + DeltaDSyM_N_Sum ,  data = REDmasterPassageNbelow90)
summary(MN90)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N + DeltaDSyM_N_Sum, data = REDmasterPassageNbelow90)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -12.6008  -4.4299  -0.8596   6.1221  17.4554 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)  
## (Intercept)      63.6176    35.1690   1.809    0.083 .
## CWPM_N            0.4009     0.2591   1.547    0.135  
## DeltaDSyM_N_Sum  -0.5779     0.3635  -1.590    0.125  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 8.182 on 24 degrees of freedom
## Multiple R-squared:  0.1544, Adjusted R-squared:  0.08397 
## F-statistic: 2.192 on 2 and 24 DF,  p-value: 0.1336
anova(PN90, MN90)

above 90 narrative

detach(REDmasterPassageNbelow90)
REDmasterPassageNabove90 <- read.csv("C:/Users/pickrese/Dropbox/Vanderbilt/Cutting Lab/Reading Error Measure (RED)/PCRC/Poster Data/Wide Data Y2/REDmasterPassageNabove90.csv")
View(REDmasterPassageNabove90)
attach(REDmasterPassageNabove90)

PNa90 <- lm(rc3_wj_br_ss_age ~ CWPM_N,  data = REDmasterPassageNabove90)
summary(PNa90)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N, data = REDmasterPassageNabove90)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -19.4390  -6.5740  -0.7089   7.6941  19.6257 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  96.1488    23.3595   4.116 0.000242 ***
## CWPM_N        0.1331     0.1363   0.977 0.335847    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 9.02 on 33 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.02809,    Adjusted R-squared:  -0.001359 
## F-statistic: 0.9538 on 1 and 33 DF,  p-value: 0.3358
MNa90 <- lm(rc3_wj_br_ss_age~ CWPM_N + DeltaDSyM_N_Sum ,  data = REDmasterPassageNabove90)
summary(MNa90)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_N + DeltaDSyM_N_Sum, data = REDmasterPassageNabove90)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -20.648  -5.102  -1.365   7.029  17.663 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     103.2535    23.3729   4.418 0.000107 ***
## CWPM_N            0.1103     0.1344   0.821 0.417849    
## DeltaDSyM_N_Sum  -0.8280     0.5429  -1.525 0.137074    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 8.844 on 32 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.09394,    Adjusted R-squared:  0.03732 
## F-statistic: 1.659 on 2 and 32 DF,  p-value: 0.2063
anova(PNa90, MNa90)

above 90 expository

detach(REDmasterPassageNabove90)
REDmasterPassage90above <- read.csv("C:/Users/pickrese/Dropbox/Vanderbilt/Cutting Lab/Reading Error Measure (RED)/PCRC/Poster Data/Wide Data Y2/REDmasterPassage90above.csv")
View(REDmasterPassage90above)
attach(REDmasterPassage90above)

PNe90 <- lm(rc3_wj_br_ss_age ~ CWPM_E,  data = REDmasterPassage90above)
summary(PNe90)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_E, data = REDmasterPassage90above)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -20.188  -7.204  -0.364   6.982  16.588 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 125.96223   30.34386   4.151 0.000315 ***
## CWPM_E       -0.03204    0.21080  -0.152 0.880356    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 9.057 on 26 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0008879,  Adjusted R-squared:  -0.03754 
## F-statistic: 0.02311 on 1 and 26 DF,  p-value: 0.8804
MNe90 <- lm(rc3_wj_br_ss_age~ CWPM_E + DeltaDSyM_E_Sum ,  data = REDmasterPassage90above)
summary(MNe90)
## 
## Call:
## lm(formula = rc3_wj_br_ss_age ~ CWPM_E + DeltaDSyM_E_Sum, data = REDmasterPassage90above)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -18.8046  -5.6736   0.1113   6.6127  14.5532 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     142.8161    35.5714   4.015 0.000477 ***
## CWPM_E           -0.1325     0.2382  -0.556 0.583083    
## DeltaDSyM_E_Sum  -0.9539     1.0418  -0.916 0.368613    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 9.085 on 25 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.03331,    Adjusted R-squared:  -0.04403 
## F-statistic: 0.4307 on 2 and 25 DF,  p-value: 0.6548
anova(PNe90, MNe90)