**0**wrote:

I have an RNA-seq dataset with 3 treatment groups, where gene g has group means mu_{g, i} for i = 1, 2, 3. Is there a way to use baySeq to calculate the posterior probabilities that mu_{g, 2} > max(mu_{g, 1}, mu_{g, 3}) and the posterior probabilities that mu_{g, 2} < min(mu_{g, 1}, mu_{g, 3})? Specifications like ** groups <- list(null = c(1,1,1,1,1,1), alternative = c(1,1,2,2,1,1))** and

**do not answer this question with the package's usual functionality demonstrated in the manual and the vignette. If not with baySeq, which packages can accomplish this easily?**

`groups <- list(null = c(1,1,1,1,1,1), alternative = c(1,1,2,2,3,3))`

**38k**• written 4.8 years ago by Will Landau •

**0**