I have both the Reaction Time (RT) and Accuracy Rate (ACC,range from 0 to 1) data for my 2*2 designed experiment. I want to know whether there is an interaction between the two factors. Two-way ANOVA could be used easily for my purpose. But I would like to calculate the bayes factor (BF) to further confirm the result. I use the same R code with MCMCregress function to do this. The RT variale seems to get compatible BF results as compared to the ANOVA. However, the BFs of ACC are pretty small (less than 1) for all my data, even when the interaction effects are very significant based on the ANOVA.The code goes as follows:

model1.ACC