Question: Proper dataset form and model for DEXeq analysis
gravatar for heikki.sarin
2.2 years ago by
heikki.sarin0 wrote:


I have a bit of problem figuring out the proper way to set dataset and model for DEXeq analysis. We have run gene-wise analysis of RNA-seq with DESeq2 but we also have exon-data that we're trying to analyse with DEXeq.

1) Is it enough for the analysis if I have:

    i) Sample data with covariates and additional exon column with "this" and "other"

    ii) Countdata with columns for samples that have counts per exon and the sum columns for sums for other exons of the gene per sample. Rows are named ENSG3923903:E001 etc. I created the count matrix manually by calculating the sum columns and then downloading the annotation from ensemble so I could combine it with my data and construct the rownames manually.

2) Couple of questions also regarding the model. In deseq2 we created model.matrices for full and reduced model.

    i) is it possible to use model.matrices with DEXeq?

    ii) We used formulas :


           reduced = CASE_CONTROL + PRE_MID_POST+ CASE_CONTROL:subject.nested

        How can we apply the  "exon" variable for the models and run the analysis on DEXeq? We're trying to test whether there DE can be detected between two condition across three time points so that we're accounting for the subject variability.

Help would be really appreciated. (Especially from Alejandro Reyes, Simon Anders, Wolfgang Huber.)



ADD COMMENTlink modified 2.2 years ago by Michael Love25k • written 2.2 years ago by heikki.sarin0

Hi Heikki, could you include the output of colData(dxd) ) ?


ADD REPLYlink written 2.2 years ago by Alejandro Reyes1.7k
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: 146 users visited in the last hour