Kernel Estimation of Rate Function for Recurrent Event Data

Recurrent event data are largely characterized by the rate function but smoothing techniques for estimating the rate function have never been rigorously developed or studied in statistical literature. This paper considers the moment and least squares methods for estimating the rate function from recurrent event data. With an independent censoring assumption on the recurrent event process, we study statistical properties of the proposed estimators and propose bootstrap procedures for the bandwidth selection and for the approximation of confidence intervals in the estimation of the occurrence rate function. It is identified that the moment method without resmoothing via a smaller bandwidth will produce curve with nicks occurring at the censoring times, whereas there is no such problem with the least squares method. Furthermore, the asymptotic variance of the least squares estimator is shown to be smaller under regularity conditions. However, in the implementation of the bootstrap procedures, the moment method is computationally more efficient than the least squares method because the former approach uses condensed bootstrap data. The performance of the proposed procedures is studied through Monte Carlo simulations and an epidemiological example on intravenous drug users.

Seasonal Analyses of Air Pollution and Mortality in 100 U.S. Cities

Time series models relating short-term changes in air pollution levels to daily mortality counts typically assume that the effects of air pollution on the log relative rate of mortality do not vary with time. However, these short-term effects might plausibly vary by season. Changes in the sources of air pollution and meteorology can result in changes in characteristics of the air pollution mixture across seasons. The authors develop Bayesian semi-parametric hierarchical models for estimating time-varying effects of pollution on mortality in multi-site time series studies. The methods are applied to the updated National Morbidity and Mortality Air Pollution Study database for the period 1987--2000, which includes data for 100 U.S. cities. At the national level, a 10 micro-gram/m3 increase in PM(10) at lag 1 is associated with a 0.15 (95% posterior interval: -0.08, 0.39),0.14 (-0.14, 0.42), 0.36 (0.11, 0.61), and 0.14 (-0.06, 0.34) percent increase in mortality for winter, spring, summer, and fall, respectively. An analysis by geographical regions finds a strong seasonal pattern in the northeast (with a peak in summer) and little seasonal variation in the southern regions of the country. These results provide useful information for understanding particle toxicity and guiding future analyses of particle constituent data.

On Time Series Analysis of Public Health and Biomedical Data

A time series is a sequence of observations made over time. Examples in public health include daily ozone concentrations, weekly admissions to an emergency department or annual expenditures on health care in the United States. Time series models are used to describe the dependence of the response at each time on predictor variables including covariates and possibly previous values in the series. Time series methods are necessary to account for the correlation among repeated responses over time. This paper gives an overview of time series ideas and methods used in public health research.

Self-Reported Memory Symptoms with Coronary Artery Disease: A Prospective of CABG Patients and Nonsurgical Controls

Background. Subjective memory complaints are common after coronary artery bypass grafting (CABG), but previous studies have concluded that such symptoms are more closely associated with depressed mood than objective cognitive dysfunction. We compared the incidence of self-reported memory symptoms at 3 and 12 months after CABG with that of a control group of patients with comparable risk factors for coronary artery disease but without surgery. Methods. Patients undergoing CABG (n = 140) and a demographically similar nonsurgical control group with coronary artery disease (n = 92) were followed prospectively at 3 and 12 months. At each follow-up time, participants were asked about changes since the previous evaluation in areas of memory, calculations, reading, and personality. A Functional Status Questionnaire (FSQ) and self-report measure of symptoms of depression (CES-D) were also completed. Results. The frequency of self-reported changes in memory, personality, and reading at 3 months was significantly higher among CABG patients than among nonsurgical controls. By contrast, there were no differences in the frequency of self-reported symptoms relating to calculations or overall rating of functional status. After adjusting for a measure of depression (CES-D rating score), the risk for self-reported memory changes remained nearly 5 times higher among the CABG patients than control subjects. The relative risk of developing new self-reported memory symptoms between 3 and 12 months was 2.5 times higher among CABG patients than among nonsurgical controls (CI 1.24 – 5.02), and the overall prevalence of memory symptoms at 12 months was also higher among CABG patients (39%) than controls (14%). Conclusions. The frequency of self-reported memory symptoms 3 and 12 months after baseline is significantly higher among CABG patients than control patients with comparable risk factors for coronary and cerebrovascular disease. These differences could not be accounted for by symptoms of depression. The self-reported cognitive symptoms appear to be relatively specific for memory, and may reflect aspects of memory functioning that are not captured by traditional measures of new verbal learning and memory. The etiology of these self-reported memory symptoms remains unclear, but our findings as well as those of others, may implicate factors other than cardiopulmonary bypass itself.

ON THE MERITS OF VOXEL-BASED MORPHOMETRIC PATH-ANALYSIS FOR INVESTIGATING VOLUMETRIC MEDIATION OF A TOXICANT'S INFLUENCE ON COGNITIVE FUNCTION

We previously showed that lifetime cumulative lead dose, measured as lead concentration in the tibia bone by X-ray fluorescence, was associated with persistent and progressive declines in cognitive function and with decreases in MRI-based brain volumes in former lead workers. Moreover, larger region-specific brain volumes were associated with better cognitive function. These findings motivated us to explore a novel application of path analysis to evaluate effect mediation. Voxel-wise path analysis, at face value, represents the natural evolution of voxel-based morphometry methods to answer questions of mediation. Application of these methods to the former lead worker data demonstrated potential limitations in this approach where there was a tendency for results to be strongly biased towards the null hypothesis (lack of mediation). Moreover, a complimentary analysis using anatomically-derived regions of interest volumes yielded opposing results, suggesting evidence of mediation. Specifically, in the ROI-based approach, there was evidence that the association of tibia lead with function in three cognitive domains was mediated through the volumes of total brain, frontal gray matter, and/or possibly cingulate. A simulation study was conducted to investigate whether the voxel-wise results arose from an absence of localized mediation, or more subtle defects in the methodology. The simulation results showed the same null bias evidenced as seen in the lead workers data. Both the lead worker data results and the simulation study suggest that a null-bias in voxel-wise path analysis limits its inferential utility for producing confirmatory results.

NONLINEAR TUBE-FITTING FOR THE ANALYSIS OF ANATOMICAL AND FUNCTIONAL STRUCTURES

We are concerned with the estimation of the exterior surface of tube-shaped anatomical structures. This interest is motivated by two distinct scientific goals, one dealing with the distribution of HIV microbicide in the colon and the other with measuring degradation in white-matter tracts in the brain. Our problem is posed as the estimation of the support of a distribution in three dimensions from a sample from that distribution, possibly measured with error. We propose a novel tube-fitting algorithm to construct such estimators. Further, we conduct a simulation study to aid in the choice of a key parameter of the algorithm, and we test our algorithm with validation study tailored to the motivating data sets. Finally, we apply the tube-fitting algorithm to a colon image produced by single photon emission computed tomography (SPECT)and to a white-matter tract image produced using diffusion tensor `imaging (DTI).

Occurrence and Origin of Placer Gold in Montana

From time immemorial man has used gold as a medium of exchange, a mea­sure of value, as jewelry and for ornamentation. Placer gold has led dir­ectly or indirectly to the settlement of lands, California and Alaska being the two best known examples. It has led the way to the discovery of other important mineral wealth, the discovery of the copper and silver deposits at Butte, Montana and the discovery of the silver deposits at Cripple Creek, Colorado being two good examples.

Spatial Ecology of the Great Basin Rattlesnake on the Snake River Plain

Dr. Charles Peterson, Herpetologist and Professor of Biological Sciences at Idaho State University. Curator of Ichthyology and Herpetology at the Idaho Museum of Natural History.