Stochastic exponential stability of singular linear systems with Markov jump parameters

This paper deals with exponential stability of discrete-time singular systems with Markov jump parameters. We propose a set of coupled generalized Lyapunov equations (CGLE) that provides sufficient conditions to check this property for this class of systems. A method for solving the obtained CGLE is also presented, based on iterations of standard singular Lyapunov equations. We present also a numerical example to illustrate the effectiveness of the approach we are proposing.

Depositional history of sedimentary linear alkylbenzenes (LABs) in a large South American industrial coastal area (Santos Estuary, Southeastern Brazil)

This paper reports the reconstruction of the contamination history of a large South American industrial coastal area (Santos Estuary, Brazil) using linear alkylbenzenes (LABs). Three sediment cores were dated by (137)Cs Concentrations in surficial layers were comparable to the midrange concentrations reported for coastal sediments worldwide LAB concentrations increased towards the surface. indicating increased waste discharges into the estuary in recent decades. The highest concentration values occurred in the early 1970s, a time of intense industrial activity and marked population growth. The decreased LAB concentration, in the late 1970s was assumed to be the result of the world oil crisis Treatment of industrial effluents, which began in 1984, was represented by decreased LAB levels Microbial degradation of LABs may be more intense in the industrial area sediments. The results show that industrial and domestic waste discharges are a historical problem in the area. (C) 2010 Elsevier Ltd. All rights reserved.

The collapse of linear polyelectrolyte chains in a poor solvent: When does a collapsing polyelectrolyte collect its counter ions?

The collapse of linear polyelectrolyte chains in a poor solvent: When does a collapsing polyelectrolyte collect its counter ions? The collapse of polyions in a poor solvent is a complex system and is an active research subject in the theoretical polyelectrolyte community. The complexity is due to the subtle interplay between hydrophobic effects, electrostatic interactions, entropy elasticity, intrinsic excluded volume as well as specific counter-ion and co-ion properties. Long range Coulomb forces can obscure single molecule properties. The here presented approach is to use just a small amount of screening salt in combination with a very high sample dilution in order to screen intermolecular interaction whereas keeping intramolecular interaction as much as possible (polyelectrolyte concentration cp ≤ 12 mg/L, salt concentration; Cs = 10^-5 mol/L). This is so far not described in literature. During collapse, the polyion is subject to a drastic change in size along with strong reduction of free counterions in solution. Therefore light scattering was utilized to obtain the size of the polyion whereas a conductivity setup was developed to monitor the proceeding of counterion collection by the polyion. Partially quaternized PVP’s below and above the Manning limit were investigated and compared to the collapse of their uncharged precursor. The collapses were induced by an isorefractive solvent/non-solvent mixture consisting of 1-propanol and 2-pentanone, with nearly constant dielectric constant. The solvent quality for the uncharged polyion could be quantified which, for the first time, allowed the experimental investigation of the effect of electrostatic interaction prior and during polyion collapse. Given that the Manning parameter M for QPVP4.3 is as low as lB / c = 0.6 (lB the Bjerrum length and c the mean contour distance between two charges), no counterion binding should occur. However the Walden product reduces with first addition of non solvent and accelerates when the structural collapse sets in. Since the dielectric constant of the solvent remains virtually constant during the chain collapse, the counterion binding is entirely caused by the reduction in the polyion chain dimension. The collapse is shifted to lower wns with higher degrees of quaternization as the samples QPVP20 and QPVP35 show (M = 2.8 respectively 4.9). The combination of light scattering and conductivity measurement revealed for the first time that polyion chains already collect their counter ions well above the theta-dimension when the dimensions start to shrink. Due to only small amounts of screening salt, strong electrostatic interactions bias dynamic as well as static light scattering measurements. An extended Zimm formula was derived to account for this interaction and to obtain the real chain dimensions. The effective degree of dissociation g could be obtained semi quantitatively using this extrapolated static in combination with conductivity measurements. One can conclude the expansion factor a and the effective degree of ionization of the polyion to be mutually dependent. In the good solvent regime g of QPVP4.3, QPVP20 and QPVP35 exhibited a decreasing value in the order 1 > g4.3 > g20 > g35. The low values of g for QPVP20 and QPVP35 are assumed to be responsible for the prior collapse of the higher quaternized samples. Collapse theory predicts dipole-dipole attraction to increase accordingly and even predicts a collapse in the good solvent regime. This could be exactly observed for the QPVP35 sample. The experimental results were compared to a newly developed theory of uniform spherical collapse induced by concomitant counterion binding developed by M. Muthukumar and A. Kundagrami. The theory agrees qualitatively with the location of the phase boundary as well as the trend of an increasing expansion with an increase of the degree of quaternization. However experimental determined g for the samples QPVP4.3, QPVP20 and QPVP35 decreases linearly with the degree of quaternization whereas this theory predicts an almost constant value.

Thin films of polythiophene: Linear and nonlinear optical characterization

Wir haben die linearen und nichtlinearen optischen Eigenschaften von dünnen Schichten und planaren Wellenleitern aus mehreren konjugierten Polymeren (MEH-PPV und P3AT) und Polymeren mit -Elektronen Systemen in der Seitenkette (PVK und PS) untersucht und verglichen. PVK und PS haben relativ kleine Werte des nichtlinearen Brechungsindex n2 bei 532 nm, nämlich (1,2 ± 0,5)10-14 cm2/W und (2,6 ± 0,5) 10-14 cm2/W.rnWir haben die linearen optischen Konstanten von mehreren P3ATs untersucht, insbesondere den Einfluss der Regioregularität und Kettenlänge der Alkylsubstituenten. Wir haben das am besten geeignete Polymere für Wellenleiter Anwendungen identifiziert, welches P3BT-ra genannt ist. Wir haben die linearen optischen Eigenschaften dünner Schichten des P3BT-ra untersucht, die mit Spincoating aus verschiedenen Lösungsmitteln mit unterschiedlichen Siedetemperaturen präparieret wurden. Wir haben festgestellt, dass P3BT-ra Filme aus Toluol-Lösungen die am besten geeigneten Wellenleiter für die intensitätsabhängigen Prismen-Kopplungs Experimente sind, weil diese geringe Wellenleiterdämpfungsverluste bei = 1064 nm haben. rnWir haben die Dispersionen des Wellenleiterdämfungsverlustes gw, des nichtlinearen Brechungsindex n2 und des nichtlinearen Absorptionskoeffizienten 2 von Wellenleitern aus P3BT-ra im Bereich von 700 - 1500 nm gemessen. Wir haben große Werte des nichtlinearen Brechungsindex bis 1,5x10-13 cm2/W bei 1150 nm beobachtet. Wir haben gefunden, dass die Gütenkriterien (“figures of merit“) für rein optische Schalter im Wellenlängebereich 1050 - 1200 nm erfüllt sind. Dieser Bereich entspricht dem niederenergetischen Ausläufer der Zwei-Photonen-Absorption. Die Gütekriterien von P3BT-ra gehören zu den besten der bisher bekannten Werte von konjugierten Polymeren.rnWir haben gefunden, dass P3BT-ra ein vielversprechender Kandidat für integriert-optische Schalter ist, weil es eine gute Kombination aus großer Nichtlinearität dritter Ordnung, geringen Wellenleiterdämpfungverlusten und ausreichender Photostabilität zeigt. rnWir haben einen Vergleich der gemessenen Dispersion von gw, n2 und 2 mit der Theorie durchgeführt. Durch Kurvenanpassung der Dispersion von gw haben wir gefunden, dass Rayleigh-Streuung der dominierende Dämpfungsmechanismus in MEH-PPV und P3BT-ra Wellenleitern ist. Ein quantenmechanischer Ansatz wurde zur Berechnung der nichtlinearen Suszeptibilität dritter Ordnung (3) verwendet, um die gemessenen Spektren von n2 und 2 von P3BT-ra und MEH-PPV zu simulieren. Dies kann erklären, dass sättigbare Absorption und Zwei-Photonen Absorption die hauptsächlichen Effekte sind, welche die Dispersion von n2 und 2 verursachen. rn

COVARIATE-ADJUSTED NONPARAMETRIC ANALYSIS OF MAGNETIC RESONANCE IMAGES USING MARKOV CHAIN MONTE CARLO

Permutation tests are useful for drawing inferences from imaging data because of their flexibility and ability to capture features of the brain that are difficult to capture parametrically. However, most implementations of permutation tests ignore important confounding covariates. To employ covariate control in a nonparametric setting we have developed a Markov chain Monte Carlo (MCMC) algorithm for conditional permutation testing using propensity scores. We present the first use of this methodology for imaging data. Our MCMC algorithm is an extension of algorithms developed to approximate exact conditional probabilities in contingency tables, logit, and log-linear models. An application of our non-parametric method to remove potential bias due to the observed covariates is presented.

Analysis of spring break-up and its effects on a biomass feedstock supply chain in northern Michigan

Demand for bio-fuels is expected to increase, due to rising prices of fossil fuels and concerns over greenhouse gas emissions and energy security. The overall cost of biomass energy generation is primarily related to biomass harvesting activity, transportation, and storage. With a commercial-scale cellulosic ethanol processing facility in Kinross Township of Chippewa County, Michigan about to be built, models including a simulation model and an optimization model have been developed to provide decision support for the facility. Both models track cost, emissions and energy consumption. While the optimization model provides guidance for a long-term strategic plan, the simulation model aims to present detailed output for specified operational scenarios over an annual period. Most importantly, the simulation model considers the uncertainty of spring break-up timing, i.e., seasonal road restrictions. Spring break-up timing is important because it will impact the feasibility of harvesting activity and the time duration of transportation restrictions, which significantly changes the availability of feedstock for the processing facility. This thesis focuses on the statistical model of spring break-up used in the simulation model. Spring break-up timing depends on various factors, including temperature, road conditions and soil type, as well as individual decision making processes at the county level. The spring break-up model, based on the historical spring break-up data from 27 counties over the period of 2002-2010, starts by specifying the probability distribution of a particular county’s spring break-up start day and end day, and then relates the spring break-up timing of the other counties in the harvesting zone to the first county. In order to estimate the dependence relationship between counties, regression analyses, including standard linear regression and reduced major axis regression, are conducted. Using realizations (scenarios) of spring break-up generated by the statistical spring breakup model, the simulation model is able to probabilistically evaluate different harvesting and transportation plans to help the bio-fuel facility select the most effective strategy. For early spring break-up, which usually indicates a longer than average break-up period, more log storage is required, total cost increases, and the probability of plant closure increases. The risk of plant closure may be partially offset through increased use of rail transportation, which is not subject to spring break-up restrictions. However, rail availability and rail yard storage may then become limiting factors in the supply chain. Rail use will impact total cost, energy consumption, system-wide CO2 emissions, and the reliability of providing feedstock to the bio-fuel processing facility.

A comparison of Markov and generalized linear models of hospital mortality

This paper reports a comparison of three modeling strategies for the analysis of hospital mortality in a sample of general medicine inpatients in a Department of Veterans Affairs medical center. Logistic regression, a Markov chain model, and longitudinal logistic regression were evaluated on predictive performance as measured by the c-index and on accuracy of expected numbers of deaths compared to observed. The logistic regression used patient information collected at admission; the Markov model was comprised of two absorbing states for discharge and death and three transient states reflecting increasing severity of illness as measured by laboratory data collected during the hospital stay; longitudinal regression employed Generalized Estimating Equations (GEE) to model covariance structure for the repeated binary outcome. Results showed that the logistic regression predicted hospital mortality as well as the alternative methods but was limited in scope of application. The Markov chain provides insights into how day to day changes of illness severity lead to discharge or death. The longitudinal logistic regression showed that increasing illness trajectory is associated with hospital mortality. The conclusion is reached that for standard applications in modeling hospital mortality, logistic regression is adequate, but for new challenges facing health services research today, alternative methods are equally predictive, practical, and can provide new insights. ^

A re-examination of petrogenesis and 40Ar/39Ar systematics in the the Chain of Ponds K-feldspar: "diffusion domain" archetype versus polyphase hygrochronology

K-feldspar (Kfs) from the Chain of Ponds Pluton (CPP) is the archetypal reference material, on which thermochronological modeling of Ar diffusion in discrete “domains” was founded. We re-examine the CPP Kfs using cathodoluminescence and back-scattered electron imaging, transmission electron microscopy, and electron probe microanalysis. 40Ar/39Ar stepwise heating experiments on different sieve fractions, and on handpicked and unpicked aliquots, are compared. Our results reproduce the staircase-shaped age spectrum and the Arrhenius trajectory of the literature sample, confirming that samples collected from the same locality have an identical Ar isotope record. Even the most pristine-looking Kfs from the CPP contains successive generations of secondary, metasomatic/retrograde mineral replacements that post-date magmatic crystallization. These chemically and chronologically distinct phases are responsible for its staircase-shaped age spectra, which are modified by handpicking. While genuine within-grain diffusion gradients are not ruled out by these data, this study demonstrates that the most important control on staircase-shaped age spectra is the simultaneous presence of heterochemical, diachronous post-magmatic mineral growth. At least five distinct mineral species were identified in the Kfs separate, three of which can be traced to external fluids interacting with the CPP in a chemically open system. Sieve fractions have size-shifted Arrhenius trajectories, negating the existence of the smallest “diffusion domains”. Heterochemical phases also play an important role in producing non-linear trajectories. In vacuo degassing rates recovered from Arrhenius plots are neither related to true Fick’s Law diffusion nor to the staircase shape of the age spectra. The CPP Kfs used to define the "diffusion domain" model demonstrates the predominance of metasomatic alteration by hydrothermal fluids and recrystallization in establishing the natural Ar distribution amongst different coexisting phases that gives rise to the staircase-shaped age spectrum. Microbeam imaging of textures is as essential for 40Ar-39Ar hygrochronology as it is for U-Pb geochronology.

Bayesian generalized linear models for meta-analysis of diagnostic tests

With the recognition of the importance of evidence-based medicine, there is an emerging need for methods to systematically synthesize available data. Specifically, methods to provide accurate estimates of test characteristics for diagnostic tests are needed to help physicians make better clinical decisions. To provide more flexible approaches for meta-analysis of diagnostic tests, we developed three Bayesian generalized linear models. Two of these models, a bivariate normal and a binomial model, analyzed pairs of sensitivity and specificity values while incorporating the correlation between these two outcome variables. Noninformative independent uniform priors were used for the variance of sensitivity, specificity and correlation. We also applied an inverse Wishart prior to check the sensitivity of the results. The third model was a multinomial model where the test results were modeled as multinomial random variables. All three models can include specific imaging techniques as covariates in order to compare performance. Vague normal priors were assigned to the coefficients of the covariates. The computations were carried out using the 'Bayesian inference using Gibbs sampling' implementation of Markov chain Monte Carlo techniques. We investigated the properties of the three proposed models through extensive simulation studies. We also applied these models to a previously published meta-analysis dataset on cervical cancer as well as to an unpublished melanoma dataset. In general, our findings show that the point estimates of sensitivity and specificity were consistent among Bayesian and frequentist bivariate normal and binomial models. However, in the simulation studies, the estimates of the correlation coefficient from Bayesian bivariate models are not as good as those obtained from frequentist estimation regardless of which prior distribution was used for the covariance matrix. The Bayesian multinomial model consistently underestimated the sensitivity and specificity regardless of the sample size and correlation coefficient. In conclusion, the Bayesian bivariate binomial model provides the most flexible framework for future applications because of its following strengths: (1) it facilitates direct comparison between different tests; (2) it captures the variability in both sensitivity and specificity simultaneously as well as the intercorrelation between the two; and (3) it can be directly applied to sparse data without ad hoc correction. ^

New methods for quantification and analysis of quantitative real-time polymerase chain reaction data

Quantitative real-time polymerase chain reaction (qPCR) is a sensitive gene quantitation method that has been widely used in the biological and biomedical fields. The currently used methods for PCR data analysis, including the threshold cycle (CT) method, linear and non-linear model fitting methods, all require subtracting background fluorescence. However, the removal of background fluorescence is usually inaccurate, and therefore can distort results. Here, we propose a new method, the taking-difference linear regression method, to overcome this limitation. Briefly, for each two consecutive PCR cycles, we subtracted the fluorescence in the former cycle from that in the later cycle, transforming the n cycle raw data into n-1 cycle data. Then linear regression was applied to the natural logarithm of the transformed data. Finally, amplification efficiencies and the initial DNA molecular numbers were calculated for each PCR run. To evaluate this new method, we compared it in terms of accuracy and precision with the original linear regression method with three background corrections, being the mean of cycles 1-3, the mean of cycles 3-7, and the minimum. Three criteria, including threshold identification, max R2, and max slope, were employed to search for target data points. Considering that PCR data are time series data, we also applied linear mixed models. Collectively, when the threshold identification criterion was applied and when the linear mixed model was adopted, the taking-difference linear regression method was superior as it gave an accurate estimation of initial DNA amount and a reasonable estimation of PCR amplification efficiencies. When the criteria of max R2 and max slope were used, the original linear regression method gave an accurate estimation of initial DNA amount. Overall, the taking-difference linear regression method avoids the error in subtracting an unknown background and thus it is theoretically more accurate and reliable. This method is easy to perform and the taking-difference strategy can be extended to all current methods for qPCR data analysis.^

Long chain 1,13- and 1,15-diols in surface sediments ans SST reconstruction for sediment core GeoB6518-1

Although commonly reported in marine and freshwater environments, little is known about the biological sources of long chain alkyl 1,13- and 1,15-diols, and factors controlling their distributions. Here we analyzed the occurrence and distribution of these lipids in a comprehensive set of marine surface sediments and compare their distributions with environmental conditions like sea surface temperature (SST), salinity and nutrient concentrations. Fractional abundances of the C28 1,13-, C30 1,13- and C30 1,15-diols show a strong correlation with SST and based on these results, we propose the Long chain Diol Index (LDI), which expresses the C30 1,15-diol abundance relative to those of C28 1,13-, C30 1,13- and C30 1,15-diols. The LDI shows a strong linear correlation with SST (LDI = 0.033 × SST + 0.095; R2 = 0.969, n = 162) over a temperature range of -3 to 27 °C. Long chain diol distributions in sediments from the South Atlantic close to the Congo River outflow (West Africa) provided a 43 kyr LDI SST record. This record reflects several known climatic events and shows similarities with an alkenone-derived SST record obtained using the same suite of sediments, both in trend and in terms of absolute SST. This confirms the potential of the LDI as a proxy for palaeo-SST reconstruction.

Sources and proxy potential of long chain alkyl diols in lacustrine environments

Long chain 1,13- and 1,15-alkyl diols form the base of a number of recently proposed proxies used for climate reconstruction. However, the sources of these lipids and environmental controls on their distribution are still poorly constrained. We have analyzed the long chain alkyl diol (LCD) composition of cultures of ten eustigmatophyte species, with three species from different families grown at various temperatures, to identify the effect of species composition and growth temperature on the LCD distribution. The results were compared with the LCD distribution of sixty-two lake surface sediments, and with previously reported LCD distributions from marine environments. The different families within the Eustigmatophyceae show distinct LCD patterns, with the freshwater family Eustigmataceae most closely resembling LCD distributions in both marine and lake environments. Unlike the other two eustigmatophyte families analyzed (Monodopsidaceae and Goniochloridaceae), C28 and C30 1,13-alkyl diols and C30 and C32 1,15-alkyl diols are all relatively abundant in the family Eustigmataceae, while the mono-unsaturated C32 1,15-alkyl diol was below detection limit. In contrast to the marine environment, LCD distributions in lakes did not show a clear relationship with temperature. The Long chain Diol Index (LDI), a proxy previously proposed for sea surface temperature reconstruction, showed a relatively weak correlation (R2 = 0.33) with mean annual air temperature used as an approximation for annual mean surface temperature of the lakes. A much-improved correlation (R2 = 0.74, p-value<0.001) was observed applying a multiple linear regression analysis between LCD distributions and lake temperatures reconstructed using branched tetraether lipid distributions. The obtained regression model provides good estimates of temperatures for cultures of the family Eustigmataceae, suggesting that algae belonging to this family have an important role as a source for LCDs in lacustrine environments, or, alternatively, that the main sources of LCDs are similarly affected by temperature as the Eustigmataceae. The results suggest that LCDs may have the potential to be applicable as a palaeotemperature proxy for lacustrine environments, although further calibration work is still required.

Three functional luciferase domains in a single polypeptide chain

We report a unique case of a gene containing three homologous and contiguous repeat sequences, each of which, after excision, cloning, and expression in Escherichia coli, is shown to code for a peptide catalyzing the same reaction as the native protein, Gonyaulax polyedra luciferase (Mr = 137). This enzyme, which catalyzes the light-emitting oxidation of a linear tetrapyrrole (dinoflagellate luciferin), exhibits no sequence similarities to other luciferases in databases. Sequence analysis also reveals an unusual evolutionary feature of this gene: synonymous substitutions are strongly constrained in the central regions of each of the repeated coding sequences.

The α chain of laminin-1 is independently secreted and drives secretion of its β- and γ-chain partners

A mammalian recombinant strategy was established to dissect rules of basement membrane laminin assembly and secretion. The α-, β-, and γ-chain subunits of laminin-1 were expressed in all combinations, transiently and/or stably, in a near-null background. In the absence of its normal partners, the α chain was secreted as intact protein and protein that had been cleaved in the coiled-coil domain. In contrast, the β and γ chains, expressed separately or together, remained intracellular with formation of ββ or βγ, but not γγ, disulfide-linked dimers. Secretion of the β and γ chains required simultaneous expression of all three chains and their assembly into αβγ heterotrimers. Epitope-tagged recombinant α subunit and recombinant laminin were affinity-purified from the conditioned medium of αγ and αβγ clones. Rotary-shadow electron microscopy revealed that the free α subunit is a linear structure containing N-terminal and included globules with a foreshortened long arm, while the trimeric species has the typical four-arm morphology of native laminin. We conclude that the α chain can be delivered to the extracellular environment as a single subunit, whereas the β and γ chains cannot, and that the α chain drives the secretion of the trimeric molecule. Such an α-chain-dependent mechanism could allow for the regulation of laminin export into a nascent basement membrane, and might serve an important role in controlling basement membrane formation.

Quantitation of Pseudomonas aeruginosa in wound biopsy samples: from bacterial culture to rapid `real-time' polymerase chain reaction

We developed a real-time detection (RTD) polymerase chain reaction (PCR) with rapid thermal cycling to detect and quantify Pseudomonas aeruginosa in wound biopsy samples. This method produced a linear quantitative detection range of 7 logs, with a lower detection limit of 103 colony-forming units (CFU)/g tissue or a few copies per reaction. The time from sample collection to result was less than 1h. RTD-PCR has potential for rapid quantitative detection of pathogens in critical care patients, enabling early and individualized treatment.