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.

Early-Mid Holocene climatic variations in Tasmania, Australia: multi-proxy records in a stalagmite from Lynds Cave

Mass spectrometric uranium-series dating and C-O isotopic analysis of a stalagmite from Lynds Cave, northern Tasmania, Australia provide a high-resolution record of regional climate change between 5100 and 9200 yr before present (BP). Combined delta(18)O, delta(13)C, growth rate, initial U-234/U-238 and physical property (color, transparency and porosity) records allow recognition of seven climatic stages: Stage I ( > 9080 yr BP) - a relatively dry period at the beginning of stalagmite growth evidenced by elevated U-234/U-238; Stage II (9080-8600 yr BP) - a period of unstable climate characterized by high-frequency variability in temperature and bio-productivity; Stage 111 (8600-8000 yr BP) - a period of stable and moderate precipitation and stable and high bio-productivity, with a continuously rising temperature; Stage IV (8000-7400 yr BP) - the warmest period with high evaporation and low effective precipitation (rainfall less evaporation); Stage V (7400-7000 yr BP) - the wettest period with highest stalagmite growth and enhanced but unstable bio-productivity; Stage VI (7000-6600 yr BP) - a period with a significantly reduced precipitation and bio-productivity without noticeable change in temperature; Stage VII (6600-5100 yr BP) - a period of lowest temperature and precipitation marking a significant climatic deterioration. Overall, the records suggest that the warmest climate occurred between 8000 and 7400 yr BP, followed by a wettest period between 7400 and 7000 yr BP. These are broadly correlated with the so-called 'Mid Holocene optimum' previously proposed using pollen and lake level records. However, the timing and resolution of the speleothem. record from Lynds Cave are significantly higher than in both the pollen and lake level records. This allows us to correlate the abrupt change in physical property, delta(18)O, delta(13)C, growth rate, and initial U-234/U-238 of the stalagmite at similar to8000 yr BP with a global climatic event at Early-Mid Holocene transition. (C) 2001 Elsevier Science B.V. All rights reserved.

Integrating high and moderate spatial resolution image data to estimate forest age structure

The collection of spatial information to quantify changes to the state and condition of the environment is a fundamental component of conservation or sustainable utilization of tropical and subtropical forests, Age is an important structural attribute of old-growth forests influencing biological diversity in Australia eucalypt forests. Aerial photograph interpretation has traditionally been used for mapping the age and structure of forest stands. However this method is subjective and is not able to accurately capture fine to landscape scale variation necessary for ecological studies. Identification and mapping of fine to landscape scale vegetative structural attributes will allow the compilation of information associated with Montreal Process indicators lb and ld, which seek to determine linkages between age structure and the diversity and abundance of forest fauna populations. This project integrated measurements of structural attributes derived from a canopy-height elevation model with results from a geometrical-optical/spectral mixture analysis model to map forest age structure at a landscape scale. The availability of multiple-scale data allows the transfer of high-resolution attributes to landscape scale monitoring. Multispectral image data were obtained from a DMSV (Digital Multi-Spectral Video) sensor over St Mary's State Forest in Southeast Queensland, Australia. Local scene variance levels for different forest tapes calculated from the DMSV data were used to optimize the tree density and canopy size output in a geometric-optical model applied to a Landsat Thematic Mapper (TU) data set. Airborne laser scanner data obtained over the project area were used to calibrate a digital filter to extract tree heights from a digital elevation model that was derived from scanned colour stereopairs. The modelled estimates of tree height, crown size, and tree density were used to produce a decision-tree classification of forest successional stage at a landscape scale. The results obtained (72% accuracy), were limited in validation, but demonstrate potential for using the multi-scale methodology to provide spatial information for forestry policy objectives (ie., monitoring forest age structure).

G2 phase cell cycle arrest in human skin following UV irradiation

The contribution of the short wavelength ultraviolet (UV) component of sunlight to the aetiology of skin cancer has been widely acknowledged, although its direct contribution to tumour initiation or progression is still poorly understood. The loss of normal cell cycle controls, particularly checkpoint controls, are a common feature of cancer. UV radiation causes both GI and G2 phase checkpoint arrest in vitro cultured cells. In this study we have investigated the cell cycle responses to suberythemal doses of UV on skin. We have utilized short-term whole organ skin cultures, and multi parameter immunohistochemical and biochemical analysis to demonstrate that basal and suprabasal layer melanocytes and keratinocytes undergo a G2 phase cell cycle arrest for up to 48 h following irradiation. The arrest is associated with increased p16 expression but no apparent p53 involvement. This type of organ culture provides a very useful model system, combining the ease of in vitro manipulation with the ability to perform detailed molecular analysis in a normal tissue environment.

Design, modelling and analysis of a six component force balance for hypervelocity wind tunnel testing

A combination of modelling and analysis techniques was used to design a six component force balance. The balance was designed specifically for the measurement of impulsive aerodynamic forces and moments characteristic of hypervelocity shock tunnel testing using the stress wave force measurement technique. Aerodynamic modelling was used to estimate the magnitude and distribution of forces and finite element modelling to determine the mechanical response of proposed balance designs. Simulation of balance performance was based on aerodynamic loads and mechanical responses using convolution techniques. Deconvolution was then used to assess balance performance and to guide further design modifications leading to the final balance design. (C) 2001 Elsevier Science Ltd. All rights reserved.

Individual and societal bases of national identity - A comparative multi-level analysis

In a previous study (Jones and Smith, 1999) we established that much the same core pattern of national identity characterizes many developed countries. Using the national identity module from the 1995 International Social Survey Programme, we identified two dimensions of national identity: an ascriptive dimension resembling the concept of ethnic identity described in the historical and theoretical literature, and a voluntarist dimension closer to the notion of civic identity. Some writers view these dimensions in terms of a historical sequence but we find that both constructs coexist in the minds of individual respondents in the nations we examine (we exclude Bulgaria and the Philippines from the present but not the earlier analysis). The dataset used for the multilevel analyses reported here consists of 28 589 respondents in the remaining 21 countries included in the national identity database for the 1995 round of surveys. The macrosociological literature on national identity does not offer well-defined predictions about what precise patterns of national identification we might expect to find among the masses of the developed countries. There are, however, recurring themes from which one can construct plausible hypotheses about how countries might differ according to their level of development, broadly conceived. Thus, we hypothesize that forces such as post-industrialism and globalization tend to favour the more open voluntaristic form of national identity over the more restrictive ascribed form. We develop different multi-level models in order to evaluate specific hypotheses pertaining to such issues, by simultaneously relating individual and societal characteristics to the relative strength of individual commitment to these different types of national identity.

Supervisory control of wastewater treatment plants by combining principal component analysis and fuzzy c-means clustering

In this paper a methodology for integrated multivariate monitoring and control of biological wastewater treatment plants during extreme events is presented. To monitor the process, on-line dynamic principal component analysis (PCA) is performed on the process data to extract the principal components that represent the underlying mechanisms of the process. Fuzzy c-means (FCM) clustering is used to classify the operational state. Performing clustering on scores from PCA solves computational problems as well as increases robustness due to noise attenuation. The class-membership information from FCM is used to derive adequate control set points for the local control loops. The methodology is illustrated by a simulation study of a biological wastewater treatment plant, on which disturbances of various types are imposed. The results show that the methodology can be used to determine and co-ordinate control actions in order to shift the control objective and improve the effluent quality.

Multi-threading and message communication in Qu-Prolog

This paper presents the multi-threading and internet message communication capabilities of Qu-Prolog. Message addresses are symbolic and the communications package provides high-level support that completely hides details of IP addresses and port numbers as well as the underlying TCP/IP transport layer. The combination of the multi-threads and the high level inter-thread message communications provide simple, powerful support for implementing internet distributed intelligent applications.

Multi-fluorescent silica colloids for encoding large combinatorial libraries

Large chemical libraries can be synthesized on solid-support beads by the combinatorial split-and-mix method. A major challenge associated with this type of library synthesis is distinguishing between the beads and their attached compounds. A new method of encoding these solid-support beads, 'colloidal bar-coding', involves attaching fluorescent silica colloids ('reporters') to the beads as they pass through the compound synthesis, thereby creating a fluorescent bar code on each bead. In order to obtain sufficient reporter varieties to bar code extremely large libraries, many of the reporters must contain multiple fluorescent dyes. We describe here the synthesis and spectroscopic analysis of various mono- and multi-fluorescent silica particles for this purpose. It was found that by increasing the amount of a single dye introduced into the particle reaction mixture, mono- fluorescent silica particles of increasing intensities could be prepared. This increase was highly reproducible and was observed for six different fluorescent dyes. Multi-fluorescent silica particles containing up to six fluorescent dyes were also prepared. The resultant emission intensity of each dye in the multi-fluorescent particles was found to be dependent upon a number of factors; the hydrolysis rate of each silane-dye conjugate, the magnitude of the inherent emission intensity of each dye within the silica matrix, and energy transfer effects between dyes. We show that by varying the relative concentration of each silane-dye conjugate in the synthesis of multi-fluorescent particles, it is possible to change and optimize the resultant emission intensity of each dye to enable viewing in a fluorescence detection instrument.

The use of liquid phase adsorption isotherms for characterization of activated carbons

The characterization of three commercial activated carbons was carried out using the adsorption of various compounds in the aqueous phase. For this purpose the generalized adsorption isotherm was employed, and a modification of the Dubinin-Radushkevich pore filling model, incorporating repulsive contributions to the pore potential as well as bulk liquid phase nonideality, was used as the local isotherm. Eight different flavor compounds were used as adsorbates, and the isotherms were jointly fitted to yield a common pore size distribution for each carbon. The bulk liquid phase nonideality was incorporated through the UNIFAC activity coefficient model, and the repulsive contribution to the pore potential was incorporated through the Steele 10-4-3 potential model. The mean micropore network coordination number for each carbon was also determined from the fitted saturation capacity based on percolation theory. Good agreement between the model and the experimental data was observed. In addition, excellent agreement between the bimodal gamma pore size distribution and density functional theory-cum-regularization-based pore size distribution obtained by argon adsorption was also observed, supporting the validity of the model. The results show that liquid phase adsorption, using adsorptive molecules of different sizes, can be an effective means of characterizing the pore size distribution as well as connectivity. Alternately, if the carbon pore size distribution is independently known, the method can be used to measure critical molecular sizes. (C) 2001 Elsevier Science.

Xmds: eXtensible multi-dimensional simulator

A software package that efficiently solves a comprehensive range of problems based on coupled complex nonlinear stochastic ODEs and PDEs is outlined. Its input and output syntax is formulated as a subset of XML, thus making a step towards a standard for specifying numerical simulations.

Genotype-by-management interactions for grain yield and grain protein concentration of wheat

The magnitude of genotype-by-management (G x M) interactions for grain yield and grain protein concentration was examined in a multi-environment trial (MET) involving a diverse set of 272 advanced breeding lines from the Queensland wheat breeding program. The MET was structured as a series of management-regimes imposed at 3 sites for 2 years. The management-regimes were generated at each site-year as separate trials in which planting time, N fertiliser application rate, cropping history, and irrigation were manipulated. irrigation was used to simulate different rainfall regimes. From the combined analysis of variance, the G x M interaction variance components were found to be the largest source of G x E interaction variation for both grain yield (0.117 +/- 0.005 t(2) ha(-2); 49% of total G x E 0.238 +/- 0.028 t(2) ha(-2)) and grain protein concentration (0.445 +/- 0.020%(2); 82% of total G x E 0.546 +/- 0.057%(2)), and in both cases this source of variation was larger than the genotypic variance component (grain yield 0.068 +/- 0.014 t(2) ha(-2) and grain protein 0.203 +/- 0.026%(2)). The genotypic correlation between the traits varied considerably with management-regime, ranging from -0.98 to -0.31, with an estimate of 0.0 for one trial. Pattern analysis identified advanced breeding lines with improved grain yield and grain protein concentration relative to the cultivars Hartog, Sunco and Meteor. It is likely that a large component of the previously documented G x E interactions for grain yield of wheat in the northern grains region are in part a result of G x M interactions. The implications of the strong influence of G x M interactions for the conduct of wheat breeding METs in the northern region are discussed. (C) 2001 Elsevier Science B.V. All rights reserved.

Some vocabulary and grammar for the analysis of multi-environment trials, as applied to the analysis of FPB and PPB trials

For the improvement of genetic material suitable for on farm use under low-input conditions, participatory and formal plant breeding strategies are frequently presented as competing options. A common frame of reference to phrase mechanisms and purposes related to breeding strategies will facilitate clearer descriptions of similarities and differences between participatory plant breeding and formal plant breeding. In this paper an attempt is made to develop such a common framework by means of a statistically inspired language that acknowledges the importance of both on farm trials and research centre trials as sources of information for on farm genetic improvement. Key concepts are the genetic correlation between environments, and the heterogeneity of phenotypic and genetic variance over environments. Classic selection response theory is taken as the starting point for the comparison of selection trials (on farm and research centre) with respect to the expected genetic improvement in a target environment (low-input farms). The variance-covariance parameters that form the input for selection response comparisons traditionally come from a mixed model fit to multi-environment trial data. In this paper we propose a recently developed class of mixed models, namely multiplicative mixed models, also called factor-analytic models, for modelling genetic variances and covariances (correlations). Mixed multiplicative models allow genetic variances and covariances to be dependent on quantitative descriptors of the environment, and confer a high flexibility in the choice of variance-covariance structure, without requiring the estimation of a prohibitively high number of parameters. As a result detailed considerations regarding selection response comparisons are facilitated. ne statistical machinery involved is illustrated on an example data set consisting of barley trials from the International Center for Agricultural Research in the Dry Areas (ICARDA). Analysis of the example data showed that participatory plant breeding and formal plant breeding are better interpreted as providing complementary rather than competing information.