Generates function that samples L from a flash fit object. The columns of L are sampled independently from their marginal posteriors (conditional on F being fixed at its expectation).
flash_l_sampler(data, f, kset = NULL, ebnm_fn = ebnm_pn)
| data | a flash data object |
|---|---|
| f | a flash fit object |
| kset | the indices of factor/loadings to include when sampling LF (defaults to all) |
| ebnm_fn | function used to solve the Empirical Bayes Normal Means problem |
A function that takes a single parameter nsamp, the number of samples of L to be produced by the sampler. This sampler returns a list of matrices.