update_x_pb: update_x_p

In sellerslab/gemini: GEMINI: Variational inference approach to infer genetic interactions from pairwise CRISPR screens

update_x_p

Description

Update values of x using data from Input and current values of other parameters.

Usage

update_x_pb(
  Model,
  mean_x = 1,
  sd_x = 1,
  mean_xx = 1,
  sd_xx = 1,
  cores = 1,
  verbose = FALSE
)

Arguments

Model	a Model object of class gemini.model
mean_x	a numeric indicating prior mean of x
sd_x	a numeric indicating prior sd of x
mean_xx	a numeric indicating prior mean of xx
sd_xx	a numeric indicating prior sd of xx
cores	a numeric indicating the number of cores to use. See gemini_parallelization for details. (default=1).
verbose	default FALSE

Value

An object of class gemini.model

Note

The structure of the screen may impede parallelization. Our ability to parallelize the updates is contingent upon the independence between guides in position 1 and 2. To account for potential dependence, we define three groups of guides: group 1 (guides only in position 1 ), group 2 (guides only in position 2), and group 3 (guides in both position 1 and position 2). Parallelization is possible for groups 1 and 2, but only serial updates are possible for group 3. As such, updates for this group will take longer.

Examples

data("Model", package = "gemini")
Model %<>% update_x_pb()

sellerslab/gemini documentation built on Dec. 5, 2022, 8:56 a.m.