plotMixDensity             package:BGmix             R Documentation

_P_l_o_t _p_r_e_d_i_c_t_i_v_e _d_e_n_s_i_t_y _o_f _d_a_t_a.

_D_e_s_c_r_i_p_t_i_o_n:

     Plot predictive density of data superimposed on histograms of
     observed data. Separate plots for ybar and sums of squares.

_U_s_a_g_e:

     plotMixDensity(res, predres, ybar, ss)

_A_r_g_u_m_e_n_t_s:

     res: list object output from 'ccParams'

 predres: list object output from 'ccPred' (need q.trace=T in 'ccPred')

    ybar: ybar data (see 'BGmix' help for details)

      ss: ss data (see 'BGmix' help for details)

_D_e_t_a_i_l_s:

     Note that this function is written for the unpaired differential
     expression design.

_A_u_t_h_o_r(_s):

     Alex Lewin

_E_x_a_m_p_l_e_s:

     ## Note this is a very short MCMC run!
     ## For good analysis need proper burn-in period.
     data(ybar,ss)
     outdir <- BGmix(ybar,ss,c(8,8),niter=100,nburn=0,nthin=1,trace.pred=1)
     pred <- ccPred(outdir,q.trace=TRUE)
     params <- ccParams(outdir)
     plotMixDensity(params,pred,ybar,ss)

