bp.dd.single: Creates control and manipulated Beta-Poisson samples for...

In RomanSchefzik/SimBPDD: Simulating differential distributions for Beta-Poisson models

Description

Creates control and manipulated Beta-Poisson (BP) samples for pre-specified BP parameters

Usage

bp.dd.single(pars,N1,N2,case,degree,seedex)

Arguments

pars	vector consisting of pre-specified values for the BP parameters alpha, beta, lambda1, lambda2 and p0 (in this order!)
N1	size of the sample drawn from the control BP model
N2	size of the sample drawn from the manipulated BP model
case	case corresponding to a specific parameter manipulation in the BP model; specifically, case can be one of "DLambda", "DAlpha", "DBeta", "DAlphaBeta" or "DPZ" following the design and the respective descriptive table in Schefzik (2021)
degree	parameter to set the degree of the created difference (low to strong), see Schefzik (2021) for details and the choice of a range of possible values
seedex	seed used for sampling from the fitted BP model to ensure reproducibility

Details

Creates control and manipulated BP samples for pre-specified BP parameters

Value

A list of two:

• sim.bp.ctrl: vector of simulated control BP sample

• sim.bp.manip: vector of simulated manipulated BP sample

References

R. Schefzik (2021). SimBPDD: Simulating differential distributions in Beta-Poisson models, in particular for single-cell RNA-sequencing data. Annales Mathematicae et Informaticae, 53:283-298. Available at https://ami.uni-eszterhazy.hu/uploads/papers/finalpdf/AMI_53_from283to298.pdf

Examples

#create vector consisting of pre-specified values for the BP parameters alpha, beta, lambda1,
#lambda2 and p0 (in this order!)
pars<-c(0.21,3.50,96.13,0.05218,0.02)
N1<-500
N2<-500
degree<-1/3
seedex<-24
dd1<-bp.dd.single(pars,N1,N2,case="DLambda",degree,seedex)
dd2<-bp.dd.single(pars,N1,N2,case="DAlpha",degree,seedex)
dd3<-bp.dd.single(pars,N1,N2,case="DBeta",degree,seedex)
dd4<-bp.dd.single(pars,N1,N2,case="DAlphaBeta",degree,seedex)
dd5<-bp.dd.single(pars,N1,N2,case="DPZ",degree,seedex)

RomanSchefzik/SimBPDD documentation built on May 24, 2021, 1:28 a.m.