Search
Question: Multiple pairwise comparisons in ballgown
0
gravatar for rachael.detar
5 months ago by
rachael.detar0 wrote:

Hi, 
I have some RNA seq data with 3 genotypes and two treatments. I'd like to do pairwise comparisons between every group. Is there an easy way to do this in ballgown, or do I need to subset my dataset into each combination of genotype/treatment, and run the ballgown pipeline on them pair by pair? What I really want as an output is a table of log2 fold change values for each pairwise comparison. 

The command where I think this may be relevent is pasted below:

> results_transcripts = stattest(bg_chrX_filt, feature="transcript", covariate="genotype_salt_mM", getFC=TRUE, meas="FPKM")
Warning message:
In stattest(bg_chrX_filt, feature = "transcript", covariate = "genotype_salt_mM",  :
  fold changes only available for 2-group comparisons

Thanks,

Rachael DeTar

ADD COMMENTlink modified 5 months ago by Alyssa Frazee170 • written 5 months ago by rachael.detar0
0
gravatar for Alyssa Frazee
5 months ago by
Alyssa Frazee170
San Francisco, CA, USA
Alyssa Frazee170 wrote:

Hi Rachael, ballgown doesn't do pairwise comparisons by default, but there are a few ways to force it to do it for you. The strategy you mention (running the pipeline on them pair by pair) is one option (you can use the "subset" function to make 6 ballgown objects, one with each genotype x pair treatment, then run "stattest" with each of the 6 objects). This is probably the easiest way to get what you want -- you'll have a two-group comparison and then you can pull the fold change values directly from the stattest function. (There are a couple other strategies I can think of as well, but they aren't as straightforward as the above and if you do want the 6 fold change values, there's some annoying data wrangling you'll have to do yourself. Ballgown doesn't automatically emit coefficients for more complex designs like yours because it's ~impossible to guess which values are important to the scientist in complex designs. But if this doesn't work, feel free to reply and folks can help brainstorm better solutions). 

ADD COMMENTlink written 5 months ago by Alyssa Frazee170
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 150 users visited in the last hour