On the Accelerated Failure Time Model for Current Status and Interval Censored Data

This paper introduces a novel approach to making inference about the regression parameters in the accelerated failure time (AFT) model for current status and interval censored data. The estimator is constructed by inverting a Wald type test for testing a null proportional hazards model. A numerically efficient Markov chain Monte Carlo (MCMC) based resampling method is proposed to simultaneously obtain the point estimator and a consistent estimator of its variance-covariance matrix. We illustrate our approach with interval censored data sets from two clinical studies. Extensive numerical studies are conducted to evaluate the finite sample performance of the new estimators.

A Pseudolikelihood Approach for Simultaneous Analysis of Array Comparative Genomic Hybridizations (aCGH)

DNA sequence copy number has been shown to be associated with cancer development and progression. Array-based Comparative Genomic Hybridization (aCGH) is a recent development that seeks to identify the copy number ratio at large numbers of markers across the genome. Due to experimental and biological variations across chromosomes and across hybridizations, current methods are limited to analyses of single chromosomes. We propose a more powerful approach that borrows strength across chromosomes and across hybridizations. We assume a Gaussian mixture model, with a hidden Markov dependence structure, and with random effects to allow for intertumoral variation, as well as intratumoral clonal variation. For ease of computation, we base estimation on a pseudolikelihood function. The method produces quantitative assessments of the likelihood of genetic alterations at each clone, along with a graphical display for simple visual interpretation. We assess the characteristics of the method through simulation studies and through analysis of a brain tumor aCGH data set. We show that the pseudolikelihood approach is superior to existing methods both in detecting small regions of copy number alteration and in accurately classifying regions of change when intratumoral clonal variation is present.

Reduced Bayesian Hierarchical Models: Estimating Health Effects of Simultaneous Exposure to Multiple Pollutants

Quantifying the health effects associated with simultaneous exposure to many air pollutants is now a research priority of the US EPA. Bayesian hierarchical models (BHM) have been extensively used in multisite time series studies of air pollution and health to estimate health effects of a single pollutant adjusted for potential confounding of other pollutants and other time-varying factors. However, when the scientific goal is to estimate the impacts of many pollutants jointly, a straightforward application of BHM is challenged by the need to specify a random-effect distribution on a high-dimensional vector of nuisance parameters, which often do not have an easy interpretation. In this paper we introduce a new BHM formulation, which we call "reduced BHM", aimed at analyzing clustered data sets in the presence of a large number of random effects that are not of primary scientific interest. At the first stage of the reduced BHM, we calculate the integrated likelihood of the parameter of interest (e.g. excess number of deaths attributed to simultaneous exposure to high levels of many pollutants). At the second stage, we specify a flexible random-effect distribution directly on the parameter of interest. The reduced BHM overcomes many of the challenges in the specification and implementation of full BHM in the context of a large number of nuisance parameters. In simulation studies we show that the reduced BHM performs comparably to the full BHM in many scenarios, and even performs better in some cases. Methods are applied to estimate location-specific and overall relative risks of cardiovascular hospital admissions associated with simultaneous exposure to elevated levels of particulate matter and ozone in 51 US counties during the period 1999-2005.

Analyzing Bivariate Survival Data with Interval Sampling and Application to Cancer Epidemiology

In medical follow-up studies, ordered bivariate survival data are frequently encountered when bivariate failure events are used as the outcomes to identify the progression of a disease. In cancer studies interest could be focused on bivariate failure times, for example, time from birth to cancer onset and time from cancer onset to death. This paper considers a sampling scheme where the first failure event (cancer onset) is identified within a calendar time interval, the time of the initiating event (birth) can be retrospectively confirmed, and the occurrence of the second event (death) is observed sub ject to right censoring. To analyze this type of bivariate failure time data, it is important to recognize the presence of bias arising due to interval sampling. In this paper, nonparametric and semiparametric methods are developed to analyze the bivariate survival data with interval sampling under stationary and semi-stationary conditions. Numerical studies demonstrate the proposed estimating approaches perform well with practical sample sizes in different simulated models. We apply the proposed methods to SEER ovarian cancer registry data for illustration of the methods and theory.

The neurophysiological time pattern of illusionary visual perceptual transitions: a simultaneous EEG and fMRI study

Several divergent cortical mechanisms generating multistability in visual perception have been suggested. Here, we investigated the neurophysiologic time pattern of multistable perceptual changes by means of a simultaneous recording with electroencephalography (EEG) and functional magnetic resonance imaging (fMRI). Volunteers responded to the subjective perception of a sudden change between stable patterns of illusionary motion (multistable transition) during a stroboscopic paradigm. We found a global deceleration of the EEG frequency prior to a transition and an occipital-accentuated acceleration after a transition, as obtained by low-resolution electromagnetic tomography analysis (LORETA) analysis. A decrease in BOLD response was found in the prefrontal cortex before, and an increase after the transitions was observed in the right anterior insula, the MT/V5 regions and the SMA. The thalamus and left superior temporal gyrus showed a pattern of decrease before and increase after transitions. No such temporal course was found in the control condition. The multimodal approach of data acquisition allows us to argue that the top-down control of illusionary visual perception depends on selective attention, and that a diminution of vigilance reduces selective attention. These are necessary conditions to allow for the occurrence of a perception discontinuity in absence of a physical change of the stimulus.

Measurement of the 1s-2s energy interval in muonium

The 1s-2s interval has been measured in the muonium (;mgr;(+)e(-)) atom by Doppler-free two-photon pulsed laser spectroscopy. The frequency separation of the states was determined to be 2 455 528 941.0(9.8) MHz, in good agreement with quantum electrodynamics. The result may be interpreted as a measurement of the muon-electron charge ratio as -1-1.1(2.1)x10(-9). We expect significantly higher accuracy at future high flux muon sources and from cw laser technology.

Simultaneous determination of 3-hydroxyanthranilic and cinnabarinic acid by high-performance liquid chromatography with photometric or electrochemical detection

A convenient and rapid method for the simultaneous determination by HPLC of 3-hydroxyanthranilic acid and the dimer derived by its oxidation, cinnabarinic acid, is described. Buffers or biological samples containing these two Trp metabolites were acidified to pH 2.0 and extracted with ethyl acetate with recoveries of 96.5 +/- 0.5 and 93.4 +/- 3.7% for 3-hydroxyanthranilic and cinnabarinic acid, respectively. The two compounds were separated on a reversed-phase (C18) column combined with ion-pair chromatography and detected photometrically or electrochemically. The method was applied successfully to biological systems in which formation of either 3-hydroxyanthranilic or cinnabarinic acid had been described previously. Thus, interferon-gamma-treated human peripheral blood mononuclear cells formed and released significant amounts of 3-hydroxyanthranilic acid into the culture medium and mouse liver nuclear fraction possessed high "cinnabarinic acid synthase" activity. In contrast, addition of 3-hydroxyanthranilic acid to human erythrocytes resulted in only marginal formation of cinnabarinic acid. We conclude that the method described is specific, sensitive, and suitable for the detection of the two Trp metabolites in biological systems.

The role of epistasis in the manifestation of heterosis: a systems-oriented approach

Heterosis is widely used in breeding, but the genetic basis of this biological phenomenon has not been elucidated. We postulate that additive and dominance genetic effects as well as two-locus interactions estimated in classical QTL analyses are not sufficient for quantifying the contributions of QTL to heterosis. A general theoretical framework for determining the contributions of different types of genetic effects to heterosis was developed. Additive x additive epistatic interactions of individual loci with the entire genetic background were identified as a major component of midparent heterosis. On the basis of these findings we defined a new type of heterotic effect denoted as augmented dominance effect di* that comprises the dominance effect at each QTL minus half the sum of additive x additive interactions with all other QTL. We demonstrate that genotypic expectations of QTL effects obtained from analyses with the design III using testcrosses of recombinant inbred lines and composite-interval mapping precisely equal genotypic expectations of midparent heterosis, thus identifying genomic regions relevant for expression of heterosis. The theory for QTL mapping of multiple traits is extended to the simultaneous mapping of newly defined genetic effects to improve the power of QTL detection and distinguish between dominance and overdominance.

The performance of computer-aided detection when analyzing prior mammograms of newly detected breast cancers with special focus on the time interval from initial imaging to detection

PURPOSE: The clinical role of CAD systems to detect breast cancer, which have not been on cancer containing mammograms not detected by the radiologist was proven retrospectively. METHODS: All patients from 1992 to 2005 with a histologically verified malignant breast lesion and a mammogram at our department, were analyzed in retrospect focussing on the time of detection of the malignant lesion. All prior mammograms were analyzed by CAD (CADx, USA). The resulting CAD printout was matched with the cancer containing images yielding to the radiological diagnosis of breast cancer. CAD performance, sensitivity as well as the association of CAD and radiological features were analyzed. RESULTS: 278 mammograms fulfilled the inclusion criteria. 111 cases showed a retrospectively visible lesion (71 masses, 23 single microcalcification clusters, 16 masses with microcalcifications, in one case two microcalcification clusters). 54/87 masses and 34/41 microcalcifications were detected by CAD. Detection rates varied from 9/20 (ACR 1) to 5/7 (ACR 4) (45% vs. 71%). The detection of microcalcifications was not influenced by breast tissue density. CONCLUSION: CAD might be useful in an earlier detection of subtle breast cancer cases, which might remain otherwise undetected.