The minb is the main function to achieve data-driven inflated values selection and identification of important predictors for multiple inflated count data modeling.
minb(X,y,pars_init=NULL,lambda1_set=NULL,lambda2_set=NULL,ntune=10,
maxiter=200,tol=1e-03,vrbs=FALSE)The minb returns a list containing the following components:
The vector of the estimated mixing proportions of the selected inflated values in the multiple-inflated negative binomial model.
The vector of selected inflated values.
The value of the dispersion parameter of Negative Binomial distribution.
The vector of estimated non-zero regression coefficients of the Negative Binomial distribution.
The design matrix, without an intercept, minb includes an intercept by default.
The response count vector.
an optional list containing the initial values for the corresponding component. See details.
A user supplied tuning sequence for inflated values. Typical usage is to have the program compute its own sequence based on ntune. Supplying a value of lambda1_set overrides this.
A user supplied tuning sequence for regression predictors. Typical usage is to have the program compute its own sequence based on ntune. Supplying a value of lambda2_set overrides this.
The number of the tuning parameter values, default to 10.
defines the maximum number of iterations, default to 200.
specifies the iteration convergency tolerance, default to 1e-03.
a boolean variable which defines whether to print the iteration details, default to FALSE.
Initialization values can be supplied or estimated by glm.fit (the default). To supply initial values of the parameters, start should be a list with elements "kappa", "omega","beta", and "phi" containing the starting values for the coefficients of the corresponding component of the model.
# This is an example of minb with simulated data
data(minb_SimuData)
X = minb_SimuData$X
y = minb_SimuData$y
result = minb(X=X,y=y,ntune=2)
result$beta
result$omega
Run the code above in your browser using DataLab