Fast Update of Conditional Simulation Ensembles

Gaussian random field (GRF) conditional simulation is a key ingredient in many spatial statistics problems for computing Monte-Carlo estimators and quantifying uncertainties on non-linear functionals of GRFs conditional on data. Conditional simulations are known to often be computer intensive, especially when appealing to matrix decomposition approaches with a large number of simulation points. This work studies settings where conditioning observations are assimilated batch sequentially, with one point or a batch of points at each stage. Assuming that conditional simulations have been performed at a previous stage, the goal is to take advantage of already available sample paths and by-products to produce updated conditional simulations at mini- mal cost. Explicit formulae are provided, which allow updating an ensemble of sample paths conditioned on n ≥ 0 observations to an ensemble conditioned on n + q observations, for arbitrary q ≥ 1. Compared to direct approaches, the proposed formulae proveto substantially reduce computational complexity. Moreover, these formulae explicitly exhibit how the q new observations are updating the old sample paths. Detailed complexity calculations highlighting the benefits of this approach with respect to state-of-the-art algorithms are provided and are complemented by numerical experiments.

A note on adding objects to an existing twoway graph

In Stata, graphs are usually generated by one call to the graph command. Sometimes, however, it would be convenient to be able to add objects to a graph after the graph has been created. In this article, I provide a command called addplot that offers such functionality for twoway graphs, capitalizing on an undocumented feature of Stata's graphics system.

Biological and protective properties of immune sera directed to the influenza virus neuraminidase

The envelope of influenza A viruses contains two large antigens, hemagglutinin (HA) and neuraminidase (NA). Conventional influenza virus vaccines induce neutralizing antibodies that are predominantly directed to the HA globular head, a domain that is subject to extensive antigenic drift. Antibodies directed to NA are induced at much lower levels, probably as a consequence of the immunodominance of the HA antigen. Although antibodies to NA may affect virus release by inhibiting the sialidase function of the glycoprotein, the antigen has been largely neglected in past vaccine design. In this study, we characterized the protective properties of monospecific immune sera that were generated by vaccination with recombinant RNA replicon particles encoding NA. These immune sera inhibited hemagglutination in an NA subtype-specific and HA subtype-independent manner and interfered with infection of MDCK cells. In addition, they inhibited the sialidase activities of various influenza viruses of the same and even different NA subtypes. With this, the anti-NA immune sera inhibited the spread of H5N1 highly pathogenic avian influenza virus and HA/NA-pseudotyped viruses in MDCK cells in a concentration-dependent manner. When chickens were immunized with NA recombinant replicon particles and subsequently infected with low-pathogenic avian influenza virus, inflammatory serum markers were significantly reduced and virus shedding was limited or eliminated. These findings suggest that NA antibodies can inhibit virus dissemination by interfering with both virus attachment and egress. Our results underline the potential of high-quality NA antibodies for controlling influenza virus replication and place emphasis on NA as a vaccine antigen. IMPORTANCE The neuraminidase of influenza A viruses is a sialidase that acts as a receptor-destroying enzyme facilitating the release of progeny virus from infected cells. Here, we demonstrate that monospecific anti-NA immune sera inhibited not only sialidase activity, but also influenza virus hemagglutination and infection of MDCK cells, suggesting that NA antibodies can interfere with virus attachment. Inhibition of both processes, virus release and virus binding, may explain why NA antibodies efficiently blocked virus dissemination in vitro and in vivo. Anti-NA immune sera showed broader reactivity than anti-HA sera in hemagglutination inhibition tests and demonstrated cross-subtype activity in sialidase inhibition tests. These remarkable features of NA antibodies highlight the importance of the NA antigen for the development of next-generation influenza virus vaccines.

MASCG: Multi-Atlas Segmentation Constrained Graph method for accurate segmentation of hip CT images

This paper addresses the issue of fully automatic segmentation of a hip CT image with the goal to preserve the joint structure for clinical applications in hip disease diagnosis and treatment. For this purpose, we propose a Multi-Atlas Segmentation Constrained Graph (MASCG) method. The MASCG method uses multi-atlas based mesh fusion results to initialize a bone sheetness based multi-label graph cut for an accurate hip CT segmentation which has the inherent advantage of automatic separation of the pelvic region from the bilateral proximal femoral regions. We then introduce a graph cut constrained graph search algorithm to further improve the segmentation accuracy around the bilateral hip joint regions. Taking manual segmentation as the ground truth, we evaluated the present approach on 30 hip CT images (60 hips) with a 15-fold cross validation. When the present approach was compared to manual segmentation, an average surface distance error of 0.30 mm, 0.29 mm, and 0.30 mm was found for the pelvis, the left proximal femur, and the right proximal femur, respectively. A further look at the bilateral hip joint regions demonstrated an average surface distance error of 0.16 mm, 0.21 mm and 0.20 mm for the acetabulum, the left femoral head, and the right femoral head, respectively.

Binary recombinase systems for high-resolution conditional mutagenesis.

Conditional mutagenesis using Cre recombinase expressed from tissue specific promoters facilitates analyses of gene function and cell lineage tracing. Here, we describe two novel dual-promoter-driven conditional mutagenesis systems designed for greater accuracy and optimal efficiency of recombination. Co-Driver employs a recombinase cascade of Dre and Dre-respondent Cre, which processes loxP-flanked alleles only when both recombinases are expressed in a predetermined temporal sequence. This unique property makes Co-Driver ideal for sequential lineage tracing studies aimed at unraveling the relationships between cellular precursors and mature cell types. Co-InCre was designed for highly efficient intersectional conditional transgenesis. It relies on highly active trans-splicing inteins and promoters with simultaneous transcriptional activity to reconstitute Cre recombinase from two inactive precursor fragments. By generating native Cre, Co-InCre attains recombination rates that exceed all other binary SSR systems evaluated in this study. Both Co-Driver and Co-InCre significantly extend the utility of existing Cre-responsive alleles.

Altered directed functional connectivity in temporal lobe epilepsy in the absence of interictal spikes: A high density EEG study.

OBJECTIVE In patients with epilepsy, seizure relapse and behavioral impairments can be observed despite the absence of interictal epileptiform discharges (IEDs). Therefore, the characterization of pathologic networks when IEDs are not present could have an important clinical value. Using Granger-causal modeling, we investigated whether directed functional connectivity was altered in electroencephalography (EEG) epochs free of IED in left and right temporal lobe epilepsy (LTLE and RTLE) compared to healthy controls. METHODS Twenty LTLE, 20 RTLE, and 20 healthy controls underwent a resting-state high-density EEG recording. Source activity was obtained for 82 regions of interest (ROIs) using an individual head model and a distributed linear inverse solution. Granger-causal modeling was applied to the source signals of all ROIs. The directed functional connectivity results were compared between groups and correlated with clinical parameters (duration of the disease, age of onset, age, and learning and mood impairments). RESULTS We found that: (1) patients had significantly reduced connectivity from regions concordant with the default-mode network; (2) there was a different network pattern in patients versus controls: the strongest connections arose from the ipsilateral hippocampus in patients and from the posterior cingulate cortex in controls; (3) longer disease duration was associated with lower driving from contralateral and ipsilateral mediolimbic regions in RTLE; (4) aging was associated with a lower driving from regions in or close to the piriform cortex only in patients; and (5) outflow from the anterior cingulate cortex was lower in patients with learning deficits or depression compared to patients without impairments and to controls. SIGNIFICANCE Resting-state network reorganization in the absence of IEDs strengthens the view of chronic and progressive network changes in TLE. These resting-state connectivity alterations could constitute an important biomarker of TLE, and hold promise for using EEG recordings without IEDs for diagnosis or prognosis of this disorder.

PSHARE: Stata module to compute and graph percentile shares

-pshare- computes and graphs percentile shares from individual level data. Percentile shares are often used in inequality research to study the distribution of income or wealth. They are defined as differences between Lorenz ordinates of the outcome variable. Technically, the observations are sorted in increasing order of the outcome variable and the specified percentiles are computed from the running sum of the outcomes. Percentile shares are then computed as differences between percentiles, divided by total outcome. pshare requires moremata to be installed on the system; see ssc describe moremata.

Analysis of BMP signalling through the receptor type IA in embryogenesis using conditional gene inactivation in mice

Although bone morphogenetic proteins (BMPs) were initially identified for their potent bone-inducing activity, their precise roles in processes of endochondral and intramembranous bone formation are far from being clear. Tissue-specific loss-of-function experiments using the BMP receptor type IA (BMPR-IA) are particularly attractive since this receptor is thought to be essential for signaling by the closely related BMPs -2, 4, and 7. To ablate signaling through this receptor during chondrogenesis, we have generated transgenic mice expressing Cre recombinase under the control of the collagen type II (Col2a1) gene regulatory sequences. Mice lacking BMPR-IA function in chondrocytes display a number of skeletal abnormalities, including defects in bones of the chondrocranium, abnormal dorsal vertebral processes, scapulae with severe hypoplasia of dorsal elements, and shortening of the long bones. Alterations in the growth plate of long bones in mutants suggest that BMPR-IA is not required for early steps of the chondrocyte specification, but is rather important in regulation of terminal differentiation. Molecular analysis revealed noticeable downregulation of the Ihh/Ptch signalling pathway, decreased chondrocyte proliferation rate and deregulation of hypertrophy. ^ In order to elucidate the role of BMP signalling in development of the limb and intramembranous ossification, we have used mice expressing Cre recombinase under control of the Prx1 (MHox) regulatory elements (M. Logan, pers comm.). Cre activity was found in those mice in the developing limb bud mesenchyme, as well as in a subset of cranial neural crest cells. Prx1-Cre-induced conditional mutants display prominent defects in distal limb outgrowth, as well as ossification defects in a number of neural crest-derived calvarial bones. Intriguingly, mutant limbs displayed alterations in patterning along all three axes. Molecular analysis revealed ectopic anterior Shh/Ptch signalling pathway activation and expression of some Hox genes. Observed loss of Msx1 and Msx2 expression in the progress zone correlates with downregulation of Cyclin D1 and decreased distal outgrowth. Abnormal ventral localization of Lmx1b-expressing cells along with observed later morphological abnormalities suggest a novel role for BMP signalling in establishment or maintaining of the dorso-ventral polarity in the limb mesoderm. ^

Conditional expression of Sry, and Wnt4 by gene-targeting: Studies in sexual and skeletal development

Sry and Wnt4 cDNAs were individually introduced into the ubiquitously-expressed Rosa26 ( R26) locus by gene targeting in embryonic stem (ES) cells to create a conditional gene expression system in mice. In the targeted alleles, expression of these cDNAs should be blocked by a neomycin resistance selection cassette that is flanked by loxP sites. Transgene expression should be activated after the blocking cassette is deleted by Cre recombinase. ^ To test this conditional expression system, I have bred R26-stop- Sry and R26-stop-Wnt4 heterozygotes with a MisRII-Cre mouse line that expresses Cre in the gonads of both sexes. Analysis of these two types of bigenic heterozygotes indicated that their gonads developed normally like those of wild types. However, one XX R26-Sry/R26-Sry; MisR2-Cre/+ showed epididymis-like structures resembling those of males. In contrast, only normal phenotypes were observed in XY R26-Wnt4/R26-Wnt4; MisR2-Cre /+ mice. To interpret these results, I have tested for Cre recombinase activity by Southern blot and transcription of the Sry and Wnt4 transgenes by RT-PCR. Results showed that bigenic mutants had insufficient activation of the transgenes in their gonads at E12.5 and E13.5. Therefore, the failure to observe mutant phenotypes may have resulted from low activity of MisR2-Cre recombination at the appropriate time. ^ Col2a1-Cre transgenic mice express Cre in differentiating chondrocytes. R26-Wnt4; Col2a1-Cre bigenic heterozygous mice were found to exhibit a dramatic alteration in growth presumably caused by Wnt4 overexpression during chondrogenesis. R26-Wnt4; Col2a1-Cre mice exhibited dwarfism beginning approximately 10 days after birth. In addition, they also had craniofacial abnormalities, and had delayed ossification of the lumbar vertebrate and pelvic bones. Histological analysis of the growth plates of R26-Wnt4; Col2a1-Cre mice revealed less structural organization and a delay in onset of the primary and secondary ossification centers. Molecular studies confirmed that overexpression of Wnt4 causes decreased proliferation and early maturation of chondrocytes. In addition, R26-Wnt4; Col2a1-Cre mice had decreased expression of vascular endothelial growth factor (VEGF), suggesting that defects in vascularization may contribute to the dwarf phenotype. Finally, 9-month-old R26-Wnt4; Col2a1-Cre mice had significantly more fat cells in the marrow cavities of their metaphysis long bones, implying that long-term overexpression of Wnt4may cause bone marrow pathologies. In conclusion, Wnt4 was activated by Col2a1-Cre recombinase and was overexpressed in the growth plate, resulting in aberrant proliferation and differentiation of chondrocytes, and ultimately leads to dwarfism in mice. ^

Empirical Analysis of Credit Risk Regime Switching and Temporal Conditional Default Correlation in Credit Default Swap Valuation: The Market liquidity effect

In this paper, we extend the debate concerning Credit Default Swap valuation to include time varying correlation and co-variances. Traditional multi-variate techniques treat the correlations between covariates as constant over time; however, this view is not supported by the data. Secondly, since financial data does not follow a normal distribution because of its heavy tails, modeling the data using a Generalized Linear model (GLM) incorporating copulas emerge as a more robust technique over traditional approaches. This paper also includes an empirical analysis of the regime switching dynamics of credit risk in the presence of liquidity by following the general practice of assuming that credit and market risk follow a Markov process. The study was based on Credit Default Swap data obtained from Bloomberg that spanned the period January 1st 2004 to August 08th 2006. The empirical examination of the regime switching tendencies provided quantitative support to the anecdotal view that liquidity decreases as credit quality deteriorates. The analysis also examined the joint probability distribution of the credit risk determinants across credit quality through the use of a copula function which disaggregates the behavior embedded in the marginal gamma distributions, so as to isolate the level of dependence which is captured in the copula function. The results suggest that the time varying joint correlation matrix performed far superior as compared to the constant correlation matrix; the centerpiece of linear regression models.

Total Factor Productivity and Monetary Policy: Evidence from Conditional Volatility

This paper empirically assesses whether monetary policy affects real economic activity through its affect on the aggregate supply side of the macroeconomy. Analysts typically argue that monetary policy either does not affect the real economy, the classical dichotomy, or only affects the real economy in the short run through aggregate demand new Keynesian or new classical theories. Real business cycle theorists try to explain the business cycle with supply-side productivity shocks. We provide some preliminary evidence about how monetary policy affects the aggregate supply side of the macroeconomy through its affect on total factor productivity, an important measure of supply-side performance. The results show that monetary policy exerts a positive and statistically significant effect on the supply-side of the macroeconomy. Moreover, the findings buttress the importance of countercyclical monetary policy as well as support the adoption of an optimal money supply rule. Our results also prove consistent with the effective role of monetary policy in the Great Moderation as well as the more recent rise in productivity growth.

Analyzing the function of myogenin in adult satellite cells through the construction of a myogenin conditional allele

Although mechanisms regulating the formation of embryonic skeletal muscle are well characterized, less is known about muscle formation in postnatal life. This disparity is unfortunate because the largest increases in skeletal muscle mass occur after birth. Adult muscle stem cells (satellite cells) appear to recapitulate the events that occur in embryonic myoblasts. In particular, the myogenic basic helix-loop-helix factors, which have crucial functions in embryonic muscle development, are assumed to have similar roles in postnatal muscle formation. Here, I test this assumption by determining the role of the myogenic regulator myogenin in postnatal life. Myogenin-null mice die at birth, necessitating the generation of floxed alleles of myogenin and the use of cre-recombinase lines to delete myogenin. Removing myogenin before embryonic muscle development resulted in myofiber deficiencies identical to those observed in myogenin-null mice. However, mice in which myogenin was deleted following embryonic muscle development had normal skeletal muscle, except for modest alterations in MRF4 and MyoD expression. Notably, myogenin-deleted mice were 30% smaller than controls, suggesting that myogenin's absence disrupted general body growth. These results suggest that skeletal muscle growth in postnatal life is controlled by mechanisms distinct from those occurring in embryonic muscle development. ^