getGibbs: Update the current gibbs sample for a 'HurdleStructure'

In amcdavid/HurdleNormal: 'Estimation and sampling for multi-dimensional Hurdle models on a Normal density with applications to single-cell co-expression'

Update the current gibbs sample for a HurdleStructure

Description

Replaces the item 'gibbs'

Usage

getGibbs(hs, Nkeep = 2000, kcell, post_fun, mvnorm, ...)

Arguments

hs	HurdleStructure
Nkeep	final number of samples to keep
kcell	Number of cells per sample (cells added on exponential scale)
post_fun	Function to apply post-sampling, ie, to add contaminating noise
mvnorm	Sample from multivariate normal (using only K matrix)–intended for applying post-selection models
...	passed to rGibbsHurdle

Value

modified HurdleStructure

amcdavid/HurdleNormal documentation built on May 14, 2022, 11:12 p.m.