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.

Fossil trees in ancient fluvial channel deposits: evidence of seasonal and longer-term climatic variability

It has been established that large numbers of certain trees can survive in the beds of rivers of northeastern Australia where a strongly seasonal distribution of precipitation causes extreme variations in flow on both a yearly and longer-term basis. In these rivers, minimal flow occurs throughout much of any year and for periods of up to several years, allowing the trees to become established and to adapt their form in order to facilitate their survival in environments that experience periodic inundation by fast-flowing, debris-laden water. Such trees (notably paperbark trees of the angiosperm genus Melaleuca) adopt a reclined to prostrate, downstream-trailing habit, have a multiple-stemmed form, modified crown with weeping foliage, development of thick, spongy bark, anchoring of roots into firm to lithified substrates beneath the channel floor, root regeneration, and develop in flow-parallel, linear groves. Individuals from within flow-parallel, linear groves are preserved in situ within the alluvial deposit of the river following burial and death. Four examples of in situ tree fossils within alluvial channel deposits in the Permian of eastern Australia demonstrate that specialised riverbed plant communities also existed at times in the geological past. These examples, from the Lower Permian Carmila Beds, Upper Permian Moranbah Coal Measures and Baralaba Coal Measures of central Queensland and the Upper Permian Newcastle Coal Measures of central New South Wales, show several of the characteristics of trees described from modern rivers in northeastern Australia, including preservation in closely-spaced groups. These properties, together with independent sedimentological evidence, suggest that the Permian trees were adapted to an environment affected by highly variable runoff, albeit in a more temperate climatic situation than the modem Australian examples. It is proposed that occurrences of fossil trees preserved in situ within alluvial channel deposits may be diagnostic of environments controlled by seasonal and longer-term variability in fluvial runoff, and hence may have value in interpreting aspects of palaeoclimate from ancient alluvial successions. (C) 2001 Elsevier Science B.V. All rights reserved.

The anatomy of two-rooted mandibular canines determined using micro-computed tomography

P>Aim To investigate the internal and external anatomy of extracted human mandibular canines with two roots and two distinct canals using micro-computed tomography (mu CT). Methodology Fourteen two-rooted human mandibular canines were scanned using a high-resolution mu CT system (SkyScan 1174v2; SkyScan N.V., Kontich, Belgium). The images were processed to evaluate the size of the roots, the furcation regions, the presence of accessory canals, the mean distances between several anatomical landmarks, the position of the apical foramina, the direction of root curvatures, the cross-sectional appearances (SMI index), the volume and surface areas of the root canals. Results Root bifurcation was located in both apical (44%, n = 6) and middle (58%, n = 8) thirds of the root. The size of the buccal and lingual roots was similar in 29% of the sample. From a buccal view, no curvature towards the lingual or buccal direction occurred in either roots. From a proximal view, no straight lingual root occurred. In both views, S-shaped roots were found in 21% of the specimens. Location of the apical foramen varied considerably, tending to the mesio-buccal aspect of both roots. Lateral and furcation canals were observed mostly in the cervical third in 29% and 65% of the sample, respectively. The structure model index (SMI) index ranged from 1.87 to 3.86, with a mean value of 2.93 +/- 0.46. Mean volume and area of the root canals were 11.52 +/- 3.44 mm3 and 71.16 +/- 11.83 mm2, respectively. Conclusions The evaluation of two-rooted mandibular canines revealed that bifurcations occurred in the apical and middle third. S-shaped roots were found in 21% of the specimens. Mean volume, surface area and SMI index of the root canals were 11.52 mm3, 71.16 mm2 and 2.93, respectively.

On surrogate methods for detecting lateral gene transfer

Surrogate methods for detecting lateral gene transfer are those that do not require inference of phylogenetic trees. Herein I apply four such methods to identify open reading frames (ORFs) in the genome of Escherichia coli K12 that may have arisen by lateral gene transfer. Only two of these methods detect the same ORFs more frequently than expected by chance, whereas several intersections contain many fewer ORFs than expected. Each of the four methods detects a different non-random set of ORFs. The methods may detect lateral ORFs of different relative ages; testing this hypothesis will require rigorous inference of trees. (C) 2001 Federation of European Microbiological Societies. Published by Elsevier Science BN. All rights reserved.

Computation of the linear elastic properties of random porous materials with a wide variety of microstructure

A finite-element method is used to study the elastic properties of random three-dimensional porous materials with highly interconnected pores. We show that Young's modulus, E, is practically independent of Poisson's ratio of the solid phase, nu(s), over the entire solid fraction range, and Poisson's ratio, nu, becomes independent of nu(s) as the percolation threshold is approached. We represent this behaviour of nu in a flow diagram. This interesting but approximate behaviour is very similar to the exactly known behaviour in two-dimensional porous materials. In addition, the behaviour of nu versus nu(s) appears to imply that information in the dilute porosity limit can affect behaviour in the percolation threshold limit. We summarize the finite-element results in terms of simple structure-property relations, instead of tables of data, to make it easier to apply the computational results. Without using accurate numerical computations, one is limited to various effective medium theories and rigorous approximations like bounds and expansions. The accuracy of these equations is unknown for general porous media. To verify a particular theory it is important to check that it predicts both isotropic elastic moduli, i.e. prediction of Young's modulus alone is necessary but not sufficient. The subtleties of Poisson's ratio behaviour actually provide a very effective method for showing differences between the theories and demonstrating their ranges of validity. We find that for moderate- to high-porosity materials, none of the analytical theories is accurate and, at present, numerical techniques must be relied upon.

Partitioned Bremer support and multiple trees

Partitioned Bremer support (PBS) is a valuable means of assessing congruence in combined data sets, but some aspects require clarification. When more than one equally parsimonious tree is found during the constrained search for trees lacking the node of interest, averaging PBS for each data set across these trees can conceal conflict, and PBS should ideally be examined for each constrained tree. Similarly, when multiple most parsimonious trees (MPTs) are generated during analysis of the combined data, PBS is usually calculated on the consensus tree. However, extra information can be obtained if PBS is calculated on each of the MPTs or even suboptimal trees. (C) 2002 The Willi Hennig Society. Published by Elsevier Science (USA). All rights reserved.

Implementing the quantum random walk

Recently, several groups have investigated quantum analogues of random walk algorithms, both on a line and on a circle. It has been found that the quantum versions have markedly different features to the classical versions. Namely, the variance on the line, and the mixing time on the circle increase quadratically faster in the quantum versions as compared to the classical versions. Here, we propose a scheme to implement the quantum random walk on a line and on a circle in an ion trap quantum computer. With current ion trap technology, the number of steps that could be experimentally implemented will be relatively small. However, we show how the enhanced features of these walks could be observed experimentally. In the limit of strong decoherence, the quantum random walk tends to the classical random walk. By measuring the degree to which the walk remains quantum, '' this algorithm could serve as an important benchmarking protocol for ion trap quantum computers.

Inferring genome trees by using a filter to eliminate phylogenetically discordant sequences and a distance matrix based on mean normalized BLASTP scores

Darwin's paradigm holds that the diversity of present-day organisms has arisen via a process of genetic descent with modification, as on a bifurcating tree. Evidence is accumulating that genes are sometimes transferred not along lineages but rather across lineages. To the extent that this is so, Darwin's paradigm can apply only imperfectly to genomes, potentially complicating or perhaps undermining attempts to reconstruct historical relationships among genomes (i.e., a genome tree). Whether most genes in a genome have arisen via treelike (vertical) descent or by lateral transfer across lineages can be tested if enough complete genome sequences are used. We define a phylogenetically discordant sequence (PDS) as an open reading frame (ORF) that exhibits patterns of similarity relationships statistically distinguishable from those of most other ORFs in the same genome. PDSs represent between 6.0 and 16.8% (mean, 10.8%) of the analyzable ORFs in the genomes of 28 bacteria, eight archaea, and one eukaryote (Saccharomyces cerevisiae). In this study we developed and assessed a distance-based approach, based on mean pairwise sequence similarity, for generating genome trees. Exclusion of PDSs improved bootstrap support for basal nodes but altered few topological features, indicating that there is little systematic bias among PDSs. Many but not all features of the genome tree from which PDSs were excluded are consistent with the 16S rRNA tree.

Conditional simulation of random fields by successive residuals

This paper presents a new approach to the LU decomposition method for the simulation of stationary and ergodic random fields. The approach overcomes the size limitations of LU and is suitable for any size simulation. The proposed approach can facilitate fast updating of generated realizations with new data, when appropriate, without repeating the full simulation process. Based on a novel column partitioning of the L matrix, expressed in terms of successive conditional covariance matrices, the approach presented here demonstrates that LU simulation is equivalent to the successive solution of kriging residual estimates plus random terms. Consequently, it can be used for the LU decomposition of matrices of any size. The simulation approach is termed conditional simulation by successive residuals as at each step, a small set (group) of random variables is simulated with a LU decomposition of a matrix of updated conditional covariance of residuals. The simulated group is then used to estimate residuals without the need to solve large systems of equations.

Nitrogen use strategies of neotropical rainforest trees in threatened Atlantic Forest

The characteristics of nitrogen acquisition, transport and assimilation were investigated in species of an Atlantic Forest succession over calcareous soil in south-eastern Brazil. Differences in behaviour were observed within the regeneration guilds. Pioneer species showed high leaf nitrogen contents, a high capacity to respond to increased soil nitrogen availability, a high capacity for leaf nitrate assimilation and were characterized by the transport of nitrate + asparagine. At the other end of the succession, late secondary species had low leaf nitrogen contents, little capacity to respond to increased soil nitrogen availability, low leaf nitrate assimilation and were active in the transport of asparagine + arginine. The characteristics of nitrogen nutrition in some early secondary species showed similarities to those of pioneer species whereas others more closely resembled late secondary species. Average leaf delta(15)N values increased along the successional gradient. The results indicate that the nitrogen metabolism characteristics of species may be an additional ecophysiological tool in classifying tropical forest tree species into ecological guilds, and may have implications for regeneration programmes in degraded areas.

Mean anisotropy of homogeneous Gaussian random fields and anisotropic norms of linear translation-invariant operators on multidimensional integer lattices

Sensitivity of output of a linear operator to its input can be quantified in various ways. In Control Theory, the input is usually interpreted as disturbance and the output is to be minimized in some sense. In stochastic worst-case design settings, the disturbance is considered random with imprecisely known probability distribution. The prior set of probability measures can be chosen so as to quantify how far the disturbance deviates from the white-noise hypothesis of Linear Quadratic Gaussian control. Such deviation can be measured by the minimal Kullback-Leibler informational divergence from the Gaussian distributions with zero mean and scalar covariance matrices. The resulting anisotropy functional is defined for finite power random vectors. Originally, anisotropy was introduced for directionally generic random vectors as the relative entropy of the normalized vector with respect to the uniform distribution on the unit sphere. The associated a-anisotropic norm of a matrix is then its maximum root mean square or average energy gain with respect to finite power or directionally generic inputs whose anisotropy is bounded above by a≥0. We give a systematic comparison of the anisotropy functionals and the associated norms. These are considered for unboundedly growing fragments of homogeneous Gaussian random fields on multidimensional integer lattice to yield mean anisotropy. Correspondingly, the anisotropic norms of finite matrices are extended to bounded linear translation invariant operators over such fields.

Carboxy-fluorescein diacetate, succinimidyl ester labeled papillomavirus virus-like particles fluoresce after internalization and interact with heparan sulfate for binding and entry

Human papillomaviruses (HPVs) infect epithelial cells and are associated with genital carcinoma. Most epithelial cell lines express cell-surface glycosaminoglycans (GAGs) usually found attached to the protein core of proteoglycans. Our aim was to study how GAGs influenced HPV entry. Using a human keratinocyte cell line (HaCaT), preincubation of HPV virus-like particles (VLPs) with GAGs showed a dose-dependent inhibition of binding. The IC50 (50% inhibition) was only 0.5 mug/ml for heparin, 1 mug/ml for dextran sulfate, and 5-10 mug/ml for heparan sulfate from mucosal origin. Mutated chinese hamster ovary (CHO) cell lines lacking heparan sulfate or all GAGs were unable to bind HPV VLPs. Here we also report a method to study internalization by using VLPs labeled with carboxy-fluorescein diacetate, succinimidyl ester, a fluorochrome that is only activated after cell entry. Pretreatment of labeled HPV VLPs with heparin inhibited uptake, suggesting a primary interaction between HPV and cell-surface heparan sulfate. (C) 2003 Elsevier Science (USA). All rights reserved.