comp_postsd
output posterior sd for beta for each component of prior mixture m,given data
comp_postsd(m, data)
m
: mixture distribution with k componentsdata
: details depend on the modelbeta = rnorm(100,0,1) betahat= beta+rnorm(100,0,1) ash.beta = ash(betahat,1,mixcompdist="normal") data= set_data(betahat,rep(1,100)) comp_postmean(get_fitted_g(ash.beta),data) comp_postsd(get_fitted_g(ash.beta),data) comp_postprob(get_fitted_g(ash.beta),data)