powered by
cdf_post
pcdf_post(m, c, data)
beta = rnorm(100,0,1) betahat= beta+rnorm(100,0,1) sebetahat=rep(1,100) ash.beta = ash(betahat,1,mixcompdist="normal") c = pcdf_post(get_fitted_g(ash.beta),beta,set_data(betahat,sebetahat))
Run the code above in your browser using DataLab