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.

Building the View Graph of a Category by Exploiting Image Realism

We propose a weakly supervised method to arrange images of a given category based on the relative pose between the camera and the object in the scene. Relative poses are points on a sphere centered at the object in a given canonical pose, which we call object viewpoints. Our method builds a graph on this sphere by assigning images with similar viewpoint to the same node and by connecting nodes if they are related by a small rotation. The key idea is to exploit a large unlabeled dataset to validate the likelihood of dominant 3D planes of the object geometry. A number of 3D plane hypotheses are evaluated by applying small 3D rotations to each hypothesis and by measuring how well the deformed images match other images in the dataset. Correct hypotheses will result in deformed images that correspond to plausible views of the object, and thus will likely match well other images in the same category. The identified 3D planes are then used to compute affinities between images related by a change of viewpoint. We then use the affinities to build a view graph via a greedy method and the maximum spanning tree.

LOGITCPRPLOT: Stata module to graph component-plus-residual plot for logistic regression

logitcprplot can be used after logistic regression for graphing a component-plus-residual plot (a.k.a. partial residual plot) for a given predictor, including a lowess, local polynomial, restricted cubic spline, fractional polynomial, penalized spline, regression spline, running line, or adaptive variable span running line smooth

CPRPLOT2: Stata module to graph component-plus-residual plots for transformed regressors

cprplot2 is a variation of official Stata's cprplot and is used for graphing component-plus-residual plots (a.k.a. partial residual plots). Additional features (compared to cprplot): (1) cprplot2 can handle variables that enter the model repeatedly via different transformations (for example, polynomials). (2) cprplot2 can display component-plus-residual plots using the original units for transformed variables in the model. (3) A wrapper is provided to quickly display several component-plus-residual plots in a single image.