Open SoniaRuiz opened 1 year ago
Hm.... http://research.libd.org/recount3/reference/transform_counts.html shows a working case where computing TPMs does work.
Also, using your example R code I couldn't reproduce the issue:
devel > rse
class: RangedSummarizedExperiment
dim: 64837 1048
metadata(8): time_created recount3_version ... annotation recount3_url
assays(3): raw_counts counts TPM
rownames(64837): ENSG00000278704.1 ENSG00000277400.1 ... ENSG00000182484.15_PAR_Y ENSG00000227159.8_PAR_Y
rowData names(10): source type ... havana_gene tag
colnames(1048): GTEX-R55C-0005-SM-3GAE9.1 GTEX-QMR6-0005-SM-32PKY.1 ... GTEX-1NSGN-0005-SM-DKPPX.1
GTEX-1QP29-0005-SM-DLZQX.1
colData names(198): rail_id external_id ... recount_seq_qc.errq BigWigURL
Description of the issue
The function "recount::getTPM(rse)" returns the unexpected error message: 'Error: subscript contains invalid names'.
Below, is a reproducible example that I have been using:
When I look at the 'rowData(rse)' dimension of the object 'rse', I cannot see the column name 'bp_length':
Maybe that's the issue?
Expected behavior
The function 'recount::getTPM(rse)' to return the gene-level TPM values across the samples of the project provided on the function _'recount3::create_rsemanual()', as it is reported on https://rdrr.io/bioc/recount/man/getTPM.html.
R Session Information
Below is the output of
sessioninfo::session_info()
.Indicate whether
BiocManager::valid()
returnsTRUE
.BiocManager::valid()
isTRUE
Is the package installed via bioconda?
NO