feeding RPKM values into DE analysis tools
0
0
Entering edit mode
8.8 years ago
zizigolu ★ 4.4k

hi,

I have two conditions A and B with 2 replicates. does it make any sense to use RPKMs for fold change or t-test analysis to identify differential expressed genes????

thank you

RNA-Seq DE RPKM t-test • 2.4k views
ADD COMMENT
1
Entering edit mode

Is there a reason you don't want to use DESeq2/edgeR/limma?

ADD REPLY
0
Entering edit mode

you all right Devon,

I back home and I am with a lap top now... I should perform miRNA DE analysis..I have remote access to linux machine in MPI but it is too slow, I have fedora virtual box on my lap top but without anything installed such as python last night Ram helped me to equipped that but actually I am too naive in computer to be able to install python or like so, I tried galaxy but even uploading my files took many hours at last failed to complete... then my last chance is something working on windows like BBmap might be have at least RPKMs to perform t-test by stand alone tools. I don't have raw counts thereby no way to use featurecounts and then DESeq2 via R

ADD REPLY
1
Entering edit mode

If you have RPKM at gene level, i.e if they are "not estimated counts", you can simply convert them to raw counts using the RPKM formula.

ADD REPLY
0
Entering edit mode

sorry you mean t-test between condition by RPKMs does not make sense?

ADD REPLY
0
Entering edit mode

sorry Devon,

I want to count the miRNA reads using the sorted .bam file containing only mapped reads (all generated by samtools from Bowtie output --SAM file) by featurecounts for Aspergillus fumigatus but in mirbase or wherever there was not a gtf for miRNA then I am using a GTF from ensembl. there is a gff3 contains miRNA but featurecounts did not recognize that. I converted that to gtf by cufflinks but again featurecounts did not recognize that. the GTF from enseble does not contain any miRNA then I am going wrong with this GTF??

thank you

ADD REPLY
1
Entering edit mode

You might need to write a little script to correctly convert the gff3 file to a standard GTF that other things can handle (GFF3 is really an annoying format).

ADD REPLY
0
Entering edit mode

thank you...scripting, my favorite!!! :)

when I converted gff3 to GTF by cufflink all miRNAs which were in gff3 disappeared in GTF

ADD REPLY

Login before adding your answer.

Traffic: 3305 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6