R: Sample from posterior

get_post_sample {ashr} | R Documentation

## 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)

[Package ashr version 2.2-54 Index]