Utility functions for calculating nucleosome banding scores from ATAC-seq data.
Adapted from Andrew Hill's R script, which contains methods used in Cusanovich et al.'s publication.
See this doc for an introduction to spectral analysis.
library(devtools)
install_github("anthony-aylward/nucband")
Assuming insert_size
is an integer vector of insert sizes:
library(nucband)
density <- insert_size_density(insert_size)
score <- nucleosome_banding_score(density)
log(score)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.