paulino            package:hyperdirichlet            R Documentation

_D_a_t_a_s_e_t _u_s_e_d _b_y _P_a_u_l_i_n_o

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

     A hyperdirichlet distribution for the dataset considered by
     Paulino and de Braganca Pereira (1995).

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

     data(paulino)

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

     Paulino and de Braganca Pereiria considered 97 subjects who each
     were one of three risk types: low, medium, and high (or 1,
     2, 3)

     Of the 97 subjects, 51 were fully categorized: 14 were type 1, 17
     type 2, and 20 type 3; the remaining 46 were partly categorized:
     28 were of types 1 or 2, and 18 were of types 2 or 3.

     The likelihood function is thus

               p1^14.p2^17.p3^20.(p1+p2)^28.(p2+p3)^18


     and object 'paulino' gives the PDF.

     *Note* the object gives the _likelihood_, not the posterior PDF. 
     Paulino and de Branganca Pereira use prior information of a type
     not readily representable in the hyperdirichlet paradigm.

_S_o_u_r_c_e:

     C. D. M. Paulino and C. A. de Braganca Pereira 1995. _Bayesian
     Methods for Categorical Data Under Informative General Censoring_,
     Biometrika, volume 82, number 2, pages 439-446

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

     data(paulino)
     maximum_likelihood(paulino)

