-
Notifications
You must be signed in to change notification settings - Fork 4
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
recount::getTPM(rse) method #30
Labels
bug
Something isn't working
Milestone
Comments
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 <- recount3::create_rse_manual(
+ project = "BLOOD",
+ project_home = "data_sources/gtex",
+ organism = "human",
+ annotation = "gencode_v29",
+ type = "gene")
2023-05-06 19:54:15.455204 downloading and reading the metadata.
2023-05-06 19:54:15.899823 caching file gtex.gtex.BLOOD.MD.gz.
adding rname 'http://duffel.rail.bio/recount3/human/data_sources/gtex/metadata/OD/BLOOD/gtex.gtex.BLOOD.MD.gz'
|=================================================================================| 100%
2023-05-06 19:54:17.301443 caching file gtex.recount_project.BLOOD.MD.gz.
adding rname 'http://duffel.rail.bio/recount3/human/data_sources/gtex/metadata/OD/BLOOD/gtex.recount_project.BLOOD.MD.gz'
|=================================================================================| 100%
2023-05-06 19:54:18.317409 caching file gtex.recount_qc.BLOOD.MD.gz.
adding rname 'http://duffel.rail.bio/recount3/human/data_sources/gtex/metadata/OD/BLOOD/gtex.recount_qc.BLOOD.MD.gz'
|=================================================================================| 100%
2023-05-06 19:54:19.555473 caching file gtex.recount_seq_qc.BLOOD.MD.gz.
adding rname 'http://duffel.rail.bio/recount3/human/data_sources/gtex/metadata/OD/BLOOD/gtex.recount_seq_qc.BLOOD.MD.gz'
|=================================================================================| 100%
2023-05-06 19:54:20.419423 downloading and reading the feature information.
2023-05-06 19:54:20.787712 caching file human.gene_sums.G029.gtf.gz.
adding rname 'http://duffel.rail.bio/recount3/human/annotations/gene_sums/human.gene_sums.G029.gtf.gz'
|=================================================================================| 100%
2023-05-06 19:54:22.512853 downloading and reading the counts: 1048 samples across 64837 features.
2023-05-06 19:54:22.94642 caching file gtex.gene_sums.BLOOD.G029.gz.
adding rname 'http://duffel.rail.bio/recount3/human/data_sources/gtex/gene_sums/OD/BLOOD/gtex.gene_sums.BLOOD.G029.gz'
|=================================================================================| 100%
2023-05-06 19:54:40.040121 constructing the RangedSummarizedExperiment (rse) object.
devel >
devel > SummarizedExperiment::assays(rse)$counts <- recount3::transform_counts(rse)
devel > SummarizedExperiment::assays(rse)$TPM <- recount::getTPM(rse)
Setting options('download.file.method.GEOquery'='auto')
Setting options('GEOquery.inmemory.gpl'=FALSE)
devel > options(width = 120)
devel > sessioninfo::session_info()
─ Session info ───────────────────────────────────────────────────────────────────────────────────────────────────────
setting value
version R version 4.3.0 (2023-04-21)
os macOS Ventura 13.0.1
system aarch64, darwin20
ui RStudio
language (EN)
collate en_US.UTF-8
ctype en_US.UTF-8
tz America/New_York
date 2023-05-06
rstudio 2023.03.0+386 Cherry Blossom (desktop)
pandoc 2.19.2 @ /Applications/RStudio.app/Contents/Resources/app/quarto/bin/tools/ (via rmarkdown)
─ Packages ───────────────────────────────────────────────────────────────────────────────────────────────────────────
package * version date (UTC) lib source
AnnotationDbi 1.62.1 2023-05-02 [1] Bioconductor
backports 1.4.1 2021-12-13 [1] CRAN (R 4.3.0)
base64enc 0.1-3 2015-07-28 [1] CRAN (R 4.3.0)
Biobase 2.60.0 2023-04-25 [1] Bioconductor
BiocFileCache 2.8.0 2023-04-25 [1] Bioconductor
BiocGenerics 0.46.0 2023-04-25 [1] Bioconductor
BiocIO 1.10.0 2023-04-25 [1] Bioconductor
BiocParallel 1.34.0 2023-04-25 [1] Bioconductor
biocthis 1.11.1 2023-05-06 [1] Github (lcolladotor/biocthis@42dc8df)
biomaRt 2.56.0 2023-04-25 [1] Bioconductor
Biostrings 2.68.0 2023-04-25 [1] Bioconductor
bit 4.0.5 2022-11-15 [1] CRAN (R 4.3.0)
bit64 4.0.5 2020-08-30 [1] CRAN (R 4.3.0)
bitops 1.0-7 2021-04-24 [1] CRAN (R 4.3.0)
blob 1.2.4 2023-03-17 [1] CRAN (R 4.3.0)
brio 1.1.3 2021-11-30 [1] CRAN (R 4.3.0)
BSgenome 1.68.0 2023-04-25 [1] Bioconductor
bumphunter 1.42.0 2023-04-25 [1] Bioconductor
cachem 1.0.8 2023-05-01 [1] CRAN (R 4.3.0)
callr 3.7.3 2022-11-02 [1] CRAN (R 4.3.0)
checkmate 2.2.0 2023-04-27 [1] CRAN (R 4.3.0)
cli 3.6.1 2023-03-23 [1] CRAN (R 4.3.0)
cluster 2.1.4 2022-08-22 [1] CRAN (R 4.3.0)
codetools 0.2-19 2023-02-01 [1] CRAN (R 4.3.0)
colorout 1.2-2 2023-05-06 [1] Github (jalvesaq/colorout@79931fd)
colorspace 2.1-0 2023-01-23 [1] CRAN (R 4.3.0)
crayon 1.5.2 2022-09-29 [1] CRAN (R 4.3.0)
curl 5.0.0 2023-01-12 [1] CRAN (R 4.3.0)
data.table 1.14.8 2023-02-17 [1] CRAN (R 4.3.0)
DBI 1.1.3 2022-06-18 [1] CRAN (R 4.3.0)
dbplyr 2.3.2 2023-03-21 [1] CRAN (R 4.3.0)
DelayedArray 0.26.1 2023-05-01 [1] Bioconductor
derfinder 1.34.0 2023-04-25 [1] Bioconductor
derfinderHelper 1.34.0 2023-04-25 [1] Bioconductor
devtools * 2.4.5 2022-10-11 [1] CRAN (R 4.3.0)
digest 0.6.31 2022-12-11 [1] CRAN (R 4.3.0)
doRNG 1.8.6 2023-01-16 [1] CRAN (R 4.3.0)
downloader 0.4 2015-07-09 [1] CRAN (R 4.3.0)
dplyr 1.1.2 2023-04-20 [1] CRAN (R 4.3.0)
ellipsis 0.3.2 2021-04-29 [1] CRAN (R 4.3.0)
evaluate 0.21 2023-05-05 [1] CRAN (R 4.3.0)
fansi 1.0.4 2023-01-22 [1] CRAN (R 4.3.0)
fastmap 1.1.1 2023-02-24 [1] CRAN (R 4.3.0)
filelock 1.0.2 2018-10-05 [1] CRAN (R 4.3.0)
foreach 1.5.2 2022-02-02 [1] CRAN (R 4.3.0)
foreign 0.8-84 2022-12-06 [1] CRAN (R 4.3.0)
Formula 1.2-5 2023-02-24 [1] CRAN (R 4.3.0)
fs 1.6.2 2023-04-25 [1] CRAN (R 4.3.0)
generics 0.1.3 2022-07-05 [1] CRAN (R 4.3.0)
GenomeInfoDb 1.36.0 2023-04-25 [1] Bioconductor
GenomeInfoDbData 1.2.10 2023-05-06 [1] Bioconductor
GenomicAlignments 1.36.0 2023-04-25 [1] Bioconductor
GenomicFeatures 1.52.0 2023-04-25 [1] Bioconductor
GenomicFiles 1.36.0 2023-04-25 [1] Bioconductor
GenomicRanges 1.52.0 2023-04-25 [1] Bioconductor
GEOquery 2.68.0 2023-04-25 [1] Bioconductor
ggplot2 3.4.2 2023-04-03 [1] CRAN (R 4.3.0)
glue 1.6.2 2022-02-24 [1] CRAN (R 4.3.0)
gridExtra 2.3 2017-09-09 [1] CRAN (R 4.3.0)
gtable 0.3.3 2023-03-21 [1] CRAN (R 4.3.0)
Hmisc 5.0-1 2023-03-08 [1] CRAN (R 4.3.0)
hms 1.1.3 2023-03-21 [1] CRAN (R 4.3.0)
htmlTable 2.4.1 2022-07-07 [1] CRAN (R 4.3.0)
htmltools 0.5.5 2023-03-23 [1] CRAN (R 4.3.0)
htmlwidgets 1.6.2 2023-03-17 [1] CRAN (R 4.3.0)
httpuv 1.6.9 2023-02-14 [1] CRAN (R 4.3.0)
httr 1.4.5 2023-02-24 [1] CRAN (R 4.3.0)
IRanges 2.34.0 2023-04-25 [1] Bioconductor
iterators 1.0.14 2022-02-05 [1] CRAN (R 4.3.0)
jsonlite 1.8.4 2022-12-06 [1] CRAN (R 4.3.0)
KEGGREST 1.40.0 2023-04-25 [1] Bioconductor
knitr 1.42 2023-01-25 [1] CRAN (R 4.3.0)
later 1.3.1 2023-05-02 [1] CRAN (R 4.3.0)
lattice 0.21-8 2023-04-05 [1] CRAN (R 4.3.0)
lifecycle 1.0.3 2022-10-07 [1] CRAN (R 4.3.0)
limma 3.56.0 2023-04-25 [1] Bioconductor
locfit 1.5-9.7 2023-01-02 [1] CRAN (R 4.3.0)
lubridate 1.9.2 2023-02-10 [1] CRAN (R 4.3.0)
magrittr 2.0.3 2022-03-30 [1] CRAN (R 4.3.0)
Matrix 1.5-4 2023-04-04 [1] CRAN (R 4.3.0)
MatrixGenerics 1.12.0 2023-04-25 [1] Bioconductor
matrixStats 0.63.0 2022-11-18 [1] CRAN (R 4.3.0)
memoise 2.0.1 2021-11-26 [1] CRAN (R 4.3.0)
mime 0.12 2021-09-28 [1] CRAN (R 4.3.0)
miniUI 0.1.1.1 2018-05-18 [1] CRAN (R 4.3.0)
munsell 0.5.0 2018-06-12 [1] CRAN (R 4.3.0)
nnet 7.3-19 2023-05-03 [1] CRAN (R 4.3.0)
pillar 1.9.0 2023-03-22 [1] CRAN (R 4.3.0)
pkgbuild 1.4.0 2022-11-27 [1] CRAN (R 4.3.0)
pkgconfig 2.0.3 2019-09-22 [1] CRAN (R 4.3.0)
pkgload 1.3.2 2022-11-16 [1] CRAN (R 4.3.0)
plyr 1.8.8 2022-11-11 [1] CRAN (R 4.3.0)
png 0.1-8 2022-11-29 [1] CRAN (R 4.3.0)
prettyunits 1.1.1 2020-01-24 [1] CRAN (R 4.3.0)
processx 3.8.1 2023-04-18 [1] CRAN (R 4.3.0)
profvis 0.3.8 2023-05-02 [1] CRAN (R 4.3.0)
progress 1.2.2 2019-05-16 [1] CRAN (R 4.3.0)
promises 1.2.0.1 2021-02-11 [1] CRAN (R 4.3.0)
prompt 1.0.1 2023-05-06 [1] Github (gaborcsardi/prompt@7ef0f2e)
ps 1.7.5 2023-04-18 [1] CRAN (R 4.3.0)
purrr 1.0.1 2023-01-10 [1] CRAN (R 4.3.0)
qvalue 2.32.0 2023-04-25 [1] Bioconductor
R.cache 0.16.0 2022-07-21 [1] CRAN (R 4.3.0)
R.methodsS3 1.8.2 2022-06-13 [1] CRAN (R 4.3.0)
R.oo 1.25.0 2022-06-12 [1] CRAN (R 4.3.0)
R.utils 2.12.2 2022-11-11 [1] CRAN (R 4.3.0)
R6 2.5.1 2021-08-19 [1] CRAN (R 4.3.0)
rappdirs 0.3.3 2021-01-31 [1] CRAN (R 4.3.0)
Rcpp 1.0.10 2023-01-22 [1] CRAN (R 4.3.0)
RCurl 1.98-1.12 2023-03-27 [1] CRAN (R 4.3.0)
readr 2.1.4 2023-02-10 [1] CRAN (R 4.3.0)
recount 1.26.0 2023-04-25 [1] Bioconductor
recount3 1.9.1 2023-01-11 [1] Bioconductor
remotes 2.4.2 2021-11-30 [1] CRAN (R 4.3.0)
rentrez 1.2.3 2020-11-10 [1] CRAN (R 4.3.0)
reshape2 1.4.4 2020-04-09 [1] CRAN (R 4.3.0)
restfulr 0.0.15 2022-06-16 [1] CRAN (R 4.3.0)
rjson 0.2.21 2022-01-09 [1] CRAN (R 4.3.0)
rlang 1.1.1 2023-04-28 [1] CRAN (R 4.3.0)
rmarkdown 2.21 2023-03-26 [1] CRAN (R 4.3.0)
rngtools 1.5.2 2021-09-20 [1] CRAN (R 4.3.0)
rpart 4.1.19 2022-10-21 [1] CRAN (R 4.3.0)
Rsamtools 2.16.0 2023-04-25 [1] Bioconductor
RSQLite 2.3.1 2023-04-03 [1] CRAN (R 4.3.0)
rsthemes 0.4.0 2023-05-06 [1] Github (gadenbuie/rsthemes@34a55a4)
rstudioapi 0.14 2022-08-22 [1] CRAN (R 4.3.0)
rtracklayer 1.60.0 2023-04-25 [1] Bioconductor
S4Arrays 1.0.1 2023-05-01 [1] Bioconductor
S4Vectors 0.38.1 2023-05-02 [1] Bioconductor
scales 1.2.1 2022-08-20 [1] CRAN (R 4.3.0)
sessioninfo 1.2.2 2021-12-06 [1] CRAN (R 4.3.0)
shiny 1.7.4 2022-12-15 [1] CRAN (R 4.3.0)
stringi 1.7.12 2023-01-11 [1] CRAN (R 4.3.0)
stringr 1.5.0 2022-12-02 [1] CRAN (R 4.3.0)
styler 1.9.1 2023-03-04 [1] CRAN (R 4.3.0)
SummarizedExperiment 1.30.1 2023-05-01 [1] Bioconductor
suncalc 0.5.1 2022-09-29 [1] CRAN (R 4.3.0)
testthat * 3.1.8 2023-05-04 [1] CRAN (R 4.3.0)
tibble 3.2.1 2023-03-20 [1] CRAN (R 4.3.0)
tidyr 1.3.0 2023-01-24 [1] CRAN (R 4.3.0)
tidyselect 1.2.0 2022-10-10 [1] CRAN (R 4.3.0)
timechange 0.2.0 2023-01-11 [1] CRAN (R 4.3.0)
tzdb 0.3.0 2022-03-28 [1] CRAN (R 4.3.0)
urlchecker 1.0.1 2021-11-30 [1] CRAN (R 4.3.0)
usethis * 2.1.6 2022-05-25 [1] CRAN (R 4.3.0)
utf8 1.2.3 2023-01-31 [1] CRAN (R 4.3.0)
VariantAnnotation 1.46.0 2023-04-25 [1] Bioconductor
vctrs 0.6.2 2023-04-19 [1] CRAN (R 4.3.0)
withr 2.5.0 2022-03-03 [1] CRAN (R 4.3.0)
xfun 0.39 2023-04-20 [1] CRAN (R 4.3.0)
XML 3.99-0.14 2023-03-19 [1] CRAN (R 4.3.0)
xml2 1.3.4 2023-04-27 [1] CRAN (R 4.3.0)
xtable 1.8-4 2019-04-21 [1] CRAN (R 4.3.0)
XVector 0.40.0 2023-04-25 [1] Bioconductor
yaml 2.3.7 2023-01-23 [1] CRAN (R 4.3.0)
zlibbioc 1.46.0 2023-04-25 [1] Bioconductor
[1] /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library
──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
devel > |
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 |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
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_rse_manual()', 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
The text was updated successfully, but these errors were encountered: