paulino {hyperdirichlet} | R Documentation

## Dataset used by Paulino

### Description

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

### Usage

    data(paulino)

### Details

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.

### Source

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

### Examples

    data(paulino)
    maximum_likelihood(paulino)

[Package hyperdirichlet version 1.3-8 Index]