On a resampling approach for tests on the number of clusters with mixture model-based clustering of tissue samples

We consider the problem of assessing the number of clusters in a limited number of tissue samples containing gene expressions for possibly several thousands of genes. It is proposed to use a normal mixture model-based approach to the clustering of the tissue samples. One advantage of this approach is that the question on the number of clusters in the data can be formulated in terms of a test on the smallest number of components in the mixture model compatible with the data. This test can be carried out on the basis of the likelihood ratio test statistic, using resampling to assess its null distribution. The effectiveness of this approach is demonstrated on simulated data and on some microarray datasets, as considered previously in the bioinformatics literature. (C) 2004 Elsevier Inc. All rights reserved.

Statistical tests of load-unload response ratio signals by lattice solid model: Implication to tidal triggering and earthquake prediction

Statistical tests of Load-Unload Response Ratio (LURR) signals are carried in order to verify statistical robustness of the previous studies using the Lattice Solid Model (MORA et al., 2002b). In each case 24 groups of samples with the same macroscopic parameters (tidal perturbation amplitude A, period T and tectonic loading rate k) but different particle arrangements are employed. Results of uni-axial compression experiments show that before the normalized time of catastrophic failure, the ensemble average LURR value rises significantly, in agreement with the observations of high LURR prior to the large earthquakes. In shearing tests, two parameters are found to control the correlation between earthquake occurrence and tidal stress. One is, A/(kT) controlling the phase shift between the peak seismicity rate and the peak amplitude of the perturbation stress. With an increase of this parameter, the phase shift is found to decrease. Another parameter, AT/k, controls the height of the probability density function (Pdf) of modeled seismicity. As this parameter increases, the Pdf becomes sharper and narrower, indicating a strong triggering. Statistical studies of LURR signals in shearing tests also suggest that except in strong triggering cases, where LURR cannot be calculated due to poor data in unloading cycles, the larger events are more likely to occur in higher LURR periods than the smaller ones, supporting the LURR hypothesis.

Parameters affecting the performance of wetting and drying in a two-dimensional finite element long wave hydrodynamic model

The marsh porosity method, a type of thin slot wetting and drying algorithm in a two-dimensional finite element long wave hydrodynamic model, is discussed and analyzed to assess model performance. Tests, including comparisons to simple examples and theoretical calculations, examine the effects of varying the marsh porosity parameters. The findings demonstrate that the wetting and drying concept of marsh porosity, often used in finite element hydrodynamic modeling, can behave in a more complex manner than initially expected.

Genetic and environmental sources of covariance between reading tests used in neuropsychological assessment and IQ subtests

In this study, we examined genetic and environmental influences on covariation among two reading tests used in neuropsychological assessment (Cambridge Contextual Reading Test [CCRT], [Beardsall, L., and Huppert, F. A. ( 1994). J. Clin. Exp. Neuropsychol. 16: 232 - 242], Schonell Graded Word Reading Test [SGWRT], [ Schonell, F. J., and Schonell, P. E. ( 1960). Diagnostic and attainment testing. Edinburgh: Oliver and Boyd.]) and among a selection of IQ subtests from the Multidimensional Aptitude Battery (MAB), [Jackson, D. N. (1984). Multidimensional aptitude battery, Ontario: Research Psychologists Press.] and the Wechsler Adult Intelligence Scale-Revised (WAIS-R) [Wechsler, D. (1981). Manual for the Wechsler Adult Intelligence Scale-Revised (WAIS-R). San Antonio: The Psychological Corporation]. Participants were 225 monozygotic and 275 dizygotic twin pairs aged from 15 years to 18 years ( mean, 16 years). For Verbal IQ subtests, phenotypic correlations with the reading tests ranged from 0.44 to 0.65. For Performance IQ subtests, phenotypic correlations with the reading tests ranged from 0.23 to 0.34. Results of Structural Equation Modeling (SEM) supported a model with one genetic General factor and three genetic group factors ( Verbal, Performance, Reading). Reading performance was influenced by the genetic General factor ( accounting for 13% and 20% of the variance for the CCRT and SGWRT, respectively), the genetic Verbal factor ( explaining 17% and 19% of variance for the CCRT and SGWRT), and the genetic Reading factor ( explaining 21% of the variance for both the CCRT and SGWRT). A common environment factor accounted for 25% and 14% of the CCRT and SGWRT variance, respectively. Genetic influences accounted for more than half of the phenotypic covariance between the reading tests and each of the IQ subtests. The heritabilities of the CCRT and SGWRT were 0.54 and 0.65, respectively. Observable covariance between reading assessments used by neuropsychologists to estimate IQ and IQ subtests appears to be largely due to genetic effects.

HPGR model verification and scale-up

This paper summarises test results that were used to validate a model and scale-up procedure of the high pressure grinding roll (HPGR) which was developed at the JKMRC by Morrell et al. [Morrell, Lim, Tondo, David,1996. Modelling the high pressure grinding rolls. In: Mining Technology Conference, pp. 169-176.]. Verification of the model is based on results from four data sets that describe the performance of three industrial scale units fitted with both studded and smooth roll surfaces. The industrial units are currently in operation within the diamond mining industry and are represented by De Beers, BHP Billiton and Rio Tinto. Ore samples from the De Beers and BHP Billiton operations were sent to the JKMRC for ore characterisation and HPGR laboratory-scale tests. Rio Tinto contributed an historical data set of tests completed during a previous research project. The results conclude that the modelling of the HPGR process has matured to a point where the model may be used to evaluate new and to optimise existing comminution circuits. The model prediction of product size distribution is good and has been found to be strongly dependent of the characteristics of the material being tested. The prediction of throughput and corresponding power draw (based on throughput) is sensitive to inconsistent gap/diameter ratios observed between laboratory-scale tests and full-scale operations. (C) 2004 Elsevier Ltd. All rights reserved.

Noninvasive tests for arterial structure, function, and compliance: Do they identify risk or diagnose disease?

Background: Brachial artery reactivity (BAR), carotid intima-media thickness (IMT), and applanation tonometry for evaluation of total arterial compliance may provide information about preclinical vascular disease. We sought to determine whether these tests could be used to identify patients with coronary artery disease (CAD) without being influenced by their ability to identify those at risk ford CAD developing. Methods: We studied 100 patients and compared 3 groups: 35 patients with known CAD; 34 patients with symptoms and risk factors but no CAD identified by stress echocardiography (risk group); and 31 control subjects. BAR and IMT were measured using standard methods, and total arterial compliance was calculated by the pulse-pressure method from simultaneous radial applanation tonometry and pulsed wave Doppler of the left ventricular outflow. Ischemia was identified as a new or worsening wall-motion abnormality induced by stress. Results: In a comparison between the control subjects and patients either at risk for developing CAD or with CAD, the predictors of risk for CAD were: age (P = .01); smoking history (P = .002); hypercholesterolemia (P = .002); and hypertension (P = .004) (model R = 0.82; P = .0001). The independent predictors of CAD were: IMT (P = .001); BAR (P = .04); sex (P = .005); and hypertension (P = .005) (model R = 0.80; P = .0001). Conclusion: IMT, BAR, and traditional cardiovascular risk factors appear to identify patients at risk for CAD developing. However, only IMT was significantly different between patients at risk for developing CAD and those with overt CAD.

Colorectal cancer screening in Australia: an economic evaluation of a potential biennial screening program using faecal occult blood tests

Objective: To evaluate whether the introduction of a national, co-ordinated screening program using the faecal occult blood test represents 'value-for-money' from the perspective of the Australian Government as third-party funder. Methods: The annual equivalent costs and consequences of a biennial screening program in 'steady-state' operation were estimated for the Australian population using 1996 as the reference year. Disability-adjusted life years (DALYs) and the years of life lost (YLLs) averted, and the health service costs were modelled, based on the epidemiology and the costs of colorectal cancer in Australia together with the mortality reduction achieved in randomised controlled trials. Uncertainty in the model was examined using Monte Carlo simulation methods. Results: We estimate a minimum or 'base program' of screening those aged 55 to 69 years could avert 250 deaths per annum (95% uncertainty interval 99-400), at a gross cost of $A55 million (95% UI $A46 million to $A96 million) and a gross incremental cost-effectiveness ratio of $A17,000/DALY (95% UI $A13,000/DALY to $A52,000/DALY). Extending the program to include 70 to 74-year-olds is a more effective option (cheaper and higher health gain) than including the 50 to 54-year-olds. Conclusions: The findings of this study support the case for a national program directed at the 55 to 69-year-old age group with extension to 70 to 74-year-olds if there are sufficient resources. The pilot tests recently announced in Australia provide an important opportunity to consider the age range for screening and the sources of uncertainty, identified in the modelled evaluation, to assist decisions on implementing a full national program.

Testing the 'photoinhibition' model of coral bleaching using chemical inhibitors

Explants of the hard coral Seriatopora hystrix were exposed to sublethal concentrations of the herbicide diuron DCMU (N'-(3,4-dichlorophenyl,-N,N-dimethylurea)) and the heavy metal copper. Pulse amplitude modulated (PAM) chlorophyll fluorescence techniques were used to assess the effects on the photosynthetic efficiency of the algal symbionts in the tissue (in Symbio), and chlorophyll fluorescence and counts of symbiotic algae (normalised to surface area) were used to assess the extent of coral bleaching. At 30 mug DCMU l(-1), there was a reduction in both the maximum effective quantum yield (DeltaF/F-m') and maximum potential quantum yield (F-v/F-m) of the algal symbionts in symbio. Corals subsequently lost their algal symbionts and discoloured (bleached), especially on their upper sunlight-exposed surfaces. At the same DCMU concentration but under low light (5% of growth irradiance), there was a marked reduction in DeltaF/F-m' but only a slight reduction in F-v/F-m and slight loss of algae. Loss of algal symbionts was also noted after a 7 d exposure to concentrations as low as 10 mug DCMU l(-1) under normal growth irradiance, and after 14 d exposure to 10 mug DCMU l(-1) under reduced irradiance. Collectively the results indicate that DCMU-induced bleaching is caused by a light-dependent photoinactivation of algal symbionts, and that bleaching occurs when F-v/F-n, (measured 2 h after sunset) is reduced to a value of less than or equal to 0.6. Elevated copper concentrations (60 mug Cu l(-1) for 10 h) also induced a rapid bleaching in S. hystrix but without affecting the quantum yield of the algae in symbio. Tests with isolated algae indicated that substantially higher concentrations (300 mug Cu l(-1) for 8 h) were needed to significantly reduce the quantum yield. Thus, copper-induced bleaching occurs without affecting the algal photosynthesis and may be related to effects on the host (animal). It is argued that warm-water bleaching of corals resembles both types of chemically induced bleaching, suggesting the need for an integrated model of coral bleaching involving the effect of temperature on both host (coral) and algal symbionts.

Mixed model formulations for multi-environment trials

When studying genotype X environment interaction in multi-environment trials, plant breeders and geneticists often consider one of the effects, environments or genotypes, to be fixed and the other to be random. However, there are two main formulations for variance component estimation for the mixed model situation, referred to as the unconstrained-parameters (UP) and constrained-parameters (CP) formulations. These formulations give different estimates of genetic correlation and heritability as well as different tests of significance for the random effects factor. The definition of main effects and interactions and the consequences of such definitions should be clearly understood, and the selected formulation should be consistent for both fixed and random effects. A discussion of the practical outcomes of using the two formulations in the analysis of balanced data from multi-environment trials is presented. It is recommended that the CP formulation be used because of the meaning of its parameters and the corresponding variance components. When managed (fixed) environments are considered, users will have more confidence in prediction for them but will not be overconfident in prediction in the target (random) environments. Genetic gain (predicted response to selection in the target environments from the managed environments) is independent of formulation.

Unidirectional inversion of ibuprofen in Caco-2 cells: Developing a suitable model for presystemic chiral inversion study

Intestinal chiral inversion of ibuprofen is still lacking direct evidence. In a preliminary experiment, ibuprofen was found to undergo inversion in Caco-2 cells. This investigation was thus conducted to determine the characteristics and influence of some biochemical factors on the chiral inversion of ibuprofen in Caco-2 cells. The effects of substrate concentration (2.5-40 mu g/ml), cell density (0.5-2 x 10(6) cells/ well), content of serum (0-20%), coexistence of S ibuprofen (corresponding doses), sodium azide (10mm), exogenous Coenzyme A (CoA) (0.1 - 0.4 mm),. and palmitic acid (5-25 mu m) on inversion were examined. A stereoselective HPLC method based on the Chromasil-CHI-TBB column was developed for quantitative analysis of the drug in cell culture medium. The inversion ratio (F-i) and elimination rate constant were calculated as the indexes of inversion extent. Inversion of ibuprofen in Caeo-2 cells was found to be both dose and cell density dependent, indicating saturable characteristics. Addition of serum significantly inhibited the inversion, to an extent of 2.7 fold decrease at 20% content. Preexistence of S enantiomer exerted a significant inhibitory effect (p < 0.01 for all tests). Sodium azide decreased the inversion ratio from 0.43 to 0.32 (p < 0.01). Exogenous CoA and palmitic acid significantly promoted the inversion at all tested doses (p < 0.01 for all tests). This research provided strong evidence to the capacity and capability of intestinal chiral inversion. Although long incubation times up to 120 h were required, Caco-2 cells should be a suitable model for chiral inversion research of 2-APAs considering the human-resourced and well-defined characteristics from the present study.

Modelling of two-stage anaerobic digestion using the IWA Anaerobic Digestion Model No. 1 (ADM1)

The aim of the study presented was to implement a process model to simulate the dynamic behaviour of a pilot-scale process for anaerobic two-stage digestion of sewage sludge. The model implemented was initiated to support experimental investigations of the anaerobic two-stage digestion process. The model concept implemented in the simulation software package MATLAB(TM)/Simulink(R) is a derivative of the IWA Anaerobic Digestion Model No.1 (ADM1) that has been developed by the IWA task group for mathematical modelling of anaerobic processes. In the present study the original model concept has been adapted and applied to replicate a two-stage digestion process. Testing procedures, including balance checks and 'benchmarking' tests were carried out to verify the accuracy of the implementation. These combined measures ensured a faultless model implementation without numerical inconsistencies. Parameters for both, the thermophilic and the mesophilic process stage, have been estimated successfully using data from lab-scale experiments described in literature. Due to the high number of parameters in the structured model, it was necessary to develop a customised procedure that limited the range of parameters to be estimated. The accuracy of the optimised parameter sets has been assessed against experimental data from pilot-scale experiments. Under these conditions, the model predicted reasonably well the dynamic behaviour of a two-stage digestion process in pilot scale. (C) 2004 Elsevier Ltd. All rights reserved.

Implementation of particle-scale rotation in the 3-d Lattice Solid Model

In this study, 3-D Lattice Solid Model (LSMearth or LSM) was extended by introducing particle-scale rotation. In the new model, for each 3-D particle, we introduce six degrees of freedom: Three for translational motion, and three for orientation. Six kinds of relative motions are permitted between two neighboring particles, and six interactions are transferred, i.e., radial, two shearing forces, twisting and two bending torques. By using quaternion algebra, relative rotation between two particles is decomposed into two sequence-independent rotations such that all interactions due to the relative motions between interactive rigid bodies can be uniquely decided. After incorporating this mechanism and introducing bond breaking under torsion and bending into the LSM, several tests on 2-D and 3-D rock failure under uni-axial compression are carried out. Compared with the simulations without the single particle rotational mechanism, the new simulation results match more closely experimental results of rock fracture and hence, are encouraging. Since more parameters are introduced, an approach for choosing the new parameters is presented.