fishpond

Differential expression and allelic analysis, nonparametric statistics

https://github.com/thelovelab/fishpond

Science Score: 33.0%

This score indicates how likely this project is to be science-related based on various indicators:

  • CITATION.cff file
  • codemeta.json file
  • .zenodo.json file
  • DOI references
    Found 3 DOI reference(s) in README
  • Academic publication links
    Links to: ncbi.nlm.nih.gov
  • Committers with academic emails
    3 of 12 committers (25.0%) from academic institutions
  • Institutional organization owner
  • JOSS paper metadata
  • Scientific vocabulary similarity
    Low similarity (12.2%) to scientific vocabulary

Keywords

bioconductor gene-expression genomics rnaseq salmon scrnaseq statistics transcriptomics

Keywords from Contributors

bioconductor-packages grna-sequence gene ontologies gsea sequencing deseq2 gene-symbols mds-dimensions pca
Last synced: 6 months ago · JSON representation

Repository

Differential expression and allelic analysis, nonparametric statistics

Basic Info
Statistics
  • Stars: 29
  • Watchers: 2
  • Forks: 10
  • Open Issues: 1
  • Releases: 1
Topics
bioconductor gene-expression genomics rnaseq salmon scrnaseq statistics transcriptomics
Created over 7 years ago · Last pushed about 1 year ago
Metadata Files
Readme Code of conduct

README.md

fishpond

R build status

Fishpond: downstream methods and tools for expression data

Fishpond contains a method, swish(), for differential transcript and gene expression analysis of RNA-seq data using inferential replicates. Also the package contains utilities for working with Salmon, alevin, and alevin-fry quantification data, including loadFry().

Quick start

The following paradigm is used for running a Swish analysis:

y <- tximeta(coldata) # reads in counts and inf reps y <- scaleInfReps(y) # scales counts y <- labelKeep(y) # labels features to keep set.seed(1) # for reproducibility y <- swish(y, x="condition") # simplest Swish case

How does Swish work

Swish accounts for inferential uncertainty in expression estimates by averaging test statistics over a number of inferential replicate datasets, either posterior samples or bootstrap samples. This is inspired by a method called SAMseq, hence we named our method Swish, for "SAMseq With Inferential Samples Helps". Averaging over inferential replicates produces a different test statistic than what one would obtain using only point estimates for expression level.

For example, one of the tests possible with swish() is a correlation test of expression level over a condition variable. We can visualize the distribution of inferential replicates with plotInfReps():

The test statistic is formed by averaging over these sets of data:

p-values and q-values are computed through permutation of samples (see vignette for details on permutation schemes).

The Swish method is described in the following publication:

Zhu, A., Srivastava, A., Ibrahim, J.G., Patro, R., Love, M.I. "Nonparametric expression analysis using inferential replicate counts" Nucleic Acids Research (2019) 47(18):e105 PMC6765120

The SEESAW method for allelic expression analysis is described in the following preprint:

Euphy Wu, Noor P. Singh, Kwangbom Choi, Mohsen Zakeri, Matthew Vincent, Gary A. Churchill, Cheryl L. Ackert-Bicknell, Rob Patro, Michael I. Love. "Detecting isoform-level allelic imbalance accounting for inferential uncertainty" bioRxiv (2022) doi: 10.1101/2022.08.12.503785

Installation

This package can be installed via Bioconductor:

BiocManager::install("fishpond")

Funding

This work was funded by NIH NHGRI R01-HG009937.

Owner

  • Name: The Love Lab
  • Login: thelovelab
  • Kind: organization
  • Email: michaelisaiahlove [at] gmail
  • Location: United States of America

Software produced by members of the Love Lab (UNC-Chapel Hill)

GitHub Events

Total
  • Issues event: 1
  • Watch event: 2
  • Issue comment event: 3
  • Push event: 6
  • Pull request event: 2
Last Year
  • Issues event: 1
  • Watch event: 2
  • Issue comment event: 3
  • Push event: 6
  • Pull request event: 2

Committers

Last synced: over 2 years ago

All Time
  • Total Commits: 461
  • Total Committers: 12
  • Avg Commits per committer: 38.417
  • Development Distribution Score (DDS): 0.167
Past Year
  • Commits: 29
  • Committers: 6
  • Avg Commits per committer: 4.833
  • Development Distribution Score (DDS): 0.379
Top Committers
Name Email Commits
mikelove m****e@g****m 384
jgilis j****s@u****e 24
Nitesh Turaga n****a@g****m 14
DongzeHe d****a@g****m 11
DongzeHE 3****E 9
Steve Lianoglou l****u@d****m 6
unknown w****s@g****m 4
Anqi Zhu a****u@l****u 2
J Wokaty j****y@s****u 2
Mike Love m****e 2
J Wokaty j****y 2
k3yavi a****a@c****u 1
Committer Domains (Top 20 + Academic)

Issues and Pull Requests

Last synced: 6 months ago

All Time
  • Total issues: 18
  • Total pull requests: 23
  • Average time to close issues: 18 days
  • Average time to close pull requests: about 23 hours
  • Total issue authors: 14
  • Total pull request authors: 7
  • Average comments per issue: 3.22
  • Average comments per pull request: 1.96
  • Merged pull requests: 20
  • Bot issues: 0
  • Bot pull requests: 0
Past Year
  • Issues: 2
  • Pull requests: 2
  • Average time to close issues: 5 days
  • Average time to close pull requests: 1 day
  • Issue authors: 2
  • Pull request authors: 1
  • Average comments per issue: 3.0
  • Average comments per pull request: 1.0
  • Merged pull requests: 2
  • Bot issues: 0
  • Bot pull requests: 0
Top Authors
Issue Authors
  • mikelove (3)
  • sunliang3361 (2)
  • zhangjy859 (2)
  • DrHogart (1)
  • DongzeHE (1)
  • qxxxd (1)
  • smlchen (1)
  • HenrikBengtsson (1)
  • rsak-384 (1)
  • ChristelKrueger (1)
  • biozzq (1)
  • AnnaAMonaco (1)
  • JosephLalli (1)
  • jeremymsimon (1)
Pull Request Authors
  • DongzeHE (8)
  • jgilis (8)
  • Varix (2)
  • an-altosian (2)
  • k3yavi (1)
  • lianos (1)
  • JosephLalli (1)
Top Labels
Issue Labels
enhancement (1)
Pull Request Labels

Packages

  • Total packages: 1
  • Total downloads:
    • bioconductor 37,033 total
  • Total dependent packages: 1
  • Total dependent repositories: 0
  • Total versions: 5
  • Total maintainers: 1
bioconductor.org: fishpond

Fishpond: downstream methods and tools for expression data

  • Versions: 5
  • Dependent Packages: 1
  • Dependent Repositories: 0
  • Downloads: 37,033 Total
Rankings
Dependent repos count: 0.0%
Dependent packages count: 0.0%
Average: 13.5%
Downloads: 40.5%
Maintainers (1)
Last synced: 6 months ago

Dependencies

DESCRIPTION cran
  • GenomicRanges * imports
  • IRanges * imports
  • Matrix * imports
  • S4Vectors * imports
  • SingleCellExperiment * imports
  • SummarizedExperiment * imports
  • abind * imports
  • graphics * imports
  • gtools * imports
  • jsonlite * imports
  • matrixStats * imports
  • methods * imports
  • qvalue * imports
  • stats * imports
  • svMisc * imports
  • utils * imports
  • AnnotationDbi * suggests
  • DESeq2 * suggests
  • EnsDb.Hsapiens.v86 * suggests
  • GenomeInfoDb * suggests
  • GenomicFeatures * suggests
  • Gviz * suggests
  • apeglm * suggests
  • data.table * suggests
  • ensembldb * suggests
  • knitr * suggests
  • limma * suggests
  • macrophage * suggests
  • org.Hs.eg.db * suggests
  • pheatmap * suggests
  • rmarkdown * suggests
  • samr * suggests
  • testthat * suggests
  • tximeta * suggests
  • tximportData * suggests
.github/workflows/check-bioc.yml actions
  • actions/cache v2 composite
  • actions/checkout v2 composite
  • actions/upload-artifact master composite
  • docker/build-push-action v1 composite
  • r-lib/actions/setup-pandoc v2 composite
  • r-lib/actions/setup-r v2 composite