dprobit <- design(Rlevels = c("left","right"),
factors=list(subjects=1,S=c("left","right")),
formula=list(mean ~ 0+S, sd ~ 1,threshold ~ 1),
matchfun=function(d)d$S==d$lR,
constants=c(sd=log(1),threshold=0),
model=SDT)
p_vector <- sampled_pars(dprobit)
Run the code above in your browser using DataLab