Using biplots to interpret gene expression patterns in plants

Sum: Plant biologists in fields of ecology, evolution, genetics and breeding frequently use multivariate methods. This paper illustrates Principal Component Analysis (PCA) and Gabriel's biplot as applied to microarray expression data from plant pathology experiments. Availability: An example program in the publicly distributed statistical language R is available from the web site (www.tpp.uq.edu.au) and by e-mail from the contact. Contact: scott.chapman@csiro.au.

Relation of Aortic Valve Calcium Detected by Cardiac Computed Tomography to All-Cause Mortality

Aortic valve calcium (AVC) can be quantified on the same computed tomographic scan as coronary artery calcium (CAC). Although CAC is an established predictor of cardiovascular events, limited evidence is available for an independent predictive value for AVC. We studied a cohort of 8,401 asymptomatic subjects (mean age 53 10 years, 69% men), who were free of known coronary heart disease and were undergoing electron beam computed tomography for assessment of subclinical atherosclerosis. The patients were followed for a median of 5 years (range 1 to 7) for the occurrence of mortality from any cause. Multivariate Cox regression models were developed to predict all-cause mortality according to the presence of AVC. A total of 517 patients (6%) had AVC on electron beam computed tomography. During follow-up, 124 patients died (1.5%), for an overall survival rate of 96.1% and 98.7% for those with and without AVC, respectively (hazard ratio 3.39, 95% confidence interval 2.09 to 5.49). After adjustment for age, gender, hypertension, dyslipidemia, diabetes mellitus, smoking, and a family history of premature coronary heart disease, AVC remained a significant predictor of mortality (hazard ratio 1.82, 95% confidence interval 1.11 to 2.98). Likelihood ratio chi-square statistics demonstrated that the addition of AVC contributed significantly to the prediction of mortality in a model adjusted for traditional risk factors (chi-square = 5.03, p = 0.03) as well as traditional risk factors plus the presence of CAC (chi-square = 3.58, p = 0.05). In conclusion, AVC was associated with increased all-cause mortality, independent of the traditional risk factors and the presence of CAC. (C) 2010 Published by Elsevier Inc. (Am J Cardiol 2010;106:1787-1791)

High-breakdown linear discriminant analysis

The classification rules of linear discriminant analysis are defined by the true mean vectors and the common covariance matrix of the populations from which the data come. Because these true parameters are generally unknown, they are commonly estimated by the sample mean vector and covariance matrix of the data in a training sample randomly drawn from each population. However, these sample statistics are notoriously susceptible to contamination by outliers, a problem compounded by the fact that the outliers may be invisible to conventional diagnostics. High-breakdown estimation is a procedure designed to remove this cause for concern by producing estimates that are immune to serious distortion by a minority of outliers, regardless of their severity. In this article we motivate and develop a high-breakdown criterion for linear discriminant analysis and give an algorithm for its implementation. The procedure is intended to supplement rather than replace the usual sample-moment methodology of discriminant analysis either by providing indications that the dataset is not seriously affected by outliers (supporting the usual analysis) or by identifying apparently aberrant points and giving resistant estimators that are not affected by them.

Standard errors of fitted component means of normal mixtures

In this paper use consider the problem of providing standard errors of the component means in normal mixture models fitted to univariate or multivariate data by maximum likelihood via the EM algorithm. Two methods of estimation of the standard errors are considered: the standard information-based method and the computationally-intensive bootstrap method. They are compared empirically by their application to three real data sets and by a small-scale Monte Carlo experiment.

Monte Carlo evidence on cointegration and causation

The small sample performance of Granger causality tests under different model dimensions, degree of cointegration, direction of causality, and system stability are presented. Two tests based on maximum likelihood estimation of error-correction models (LR and WALD) are compared to a Wald test based on multivariate least squares estimation of a modified VAR (MWALD). In large samples all test statistics perform well in terms of size and power. For smaller samples, the LR and WALD tests perform better than the MWALD test. Overall, the LR test outperforms the other two in terms of size and power in small samples.

Survival and Recurrence in Nonmycosis Fungoides Primary Cutaneous Lymphoma

Purpose: To evaluate overall and relapse-free survival (RFS) in patients with nonmycosis fungoides (non-MF) primary cutaneous lymphoma (PCL). Methods: Thirty-eight patients with PCL excluding cases of MF treated between 1993 and 2006 were analyzed retrospectively. Survival statistics were estimated by the methods of Kaplan and Meier, and univariate and multivariate significance testing were performed by Cox regression analysis. Results: The median follow-up was 34.6 months (range, 2-138.3 months). The overall survival for the entire study population, at 5 and 10 years, was 97% and 78%, respectively. The RFS for the entire study population, at 5 and 10 years, was 30% and 22%, respectively. For those who received radiotherapy (RT) as a component of their initial therapy, the RFS at 5 and 10 years was 48% and 36%, respectively. Among those receiving RT who relapsed, the site of relapse was out-of-field in 82% of the cases. In our multivariate analysis, only RT as a component of the initial therapy and the absence of bulky disease had a statistically significant improvement in RFS (P = 0.01 and < 0.01, respectively). Conclusion: RT improves the local control and RFS of patients with non-MF PCL.

Elevated left and reduced right orbitomedial prefrontal fractional anisotropy in adults with bipolar disorder revealed by tract-based spatial statistics

Context Diffusion tensor imaging (DTI) studies in adults with bipolar disorder (BD) indicate altered white matter (WM) in the orbitomedial prefrontal cortex (OMPFC), potentially underlying abnormal prefrontal corticolimbic connectivity and mood dysregulatioin in BD. Objective: To use tract-based spatial statistics (TBSS) to examine VVM skeleton (ie, the most compact whole-brain WM) in subjects with BD vs healthy control subjects. Design: Cross-sectional, case-control, whole-brain DTI using TBSS. Setting: University research institute. Participants: Fifty-six individuals, 31 having a DSM-IV diagnosis of BD type 1 (mean age, 35.9 years [age range, 24-52 years]) and 25 controls (mean age, 29.5 years [age range, 19-52 years]). Main Outcome Measures: Fractional anisotropy (FA) longitudinal and radial diffusivities in subjects with BD vs controls (covarying for age) and their relationships with clinical and demographic variables. Results: Subjects with BD vs controls had significantly greater FA (t > 3.0, P <=.05 corrected) in the left uncinate fasciculus (reduced radial diffusivity distally and increased longitudinal diffusivity centrally), left optic radiation (increased longitudinal diffusivity), and right anterothalamic radiation (no significant diffusivity change). Subjects with BD vs controls had significantly reduced FA (t > 3.0, P <=.05 corrected) in the right uncinate fasciculus (greater radial diffusivity). Among subjects with BD, significant negative correlations (P <.01) were found between age and FA in bilateral uncinate fasciculi and in the right anterothalamic radiation, as well as between medication load and FA in the left optic radiation. Decreased FA (P <.01) was observed in the left optic radiation and in the right anterothalamic radiation among subjects with BD taking vs those not taking mood stabilizers, as well as in the left optic radiation among depressed vs remitted subjects with BD. Subjects having BD with vs without lifetime alcohol or other drug abuse had significantly decreased FA in the left uncinate fasciculus. Conclusions: To our knowledge, this is the first study to use TBSS to examine WM in subjects with BD. Subjects with BD vs controls showed greater WM FA in the left OMPFC that diminished with age and with alcohol or other drug abuse, as well as reduced WM FA in the right OMPFC. Mood stabilizers and depressed episode reduced WM FA in left-sided sensory visual processing regions among subjects with BD. Abnormal right vs left asymmetry in FA in OMPFC WM among subjects with BD, likely reflecting increased proportions of left-sided longitudinally aligned and right-sided obliquely aligned myelinated fibers, may represent a biologic mechanism for mood dysregulation in BD.

Long-term survivor mixture model with random effects: application to a multi-centre clinical trial of carcinoma

A mixture model incorporating long-term survivors has been adopted in the field of biostatistics where some individuals may never experience the failure event under study. The surviving fractions may be considered as cured. In most applications, the survival times are assumed to be independent. However, when the survival data are obtained from a multi-centre clinical trial, it is conceived that the environ mental conditions and facilities shared within clinic affects the proportion cured as well as the failure risk for the uncured individuals. It necessitates a long-term survivor mixture model with random effects. In this paper, the long-term survivor mixture model is extended for the analysis of multivariate failure time data using the generalized linear mixed model (GLMM) approach. The proposed model is applied to analyse a numerical data set from a multi-centre clinical trial of carcinoma as an illustration. Some simulation experiments are performed to assess the applicability of the model based on the average biases of the estimates formed. Copyright (C) 2001 John Wiley & Sons, Ltd.

Recent trends in pay equity: Beyond the aggregate statistics

The earnings gap between men and women has remained comparatively stable at an aggregate level over the 1990s in Australia. From one perspective, this is a reminder of the considerable difficulty of addressing wage differentials once the most overt forms of wage discrimination have been removed, and of the limited impact of most policy initiatives. From another, it may be seen as evidence that dire predictions about the effects of decentralisation on the earnings gap have failed to materialise. In this paper, I use Australian Bureau of Statistics data to show that a number of different trends are evident underneath the relatively static picture shown by the aggregate statistics, particularly as wage dispersion has increased. The data suggest not only that the prospects for pay equity are far from benign, but also that in the current labour market the issue of gender pay inequality cannot be effectively addressed separately from wage inequality more generally.