Corrected kriging update formulae for batch-sequential data assimilation

Recently, a lot of effort has been spent in the efficient computation of kriging predictors when observations are assimilated sequentially. In particular, kriging update formulae enabling significant computational savings were derived. Taking advantage of the previous kriging mean and variance computations helps avoiding a costly matrix inversion when adding one observation to the TeX already available ones. In addition to traditional update formulae taking into account a single new observation, Emery (2009) proposed formulae for the batch-sequential case, i.e. when TeX new observations are simultaneously assimilated. However, the kriging variance and covariance formulae given in Emery (2009) for the batch-sequential case are not correct. In this paper, we fix this issue and establish correct expressions for updated kriging variances and covariances when assimilating observations in parallel. An application in sequential conditional simulation finally shows that coupling update and residual substitution approaches may enable significant speed-ups.

An analysis of variance type method to describe and compare steady states in clinical data

Identifying and comparing different steady states is an important task for clinical decision making. Data from unequal sources, comprising diverse patient status information, have to be interpreted. In order to compare results an expressive representation is the key. In this contribution we suggest a criterion to calculate a context-sensitive value based on variance analysis and discuss its advantages and limitations referring to a clinical data example obtained during anesthesia. Different drug plasma target levels of the anesthetic propofol were preset to reach and maintain clinically desirable steady state conditions with target controlled infusion (TCI). At the same time systolic blood pressure was monitored, depth of anesthesia was recorded using the bispectral index (BIS) and propofol plasma concentrations were determined in venous blood samples. The presented analysis of variance (ANOVA) is used to quantify how accurately steady states can be monitored and compared using the three methods of measurement.

Landscape connectivity, habitat structure and activity of bat guilds in farmland-dominated matrices

Agricultural intensification has caused a decline in structural elements in European farmland, where natural habitats are increasingly fragmented. The loss of habitat structures has a detrimental effect on biodiversity and affects bat species that depend on vegetation structures for foraging and commuting. We investigated the impact of connectivity and configuration of structural landscape elements on flight activity, species richness and diversity of insectivorous bats and distinguished three bat guilds according to species-specific bioacoustic characteristics. We tested whether bats with shorter-range echolocation were more sensitive to habitat fragmentation than bats with longer-range echolocation. We expected to find different connectivity thresholds for the three guilds and hypothesized that bats prefer linear over patchy landscape elements. Bat activity was quantified using repeated acoustic monitoring in 225 locations at 15 study plots distributed across the Swiss Central Plateau, where connectivity and the shape of landscape elements were determined by spatial analysis (GIS). Spectrograms of bat calls were assigned to species with the software batit by means of image recognition and statistical classification algorithms. Bat activity was significantly higher around landscape elements compared to open control areas. Short- and long-range echolocating bats were more active in well-connected landscapes, but optimal connectivity levels differed between the guilds. Species richness increased significantly with connectivity, while species diversity did not (Shannon's diversity index). Total bat activity was unaffected by the shape of landscape elements. Synthesis and applications. This study highlights the importance of connectivity in farmland landscapes for bats, with shorter-range echolocating bats being particularly sensitive to habitat fragmentation. More structurally diverse landscape elements are likely to reduce population declines of bats and could improve conditions for other declining species, including birds. Activity was highest around optimal values of connectivity, which must be evaluated for the different guilds and spatially targeted for a region's habitat configuration. In a multi-species approach, we recommend the reintroduction of structural elements to increase habitat heterogeneity should become part of agri-environment schemes.

On the use of Cramér-Rao minimum variance bounds for the design of magnetic resonance spectroscopy experiments

Localized Magnetic Resonance Spectroscopy (MRS) is in widespread use for clinical brain research. Standard acquisition sequences to obtain one-dimensional spectra suffer from substantial overlap of spectral contributions from many metabolites. Therefore, specially tuned editing sequences or two-dimensional acquisition schemes are applied to extend the information content. Tuning specific acquisition parameters allows to make the sequences more efficient or more specific for certain target metabolites. Cramér-Rao bounds have been used in other fields for optimization of experiments and are now shown to be very useful as design criteria for localized MRS sequence optimization. The principle is illustrated for one- and two-dimensional MRS, in particular the 2D separation experiment, where the usual restriction to equidistant echo time spacings and equal acquisition times per echo time can be abolished. Particular emphasis is placed on optimizing experiments for quantification of GABA and glutamate. The basic principles are verified by Monte Carlo simulations and in vivo for repeated acquisitions of generalized two-dimensional separation brain spectra obtained from healthy subjects and expanded by bootstrapping for better definition of the quantification uncertainties.

How large are actor and partner effects of personality on relationship satisfaction? The importance of controlling for shared method variance

Previous research suggests that the personality of a relationship partner predicts not only the individual’s own satisfaction with the relationship but also the partner’s satisfaction. Based on the actor–partner interdependence model, the present research tested whether actor and partner effects of personality are biased when the same method (e.g., self-report) is used for the assessment of personality and relationship satisfaction and, consequently, shared method variance is not controlled for. Data came from 186 couples, of whom both partners provided self- and partner reports on the Big Five personality traits. Depending on the research design, actor effects were larger than partner effects (when using only self-reports), smaller than partner effects (when using only partner reports), or of about the same size as partner effects (when using self- and partner reports). The findings attest to the importance of controlling for shared method variance in dyadic data analysis.

Relaxed trait covariance in interspecific cichlid hybrids predicts morphological diversity in adaptive radiations

The process of adaptive radiation involves multiple events of speciation in short succession, associated with ecological diversification. Understanding this process requires identifying the origins of heritable phenotypic variation that allows adaptive radiation to progress. Hybridization is one source of genetic and morphological variation that may spur adaptive radiation. We experimentally explored the potential role of hybridization in facilitating the onset of adaptive radiation. We generated first- and second-generation hybrids of four species of African cichlid fish, extant relatives of the putative ancestors of the adaptive radiations of Lakes Victoria and Malawi. We com- pared patterns in hybrid morphological variation with the variation in the lake radiations. We show that significant fractions of the interspecific mor- phological variation and the major trajectories in morphospace that charac- terize whole radiations can be generated in second-generation hybrids. Furthermore, we show that covariation between traits is relaxed in second- generation hybrids, which may facilitate adaptive diversification. These results support the idea that hybridization can provide the heritable pheno- typic diversity necessary to initiate adaptive radiation.

Thermal release rate studies of nuclear reaction products from polycrystalline metal matrices

The thermal release rate of nuclear reaction products was investigated in offline annealing experiments. This work was motivated by the search for a high melting catcher material for recoiling products from heavy ion induced nuclear fusion reactions. Polycrystalline refractory metal foils of Ni, Y, Zr, Nb, Mo, Hf, W, and Re were investigated as catcher metals. Diffusion data for various tracer/host combinations were deduced from the measured release rates. This work focuses on the diffusion and the release rate of volatile p-elements from row 5 and 6 of the periodic table as lighter homologues of the superheavy elements with Z ≥ 113 to be studied in future experiments. A massive radiation damage enhancement of the diffusion velocity was observed. Diffusion trends have been established along the groups and rows of the periodic table based on the dependence of diffusion velocity on atomic sizes.

PREDICTION OF DNA METHYLATION BASED ON GENOMIC ARCHITECTURE AND APPLICATIONS OF POSITIONAL WEIGHT MATRICES

Gene silencing due to epigenetic mechanisms shows evidence of significant contributions to cancer development. We hypothesis that the genetic architecture based on retrotransposon elements surrounding the transcription start site, plays an important role in the suppression and promotion of DNA methylation. In our investigation we found a high rate of SINE and LINEs retrotransposon elements near the transcription start site of unmethylated genes when compared to methylated genes. The presence of these elements were positively associated with promoter methylation, contrary to logical expectations, due to the malicious effects of retrotransposon elements which insert themselves randomly into the genome causing possible loss of gene function. In our genome wide analysis of human genes, results suggested that 22% of the genes in cancer were predicted to be methylation-prone; in cancer these genes are generally down-regulated and function in the development process. In summary, our investigation validated our hypothesis and showed that these widespread genomic elements in cancer are highly associated with promoter DNA methylation and may further participate in influencing epigenetic regulation.

A nonparametric method of comparing the partial areas under two correlated receiver operating characteristic curves

A non-parametric method was developed and tested to compare the partial areas under two correlated Receiver Operating Characteristic curves. Based on the theory of generalized U-statistics the mathematical formulas have been derived for computing ROC area, and the variance and covariance between the portions of two ROC curves. A practical SAS application also has been developed to facilitate the calculations. The accuracy of the non-parametric method was evaluated by comparing it to other methods. By applying our method to the data from a published ROC analysis of CT image, our results are very close to theirs. A hypothetical example was used to demonstrate the effects of two crossed ROC curves. The two ROC areas are the same. However each portion of the area between two ROC curves were found to be significantly different by the partial ROC curve analysis. For computation of ROC curves with large scales, such as a logistic regression model, we applied our method to the breast cancer study with Medicare claims data. It yielded the same ROC area computation as the SAS Logistic procedure. Our method also provides an alternative to the global summary of ROC area comparison by directly comparing the true-positive rates for two regression models and by determining the range of false-positive values where the models differ. ^

Protein cross-linking mediated by tissue transglutaminase correlates with the maturation of extracellular matrices during lung development

At birth, the mammalian lung is still immature. The alveoli are not yet formed and the interairspace walls contain two capillary layers which are separated by an interstitial core. After alveolarization (first 2 postnatal weeks in rats) the alveolar septa mature: their capillary layers merge, the amount of connective tissue decreases, and the mature lung parenchyma is formed (second and third week). During the first 3 wk of life the role of tissue transglutaminase (tTG) was studied in rat lung by immunostaining of cryostat and paraffin sections, by Northern and Western blotting, and by a quantitative determination of gamma-glutamyl-epsilon-lysine. While enzyme activity and intracellular tTG were already present before term, the enzyme product (gamma-glutamyl-epsilon-lysine-crosslink) and extracellular tTG appeared between postnatal days 10 and 19 in the lung parenchyma. In large blood vessels and large airways, which mature earlier than the parenchyma, both the enzyme product and extracellular tTG had already appeared at the end of the first postnatal week. We conclude that tTG is expressed and externalized into the extracellular matrix of lung shortly before maturation of an organ area. Because tTG covalently and irreversibly crosslinks extracellular matrix proteins, we hypothesize that it may prevent or delay further remodeling of basement membranes and may stabilize other extracellular components, such as microfibrils.