**************************************************************************
*                  DYALOG APL WORKSPACE - ZIPGLMM.DWS                    *
**************************************************************************
*                            INSTRUCTIONS                                *
**************************************************************************

The Dyalog APL workspace "ZIPGLMM.DWS" obtains REML estimators in both 
logistic and truncated poisson parts for Zero-inflated Poisson regression
with random effects by using the ZIPREG function.  GLMM approach is 
adopted to obtain estimates and standard errors.

This file gives a brief instruction on the data format, the interaction
with the program, and the output format.

__________________________________________________________________________
Data File 
~~~~~~~~~
In order to activate the program, two data sets have to be inputed (one for
the logistic part [DATA1] and the other for the truncated Poisson part 
[DATA2]). 

Format of DATA1:
Column 1 :  subject number
Column 2 :  0-1 indicator for event occurrence
Column 3 :  Exposure time which act as offset (set all to 1 for same exposure)
Column 4+:  Covariates (need to specify an intercept column)

Format of DATA2:
Column 1 :  subject number for those non-zero events (has to be recoded)
Column 2 :  number of event occurence
Column 3 :  Exposure time which act as offset (set all to 1 for same exposure)
Column 4+:  Covariates (need to specify an intercept column)

__________________________________________________________________________
Run the Program
~~~~~~~~~~~~~~~
For example, the data files are named respectively "DATA1" and "DATA2". 
To run the program, type:
      
      DATA1 ZIPREG DATA2

__________________________________________________________________________
Outputs of the Program
~~~~~~~~~~~~~~~~~~~~~~

The following outputs will be given once the program is started.

***LOGISTIC PART***
Coefficients in logistic part
variance in logistic part
SE of coefficents
SE of variance
Standardised random effects
Deviance
Degree of freedom for goodness of fit
Deviance residuals

***TRUNCATED POISSON PART***
Coefficients in truncated Poisson part
variance in truncated Poisson part
SE of coefficents
SE of variance
Standardised random effects
Deviance 
Degree of freedom for goodness of fit
Deviance residuals