statistical analyses for transcriptomes
hi there,
I need to do some statistical analyses of my transcriptomes.
I have a database with 4 columns (gene ID (categoric), expression level (numeric), individual, species). I have 2 different species and 5 ind per species. Per each individual I have more than 20000 genes (some of them are more expressed than others).
What I want to know is whether is there differences between the expression level between species. The distribution of my data doesn't follow a Gaussian distribution.
For analysing my data I run:
wilcox.test(Exp~Species, data =data)
and then,
Wilcoxon rank sum test with continuity correction
data: Exp by Species
W = 8573700000, p-value < 2.2e-16
alternative hypothesis: true location shift is not equal to 0
According to these result there should be a significan difference in the expression level between species.BUT:
- I am not sure if the analyses are appropiate for this dataset.
- Is there any way where I can take into account (as a random factor) the ID gene?
Thank you so much in advance
• 1,554 views
•
link
updated
by
ATpoint
•
written
by
luzglongoria •
•
written
by
luzglongoria •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Permutation test for convergent gene expression changes
written by Iulia •
-
How can one determine if a sample size is representative in scRNA-seq analysis?
written by fifty_fifty •
-
Approach for multi species comparative transcriptomics and differential expression analyses
written by fish_enthusiast •
-
Can one sample t.test be used to determine statistical significance for log2FC values?
written by aUser •
-
Differential expression analysis using Lung cancer CPTAC RNA SEQ Data
written by Ezequiel •
-
Select and keep only some results from DESEq2
written by luzglongoria •
-
Feature Counts how to identify each sample
written by luzglongoria •
-
Differential Expression Analysis of Two Groups Time Series Data with Limma
written by hkarakurt
-
Evaluate nucleotide and haplotype diversity using sequences from single-pool RAD data
written by chrys.delord •
-
Which statistical test to compare expression of subset of genes from two cell types
written by florian.noack •
which type of expression data you have (rna-seq, microarray, etc...) ?
Sorry, I didn't say. It is RNA-seq
I think you are looking for differential expression analysis. Check the Bioconductor 2018 Workshop chapter 6 and 7 for more details.
It probably follows the negative binomial. This is normal and expected. Check out the common differential analysis pipelines, such as DESeq2, edgeR or limma/voom. All are well-documented.