title: Normal Distribution bayesAB - negative mean values
body: <pre><code>A_norm &lt;- rnorm(300, 5, 2)
  B_norm &lt;- rnorm(300, -5, 2)
  
  AB2 &lt;- bayesTest(A_norm, B_norm,
               priors = c('mu' = 5, 'lambda' = 1, 'alpha' = 3, 'beta' = 1), 
               distribution = 'normal')
  
  plot(AB2)
  summary(AB2)
  </code></pre>
  
  <p>Using the above code, <code>summary(AB2)</code> shows that the probability of A > B is equal to 0, however when I look at the posterior plots of mu for A and B (<code>plot(AB2)</code>) - the distribution of A looks to be always greater than B with no overlap. I think it has to do with comparing two normal distributions with both a positive and negative mean respectively. However, I cant understand why this would affect it. Can someone help me understand?</p>