estimateTheta: Estimate theta values based on hard clustering

In daewoooo/SaaRclust: R Package to cluster long sequencing reads into chromosomes to facilitate de novo genome assembly.

Description

This function takes results of hard clustering and estimates majority cell types for each Strand-seq library

Usage

estimateTheta(counts.l = NULL, hard.clust = NULL, alpha = 0.1)

Arguments

counts.l	A list of directional read counts per PB read per library.
hard.clust	A integer of cluster assignments for each PacBio read.
alpha	Estimated level of background in Strand-seq reads.

Value

A list of estimated theta values for every cluster and cell.

Author(s)

David Porubsky

daewoooo/SaaRclust documentation built on May 28, 2019, 7:50 p.m.