Generalized estimating equations will be used to estimate the parameters. Observations from the same person will be treated as repeated measurements. A first-order autoregressive correlation structure is assumed. Should the estimated covariance matrix fail to converge then we will resort to a variance component correlation structure. The expected difference between the exposed and the nonexposed will be estimated and presented with a 95% CI, based on the empiric SE.