# Figure 3
# assumes data Posner2006_2.rda have been loaded
# jochen

quartz(width=11.13,  height=7.35)
theme_set(theme_bw())
theme_update(panel.grid.major = theme_line(colour="white"),  panel.grid.minor = theme_line(colour="white"), 
plot.margin = unit(c(0, 0, 0, 0)+.01, "lines"))
CV1<-relevel(d$CV, ref="neutral")
d$CV2<-relevel(CV1, ref="valid")
# ... main effect: MTC x TYP
m <- ggplot(d, aes(x=mt, group=MTC, colour=MTC, subset=EXP=="manual RT"))  + scale_colour_grey(end=0.7,start=0,name='')
m <- m + facet_grid( CV2 ~ TYP )  + scale_x_continuous(limits=c(0,1000), breaks=seq(from=0,to=1000,200)) 
m <- m + geom_density(aes(x=mt), kernel = "gaussian", adjust=1) + xlab("Time from cue [ms]") + opts(title = "manual responses")
# m <- m + geom_density(kernel="rectangular",  adjust=.5,  aes(x=mt)) 
# m <- m + stat_bin(aes(x=mt, y = ..density..), binwidth=30,  geom="line", position="identity")
m


quartz(width=11.13,  height=7.35)
theme_set(theme_bw())
theme_update(panel.grid.major = theme_line(colour="white"),  panel.grid.minor = theme_line(colour="white"), 
plot.margin = unit(c(0, 0, 0, 0)+.01, "lines"))
CV1<-relevel(d$CV, ref="neutral")
d$CV2<-relevel(CV1, ref="valid")
# ... main effect: MTC x TYP
m <- ggplot(d, aes(x=mt, group=MTC, colour=MTC, subset=EXP=="saccadic RT"))  + scale_colour_grey(end=0.7,start=0,name='')
m <- m + facet_grid( CV2 ~ TYP )  + scale_x_continuous(limits=c(0,1000), breaks=seq(from=0,to=1000,200)) 
m <- m + geom_density(aes(x=mt), kernel = "gaussian", adjust=1) + xlab("Time from cue [ms]") + opts(title = "saccadic responses")
# m <- m + geom_density(kernel="rectangular",  adjust=.5,  aes(x=mt)) 
# m <- m + stat_bin(aes(x=mt, y = ..density..), binwidth=30,  geom="line", position="identity")
m


quartz(width=11.13,  height=7.35)
theme_set(theme_bw())
theme_update(panel.grid.major = theme_line(colour="white"),  panel.grid.minor = theme_line(colour="white"), 
plot.margin = unit(c(0, 0, 0, 0)+.01, "lines"))
CV1<-relevel(d$CV, ref="neutral")
d$CV2<-relevel(CV1, ref="valid")
d$EXPMTC <- as.factor(10*as.numeric(d$EXP) + as.numeric(d$MTC))
levels(d$EXPMTC) <- c("manual incongruent", "manual congruent", "saccadic incongruent", "saccadic congruent")
# ... main effect: MTC x TYP
m <- ggplot(d, aes(x=mt, group=EXPMTC, colour=MTC, size=EXP)) 
m <- m + facet_grid( CV2 ~ TYP )  + scale_x_continuous(limits=c(0,1000), breaks=seq(from=0,to=1000,200)) 
# m <- m + geom_density(aes(x=mt), kernel = "gaussian", adjust=1, size=1)
m <- m + geom_density(aes(x=mt), kernel = "gaussian", adjust=1) + scale_y_continuous('Microsaccade rate') 
m <- m + scale_size_manual(name="Response", values = c(1.3, 1.7))
m <- m + scale_colour_grey(end=0.5,start=0,name='MS-Target\nCongruency') 
m <- m + xlab("Time from cue [ms]") + ylab("Microsaccade rate")
# m <- m + geom_density(kernel="rectangular",  adjust=.5,  aes(x=mt)) 
# m <- m + stat_bin(aes(x=mt, y = ..density..), binwidth=30,  geom="line", position="identity")
m