weightedScores (version 0.9.5.3)

iee: INDEPENDENT ESTIMATING EQUATIONS FOR BINARY AND COUNT REGRESSION

Description

Independent estimating equations for binary and count regression.

Usage

iee(xdat,ydat,margmodel,link)

Arguments

xdat

\((\mathbf{x}_1 , \mathbf{x}_2 , \ldots , \mathbf{x}_n )^\top\), where the matrix \(\mathbf{x}_i,\,i=1,\ldots,n\) for a given unit will depend on the times of observation for that unit (\(j_i\)) and will have number of rows \(j_i\), each row corresponding to one of the \(j_i\) elements of \(y_i\) and \(p\) columns where \(p\) is the number of covariates including the unit first column to account for the intercept. This xdat matrix is of dimension \((N\times p),\) where \(N =\sum_{i=1}^n j_i\) is the total number of observations from all units.

ydat

\((y_1 , y_2 , \ldots , y_n )^\top\), where the response data vectors \(y_i,\,i=1,\ldots,n\) are of possibly different lengths for different units. In particular, we now have that \(y_i\) is (\(j_i \times 1\)), where \(j_i\) is the number of observations on unit \(i\). The total number of observations from all units is \(N =\sum_{i=1}^n j_i\). The ydat are the collection of data vectors \(y_i, i = 1,\ldots,n\) one from each unit which summarize all the data together in a single, long vector of length \(N\).

margmodel

Indicates the marginal model. Choices are “poisson” for Poisson, “bernoulli” for Bernoulli, and “nb1” , “nb2” for the NB1 and NB2 parametrization of negative binomial in Cameron and Trivedi (1998).

link

The link function. Choices are “log” for the log link function, “logit” for the logit link function, and “probit” for the probit link function.

Value

A list containing the following components:

coef

The vector with the estimated regression parameters.

gam

The vector with the estimated parameters that are not regression parameters. This is NULL for Poisson and binary regression.

Details

The univariate parameters are estimated from the sum of univariate marginal log-likelihoods.

References

Cameron, A. C. and Trivedi, P. K. (1998) Regression Analysis of Count Data. Cambridge: Cambridge University Press.

See Also

marglik

Examples

# NOT RUN {
################################################################################
#                      read and set up data set
################################################################################
data(toenail)
# covariates
xdat<-cbind(1,toenail$treat,toenail$time,toenail$treat*toenail$time)
# response
ydat<-toenail$y
#id
id<-toenail$id
#time
tvec<-toenail$time
################################################################################
#                      select the marginal model
################################################################################
margmodel="bernoulli"
################################################################################
#                      perform the IEE method
################################################################################
i.est<-iee(xdat,ydat,margmodel)
cat("\niest: IEE estimates\n")
print(c(i.est$reg,i.est$gam))
# }