Closed Chloe-Thangavelu closed 1 year ago
In addition, I would like to confirm. I have an output file from feature counts that looks like this. Will rnanorm be able convert this into TPM? I do not have access to the original BAM file.
geneID geneLength read counts
ENSG00000121410 4006 3
ENSG00000268895 2793 3
ENSG00000148584 9603 11
ENSG00000175899 6384 1
Hi @CthangavUCI
Which OS are you using? Windows? If so, issues with installing of (py)bedtools are well known for Windows users.
Are you familiar with running things in Docker?
The other option would be to install the rnanorm 2.0.0, which is not yet published o Pypi. 2.0.0 does not depend on pybedtools
.
In addition, I would like to confirm. I have an output file from feature counts that looks like this. Will rnanorm be able convert this into TPM? I do not have access to the original BAM file.
This is fine, do you have access to the GTF file used to obtain the expressions counts?
Thank you very much for your reply. Could you give me instructions on how to access/install rnanorm 2.0.0? I can try to download docker if you have advice on how to do that- I may try both methods and see which one works. Also, just to confirm, I only have a featurecounts output file containing the counts and feature lengths- is that sufficient for RNAnorm to calculate TPMs? THank you again for your help!
On Tue, Jun 13, 2023 at 1:14 AM Jure Zmrzlikar @.***> wrote:
Which OS are you using? Windows? If so, issues with installing of (py)bedtools are well known https://urldefense.com/v3/__https://www.biostars.org/p/9561321/__;!!CzAuKJ42GuquVTTmVmPViYEvSg!ODBHCXrzVPlyoZvj6yfBa6O5zDD9vT6WTANEoQYosjyt-Kz_gQMMrG8hZ91IRnwKwktEuS15_9I4BEASuzZoJExh$ for Windows users.
Are you familiar with running things in Docker?
The other option would be to install the rnanorm 2.0.0, which is not yet published o Pypi. 2.0.0 does not depend on pybedtools.
— Reply to this email directly, view it on GitHub https://urldefense.com/v3/__https://github.com/genialis/RNAnorm/issues/18*issuecomment-1588780728__;Iw!!CzAuKJ42GuquVTTmVmPViYEvSg!ODBHCXrzVPlyoZvj6yfBa6O5zDD9vT6WTANEoQYosjyt-Kz_gQMMrG8hZ91IRnwKwktEuS15_9I4BEASu-ivbrgy$, or unsubscribe https://urldefense.com/v3/__https://github.com/notifications/unsubscribe-auth/ATZSYG6LJ655R7GZPVLC5XDXLAOPDANCNFSM6AAAAAAZD6JM74__;!!CzAuKJ42GuquVTTmVmPViYEvSg!ODBHCXrzVPlyoZvj6yfBa6O5zDD9vT6WTANEoQYosjyt-Kz_gQMMrG8hZ91IRnwKwktEuS15_9I4BEASu7kA47Kn$ . You are receiving this because you were mentioned.Message ID: @.***>
-- Chloe Thangavelu PhD Candidate, University of California, Irvine Department of Biological Chemistry
Hi @CthangavUCI
If you only have is FeatureCounts file than I suggest the following
Install rnanorm 2.0.0 release candidate with: pip install rnanorm==2.0.0.rc2
Warning: this is development version. It should work fine, but not 100% bulletproof.
If this works you can use the (1)command line or (2) Python
rnanorm tpm exp.csv --gene-lengths lengths.csv
The format of the exp.csv must be samples in rows and genes in columns, e.g.
,ENSG00000121410,ENSG00000268895,...
SAMPLE_1,31,54, ...
SAMPLE_2,26,34, ...
...
The format of the lengths.csv must be two columns, first gene IDs, second gene lengths, e.g.
gene_id,gene_length
ENSG00000121410,4006
ENSG00000268895,2793
ENSG00000148584,9603
ENSG00000175899,6384
...
import pandas as pd
from rnanorm import TPM
fc = pd.read_csv(
exp_tpm = TPM(gene_lengths=gene_lengths).fit_transform(exp)
Thank you very much- option 2 worked in python for me. Is the code similar to get TMM normalization? Something like: exp_tmm = TMM().fit_transform(exp)
Is the code similar to get TMM normalization? Something like: exp_tmm = TMM().fit_transform(exp)
Exactly. All the methods have as similar interface as possible.
Hello, I am struggling to install rnanorm. When I type "pip install rnanorm" in jupyter notebooks I get the following error. I have tried to install with python versions 3.8.3 and 3.10.9, but get the same error.