
\newpage


\section{Plots of RTs by condition, subject, and item}

\subsection{All the subjects}

\begin{Schunk}
\begin{Sinput}
> plot.data <- d.melt[d.melt$roi <= 11, ]
> reorder.frame <- data.frame(condition = c(rep("a", 11), rep("b", 11), rep("c", 11), rep("d", 
+     11), rep("e", 11), rep("f", 11)), roi = rep(1:11, 6), plot.order = c(1, 5, 11, NA, NA, 
+     NA, NA, NA, NA, NA, NA, 1, 5, 8, 9, 10, 11, NA, NA, NA, NA, NA, 1, 5, 6, 7, 8, 9, 10, 
+     11, NA, NA, NA, 1, 2, 3, 4, 5, 11, NA, NA, NA, NA, NA, 1, 2, 3, 4, 5, 8, 9, 10, 11, 
+     NA, NA, 1:11))
> plot.data <- merge(plot.data, reorder.frame)
\end{Sinput}
\end{Schunk}

\begin{Schunk}
\begin{Sinput}
> trim.value <- 60000
> rt.mean <- cast(plot.data, condition + plot.order ~ ., function(x) c(m = round(mean(x)), 
+     CI = round(ci(x))), subset = ((value < trim.value)))
> rt.mean <- rt.mean[!(rt.mean$plot.order == "NA"), ]
> rt.mean <- na.delete(rt.mean)
> dummy <- data.frame(region = c(1:11), m = rep(-999, 11))
> subtitle <- "Error bars are one standard error."
\end{Sinput}
\end{Schunk}



\begin{Schunk}
\begin{Sinput}
> color1 <- gray(0.7)
> color2 <- gray(0.3)
> color3 <- "black"
> lty1 <- "dashed"
> lty2 <- "solid"
\end{Sinput}
\end{Schunk}



\begin{center}
\begin{Schunk}
\begin{Sinput}
> plot(m ~ region, data = dummy, xaxt = "n", type = "n", cex.axis = 1.2, cex.lab = 1.5, cex.main = 1.5, 
+     cex.sub = 0.7, xlab = "Region", ylab = "Reading time (ms)", main = "Experiment 3 self-paced reading times", 
+     sub = subtitle, ylim = rtlimits)
> axis(side = 1, at = 1:11, labels = FALSE)
> mtext(side = 1, at = c(1:11), cex = 0.7, line = 1, c("The", "(sports", "that", "the)", "child", 
+     "(who", "was", "(from", "the", "school))", "played"))
> legend(x = rt.mean$plot.order[1], y = 1200, legend = c("matrix - unmodified", "PP", "RC", 
+     "embedded - unmodified", "PP", "RC"), lwd = c(2), lty = c(lty1, lty1, lty1, lty2, lty2, 
+     lty2), col = c(color1, color2, color3, color1, color2, color3), pch = c(1, 1, 1, 15, 
+     15, 15), bty = "n", cex = 0.9)
> plot.intervals(xrange = rep(11, 6), lower = rt.mean$CI.lower[rt.mean$plot.order == 11], 
+     upper = rt.mean$CI.upper[rt.mean$plot.order == 11], width = 0.05, lwd = 1, col = "black", 
+     lty = "solid")
> condition <- "a"
> point.type <- c(1)
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, cex = 1.2, col = color1)
> condition <- "b"
> dom1 <- 8:11
> range1 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order > 7]
> color <- color2
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, col = color, cex = 1.2)
> lines(x = dom1, y = range1, lty = lty1, col = color, pch = point.type, lwd = 1, type = "o")
> condition <- "c"
> dom1 <- 5:11
> range1 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order > 4]
> color <- color3
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, col = color, cex = 1.2)
> lines(x = dom1, y = range1, lty = lty1, col = color, pch = point.type, lwd = 1, type = "o")
> condition <- "d"
> dom1 <- 1:5
> range1 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order < 6]
> point.type <- c(15)
> color <- color1
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, col = color, cex = 1.2)
> lines(x = dom1, y = range1, lty = lty2, col = color, pch = point.type, lwd = 1, type = "o")
> condition <- "e"
> dom1 <- 1:5
> dom2 <- 8:11
> range1 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order < 6]
> range2 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order > 7]
> color <- color2
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, col = color, cex = 1.2)
> lines(x = dom1, y = range1, lty = lty2, col = color, pch = point.type, lwd = 1, type = "o")
> lines(x = dom2, y = range2, lty = lty2, col = color, pch = point.type, lwd = 1, type = "o")
> condition <- "f"
> dom1 <- 1:11
> range1 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order]
> color <- color3
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, col = color, cex = 1.2)
> lines(x = dom1, y = range1, lty = lty2, col = color, pch = point.type, lwd = 1, type = "o")
\end{Sinput}
\end{Schunk}
\includegraphics{all-roi-bb-plots-all-roi-bb-spr-by-condition}

\end{center}




\newpage

\subsection{Accurate subjects only}

Now lets look at the RT pattern for the accurate subjects (defined above).

\begin{Schunk}
\begin{Sinput}
> trim.value <- 60000
> rt.mean <- cast(plot.data, condition + plot.order ~ ., function(x) c(m = round(mean(x)), 
+     CI = round(ci(x))), subset = ((value < trim.value) & (good.subj == TRUE)))
> rt.mean <- rt.mean[!(rt.mean$plot.order == "NA"), ]
> rt.mean <- rt.mean[-c(4, 11, 20, 27, 37), ]
\end{Sinput}
\end{Schunk}



\begin{center}
\begin{Schunk}
\begin{Sinput}
> plot(m ~ region, data = dummy, xaxt = "n", type = "n", cex.axis = 1.2, cex.lab = 1.5, cex.main = 1.5, 
+     cex.sub = 0.7, xlab = "Region", ylab = "Reading time (ms)", main = "Experiment 3 self-paced reading times", 
+     sub = subtitle, ylim = rtlimits)
> axis(side = 1, at = 1:11, labels = FALSE)
> mtext(side = 1, at = c(1:11), cex = 0.7, line = 1, c("The", "(sports", "that", "the)", "child", 
+     "(who", "was", "(from", "the", "school))", "played"))
> legend(x = rt.mean$plot.order[1], y = 1200, legend = c("matrix - unmodified", "PP", "RC", 
+     "embedded - unmodified", "PP", "RC"), lwd = c(2), lty = c(lty1, lty1, lty1, lty2, lty2, 
+     lty2), col = c(color1, color2, color3, color1, color2, color3), pch = c(1, 1, 1, 15, 
+     15, 15), bty = "n", cex = 0.9)
> plot.intervals(xrange = rep(11, 6), lower = rt.mean$CI.lower[rt.mean$plot.order == 11], 
+     upper = rt.mean$CI.upper[rt.mean$plot.order == 11], width = 0.05, lwd = 1, col = "black", 
+     lty = "solid")
> condition <- "a"
> point.type <- c(1)
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, cex = 1.2, col = color1)
> condition <- "b"
> dom1 <- 8:11
> range1 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order > 7]
> color <- color2
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, col = color, cex = 1.2)
> lines(x = dom1, y = range1, lty = lty1, col = color, pch = point.type, lwd = 1, type = "o")
> condition <- "c"
> dom1 <- 5:11
> range1 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order > 4]
> color <- color3
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, col = color, cex = 1.2)
> lines(x = dom1, y = range1, lty = lty1, col = color, pch = point.type, lwd = 1, type = "o")
> condition <- "d"
> dom1 <- 1:5
> range1 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order < 6]
> point.type <- c(15)
> color <- color1
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, col = color, cex = 1.2)
> lines(x = dom1, y = range1, lty = lty2, col = color, pch = point.type, lwd = 1, type = "o")
> condition <- "e"
> dom1 <- 1:5
> dom2 <- 8:11
> range1 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order < 6]
> range2 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order > 7]
> color <- color2
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, col = color, cex = 1.2)
> lines(x = dom1, y = range1, lty = lty2, col = color, pch = point.type, lwd = 1, type = "o")
> lines(x = dom2, y = range2, lty = "solid", col = color, pch = point.type, lwd = 1, type = "o")
> condition <- "f"
> dom1 <- 1:11
> range1 <- rt.mean$m[rt.mean$condition == condition & rt.mean$plot.order]
> color <- color3
> points(rt.mean$plot.order[rt.mean$condition == condition], rt.mean$m[rt.mean$condition == 
+     condition], pch = point.type, col = color, cex = 1.2)
> lines(x = dom1, y = range1, lty = lty2, col = color, pch = point.type, lwd = 1, type = "o")
\end{Sinput}
\end{Schunk}
\includegraphics{all-roi-bb-plots-all-roi-bb-spr-by-condition-good-subj-only}

\end{center}
