Learn R Programming

SnapATAC (Latest Updates: 2019-09-19)

SnapATAC (Single Nucleus Analysis Pipeline for ATAC-seq) is a fast, accurate and comprehensive method for analyzing single cell ATAC-seq datasets.

Latest News

FAQs

Requirements

  • Linux/Unix
  • Python (>= 2.7 & < 3.0) (SnapTools) (highly recommanded for 2.7);
  • R (>= 3.4.0 & < 3.6.0) (SnapATAC) (3.6 does not work for rhdf5 package);

Pre-print

Rongxin Fang, Sebastian Preissl, Xiaomeng Hou, Jacinta Lucero, Xinxin Wang, Amir Motamedi, Andrew K. Shiau, Eran A. Mukamel, Yanxiao Zhang, M. Margarita Behrens, Joseph Ecker, Bing Ren. Fast and Accurate Clustering of Single Cell Epigenomes Reveals Cis-Regulatory Elements in Rare Cell Types. bioRxiv 615179; doi: https://doi.org/10.1101/615179

Installation

SnapATAC has two components: Snaptools and SnapATAC.

  • SnapTools - a python module for pre-processing and working with snap file.
  • SnapATAC - a R package for the clustering, annotation, motif discovery and downstream analysis.

Install snaptools from PyPI. See how to install snaptools on FAQs. NOTE: Please use python 2.7 if possible.

$ pip install snaptools

Install SnapATAC R pakcage (development version).

$ R
> library(devtools)
> install_github("r3fang/SnapATAC")

Galleries & Tutorials (click on the image for details)

Copy Link

Version

Version

1.0.0

License

GPL-3 | file LICENSE

Issues

Pull Requests

Stars

Forks

Maintainer

Rongxin Fang

Last Published

March 4th, 2021

Functions in SnapATAC (1.0.0)

addBmatToSnap

Add cell-by-bin matrix
addGmatToSnap

Add cell-by-gene matrix
createPmat

Create Cell-by-Peak Matrix
readMetaData

Read meta data from a snap file
createGmatFromMat

Create Cell-by-Gene Matrix
demo.sp

Single Nucleus ATAC-seq Object
rmBmatFromSnap

Remove cell-by-bin matrix
isSnapFile

Check a snap-format file
jaccard-class

An S4 class jaccard to represent a jaccard object.
dim.reduct-class

An S4 class to represent dimentionality reduction object.
runHomer

Motif Analysis Using Homer
runJDA

Jaccard Distance Analysis
plotFeatureSingle

plotFeatureSingle
snapToSeurat

Convert a snap object to a seurat object
[,dim.reduct,ANY,ANY,ANY-method

subsetting for dim.reduct objects
boxPlotFeature

Feature Enrichment Boxplot
plotGene

Plot gene-body accessibility level
createSnap

Create a snap object from a snap file
calBmatCor

Check correlation of cell-by-bin matrix
filterBins

Feature filtration
filterCells

Cell filtration
createSnapFromBmat

Create a snap object from cell-by-bin matrix
plotBarcode

Plot Barocde Quality Control Distribution
kgraph-class

An S4 class jaccard to represent a knn graph object.
rowMeans

rowMeans for snap objects
runSnapAddPmat

Add cell-by-peak matrix to a snap file
makeBinary

Convert the count matrix to binary matrix
rowSums

rowSums for snap objects
plotBinCoverage

Plot Bin Coverage Distribution
colMeans

colMeans for snap objects
colSums

colSums for snap objects
createSnapFromGmat

Create a snap object from cell-by-gene matrix
createSnapFromPmat

Create a snap object from cell-by-peak matrix
runLDA

Latent Dirichlet Allocation
is.snap

Check snap object
plotDimReductElbow

Elbow Plot for Dimentionality Reduction Result
findDAR

Identifiy Differential Accessible Regions (DARs)
addPmatToSnap

Add cell-by-peak matrix
barcodeInSnapFile

Check barcode existance in snap file
plotDimReductPW

Pairwise plot for Dimentionality Reduction Result
plotViz

Visulization
runDimReduct

Linear Dimentionality Reduction
predictGenePeakPair

Link Distal Regulatory Elements (peaks) to Putative Target Genes
exportMetaData

Export barcode meta data
extractReads

Extract Reads By Barcodes
runMACS

Call Peaks Using MACS2
runHarmony

Perform batch effect using harmony
newSnap

Create an empty snap object
rmGmatFromSnap

Remove cell-by-gene matrix
nrow

nrow for snap object.
runDiffusionMaps

Dimentionality Reduction by Diffusion Maps Algorithm
runLSA

Latent Semantic Analysis
rmPmatFromSnap

Remove cell-by-peak matrix
runDiffusionMapsExtension

Diffusion Maps Extension
runMagic

Affinity Graph-based Smoothing
runNormJaccard

Normalize Jaccard Index Matrix
[,jaccard,ANY,ANY,ANY-method

subsetting for jaccard objects
scaleCountMatrix

Normlaize Count Matrix
[,kgraph,ANY,ANY,ANY-method

subsetting for kgraph objects
snap-class

An S4 class snap to represent single-nucleus accessibility object.
snapRbind

Combine snap objects
runViz

Dimentionality Reduction for Visualization
showBinSizes

Show bin sizes in a snap file
[,snap,ANY,ANY,ANY-method

subsetting for snap objects
runMACSForAll

Call Peaks Using MACS2 For All Clusters
summarySnap

summarySnap for snap object.
runCluster

Find Clusters Using Louvain/Leiden Algorithm
runChromVAR

Estimate Motif Variability Using ChromVAR
runJaccard

Calcualte Jaccard Index Matrix
runKNN

K Nearest Neighbour Graph Construction