Diffusion Model Predictions for Accuracy and Reaction Time Quantiles
       and a Chi-square Goodness of Fit Value 
       for that One Condition (One Drift Rate)

     
     See how the reaction time quantiles and accuracy changes with
     diffusion model parameters and how goodness of fit changes
     for the predictions for the data entered.

	
     Parameters:

        v     a    z    Ter    η    Sz   St    P0    an
         

     Correct and Error Reaction Time Quantiles from Data:

               Prob.      Q.1     Q.3     Q.5     Q.7     Q.9    Min RT    Max RT
     Correct:         
       Error:         

    


     

        dyn.load("/var/www/html/cogsys-web/Rpad/corqout.so")
        corqout<-function(par,q,chi,out,n)
        .Fortran("corqout",as.double(par),as.double(q),numeric(1),numeric(12),as.double(n))

#        par<-cbind(.1,.4,.08,.02,.3,.01,.1,.05)
        par<-cbind(a,Ter,eta,Sz,v,P0,St,z)

        for (i in 1:8) {
          if (par[i]<.001) {
            print("  ALERT: values cannot be less than .001  ");
            par[i]=.001;
          }
        }
	
#        q<-cbind(0.936,438.,482.,522.,580.,710.,300.,2000.,0.064,440.,487.,532.,599.,753.,300.,2000.)
        q<-cbind(q1,q2,q3,q4,q5,q6,q7,q8,q9,q10,q11,q12,q13,q14,q15,q16)
        z<-corqout(par,q,chi,out,an)
        chi<-format(z[[3]],digits=3)
#        out<-format(z[[4]],digits=2)
        out<-round(z[[4]]*1000)/1000

HTMLon()
H("div",BR)
H("pre","               Cumulative probabilities at the RT quantiles")
H("pre","          Prob     Q.1     Q.3     Q.5     Q.7     Q.9")
H("pre","Correct:",out[1]," ",out[2]," ",out[3]," ",out[4]," ",out[5]," ",out[6])
H("pre","Error:  ",out[7]," ",out[8]," ",out[9]," ",out[10]," ",out[11]," ",out[12])
H("pre","chi-square (from both correct and error RTs):",chi)
HTMLoff()


v = drift rate; a = boundary separation; z = starting point; Ter = nondecision component of response time;
η = standard deviation in drift across trials; Sz range of the distribution of starting point (z);
St = range of the distribution of nondecision times; P0 = proportion of contaminants;
an = number of observations.
             Copyright & Programmed 2006 Ratcliff & McKoon Lab | Designed by Andrew Webb


// this makes all Rpad updates to output boxes "yellowfade" in
dojo.event.connect(rpad, "updateResults", 
function (rpadResults, data) {
       	dojo.fx.html.colorFade(rpadResults, 1500,
       	dojo.graphics.color.extractRGB("#FFFF70"),  // from color - bright yellow
       	dojo.graphics.color.extractRGB("#FFFFD0")); // to color
});