Estimating Lorenz and concentration curves in Stata

Lorenz and concentration curves are widely used tools in inequality research. In this paper I present a new Stata command called -lorenz- that estimates Lorenz and concentration curves from individual-level data and, optionally, displays the results in a graph. The -lorenz- command supports relative as well as generalized, absolute, unnormalized, or custom-normalized Lorenz or concentration curves, and provides tools for computing contrasts between different subpopulations or outcome variables. Variance estimation for complex samples is fully supported.

Quick time planning using "s" curves and cost based durations.

ABSTRACT: The comparison of the different bids in the tender for a project, with the traditional contract system based on unit rates open to and re-measurement, requires analysis tools that are able to discriminate proposals having a similar overall economic impact, but that might show a very different behaviour during the execution of the works. RESUMEN: La estimación rápida de costes en fases iniciales del proyecto por métodos paramétricos y referencias estadísticas es un tema bien estudiado, divulgado y aplicado en el sector de la construcción. Sin embargo, existe poca literatura técnica sobre sistemas de predimensionado de tiempos, que permitan realizar rápidamente una planificación con un grado de aproximación razonable. Este texto reúne dos aspectos ya conocidos, pero hasta ahora independientes, y una aportación propia:  -La estimación del plazo final por referencias estadísticas (BCIS, 2000)  - La estimación del reparto del coste total a lo largo de la ejecución mediante curvas "S" (diversos autores)  La estimación de la duración de la ejecución de las actividades en función de su coste. El conjunto de estas tres técnicas, aplicadas a un proyecto, permite obtener una planificación con el suficiente grado de detalle y fiabilidad para tomar decisiones en fases iniciales del proyecto.

Bifurcation diagrams for polymer blends with diffuse interfaces in confined and adaptive geometries

Dynamics of binary mixtures such as polymer blends, and fluids near the critical point, is described by the model-H, which couples momentum transport and diffusion of the components [1]. We present an extended version of the model-H that allows to study the combined effect of phase separation in a polymer blend and surface structuring of the film itself [2]. We apply it to analyze the stability of vertically stratified base states on extended films of polymer blends and show that convective transport leads to new mechanisms of instability as compared to the simpler diffusive case described by the Cahn- Hilliard model [3, 4]. We carry out this analysis for realistic parameters of polymer blends used in experimental setups such as PS/PVME. However, geometrically more complicated states involving lateral structuring, strong deflections of the free surface, oblique diffuse interfaces, checkerboard modes, or droplets of a component above of the other are possible at critical composition solving the Cahn Hilliard equation in the static limit for rectangular domains [5, 6] or with deformable free surfaces [6]. We extend these results for off-critical compositions, since balanced overall composition in experiments are unusual. In particular, we study steady nonlinear solutions of the Cahn-Hilliard equation for bidimensional layers with fixed geometry and deformable free surface. Furthermore we distinguished the cases with and without energetic bias at the free surface. We present bifurcation diagrams for off-critical films of polymer blends with free surfaces, showing their free energy, and the L2-norms of surface deflection and the concentration field, as a function of lateral domain size and mean composition. Simultaneously, we look at spatial dependent profiles of the height and concentration. To treat the problem of films with arbitrary surface deflections our calculations are based on minimizing the free energy functional at given composition and geometric constraints using a variational approach based on the Cahn-Hilliard equation. The problem is solved numerically using the finite element method (FEM).

Structure and function of the Bacillus hybrid enzyme GluXyn-1: Native-like jellyroll fold preserved after insertion of autonomous globular domain

The 1,3–1,4-β-glucanase from Bacillus macerans (wtGLU) and the 1,4-β-xylanase from Bacillus subtilis (wtXYN) are both single-domain jellyroll proteins catalyzing similar enzymatic reactions. In the fusion protein GluXyn-1, the two proteins are joined by insertion of the entire XYN domain into a surface loop of cpMAC-57, a circularly permuted variant of wtGLU. GluXyn-1 was generated by protein engineering methods, produced in Escherichia coli and shown to fold spontaneously and have both enzymatic activities at wild-type level. The crystal structure of GluXyn-1 was determined at 2.1 Å resolution and refined to R = 17.7% and R(free) = 22.4%. It shows nearly ideal, native-like folding of both protein domains and a small, but significant hinge bending between the domains. The active sites are independent and accessible explaining the observed enzymatic activity. Because in GluXyn-1 the complete XYN domain is inserted into the compact folding unit of GLU, the wild-type-like activity and tertiary structure of the latter proves that the folding process of GLU does not depend on intramolecular interactions that are short-ranged in the sequence. Insertion fusions of the GluXyn-1 type may prove to be an easy route toward more stable bifunctional proteins in which the two parts are more closely associated than in linear end-to-end protein fusions.

Tumor necrosis factor (TNF)-mediated kinase cascades: Bifurcation of Nuclear Factor-κB and c-jun N-terminal kinase (JNK/SAPK) pathways at TNF receptor-associated factor 2

TNF-induced activation of the transcription factor NF-κB and the c-jun N-terminal kinase (JNK/SAPK) requires TNF receptor-associated factor 2 (TRAF2). The NF-κB-inducing kinase (NIK) associates with TRAF2 and mediates TNF activation of NF-κB. Herein we show that NIK interacts with additional members of the TRAF family and that this interaction requires the conserved “WKI” motif within the TRAF domain. We also investigated the role of NIK in JNK activation by TNF. Whereas overexpression of NIK potently induced NF-κB activation, it failed to stimulate JNK activation. A kinase-inactive mutant of NIK was a dominant negative inhibitor of NF-κB activation but did not suppress TNF- or TRAF2-induced JNK activation. Thus, TRAF2 is the bifurcation point of two kinase cascades leading to activation of NF-κB and JNK, respectively.

In vitro and in vivo differentiation into B cells, T cells, and myeloid cells of primitive yolk sac hematopoietic precursor cells expanded >100-fold by coculture with a clonal yolk sac endothelial cell line

The yolk sac, first site of hematopoiesis during mammalian development, contains not only hematopoietic stem cells but also the earliest precursors of endothelial cells. We have previously shown that a nonadherent yolk sac cell population (WGA+, density <1.077, AA4.1+) can give rise to B cells, T cells, and myeloid cells both in vitro and in vivo. We now report on the ability of a yolk sac-derived cloned endothelial cell line (C166) to provide a suitable microenvironment for expansion of these early precursor cells. Single day 10 embryonic mouse yolk sac hematopoietic stem cells were expanded >100 fold within 8 days by coculture with irradiated C166 cells. Colony-forming ability was retained for at least three passages in vitro, with retention of the ability to differentiate into T-cell, B-cell, and myeloid lineages. Stem cell properties were maintained by a significant fraction of nonadherent cells in the third passage, although these stem cells expressed a somewhat more mature cell surface phenotype than the initial yolk sac stem cells. When reintroduced into adult allogeneic immunocompromised (scid) hosts, they were able to give rise to all of the leukocyte lineages, including T cells, B cells, and myeloid cells. We conclude that yolk sac endothelial cells can support the stable proliferation of multipotential hematopoietic stem cells, thus generating adequate numbers of cells for study of the mechanisms involved in their subsequent development and differentiation, for in vivo hematopoietic restitution, and for potential use as a vehicle for gene transfer.

How does a β-hairpin fold/unfold? Competition between topology and heterogeneity in a solvable model

We study the competition between topological effects and sequence inhomogeneities in determining the thermodynamics and the un/folding kinetics of a β-hairpin. Our work utilizes a new exactly solvable model that allows for arbitrary configurations of native contacts. In general, the competition between heterogeneity and topology results in a crossover of the dominant transition state. Interestingly, near this crossover, the single reaction coordinate picture can be seriously misleading. Our results also suggest that inferring the folding pathway from unfolding simulations is not always justified.

“Coarse” stability and bifurcation analysis using time-steppers: A reaction-diffusion example

Evolutionary, pattern forming partial differential equations (PDEs) are often derived as limiting descriptions of microscopic, kinetic theory-based models of molecular processes (e.g., reaction and diffusion). The PDE dynamic behavior can be probed through direct simulation (time integration) or, more systematically, through stability/bifurcation calculations; time-stepper-based approaches, like the Recursive Projection Method [Shroff, G. M. & Keller, H. B. (1993) SIAM J. Numer. Anal. 30, 1099–1120] provide an attractive framework for the latter. We demonstrate an adaptation of this approach that allows for a direct, effective (“coarse”) bifurcation analysis of microscopic, kinetic-based models; this is illustrated through a comparative study of the FitzHugh-Nagumo PDE and of a corresponding Lattice–Boltzmann model.

Reversibly locking a protein fold in an active conformation with a disulfide bond: Integrin αL I domains with high affinity and antagonist activity in vivo

The integrin αLβ2 has three different domains in its headpiece that have been suggested to either bind ligand or to regulate ligand binding. One of these, the inserted or I domain, has a fold similar to that of small G proteins. The I domain of the αM and α2 subunits has been crystallized in both open and closed conformations; however, the αL I domain has been crystallized in only the closed conformation. We hypothesized that the αL domain also would have an open conformation, and that this would be the ligand binding conformation. Therefore, we introduced pairs of cysteine residues to form disulfides that would lock the αL I domain in either the open or closed conformation. Locking the I domain open resulted in a 9,000-fold increase in affinity to intercellular adhesion molecule-1 (ICAM-1), which was reversed by disulfide reduction. By contrast, the affinity of the locked closed conformer was similar to wild type. Binding completely depended on Mg2+. Orders of affinity were ICAM-1 > ICAM-2 > ICAM-3. The kon, koff, and KD values for the locked open I domain were within 1.5-fold of values previously determined for the αLβ2 complex, showing that the I domain is sufficient for full affinity binding to ICAM-1. The locked open I domain antagonized αLβ2-dependent adhesion in vitro, lymphocyte homing in vivo, and firm adhesion but not rolling on high endothelial venules. The ability to reversibly lock a protein fold in an active conformation with dramatically increased affinity opens vistas in therapeutics and proteomics.

Parametrizations of elliptic curves by Shimura curves and by classical modular curves

Fix an isogeny class

Euler characteristics and elliptic curves

Let E be a modular elliptic curve over ℚ, without complex multiplication; let p be a prime number where E has good ordinary reduction; and let F∞ be the field obtained by adjoining to ℚ all p-power division points on E. Write G∞ for the Galois group of F∞ over ℚ. Assume that the complex L-series of E over ℚ does not vanish at s = 1. If p ⩾ 5, we make a precise conjecture about the value of the G∞-Euler characteristic of the Selmer group of E over F∞. If one makes a standard conjecture about the behavior of this Selmer group as a module over the Iwasawa algebra, we are able to prove our conjecture. The crucial local calculations in the proof depend on recent joint work of the first author with R. Greenberg.

A joint hazard and time scaling model to compare survival curves.

To provide a more general method for comparing survival experience, we propose a model that independently scales both hazard and time dimensions. To test the curve shape similarity of two time-dependent hazards, h1(t) and h2(t), we apply the proposed hazard relationship, h12(tKt)/ h1(t) = Kh, to h1. This relationship doubly scales h1 by the constant hazard and time scale factors, Kh and Kt, producing a transformed hazard, h12, with the same underlying curve shape as h1. We optimize the match of h12 to h2 by adjusting Kh and Kt. The corresponding survival relationship S12(tKt) = [S1(t)]KtKh transforms S1 into a new curve S12 of the same underlying shape that can be matched to the original S2. We apply this model to the curves for regional and local breast cancer contained in the National Cancer Institute's End Results Registry (1950-1973). Scaling the original regional curves, h1 and S1 with Kt = 1.769 and Kh = 0.263 produces transformed curves h12 and S12 that display congruence with the respective local curves, h2 and S2. This similarity of curve shapes suggests the application of the more complete curve shapes for regional disease as templates to predict the long-term survival pattern for local disease. By extension, this similarity raises the possibility of scaling early data for clinical trial curves according to templates of registry or previous trial curves, projecting long-term outcomes and reducing costs. The proposed model includes as special cases the widely used proportional hazards (Kt = 1) and accelerated life (KtKh = 1) models.