Learn R Programming

HRQoL (version 1.0)

summary.BBreg:

Description

summary.BBreg si the BBreg specific method for the generic function summary which summarizes objects returned by modelling functions.

Usage

# S3 method for BBreg
summary(object, ...)

Arguments

object
a BBreg class model.
...
for extra arguments.

Value

summary.BBreg returns an object of class "summary.BBreg".
coefficients
a table with all the relevant information about the significance of the regression coefficients of the model. It includes the estimations, the standard errors of the estimations, the test-statistics and the p-values.
psi.table
a table which includes the estimation and standard errors of the logarithm of the dispersion parameter of the conditional beta-binomial distribution.
deviance
deviance of the model.
df
degrees of freedom of the model.
null.deviance
null-deviance, deviance of the null model.
null.df
degrees of freedom of the null model.
Goodness.of.fit
p-value of the goodness of fit test.
iter
number of iterations in the estimation method.
X
the model matrix.
y
the dependent variable in the model.
balanced
if the response variable is balanced it returns "yes", otherwise "no".
m
number of trials in each binomial observation.
nObs
number of observations.
m
number of trials in each observation.
balanced
if the response beta-binomial variable is balanced it returns "yes", otherwise "no".
conv
convergence of the methodology. If the algorithm has converged it returns "yes", otherwise "no".

Details

summary.BBreg summarizes all the relevant information about the estimation of the parameters in a BBreg class model. The function performs statistical significance hypothesis about the estimated regression parameters based on the normal distribution of the estimates. It also performs a goodness of fit test based on the difference between the calculated deviance of the model and the null deviance or deviance of the null model, which it is suppose to follow a Chi-square distribution with degrees of freedom equal to the difference in degrees of freedom of the models.

References

Forcina A. & Franconi L. (1988): Regression analysis with Beta-Binomial distribution, Revista di Statistica Applicata, 21, 7-12 Najera-Zuloaga J., Lee D.-J. & Arostegui I. (2017): Comparison of beta-binomial regression model approaches to analyze health related quality of life data, Statistical Methods in Medical Research, DOI: 10.1177/0962280217690413

Examples

Run this code
# We simulate a covariate, fix the paramters of the beta-binomial 
# distribution and simulate a response variable.

# Then we apply the model, and try to get the same values.
set.seed(18)
k <- 1000
m <- 10
x <- rnorm(k,5,3)

beta <- c(-10,2)
p <- 1/(1+exp(-(beta[1]+beta[2]*x)))
phi <- 1.2

y <- rBB(k,m,p,phi)

model <- BBreg(y~x,m)
sum.model <- summary(model)

Run the code above in your browser using DataLab