ashr: get_post_sample – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

get_post_sample

Sample from posterior

Description

Returns random samples from the posterior distribution for each observation in an ash object. A matrix is returned, with columns corresponding to observations and rows corresponding to samples.

Usage

get_post_sample(a, nsamp)

Arguments

`a`	the fitted ash object
`nsamp`	number of samples to return (for each observation)

Examples

beta = rnorm(100,0,1)
betahat= beta+rnorm(100,0,1)
ash.beta = ash(betahat,1,mixcompdist="normal")
post.beta = get_post_sample(ash.beta,1000)

ashr

Methods for Adaptive Shrinkage, using Empirical Bayes

v2.2-47

GPL (>= 3)

Authors

Matthew Stephens [aut], Peter Carbonetto [aut, cre], Chaoxing Dai [ctb], David Gerard [aut], Mengyin Lu [aut], Lei Sun [aut], Jason Willwerscheid [aut], Nan Xiao [aut], Mazon Zeng [ctb]

Initial release

2020-02-19