pcdf_post
``parallel" vector version of cdf_post
where c is a vector, of same length as betahat and sebetahat
pcdf_post(m, c, data)
m
: mixture distribution with k componentsc
: a numeric vector with n elementsdata
: depends on contextan n vector, whose ith element is the cdf for beta_i at c_i
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))