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 components |
data |
details depend on the model |
beta = 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)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.