Mixture modeling for joint analysis of survival, discrete, and continuous data

Mixture modeling is commonly used to model categorical latent variables that represent subpopulations in which population membership is unknown but can be inferred from the data. In relatively recent years, the potential of finite mixture models has been applied in time-to-event data. However, the commonly used survival mixture model assumes that the effects of the covariates involved in failure times differ across latent classes, but the covariate distribution is homogeneous. The aim of this dissertation is to develop a method to examine time-to-event data in the presence of unobserved heterogeneity under a framework of mixture modeling. A joint model is developed to incorporate the latent survival trajectory along with the observed information for the joint analysis of a time-to-event variable, its discrete and continuous covariates, and a latent class variable. It is assumed that the effects of covariates on survival times and the distribution of covariates vary across different latent classes. The unobservable survival trajectories are identified through estimating the probability that a subject belongs to a particular class based on observed information. We applied this method to a Hodgkin lymphoma study with long-term follow-up and observed four distinct latent classes in terms of long-term survival and distributions of prognostic factors. Our results from simulation studies and from the Hodgkin lymphoma study demonstrated the superiority of our joint model compared with the conventional survival model. This flexible inference method provides more accurate estimation and accommodates unobservable heterogeneity among individuals while taking involved interactions between covariates into consideration.^

Comparisons between analogue and numerical models of thrust wedge development

Analogue and finite element numerical models with frictional and viscous properties are used to model thrust wedge development. Comparison between model types yields valuable information about analogue model evolution, scaling laws and the relative strengths and limitations of the techniques. Both model types show a marked contrast in structural style between ‘frictional-viscous domains’ underlain by a thin viscous layer and purely ‘frictional domains’. Closely spaced thrusts form a narrow and highly asymmetric fold-and-thrust belt in the frictional domain, characterized by in-sequence propagation of forward thrusts. In contrast, the frictional-viscous domain shows a wide and low taper wedge and a thrust belt with a more symmetrical vergence, with both forward and back thrusts. The frictional-viscous domain numerical models show that the viscous layer initially simple shears as deformation propagates along it, while localized deformation resulting in the formation of a pop-up structure occurs in the overlying frictional layers. In both domains, thrust shear zones in the numerical model are generally steeper than the equivalent faults in the analogue model, because the finite element code uses a non-associated plasticity flow law. Nevertheless, the qualitative agreement between analogue and numerical models is encouraging. It shows that the continuum approximation used in numerical models can be used to model frictional materials, such as sand, provided caution is taken to properly scale the experiments, and some of the limitations are taken into account.

Modelling pre-clearing vegetation distribution using GIS-integrated statistical, ecological and data models: A case study from the wet tropics of Northeastern Australia

Traditional vegetation mapping methods use high cost, labour-intensive aerial photography interpretation. This approach can be subjective and is limited by factors such as the extent of remnant vegetation, and the differing scale and quality of aerial photography over time. An alternative approach is proposed which integrates a data model, a statistical model and an ecological model using sophisticated Geographic Information Systems (GIS) techniques and rule-based systems to support fine-scale vegetation community modelling. This approach is based on a more realistic representation of vegetation patterns with transitional gradients from one vegetation community to another. Arbitrary, though often unrealistic, sharp boundaries can be imposed on the model by the application of statistical methods. This GIS-integrated multivariate approach is applied to the problem of vegetation mapping in the complex vegetation communities of the Innisfail Lowlands in the Wet Tropics bioregion of Northeastern Australia. The paper presents the full cycle of this vegetation modelling approach including sampling sites, variable selection, model selection, model implementation, internal model assessment, model prediction assessments, models integration of discrete vegetation community models to generate a composite pre-clearing vegetation map, independent data set model validation and model prediction's scale assessments. An accurate pre-clearing vegetation map of the Innisfail Lowlands was generated (0.83r(2)) through GIS integration of 28 separate statistical models. This modelling approach has good potential for wider application, including provision of. vital information for conservation planning and management; a scientific basis for rehabilitation of disturbed and cleared areas; a viable method for the production of adequate vegetation maps for conservation and forestry planning of poorly-studied areas. (c) 2006 Elsevier B.V. All rights reserved.

Bayesian online algorithms for learning in discrete Hidden Markov Models

We propose and analyze two different Bayesian online algorithms for learning in discrete Hidden Markov Models and compare their performance with the already known Baldi-Chauvin Algorithm. Using the Kullback-Leibler divergence as a measure of generalization we draw learning curves in simplified situations for these algorithms and compare their performances.

Online learning in discrete hidden Markov models

We present and analyze three different online algorithms for learning in discrete Hidden Markov Models (HMMs) and compare their performance with the Baldi-Chauvin Algorithm. Using the Kullback-Leibler divergence as a measure of the generalization error we draw learning curves in simplified situations and compare the results. The performance for learning drifting concepts of one of the presented algorithms is analyzed and compared with the Baldi-Chauvin algorithm in the same situations. A brief discussion about learning and symmetry breaking based on our results is also presented. © 2006 American Institute of Physics.

Notes on the Solution of Dynamic Optimization Problems with Explicit Controls in Discrete-Time Economic Models

Йордан Йорданов, Андрей Василев - В работата се изследват методи за решаването на задачи на оптималното управление в дискретно време с безкраен хоризонт и явни управления. Дадена е обосновка на една процедура за решаване на такива задачи, базирана на множители на Лагранж, коята често се употребява в икономическата литература. Извеждени са необходимите условия за оптималност на базата на уравнения на Белман и са приведени достатъчни условия за оптималност при допускания, които често се използват в икономиката.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

Dermaseptin 01 as antimicrobial peptide with rich biotechnological potential: study of peptide interaction with membranes containing Leishmania amazonensis lipid-rich extract and membrane models

This article addresses the interactions of the synthetic antimicrobial peptide dermaseptin 01 (GLWSTIKQKGKEAAIAAA-KAAGQAALGAL-NH(2), DS 01) with phospholipid (PL) monolayers comprising (i) a lipid-rich extract of Leishmania amazonensis (LRE-La), (ii) zwitterionic PL (dipalmitoylphosphatidylcholine, DPPC), and (iii) negatively charged PL (dipalmitoylphosphatidylglycerol, DPPG). The degree of interaction of DS 01 with the different biomembrane models was quantified from equilibrium and dynamic liquid-air interface parameters. At low peptide concentrations, interactions between DS 01 and zwitterionic PL, as well as with the LRE-La monolayers were very weak, whereas with negatively charged PLs the interactions were stronger. For peptide concentrations above 1 mu g/ml, a considerable expansion of negatively charged monolayers occurred. In the case of DPPC, it was possible to return to the original lipid area in the condensed phase, suggesting that the peptide was expelled from the monolayer. However, in the case of DPPG, the average area per lipid molecule in the presence of DS 01 was higher than pure PLs even at high surface pressures, suggesting that at least part of DS 01 remained incorporated in the monolayer. For the LRE-La monolayers, DS 01 also remained in the monolayer. This is the first report on the antiparasitic activity of AMPs using Langmuir monolayers of a natural lipid extract from L. amazonensis. Copyright (C) 2011 European Peptide Society and John Wiley & Sons, Ltd.

Virtual plants - integrating architectural and physiological models

Recent advances in computer technology have made it possible to create virtual plants by simulating the details of structural development of individual plants. Software has been developed that processes plant models expressed in a special purpose mini-language based on the Lindenmayer system formalism. These models can be extended from their architectural basis to capture plant physiology by integrating them with crop models, which estimate biomass production as a consequence of environmental inputs. Through this process, virtual plants will gain the ability to react to broad environmental conditions, while crop models will gain a visualisation component. This integration requires the resolution of the fundamentally different time scales underlying the approaches. Architectural models are usually based on physiological time; each time step encompasses the same amount of development in the plant, without regard to the passage of real time. In contrast, physiological models are based in real time; the amount of development in a time step is dependent on environmental conditions during the period. This paper provides a background on the plant modelling language, then describes how widely-used concepts of thermal time can be implemented to resolve these time scale differences. The process is illustrated using a case study. (C) 1997 Elsevier Science Ltd.

Systemic lupus erythematosus exhibits a dynamic and continuum spectrum of effector/regulatory T cells

Objective: The identification of regulatory T cells (Treg cells) as CD4(+)CD25(high) cells may be upset by the increased frequency of activated effector T cells (Teff cells) in inflammatory diseases such as systemic lupus erythematosus (SLE). This study aimed to evaluate the frequency of T-cell subsets according to the expression of CD25 and CD127 in active (A-SLE) and inactive SLE (I-SLE). Methods: Peripheral blood mononuclear cells (PBMCs) from 26 A-SLE patients (SLE Disease Activity Index (SLEDAI) = 10.17 +/- 3.7), 31 I-SLE patients (SLEDAI = 0), and 26 healthy controls (HC) were analysed by multicolour flow. cytometry. Results: CD25(high) cell frequency was increased in A-SLE (5.2 +/- 5.7%) compared to I-SLE (3.4 +/- 3.4%) and HC (1.73 +/- 0.8%) (p < 0.01). However, the percentage of FoxP3(+) cells in the CD25(high) subset was decreased in A-SLE (24.6 +/- 16.4%) compared to I-SLE (33.7 +/- 16) and HC (45 +/- 25.1%) (p < 0.01). This was partly due to the increased frequency of Teff cells (CD25(high)CD127(+)FoxP3(empty set)) in A-SLE (10.7 +/- 7.3%) compared to I-SLE (8.5 +/- 6.5) and HC (6.1 +/- 1.8%) (p = 0.02). Hence the frequency of Treg cells (CD25(+/high)CD127(low/empty set)FoxP3(+)) was equivalent in A-SLE (1.4 +/- 0.8%), I-SLE (1.37 +/- 1.0%), and HC (1.13 +/- 0.59%) (p = 0.42). A-SLE presented an increased frequency of CD25(+)CD127(+)FoxP3(+) and CD25(empty set)FoxP3(+)CD127(low/empty set) T cells, which may represent intermediate phenotypes between Treg and Teff cells. Conclusions: The present study has provided data supporting normal Treg cell frequency in A-SLE and I-SLE as well as increased frequency of Teff cells in A-SLE. This scenario reflects a Treg/Teff ratio imbalance that may favour the inflammatory phenotype of the disease. In addition, the increased frequency of T cells with putative intermediate phenotypes may be compatible with a highly dynamic immune system in SLE.

Immune responses to ISCOM (R) formulations in animal and primate models

ISCOMs(R) are typically 40 nm cage-like structures comprising antigen, saponin, cholesterol and phospholipid. ISCOMs(R) have been shown to induce antibody responses and activate T helper cells and cyrolytic T lymphocytes in a number of animal species, including non-human primates. Recent clinical studies have demonstrated that ISCOMs(R) are also able to induce antibody and cellular immune responses in humans. This review describes the current understanding of the ability of ISCOMs(R) to induce immune responses and the mechanisms underlying this property. Recent progress in the characterisation and manufacture of ISCOMs(R) will also be discussed. (C) 2001 Elsevier Science Ltd. All rights reserved.

Linking physiological and architectural models of cotton

Despite the strong influence of plant architecture on crop yield, most crop models either ignore it or deal with it in a very rudimentary way. This paper demonstrates the feasibility of linking a model that simulates the morphogenesis and resultant architecture of individual cotton plants with a crop model that simulates the effects of environmental factors on critical physiological processes and resulting yield in cotton. First the varietal parameters of the models were made concordant. Then routines were developed to allocate the flower buds produced each day by the crop model amongst the potential positions generated by the architectural model. This allocation is done according to a set of heuristic rules. The final weight of individual bolls and the shedding of buds and fruit caused by water, N, and C stresses are processed in a similar manner. Observations of the positions of harvestable fruits, both within and between plants, made under a variety of agronomic conditions that had resulted in a broad range of plant architectures were compared to those predicted by the model with the same environmental inputs. As illustrated by comparisons of plant maps, the linked models performed reasonably well, though performance of the fruiting point allocation and shedding algorithms could probably be improved by further analysis of the spatial relationships of retained fruit. (C) 2002 Elsevier Science Ltd. All rights reserved.