NegBinBetaBinreg: NegBinBetaBinreg

Description Usage Arguments Details Value Author(s) References Examples

Description

Function to estimate a Negative Binomial regression models with mean and shape (or variance) regression structures, and Beta Binomial regression with mean and dispersion regression structures.

Usage

1
2
NegBinBetaBinreg(y,x,z,nsim,bpri,Bpri,
gpri,Gpri,burn,jump,bini,gini,model,m,ni,graph1,graph2)

Arguments

y

object of class matrix or vector, with the dependent variable.

x

object of class matrix, with the variables for modelling the mean.

z

object of class matrix, with the variables for modelling the shape, variance or dispersion.

nsim

a number that indicate the number of iterations.

bpri

a vector with the prior values of beta.

Bpri

a matrix with the prior values of the variance of beta.

gpri

a vector with the prior values of gamma.

Gpri

a matrix with the prior values of the variance of gamma.

burn

a proportion that indicate the number of iterations to be burn at the beginning of the chain.

jump

a number that indicate the distance between samples of the autocorrelated the chain, to be excluded from the final chain.

bini

a vector with the initial values of beta.

gini

a vector with the initial values of gamma.

model

it indicates the model that will be used. By default, is the Beta Binomial model (BB), but it could also be the Negative Binomial with mean and shape (NB1) or the Negative Binomial with mean and variance (NB2).

m

Is positive integer that In the Beta Binomial model indicates the number of trials. By default, is the number of data

ni

Is a vector of positive integer that In the Beta Binomial model indicates the number of trials to each individual. By default, is a vector of m

graph1

if it is TRUE present the graph of the chains without jump and burn.

graph2

if it is TRUE present the graph of the chains with jump and burn.

Details

The Bayesian Negative Binomial regression allow the joint modelling of mean and shape or variance of a negative binomial distributed variable, as is proposed in Cepeda (2001), with exponential link for the mean and the shape or variance. The Bayesian Beta Binomial regression allow the joint modelling of mean and precision of a beta binomial distributed variable, as is proposed in Cepeda (2001), with logit link for the mean and exponential link for the precision.

Value

object of class NegBinBetaBinreg with:

coefficients

object of class matrix with the estimated coefficients of beta and gamma.

desv

object of class matrix with the estimated desviations of beta and gamma.

interv

object of class matrix with the estimated confidence intervals of beta and gamma.

fitted.values

object of class matrix with the fitted values of y.

residuals

object of class matrix with the residuals of the regression.

estresiduals

object of class matrix with the standardized residuals of the regression.

beta.mcmc

object of class matrix with the complete chains for beta.

gamma.mcmc

object of class matrix with the complete chains for gamma.

beta.mcmc.short

object of class matrix with the chains for beta after the burned process.

gamma.mcmc.short

object of class matrix with the chains for gamma after the burned process.

aceptbeta

object of class integer with the acceptance rate for the beta values.

aceptgamma

object of class integer with the acceptance rate for the gamma values.

call

Call.

Author(s)

Edilberto Cepeda-Cuervo ecepedac@unal.edu.co, Maria Victoria Cifuentes-Amado mvcifuentesa@unal.edu.co, Margarita Marin mmarinj@unal.edu.co

References

1. Cepeda C. E. (2001). Modelagem da variabilidade em modelos lineares generalizados. Unpublished Ph.D. tesis. Instituto de Matematicas. Universidade Federal do Rio do Janeiro. //http://www.docentes.unal.edu.co/ecepedac/docs/MODELAGEM20DA20VARIABILIDADE.pdf. http://www.bdigital.unal.edu.co/9394/. 2.Cepeda, E. C. and Gamerman D. (2005). Bayesian Methodology for modeling parameters in the two-parameter exponential family. Estadistica 57, 93 105. // 3.Cepeda, E. and Garrido, L. (2011). Bayesian beta regression models: joint mean and precision modeling. Universidad Nacional // 4.Cepeda, E. and Migon, H. and Garrido, L. and Achcar, J. (2012) Generalized Linear models with random effects in the two parameter exponential family. Journal of Statistical Computation and Simulation. 1, 1 13. // 5.Cepeda-Cuervo, E. and Cifuentes-Amado, V. (2016) Double generalized beta-binomial and negative binomial regression. To appear.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
rm(list=ls(all=TRUE))

Y<-c(6,6,9,13,23,25,32,53,54,5,5,11,17,19,2,8,13,14,20,47,
     48,60,81,6,17,67,0,0,2,7,11,12,0,0,5,5,5,11,17,3,4,22,
     30,36,0,1,5,7,8,16,27,25,10,11,20,33,0,1,5,5,5,5,5,7,7,11,15,5,6,6,7,14
)
y <- Y <- Y[1:68]

x0<-rep(1,times=68)
x2<-c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,
      1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
      1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1
)
x3<-c(0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,
      0,1,1,1,1,1,1,0,0,0,0,0,0,0,1,1,0,0,0,1,1,1,1,1,1,1,
      0,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,1
)
x<-cbind(x0,x2,x3)
z0<-rep(1,times=68)
z<-cbind(z0,x2)

Bpri=diag(rep(1,3))
bpri=rep(0,3)
Gpri=diag(rep(1,2))
gpri=rep(0,2)

Bini=diag(rep(1,3))
bini=c(3,-1,-0.5)
Gini=diag(rep(1,2))
gini=c(3,-1)

nsim = 300
burn <- 0.1
jump <- 5
model <- "NB1"
m <- 360
ni <- NULL
re<- NegBinBetaBinregEst (y,x,z,nsim,bpri,
	Bpri,gpri,Gpri,burn,jump,bini,gini,
	model,m,ni,graph1=FALSE,graph2=FALSE)
summary(re)

NegBinBetaBinreg documentation built on May 2, 2019, 10:52 a.m.