Using Expression of a Gene as a Linear Regression Variable
Entering edit mode
4 months ago
jamie3355 • 0

Hi there,

I am very new to bioinformatics, so apologies in advance for any lack of clarity in my question.

I am attempting to investigate genes that may influence the expression of CCR9, in an RNAseq dataset obtained from a large cohort of children newly diagnosed with T-ALL.

The data has been combined into a unified count matrix using RSEM and I plan on using RStudio to process these data.

My initial thoughts were to create a condition within this group based on high and low expression of CCR9 and carry out gene expression analysis using DESeq2. I am now wondering whether it may be better to use another package, where I could use CCR9 expression as a regression analysis variable? I was wondering whether Limma would be a suitable tool for this use case.

Many thanks in advance for taking the time to read my question and thank you in advance for any suggestions and advice you can provide.

rna • 348 views
Entering edit mode
4 months ago
dsull ★ 6.1k

DESeq2 can do regression with continuous variables.

Another simpler option if you have a large cohort is simply to do the pearson or spearman correlation between CCR9 and every other gene.

Entering edit mode

I don't know why, but people often say that limma-voom is more suitable for continuous variables than DESeq2. Not sure of the reason.


Login before adding your answer.

Traffic: 2108 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6