MCbinomialbeta           package:MCMCpack           R Documentation

_M_o_n_t_e _C_a_r_l_o _S_i_m_u_l_a_t_i_o_n _f_r_o_m _a _B_i_n_o_m_i_a_l _L_i_k_e_l_i_h_o_o_d _w_i_t_h _a _B_e_t_a _P_r_i_o_r

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

     This function generates a sample from the posterior distribution
     of a binomial likelihood with a Beta prior.

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

     MCbinomialbeta(y, n, alpha=1, beta=1, mc=1000, ...)

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

       y: The number of successes in the independent Bernoulli trials.

       n: The number of independent Bernoulli trials.

   alpha: Beta prior distribution alpha parameter.

    beta: Beta prior distribution beta parameter.

      mc: The number of Monte Carlo draws to make.

     ...: further arguments to be passed

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

     'MCbinomialbeta' directly simulates from the posterior
     distribution.  This model is designed primarily for instructional
     use.  pi is the probability of success for each independent
     Bernoulli trial.  We assume a conjugate Beta prior:

                        pi ~ Beta(alpha, beta)

     y is the number of successes in n trials. By default, a uniform
     prior is used.

_V_a_l_u_e:

     An mcmc object that contains the posterior sample.  This  object
     can be summarized by functions provided by the coda package.

_S_e_e _A_l_s_o:

     'plot.mcmc', 'summary.mcmc'

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

     ## Not run: 
     posterior <- MCbinomialbeta(3,12,mc=5000)
     summary(posterior)
     plot(posterior)
     grid <- seq(0,1,0.01)
     plot(grid, dbeta(grid, 1, 1), type="l", col="red", lwd=3, ylim=c(0,3.6),
       xlab="pi", ylab="density")
     lines(density(posterior), col="blue", lwd=3)
     legend(.75, 3.6, c("prior", "posterior"), lwd=3, col=c("red", "blue"))
     ## End(Not run)

