Search
Question: DESEQ2 for multi-factor analysis with age(numeric), age and group
0
gravatar for songeric1107
12 days ago by
United States
songeric11070 wrote:

 

I am trying to use DESEQ2 to estimate the correlation coefficient of gene expression levels with different phenotypes.

dds = DESeqDataSetFromMatrix(
  countData = round(count),
  colData = pheno1,
  design = ~ Age+Gender+GROUP)

 

ddm <- DESeq(dds)

my questions is how could I get the coefficient for Age, gender and group, especially, age is a continuous variable

Thank you

ADD COMMENTlink modified 11 days ago • written 12 days ago by songeric11070

Great, I got it now.

One more question, in order to get coefficient for gender and group, I have to convert them to 0,1, instead of F,M, how to interpret the coefficient of gender then?

Thank you

 

ADD REPLYlink written 11 days ago by songeric11070

I moved this to a comment.

If you print the results table for sex, it will tell you the LFC direction. It’s also determined by the factor levels. See the vignette section on factor levels.

ADD REPLYlink written 11 days ago by Michael Love19k
0
gravatar for Michael Love
12 days ago by
Michael Love19k
United States
Michael Love19k wrote:

You can get the coefficients one by one with results() and the "name" argument. Or you can get all the coefficients with the coef() approach described here:

https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#access-to-all-calculated-values

ADD COMMENTlink modified 12 days ago • written 12 days ago by Michael Love19k
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: 307 users visited in the last hour