Question: Best way to compare expression values between two groups
0
gravatar for Lila M
3.2 years ago by
Lila M 860
UK
Lila M 860 wrote:

Hi guys, I have two matrix with expression values for near 19,000 transcripts for two different groups in two different matrix: group1 (n=30) and group2 (n=60). So for each transcript, I have around 90 different value.

data format for Group 1 but is the same for Group 2

GENE                Group1           Group1          Group1         Group1
ENSG00000183154.1   0.925443037     3.0369279927    2.7557872516    2.2384197806
ENSG00000227210.1   2.0079999555    3.2941987268    1.7373249132    1.8805827534
ENSG00000240591.1   1.833712588     2.8881138203    2.8422437594    1.654280957
ENSG00000279342.1   2.1612707148    2.3112540357    3.4992176678    2.3862284068
ENSG00000248383.4   2.2085886874    2.7214426016    0                0
ENSG00000253837.1   1.384282608     3.4071437949    1.7373249132    1.3734792256
ENSG00000107165.11  1.6305563       1.2443063796    2.3422637913    3.8709988618

So for each gene I want to know if there are differences between groups and also if there are differences in the expression average. I'm planing to do a chisq using R, but not sure which could be the best method, any suggestion?

Thanks!

R p-value comparation • 1.9k views
ADD COMMENTlink modified 3.1 years ago by theobroma221.1k • written 3.2 years ago by Lila M 860

Do you have access to the raw read counts? If so, please read https://f1000research.com/articles/4-1070/v2

ADD REPLYlink modified 3.2 years ago • written 3.2 years ago by James Ashmore3.0k

Yes, this is exactly what I did, in more detail, the matrix is derived from varianceStabilizingTransformation, but I can't find a method to get the differences among each transcript per group and neither if there are differences for the average expression. I have the PCA, the heatmap... but I need those values also, any idea?

ADD REPLYlink written 3.2 years ago by Lila M 860

Okay I'm slightly confused, it sounds like you want to do a differential expression analysis (an example of which you can find in the workflow I linked). Is that correct?

ADD REPLYlink written 3.2 years ago by James Ashmore3.0k

Yes it is! But I couldn't find a way to get the differential expression for each gene (" to say if ENSG00000183154.1 is more expressed in group1 and group2). So is there a function to do that?

ADD REPLYlink written 3.2 years ago by Lila M 860

If you look at the "Building the results table" section of that workflow you'll see that the "results" function will give you a table of genes which are either up- or down- regulated in either group.

ADD REPLYlink written 3.2 years ago by James Ashmore3.0k

So the function res <- results(dds) that report a pvalue for each gene is giving me that information? or should I use instead sum(res$pvalue < 0.05, na.rm=TRUE) for statically diferences?

ADD REPLYlink written 3.2 years ago by Lila M 860
1

sum(res$pvalue < 0.05, na.rm=TRUE) will give you the number of DEGs. results(dds) will give you the list.

You also might find useful to: DESeq2 proper design setting

ADD REPLYlink written 3.2 years ago by Macspider3.2k

And how can I export the transcripts ID that are deferentially expressed in both groups? is the output a data frame?

ADD REPLYlink written 3.2 years ago by Lila M 860
1

I don't remember what class the output is, but typing class(results(dds)) will tell you so. Anyway, you just need to write.table(results(dds), file="yourfilenamechoice.table")to export them into a file.

ADD REPLYlink written 3.2 years ago by Macspider3.2k

many thanks for the information!

ADD REPLYlink written 3.2 years ago by Lila M 860
0
gravatar for theobroma22
3.1 years ago by
theobroma221.1k
theobroma221.1k wrote:

If your comparing two groups you could draw up a linear model like in the Bioconductor package Limma, or use ANOVA.

ADD COMMENTlink written 3.1 years ago by theobroma221.1k
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.3.0
Traffic: 936 users visited in the last hour