BBest: Estimation of the parameters of a beta-binomial distribution

In PROreg: Patient Reported Outcomes Regression Analysis

Estimation of the parameters of a beta-binomial distribution

Description

This function performs the estimation of the parameters of a beta-binomial distribution for the given data and maximum score number in each observation.

There are two different approaches available for performing the estimation of the parameters: (i) Method of moments, and, (ii) maximum likelihood approach.

Usage

BBest(y,m,method="MM")

Arguments

y	response variable which folloes a beta-binomial distribution.
m	maximum score number in each beta-binomial observation.
method	the method used for performing the estimation of the probability and dispersion parameters of a beta-binomial distribution, "MM" for the method of moments and "MLE" for maximum likelihood estimation. Default "MM".

Details

BBest function performs estimations in the parameters of a beta-binomial distribution for the given data. The estimations can be performed using two different approaches, the methods of moments and the maximum likelihood estimation approach.

The density function of a given observation y that follows a beta-binomial distribution with paramters m, p and phi is defined as

f(y)=[\Gamma(1/phi)*\Gamma(p/phi+m)*\Gamma((1-p)/phi+m)]/[\Gamma(1/phi+m)*\Gamma(p/phi)*\Gamma((1-p)/phi)].

The first and second order moments are defined as

E[y]=mp

Var[y]=mp(1-p)[1+(n-1)phi/(1+phi)].

Hence, if y=(y_1,...,y_n) is the given data, we can conclude the method of moments from the previous as

p=E/m,

phi=[V-mp(1-p)]/[mp(1-p)m-V],

where E is the sample mean and V is the sample variance.

On the other hand, the maximum likelihood estimation of both parameters consits of solving the derivative of the log-likelihood defined by the density function with respect to each parameter and equaling them to zero. An iterative algorithm is needed for both parameter estimation as the score equations the parameters depend each other.

The variance of the estimation of the probability parameter of the beta-binomial distribution for the given data set is computed by the inverse of the Fisher information, i.e., the inverse of the negative second derivate of the log-likelihodd remplacing p by its estimation.

Value

BBest returns an object of class "BBest".

The function summary (i.e., summary.BBest) can be used to obtain or print a summary of the results.

p	estimated probability parameter of the beta-binomial distribution.
phi	estimated dispersion parameter of the beta-binomial distribution.
pVar	variance of the estimation of the probability parameter p.
psi	estimation of the logarithm of the dispersion parameter phi.
psiVar	variance of the estimation of the logarithm of the dispersion parameter psi.
m	maximum score number in each beta-binomial observation.
balanced	if the response variable is balanced it returns "yes", otherwise "no".
method	the used approach for performing the estimations.

Author(s)

J. Najera-Zuloaga

D.-J. Lee

I. Arostegui

References

Arostegui I., Nunez-Anton V. & Quintana J. M. (2006): Analysis of short-form-36 (SF-36): The beta-binomial distribution approach, Statistics in Medicine, 26, 1318-1342

Examples

# We simulate 1000 observations of a beta-binomial distribution
# for the fixed paramters.
m <- 10
k <- 1000
p <- 0.7
phi <- 1.6

set.seed(5)
y <- rBB(k,m,p,phi)

# Performing the estimation of the parameters

# Method of moments:
MM <- BBest(y,m)
MM

# Maximum likelihood approach
MLE <- BBest(y,m,method="MLE")
MLE

PROreg documentation built on May 29, 2024, 11:29 a.m.