cdf_post | R Documentation |
evaluate cdf of posterior distribution of beta at c. m is the prior on beta, a mixture; c is location of evaluation assumption is betahat | beta ~ t_v(beta,sebetahat)
cdf_post(m, c, data)
m |
mixture distribution with k components |
c |
a scalar |
data |
details depend on model |
an n vector containing the cdf for beta_i at c
beta = rnorm(100,0,1)
betahat= beta+rnorm(100,0,1)
sebetahat=rep(1,100)
ash.beta = ash(betahat,1,mixcompdist="normal")
cdf0 = cdf_post(ash.beta$fitted_g,0,set_data(betahat,sebetahat))
graphics::plot(cdf0,1-get_pp(ash.beta))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.