Hi Gang,
Thank you so much. I coded the model as you specified with the 3 groups and was able to a pillot test on a subset of my sample.
Now I'm working on extracting effect estimates, summary info, plots, etc.
As I was going through the MBA code, I noticed that when summarizing the information for the posterior samples of Subject pairs, you apply some sort of standardization to all values based on the mean of all samples for that subject and the square root of 2, per below:
for(ii in 1:nSubjs) for(jj in 1:nSubjs) ps[,ii,jj] <- sqrt(2)*(ps[,ii,jj] - mm[ii,jj]) + mm[ii,jj]
I am wondering why there was a need to do this, as I have never come across this formula.
Best,
Vasco