# In order to provide some perspective regarding the ranges of
# parameter values that may be expected when fitting the full diffusion
# model to data, this file provides an overview of the best--fitting 
# parameter values encountered in
# previous research

# We use "R" to plot the results.

## Criteria for incorporating results in review ##
# 1. Whenever there is a choice, use the fit to average data,
#    instead of individual Ss.
# 2. No explicit manipulation of starting point.
# 3. Negative drift rates multiplied by -1.
##################################################
## Studies ##
# 1. Ratcliff & Smith 2004:	Table 1, 2, 5, 6, and 8.
# 2. Van Zandt et al., 2000: Table 8 in Ratcliff & Smith, 2004.
# 3. Ratcliff, Thapar, Gomez, & McKoon, 2004: Table 3 and 4.
# 4. Ratcliff, Gomez, & McKoon, 2004: Table 6 and 10.
# 5. Ratcliff, Thapar, & McKoon, 2001: Table 2.
# 6. Ratcliff & Rouder, 2000: Table 2 and 3
# 7. Ratcliff, Van Zandt, & McKoon, 1999, Table 1.
# 8. Voss et al., 2004: Table 2 and 4.
# 9. Ratcliff, Thapar, & McKoon, 2004: Table 3 and 4.
#############

v.1.1   = c(0.0391, 0.1320, 0.1944, 0.3208)
a.1.2   = c(0.0821, 0.1440, 0.0855, 0.1593, 0.1139)
z.1.2   = c(0.0821/2, 0.1440/2, 0.0410, 0.0774, 0.0550)
Ter.1.2 = c(0.3109, 0.4082, 0.5222)
eta.1.2 = c(0.1475, 0.1135, 0.1658)
sz.1.2  = c(0.0324, 0.0531, 0.0393)
st.1.2  = c(0.1000, 0.1430, 0.2140)

v.1.5   = c(0.4316, 0.2749, 0.1818, 0.2566)
v.1.6   = c(0.2132, 0.3015, 0.3228, 0.1009, 0.1430, 0.1509, 0.0890, 0.1484, 0.1675)

v.1.8   = c(0.2161, 0.3659)
a.1.8   = c(0.1017)
z.1.8   = c(0.0589)
Ter.1.8 = c(0.3646)
eta.1.8 = c(0.0252)
sz.1.8  = c(0.0168)
st.1.8  = c(NA)

v.3.4   = c(0.481, 0.281, 0.180, 0.279, 0.435, 0.287, 0.200, 0.268, 0.437, 0.356, 0.300, 0.321, 0.465, 0.367, 0.296, 0.363)
a.3.3   = c(0.126, 0.184, 0.124, 0.166)
z.3.3   = c(0.062, 0.091, 0.069, 0.090)
Ter.3.3 = c(0.447, 0.548, 0.398, 0.480)
eta.3.3 = c(0.120, 0.097, 0.084, 0.111)
sz.3.3  = c(0.069, 0.064, 0.065, 0.042)
st.3.3  = c(0.154, 0.177, 0.141, 0.150)

v.4.6   = c(0.213, 0.396, 0.216, 0.128, 0.358, 0.477, 0.368, 0.312, 0.235, 0.337, 0.151, 0.350, 0.392, 0.274, 0.198, 0.368, 0.183, 0.104, 0.215, 0.360, 0.138, 0.401, 0.164, 0.390, 0.255, 0.158)
a.4.6   = c(0.110, 0.125, 0.112, 0.120, 0.127, 0.130, 0.131)
z.4.6   = c(0.056, 0.066, 0.056, 0.063, 0.063, 0.060, 0.063)
Ter.4.6 = c(0.435, 0.407, 0.433, 0.394, 0.438, 0.466, 0.402)
eta.4.6 = c(0.070, 0.123, 0.014, 0.101, 0.035, 0.081, 0.086)
sz.4.6  = c(0.004, 0.076, 0.054, 0.072, 0.046, 0.004, 0.025)
st.4.6  = c(0.159, 0.133, 0.164, 0.110, 0.179, 0.171, 0.146)

v.4.10   = c(0.218, 0.334, 0.181, 0.109, 0.203, 0.384, 0.260, 0.195, 0.320, 0.394, 0.325, 0.267, 0.320, 0.408, 0.355, 0.346)
a.4.10   = c(0.115, 0.115, 0.117, 0.117)
z.4.10   = c(0.057, 0.057, 0.058, 0.058)
Ter.4.10 = c(0.441, 0.441, 0.389, 0.389)
eta.4.10 = c(0.066, 0.066, 0.076, 0.076)
sz.4.10  = c(0.024, 0.024, 0.061, 0.061)
st.4.10  = c(0.003, 0.003, 0.122, 0.122)

v.5.2   = c(0.3507, 0.2727, 0.3152, 0.3008, 0.3302, 0.3302, 0.3300, 0.3300, 0.2040, 0.1472, 0.1950, 0.1801, 0.2063, 0.2063, 0.2050, 0.2050, 0.1152, 0.0861, 0.1188, 0.1163, 0.1357, 0.1357, 0.1077, 0.1077, 0.0425, 0.0037, 0.0399, 0.0410, 0.0364, 0.0364, 0.0425, 0.0425)
a.5.2   = c(0.1705, 0.1413, 0.1475, 0.1416, 0.1805, 0.1096, 0.1483, 0.0914)
z.5.2   = a.5.2/2
Ter.5.2 = c(0.3938, 0.3409, 0.3459, 0.3549, 0.3234, 0.3234, 0.2846, 0.2846)
eta.5.2 = c(0.0928, 0.0846, 0.0992, 0.0959, 0.1081, 0.1081, 0.1040, 0.1040)
sz.5.2  = c(0.0365, 0.0513, 0.0458, 0.0561, 0.0453, 0.0453, 0.0396, 0.0396)
st.5.2  = c(NA, NA, NA, NA, NA, NA, NA, NA)

v.6.3   = c(.026, .067, .124, .216, .298, .024, .080, .224, .274, .311)
a.6.2   = c(.108, .104)
z.6.2   = c(.054, .052)
Ter.6.2 = c(0.316, 0.319, 0.303)
eta.6.2 = c(.111, .071, .170)
sz.6.2  = c(.018, .014)
st.6.2  = c(NA, NA)

v.7.1   = c(NA, NA, NA, NA, NA)
a.7.1   = c(.115, .151, .141, .150, .065)
z.7.1   = a.7.1/2
Ter.7.1 = c(.256, .323, .335, .313, .266)
eta.7.1 = c(.112, .088, .174, .142, .055)
sz.7.1  = c(NA, NA, NA, NA, NA)
st.7.1  = c(NA, NA)

v.8.2   = c(0.089, 0.046, 0.091, 0.102, 0.103, 0.057, 0.118, 0.141)
a.8.2   = c(0.118, 0.119, 0.137, 0.116)
z.8.2   = c(0.49*0.118, 0.47*0.119, 0.52*0.137, 0.56*0.116)
Ter.8.2 = c(0.470, 0.470, 0.520, 0.730)
eta.8.2 = c(.080, .090, .085, .089)
sz.8.2  = c(0.018, 0.018, 0.018, 0.016)
st.8.2  = c(NA, NA, NA, NA)

v.9.4   = c(0.206, 0.215, 0.243, 0.015, 0.095, 0.109, 0.166, 0.256, 0.297, 0.163, 0.252, 0.246, 0.003, 0.114, 0.150, 0.125, 0.246, 0.274)
a.9.3   = c(.1, .138, .069, .132)
z.9.3   = c(.045, .063, .035, .058)
Ter.9.3 = c(0.594, 0.486)
eta.9.3 = c(0.218, 0.164)
sz.9.3  = c(0.006, 0.011)
st.9.3  = c(0.232, 0.203)

v   = c(v.1.1,v.1.5,v.1.6,v.1.8,v.3.4,v.4.6,v.4.10,v.5.2,v.6.3,v.7.1,v.8.2,v.9.4)
a   = c(a.1.2,a.1.8,a.3.3,a.4.6,a.4.10,a.5.2,a.6.2,a.7.1,a.8.2,a.9.3)
z   = c(z.1.2,z.1.8,z.3.3,z.4.6,z.4.10,z.5.2,z.6.2,z.7.1,z.8.2,z.9.3)
Ter = c(Ter.1.2,Ter.1.8,Ter.3.3,Ter.4.6,Ter.4.10,Ter.5.2,Ter.6.2,Ter.7.1,Ter.8.2, Ter.9.3)
eta = c(eta.1.2,eta.1.8,eta.3.3,eta.4.6,eta.4.10,eta.5.2,eta.6.2,eta.7.1,eta.8.2,eta.9.3)
sz  = c(sz.1.2,sz.1.8,sz.3.3,sz.4.6,sz.4.10,sz.5.2,sz.6.2,sz.7.1,sz.8.2,sz.9.3)
st  = c(st.1.2,st.1.8,st.3.3,st.4.6,st.4.10,st.5.2,st.6.2,st.7.1,st.8.2,st.9.3)

layout(matrix(1:6, 2, byrow = T))
################ v ############################################################
par(cex.main = 1.5, mar = c(5, 5, 4, 2) + 0.1, mgp = c(3, 1, 0), cex.lab = 1.5,
    font.lab = 2, cex.axis = 1.3, bty = "O", las=1)
# mpg = c(3, 1, 0) is default. first = axis labels!; middle = tick labels  
# mar = c(5, 4, 4, 2) + 0.1 is default

hist(v, freq=T, lwd = 2,
	  main = "Drift Rate", ylab = "Frequency", xlab = "v")
rug(v, side = 1, quiet = 1, col = "black")
#length(which(is.finite(v)))
text(0.20, 23, lab="N=145", cex=2, pos=4)
#################################################################################

################ a ############################################################
par(cex.main = 1.5, mar = c(5, 5, 4, 2) + 0.1, mgp = c(3, 1, 0), cex.lab = 1.5,
    font.lab = 2, cex.axis = 1.3, bty = "O", las=1)
# mpg = c(3, 1, 0) is default. first = axis labels!; middle = tick labels  
# mar = c(5, 4, 4, 2) + 0.1 is default

hist(a, freq=T, lwd = 2,
	  main = "Boundary Separation", ylab = "Frequency", xlab = "a")
rug(a, side = 1, quiet = 1, col = "black")
#length(which(is.finite(a)))
text(0.14, 15, lab="N=44", cex=2, pos=4)
#################################################################################

################ Ter ############################################################
par(cex.main = 1.5, mar = c(5, 5, 4, 2) + 0.1, mgp = c(3, 1, 0), cex.lab = 1.5,
    font.lab = 2, cex.axis = 1.3, bty = "O", las=1)
# mpg = c(3, 1, 0) is default. first = axis labels!; middle = tick labels  
# mar = c(5, 4, 4, 2) + 0.1 is default

hist(Ter, freq=T, lwd = 2,
	  main = "Mean of Nondecision Time", ylab = "Frequency", xlab = expression(T[er]))
rug(Ter, side = 1, quiet = 1, col = "black")
#length(which(is.finite(Ter)))
text(0.5, 10, lab="N=41", cex=2, pos=4)
#################################################################################

################ eta ############################################################
par(cex.main = 1.5, mar = c(5, 5, 4, 2) + 0.1, mgp = c(3, 1, 0), cex.lab = 1.5,
    font.lab = 2, cex.axis = 1.3, bty = "O", las=1)
# mpg = c(3, 1, 0) is default. first = axis labels!; middle = tick labels  
# mar = c(5, 4, 4, 2) + 0.1 is default

hist(eta, freq=T, lwd = 2,
	  main = "Trial-to-Trial Variability \n in Drift Rate", ylab = "Frequency", xlab = expression(eta))
rug(eta, side = 1, quiet = 1, col = "black")
#length(which(is.finite(eta)))
text(0.11, 15, lab="N=41", cex=2, pos=4)
#################################################################################

################ sz ############################################################
par(cex.main = 1.5, mar = c(5, 5, 4, 2) + 0.1, mgp = c(3, 1, 0), cex.lab = 1.5,
    font.lab = 2, cex.axis = 1.3, bty = "O", las=1)
# mpg = c(3, 1, 0) is default. first = axis labels!; middle = tick labels  
# mar = c(5, 4, 4, 2) + 0.1 is default

hist(sz, freq=T, lwd = 2,
	  main = "Trial-to-Trial Variability \n in Starting Point", ylab = "Frequency", xlab = expression(s[z]))
rug(sz, side = 1, quiet = 1, col = "black")
#length(which(is.finite(sz)))
text(0.04, 6, lab="N=35", cex=2, pos=4)
#################################################################################

################ st ############################################################
par(cex.main = 1.5, mar = c(5, 5, 4, 2) + 0.1, mgp = c(3, 1, 0), cex.lab = 1.5,
    font.lab = 2, cex.axis = 1.3, bty = "O", las=1)
# mpg = c(3, 1, 0) is default. first = axis labels!; middle = tick labels  
# mar = c(5, 4, 4, 2) + 0.1 is default

hist(st, freq=T, lwd = 2,
	  main = "Trial-to-Trial Variability \n in Nondecision Time", ylab = "Frequency", xlab = expression(s[t]))
rug(st, side = 1, quiet = 1, col = "black")
#length(which(is.finite(st)))
text(0.15, 7, lab="N=20", cex=2, pos=4)
#################################################################################

#### Figure showing how z is close to a/2:
par(cex.main = 1.5, mar = c(5, 6, 4, 2) + 0.1, mgp = c(3, 1, 0), cex.lab = 1.5,
    font.lab = 2, cex.axis = 1.3, bty = "O", las=1)
# mpg = c(3, 1, 0) is default. first = axis labels!; middle = tick labels  
# mar = c(5, 4, 4, 2) + 0.1 is default

plot(z,a, lwd = 2, ylim = c(.06,.20), xlim=c(.03, .10),
	  main = " ", ylab = " ", xlab = " ")
mtext("Boundary Separation", side=2, padj=-3.2, cex=1.8, las=0)
mtext("Starting Point", side=1, padj=2.2, cex=1.8, las=0)

abline(0, 2, lwd=2)
#length(which(is.finite(a)))
text(0.05, .18, lab="N=44", cex=2, pos=4)