bp.sim.manipulated: Creates manipulated samples for simulations involving...

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

Description

Creates manipulated samples based on the control samples for simulations involving Beta-Poisson (BP) models

Usage

bp.sim.manipulated(Res.par,N2,case,degree,seedex)

Arguments

Res.par	matrix including the fitted parameters for the well-fitted control instances that are to be manipulated (e.g. the output matrix parameters.bp.wellfit from the bp.sim.control function)
N2	size of the samples drawn from the manipulated BP models
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 models to ensure reproducibility

Details

Creates manipulated samples based on the control samples for simulations involving BP models. Details regarding the design of the manipulations can be found in Schefzik (2021).

Value

matrix with simulated manipulated single-cell RNA-sequencing expression data of dimension KxN2, where K is the number of well-fitted genes by the BP models for which the desired manipulation is feasible, see Schefzik (2021) for details

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

N1<-500
seedex<-24
ctrl<-bp.sim.control(DATA.EX,N1,seedex)
N2<-500
degree<-1/3
man1<-bp.sim.manipulated(ctrl[[2]],N2,case="DLambda",degree,seedex)
man2<-bp.sim.manipulated(ctrl[[2]],N2,case="DAlpha",degree,seedex)
man3<-bp.sim.manipulated(ctrl[[2]],N2,case="DBeta",degree,seedex)
man4<-bp.sim.manipulated(ctrl[[2]],N2,case="DAlphaBeta",degree,seedex)
man5<-bp.sim.manipulated(ctrl[[2]],N2,case="DPZ",degree,seedex)

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