Before running any RNA-seq data analysis packages, I want to do statistical test for my gene expression data. I have 4 replicates for two treatment (r and s). As RNA-seq data have -ve bionomial distribution. Does it make any sense if i do gene wise paired t-test, so that i can get p-values for each gene.
r1 r2 r2 r4 s1 s2 s3 s4 gene1 569 257 301 446 259 289 835 414 gene2 24 1 28 26 15 1 51 36 . .
The following lines (this is not a code) is what i would like to execute :
gene1 :paired t test[(569 257 301 446) Vs (259 289 835 414)] gene2 :paired t test [(24 1 28 26) Vs (15 1 51 36)]
This is just for getting an idea about the data values per conditions not for actual analysis. Can i do paired t-test in this way? also does the above calculation make any sense on getting p-value per gene ? Seeking suggetions in this context.