Question: JunctionSeq pairwise contrast with a multi-level treatment
gravatar for stuart
3.1 years ago by
stuart0 wrote:


Does anyone know how to contrast only a single pair of conditions in JunctionSeq (which is based on DEXseq AFAIK), if there are > 2 levels of the condition? I would like to test the null hypothesis that there is no condition:bin effect between one pair of conditions, but still using information from the untested conditions to better estimate the dispersion.

For example, in the JunctionSeq manual pg 40 they give code for testing a single, 3-level condition in the example data:

threeLevelVariable <- c("GroupA","GroupA",

jscs <- runJunctionSeqAnalyses(sample.files = countFiles,
sample.names = decoder$sample.ID,
flat.gff.file = gff.file,
nCores = 1,
analysis.type = "junctionsAndExons"

How could I request a Group A vs Group B contrast, without throwing away the information on dispersion contained in the data for Group C?

Thank you


ADD COMMENTlink written 3.1 years ago by stuart0
Please log in to add an answer.


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