Variance reduction methods for simulation of densities on Wiener space

We develop a general error analysis framework for the Monte Carlo simulationof densities for functionals in Wiener space. We also study variancereduction methods with the help of Malliavin derivatives. For this, wegive some general heuristic principles which are applied to diffusionprocesses. A comparison with kernel density estimates is made.

Heidelberg Retina Tomograph in Glaucoma Case-finding

Purpose: To compare the disease classification using optic nerve head imaging against clinical examination combined with perimetry in the case-finding environment.Methods: During a glaucoma screening event at the University Hospital in Lausanne 148 attendees opted to undergo an extended screening exam consisting of perimetry with Octopus (Haag-Streit, Koeniz, Switzerland), tonometry, gonioscopy, slit-lamp optic nerve head examination and Heidelberg Retinal Tomorgraph (HRT, Heidelberg engineering). Classification criteria of an abnormal clinical exam consisted of a cup to disc ratio > 0.6, other disc abnormality (e.g. notching), closed angles, Van Herrick < ¼ and intraocular pressure >21mmHg. In perimetry a square root loss variance (sLV) >3.0 dB was considered abnormal. Those cases with both an abnormal clinical exam and abnormal perimetric results were considered to be glaucoma suspects. An abnormal HRT exam was a global Moorfields Regression Analysis (MRA) result of "outside normal limits". The results from the worse eye for each instrument are reported.Results: The mean age of the patients was 59.9 years (SD ±14.8years). 46 subjects failed the clinical exam. 55 subjects had sLV>3.0dB. 86 subjects failed one or more of the testing components, 21 subjects failed both the clinical exam and perimetry. Two cases of advanced glaucoma were diagnosed on the day of the event. 20 subjects were referred due to other ocular pathology (10 cataract; 2 AMD; 8 other). 16 subjects were outside normal limits with the HRT, of which 11 failed the clinical exam, 9 failed the perimetric exam and 7 failed all three. HRT had a sensitivity of 33% CI [15%,57%] with a specificity of 93% CI [87%,97%].Conclusions: HRT shows good specificity, however the low sensitivity makes it of limited use in the proposed case finding scenario.

Duality in mean-variance frontiers with conditioning information

Portfolio and stochastic discount factor (SDF) frontiers are usually regarded as dual objects, and researchers sometimes use one to answer questions about the other. However, the introduction of conditioning information and active portfolio strategies alters this relationship. For instance, the unconditional portfolio frontier in Hansen and Richard (1987) is not dual to the unconditional SDF frontier in Gallant, Hansen and Tauchen (1990). We characterise the dual objects to those frontiers, and relate them to the frontiers generated with managed portfolios, which are commonly used in empirical work. We also study the implications of a safe asset and other special cases.

Weighted Euclidean biplots

We construct a weighted Euclidean distance that approximates any distance or dissimilarity measure between individuals that is based on a rectangular cases-by-variables data matrix. In contrast to regular multidimensional scaling methods for dissimilarity data, the method leads to biplots of individuals and variables while preserving all the good properties of dimension-reduction methods that are based on the singular-value decomposition. The main benefits are the decomposition of variance into components along principal axes, which provide the numerical diagnostics known as contributions, and the estimation of nonnegative weights for each variable. The idea is inspired by the distance functions used in correspondence analysis and in principal component analysis of standardized data, where the normalizations inherent in the distances can be considered as differential weighting of the variables. In weighted Euclidean biplots we allow these weights to be unknown parameters, which are estimated from the data to maximize the fit to the chosen distances or dissimilarities. These weights are estimated using a majorization algorithm. Once this extra weight-estimation step is accomplished, the procedure follows the classical path in decomposing the matrix and displaying its rows and columns in biplots.

NOTCH1 can initiate NF-κB activation via cytosolic interactions with components of the T cell Signalosome.

T cell stimulation requires the input and integration of external signals. Signaling through the T cell receptor (TCR) is known to induce formation of the membrane-tethered CBM complex, comprising CARMA1, BCL10, and MALT1, which is required for TCR-mediated NF-κB activation. TCR signaling has been shown to activate NOTCH proteins, transmembrane receptors also implicated in NF-κB activation. However, the link between TCR-mediated NOTCH signaling and early events leading to induction of NF-κB activity remains unclear. In this report, we demonstrate a novel cytosolic function for NOTCH1 and show that it is essential to CBM complex formation. Using a model of skin allograft rejection, we show in vivo that NOTCH1 acts in the same functional pathway as PKCθ, a T cell-specific kinase important for CBM assembly and classical NF-κB activation. We further demonstrate in vitro NOTCH1 associates physically with PKCθ and CARMA1 in the cytosol. Unexpectedly, when NOTCH1 expression was abrogated using RNAi approaches, interactions between CARMA1, BCL10, and MALT1 were lost. This failure in CBM assembly reduced inhibitor of kappa B alpha phosphorylation and diminished NF-κB-DNA binding. Finally, using a luciferase gene reporter assay, we show the intracellular domain of NOTCH1 can initiate robust NF-κB activity in stimulated T cells, even when NOTCH1 is excluded from the nucleus through modifications that restrict it to the cytoplasm or hold it tethered to the membrane. Collectively, these observations provide evidence that NOTCH1 may facilitate early events during T cell activation by nucleating the CBM complex and initiating NF-κB signaling.

Sample size and statistical power in the hierarchical analysis of variance: applications in morphometry of the nervous system.

Analysis of variance is commonly used in morphometry in order to ascertain differences in parameters between several populations. Failure to detect significant differences between populations (type II error) may be due to suboptimal sampling and lead to erroneous conclusions; the concept of statistical power allows one to avoid such failures by means of an adequate sampling. Several examples are given in the morphometry of the nervous system, showing the use of the power of a hierarchical analysis of variance test for the choice of appropriate sample and subsample sizes. In the first case chosen, neuronal densities in the human visual cortex, we find the number of observations to be of little effect. For dendritic spine densities in the visual cortex of mice and humans, the effect is somewhat larger. A substantial effect is shown in our last example, dendritic segmental lengths in monkey lateral geniculate nucleus. It is in the nature of the hierarchical model that sample size is always more important than subsample size. The relative weight to be attributed to subsample size thus depends on the relative magnitude of the between observations variance compared to the between individuals variance.

Cognitive efficacy of quetiapine in early-onset first-episode psychosis: a 12-week open label trial.

Twenty-three adolescents with psychotic disorders, aged from 13 to 18 years, participated in a 12-week open label trial (17 adolescents completed the study) in order to examine the impact of quetiapine on clinical status and cognitive functions (encompassing processing speed, attention, short-term memory, long-term memory and executive function). An improvement in Clinical Global Impression and Positive and Negative Symptom Scale (P's ≤ 0.001) was observed. In addition, after controlling for amelioration of symptoms, a significant improvement was observed on one executive function (P = 0.044; Trail Making Part B). The remaining cognitive abilities showed stability. In addition, we observed an interaction between quetiapine doses (>300 mg/day or <300 mg/day) and time, where lower doses showed more improvement in verbal short-term memory (P = 0.048), inhibition abilities (P = 0.038) and positive symptoms (P = 0.020). The neuropsychological functioning of adolescents with psychotic disorders remained mainly stable after 12 weeks of treatment with quetiapine. However, lower doses seemed to have a better impact on two components of cognition (inhibition abilities and verbal short-term memory) and on positive symptoms.

Role of NINJA in root jasmonate signaling.

Wound responses in plants have to be coordinated between organs so that locally reduced growth in a wounded tissue is balanced by appropriate growth elsewhere in the body. We used a JASMONATE ZIM DOMAIN 10 (JAZ10) reporter to screen for mutants affected in the organ-specific activation of jasmonate (JA) signaling in Arabidopsis thaliana seedlings. Wounding one cotyledon activated the reporter in both aerial and root tissues, and this was either disrupted or restricted to certain organs in mutant alleles of core components of the JA pathway including COI1, OPR3, and JAR1. In contrast, three other mutants showed constitutive activation of the reporter in the roots and hypocotyls of unwounded seedlings. All three lines harbored mutations in Novel Interactor of JAZ (NINJA), which encodes part of a repressor complex that negatively regulates JA signaling. These ninja mutants displayed shorter roots mimicking JA-mediated growth inhibition, and this was due to reduced cell elongation. Remarkably, this phenotype and the constitutive JAZ10 expression were still observed in backgrounds lacking the ability to synthesize JA or the key transcriptional activator MYC2. Therefore, JA-like responses can be recapitulated in specific tissues without changing a plant's ability to make or perceive JA, and MYC2 either has no role or is not the only derepressed transcription factor in ninja mutants. Our results show that the role of NINJA in the root is to repress JA signaling and allow normal cell elongation. Furthermore, the regulation of the JA pathway differs between roots and aerial tissues at all levels, from JA biosynthesis to transcriptional activation.

Adaptive frequency tracking of the baseline ECG identifies the site of atrial fibrillation termination by catheter ablation

Multiple organization indices have been used to predict the outcome of stepwise catheter ablation in long-standing persistent atrial fibrillation (AF), however with limited success. Our study aims at developinginnovative organization indices from baseline ECG (i.e. during the procedure, before ablation) in orderto identify the site of AF termination by catheter ablation. Seventeen consecutive male patients (age60 ± 5 years, AF duration 7 ± 5 years) underwent a stepwise catheter ablation. Chest lead V6 was placedin the back (V6b). QRST cancelation was performed from chest leads V1 to V6b. Using an innovativeadaptive harmonic frequency tracking, two measures of AF organization were computed to quantify theharmonics components of ECG activity: (1) the adaptive phase difference variance (APD) between theAF harmonic components as a measure of AF regularity, and (2) and adaptive organization index (AOI)evaluating the cyclicity of the AF oscillations. Both adaptive indices were compared to indices computedusing a time-invariant approach: (1) ECG AF cycle length (AFCL), (2) the spectrum based organizationindex (OI), and (3) the time-invariant phase difference TIPD. Long-standing persistent AF was terminatedinto sinus rhythm or atrial tachycardia in 13/17 patients during stepwise ablation, 11 during left atriumablation (left terminated patients - LT), 2 during the right atrium ablation (right terminated patients -RT), and 4 were non terminated (NT) and required electrical cardioversion. Our findings showed that LTpatients were best separated from RT/NT before ablation by the duration of sustained AF and by AOI onchest lead V1 and APD from the dorsal lead V6b as compared to ECG AFCL, OI and TIPD, respectively. Ourresults suggest that adaptive measures of AF organization computed before ablation perform better thantime-invariant based indices for identifying patients whose AF will terminate during ablation within theleft atrium. These findings are indicative of a higher baseline organization in these patients that could beused to select candidates for the termination of AF by stepwise catheter ablation.© 2013 Elsevier Ltd. All rights reserved.

Statistical analysis of maximum likelihood estimator images of human brain FDG PET studies

The work presented evaluates the statistical characteristics of regional bias and expected error in reconstructions of real positron emission tomography (PET) data of human brain fluoro-deoxiglucose (FDG) studies carried out by the maximum likelihood estimator (MLE) method with a robust stopping rule, and compares them with the results of filtered backprojection (FBP) reconstructions and with the method of sieves. The task of evaluating radioisotope uptake in regions-of-interest (ROIs) is investigated. An assessment of bias and variance in uptake measurements is carried out with simulated data. Then, by using three different transition matrices with different degrees of accuracy and a components of variance model for statistical analysis, it is shown that the characteristics obtained from real human FDG brain data are consistent with the results of the simulation studies.

Termination of Atrial Fibrillation by Catheter Ablation can be Successfully Predicted from Baseline ECG

Purpose Multiple organization indices (OIs) have been used to predict the outcome of stepwise catheter ablation (step-CA) in long-standing persistent AF (pers-AF), however with limited success. Our study aims at developing innovative OIs from baseline (BL, before ablation) ECG in order to predict the outcome of step-CA. Methods Fourteen consecutive patients (pts) (60±5 y, AF duration 21±9 m) underwent a step-CA consisting in pulmonary veins isolation, left atrial (LA) defragmentation and linear ablations, and right atrial (RA) ablations if non terminated. Chest lead V6 was placed in the back (V6b). After QRST cancellation from chest leads V1 to V6b, two OIs were computed to quantify the harmonic components of ECG atrial activity: 1) phase difference variance (PD) between the AF harmonic components as a measure of AF regularity, and 2) adaptive OI (AOI) evaluating the time evolution of the AF harmonic components. Both indices were compared to classical ones: a spectrum-based OI (SOI) and ECG AF cycle length (AFCL). Results Pers-AF was terminated into sinus rhythm or atrial tachycardia in 10/14 pts during step-CA, 8 during LA (LT), 2 during RA (RT) ablation, and 4 were non terminated (NT). The figure shows that LT was best separated from RT/NT before ablation by AOI computed on lead V1 (A) and PD from lead V6b (B) as compared to SOI and AFCL respectively. Conclusion Our results suggest that adaptive OIs computed before ablation perform better than classical OIs for separating LT from RT/NT pts. These findings are indicative of a higher baseline organization in LT pts that could be used to select candidates for the restoration of sinus rhythm by step-CA.

Drift compensation of gas sensor array data by common principal component analysis

A new drift compensation method based on Common Principal Component Analysis (CPCA) is proposed. The drift variance in data is found as the principal components computed by CPCA. This method finds components that are common for all gasses in feature space. The method is compared in classification task with respect to the other approaches published where the drift direction is estimated through a Principal Component Analysis (PCA) of a reference gas. The proposed new method ¿ employing no specific reference gas, but information from all gases ¿has shown the same performance as the traditional approach with the best-fitted reference gas. Results are shown with data lasting 7-months including three gases at different concentrations for an array of 17 polymeric sensors.