paulino {hyperdirichlet} | R Documentation |
A hyperdirichlet distribution for the dataset considered by Paulino and de Braganca Pereira (1995).
data(paulino)
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.
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
data(paulino) max.like(paulino)