AdjustSIZE_fun: Adjust MME size estimate

In WT215/bayNorm: Single-cell RNA sequencing data normalization

Adjust MME size estimate

Description

This function adjusts MME estimated size parameter of prior, which is a negative binomial distribution, using estimates from maximizing marginal distirbution (BB_SIZE). Simulation studies has shown this hybrid method of using adjusted MME size estimates is the most robust (see bayNorm paper). Hence, this is the default option for estimating size in bayNorm.

Usage

AdjustSIZE_fun(BB_SIZE, MME_MU, MME_SIZE)

Arguments

BB_SIZE	size estimated from BB_Fun.
MME_MU	mu estimated from EstPrior.
MME_SIZE	size estimated from EstPrior.

Value

MME_SIZE_adjust: A vector of estimated size. Adjusted MME_SIZE based on BB_SIZE (size estimated by maximizing marginal distribution)

Examples

data('EXAMPLE_DATA_list')
MME_MU<-rlnorm(100,meanlog=5,sdlog=1)
MME_SIZE<-rlnorm(100,meanlog=1,sdlog=1)
BB_SIZE<-rlnorm(100,meanlog=0.5,sdlog=1)
adjustt<-AdjustSIZE_fun(BB_SIZE, MME_MU, MME_SIZE)

WT215/bayNorm documentation built on Sept. 2, 2022, 1:46 a.m.