Comparative Analysis of Evolving Software Systems Using the Gini Coefficient

Software metrics offer us the promise of distilling useful information from vast amounts of software in order to track development progress, to gain insights into the nature of the software, and to identify potential problems. Unfortunately, however, many software metrics exhibit highly skewed, non-Gaussian distributions. As a consequence, usual ways of interpreting these metrics --- for example, in terms of "average" values --- can be highly misleading. Many metrics, it turns out, are distributed like wealth --- with high concentrations of values in selected locations. We propose to analyze software metrics using the Gini coefficient, a higher-order statistic widely used in economics to study the distribution of wealth. Our approach allows us not only to observe changes in software systems efficiently, but also to assess project risks and monitor the development process itself. We apply the Gini coefficient to numerous metrics over a range of software projects, and we show that many metrics not only display remarkably high Gini values, but that these values are remarkably consistent as a project evolves over time.

HIV-1 Transmission During Recent Infection and During Treatment Interruptions as Major Drivers of New Infections in the Swiss HIV Cohort Study.

BACKGROUND  Reducing the fraction of transmissions during recent human immunodeficiency virus (HIV) infection is essential for the population-level success of "treatment as prevention". METHODS  A phylogenetic tree was constructed with 19 604 Swiss sequences and 90 994 non-Swiss background sequences. Swiss transmission pairs were identified using 104 combinations of genetic distance (1%-2.5%) and bootstrap (50%-100%) thresholds, to examine the effect of those criteria. Monophyletic pairs were classified as recent or chronic transmission based on the time interval between estimated seroconversion dates. Logistic regression with adjustment for clinical and demographic characteristics was used to identify risk factors associated with transmission during recent or chronic infection. FINDINGS  Seroconversion dates were estimated for 4079 patients on the phylogeny, and comprised between 71 (distance, 1%; bootstrap, 100%) to 378 transmission pairs (distance, 2.5%; bootstrap, 50%). We found that 43.7% (range, 41%-56%) of the transmissions occurred during the first year of infection. Stricter phylogenetic definition of transmission pairs was associated with higher recent-phase transmission fraction. Chronic-phase viral load area under the curve (adjusted odds ratio, 3; 95% confidence interval, 1.64-5.48) and time to antiretroviral therapy (ART) start (adjusted odds ratio 1.4/y; 1.11-1.77) were associated with chronic-phase transmission as opposed to recent transmission. Importantly, at least 14% of the chronic-phase transmission events occurred after the transmitter had interrupted ART. CONCLUSIONS  We demonstrate a high fraction of transmission during recent HIV infection but also chronic transmissions after interruption of ART in Switzerland. Both represent key issues for treatment as prevention and underline the importance of early diagnosis and of early and continuous treatment.

Highly conservative reciprocal translocations formed by apparent joining of exchanged DNA double-strand break ends

Chromosomal translocations induced by ionizing radiation and radiomimetic drugs are thought to arise by incorrect joining of DNA double-strand breaks. To dissect such misrepair events at a molecular level, large-scale, bleomycin-induced rearrangements in the aprt gene of Chinese hamster ovary D422 cells were mapped, the breakpoints were sequenced, and the original non-aprt parental sequences involved in each rearrangement were recovered from nonmutant cells. Of seven rearrangements characterized, six were reciprocal exchanges between aprt and unrelated sequences. Consistent with a mechanism involving joining of exchanged double-strand break ends, there was, in most cases, no homology between the two parental sequences, no overlap in sequences retained at the two newly formed junctions, and little or no loss of parental sequences (usually ≤2 bp) at the breakpoints. The breakpoints were strongly correlated (P < 0.0001) with expected sites of bleomycin-induced, double-strand breaks. Fluorescence in situ hybridization indicated that, in six of the mutants, the rearrangement was accompanied by a chromosomal translocation at the aprt locus, because upstream and downstream flanking sequences were detected on separate chromosomes. The results suggest that repair of free radical-mediated, double-strand breaks in confluence-arrested cells is effected by a conservative, homology-independent, end-joining pathway that does not involve single-strand intermediate and that misjoining of exchanged ends by this pathway can directly result in chromosomal translocations.

A small asparagine-rich protein required for S-allele-specific pollen rejection in Nicotiana

Although S-locus RNases (S-RNases) determine the specificity of pollen rejection in self-incompatible (SI) solanaceous plants, they alone are not sufficient to cause S-allele-specific pollen rejection. To identify non-S-RNase sequences that are required for pollen rejection, a Nicotiana alata cDNA library was screened by differential hybridization. One clone, designated HT, hybridized strongly to RNA from N. alata styles but not to RNA from Nicotiana plumbaginifolia, a species known to lack one or more factors necessary for S-allele-specific pollen rejection. Sequence analysis revealed a 101-residue ORF including a putative secretion signal and an asparagine-rich domain near the C terminus. RNA blot analysis showed that the HT-transcript accumulates in the stigma and style before anthesis. The timing of HT-expression lags slightly behind SC10-RNase in SI N. alata SC10SC10 and is well correlated with the onset of S-allele-specific pollen rejection in the style. An antisense-HT construct was prepared to test for a role in pollen rejection. Transformed (N. plumbaginifolia × SI N. alata SC10SC10) hybrids with reduced levels of HT-protein continued to express SC10-RNase but failed to reject SC10-pollen. Control hybrids expressing both SC10-RNase and HT-protein showed a normal S-allele-specific pollen rejection response. We conclude that HT-protein is directly implicated in pollen rejection.

Protein Information Resource: a community resource for expert annotation of protein data

The Protein Information Resource, in collaboration with the Munich Information Center for Protein Sequences (MIPS) and the Japan International Protein Information Database (JIPID), produces the most comprehensive and expertly annotated protein sequence database in the public domain, the PIR-International Protein Sequence Database. To provide timely and high quality annotation and promote database interoperability, the PIR-International employs rule-based and classification-driven procedures based on controlled vocabulary and standard nomenclature and includes status tags to distinguish experimentally determined from predicted protein features. The database contains about 200 000 non-redundant protein sequences, which are classified into families and superfamilies and their domains and motifs identified. Entries are extensively cross-referenced to other sequence, classification, genome, structure and activity databases. The PIR web site features search engines that use sequence similarity and database annotation to facilitate the analysis and functional identification of proteins. The PIR-Inter­national databases and search tools are accessible on the PIR web site at http://pir.georgetown.edu/ and at the MIPS web site at http://www.mips.biochem.mpg.de. The PIR-International Protein Sequence Database and other files are also available by FTP.

Vasopressin mRNA localization in nerve cells: Characterization of cis-acting elements and trans-acting factors

mRNA localization is a complex pathway. Besides mRNA sorting per se, this process includes aspects of regulated translation. It requires protein factors that interact with defined sequences (or sequence motifs) of the transcript, and the protein/RNA complexes are finally guided along the cytoskeleton to their ultimate destinations. The mRNA encoding the vasopressin (VP) precursor protein is localized to the nerve cell processes in vivo and in primary cultured nerve cells. Sorting of VP transcripts to dendrites is mediated by the last 395 nucleotides of the mRNA, the dendritic localizer sequence, and it depends on intact microtubules. In vitro interaction studies with cytosolic extracts demonstrated specific binding of a protein, enriched in nerve cell tissues, to the radiolabeled dendritic localizer sequence probe. Biochemical purification revealed that this protein is the multifunctional poly(A)-binding protein (PABP). It is well known for its ability to bind with high affinity to poly(A) tails of mRNAs, prerequisite for mRNA stabilization and stimulation of translational initiation, respectively. With lower affinities, PABP can also associate with non-poly(A) sequences. The physiological consequences of these PABP/RNA interactions are far from clear but may include functions such as translational silencing. Presumably, the translational state of mRNAs subject to dendritic sorting is influenced by external stimuli. PABP thus could be a component required to regulate local synthesis of the VP precursor and possibly of other proteins.

Expression of error-prone polymerases in BL2 cells activated for Ig somatic hypermutation

High affinity antibodies are generated in mice and humans by means of somatic hypermutation (SHM) of variable (V) regions of Ig genes. Mutations with rates of 10−5–10−3 per base pair per generation, about 106-fold above normal, are targeted primarily at V-region hot spots by unknown mechanisms. We have measured mRNA expression of DNA polymerases ι, η, and ζ by using cultured Burkitt's lymphoma (BL)2 cells. These cells exhibit 5–10-fold increases in heavy-chain V-region mutations targeted only predominantly to RGYW (R = A or G, Y = C or T, W = T or A) hot spots if costimulated with T cells and IgM crosslinking, the presumed in vivo requirements for SHM. An ∼4-fold increase pol ι mRNA occurs within 12 h when cocultured with T cells and surface IgM crosslinking. Induction of pols η and ζ occur with T cells, IgM crosslinking, or both stimuli. The fidelity of pol ι was measured at RGYW hot- and non-hot-spot sequences situated at nicks, gaps, and double-strand breaks. Pol ι formed T⋅G mispairs at a frequency of 10−2, consistent with SHM-generated C to T transitions, with a 3-fold increased error rate in hot- vs. non-hot-spot sequences for the single-nucleotide overhang. The T cell and IgM crosslinking-dependent induction of pol ι at 12 h may indicate an SHM “triggering” event has occurred. However, pols ι, η, and ζ are present under all conditions, suggesting that their presence is not sufficient to generate mutations because both T cell and IgM stimuli are required for SHM induction.

Análise transcriptômica de miRNAs em pacientes sob dupla terapia de antiagregação

A terapia antiagregante é comumente indicada na prevenção e tratamento de doenças cardiovasculares. A dupla antiagregação com clopidrogrel e ácido acetilsalicílico (AAS) tem sido frequentemente adotada em pacientes com Doença Arterial Coronariana (DAC), mas apresenta ineficácia em uma parcela significativa da população com genótipo de respondedores. Essa falha terapêutica nos leva a questionar se outros mecanismos moleculares podem estar influenciando na resposta a esses fármacos. Recentes estudos sugerem que pequenas sequências de RNA não codificantes denominadas microRNAs (miRNAs) podem estar fortemente relacionadas com resposta ao tratamento fármaco-terapêutico, controlando as proteínas envolvidas na farmacocinética e farmacodinâmica. Entretanto, os principais miRNAs que atuam na dinâmica da resposta medicamentosa ainda não foram bem definidos. O objetivo deste estudo foi avaliar o perfil de miRNAs no sangue total periférico, procurando melhor esclarecer os mecanismos envolvidos na resposta aos antiagregantes plaquetários AAS e clopidogrel. Para isso, selecionou-se pacientes com DAC, os quais apresentavam diferentes respostas à dupla terapia de antiagregação determinadas pelo teste de agregação plaquetária. Baseados nos fenótipos, os perfis de expressão de miRNAs foram comparados entre os valores da taxa de agregação categorizados em tercis (T) de resposta. O grupo T1 foi constituído de pacientes respondedores, o T2 de respondedores intermediários e o T3 de não respondedores. Os perfis de miRNAs foram obtidos após sequenciamento de última geração e os dados obtidos foram analisados pelo pacote Deseq2. Os resultados mostraram 18 miRNAs diferentemente expressos entre os dois tercis extremos. Dentre esses miRNAs, 10 deles apresentaram importantes alvos relacionados com vias de ativação e agregação plaquetária quando analisados pelo software Ingenuity®. Dos 10 miRNAs, 4 deles, os quais apresentaram-se menos expressos no sequenciamento, demonstraram os mesmos perfis de expressão quando analisados pela reação em cadeia pela polimerase quantitativa (qPCR): hsa-miR-423-3p, hsa-miR-744-5p, hsa-miR- 30a-5p e hsa-let-7g-5p. A partir das análises de predição de alvos, pôde-se observar que os quatro miRNAs, quando menos expressos simultaneamente, predizem ativação da agregação plaquetária. Além disso, os miRNAs hsa-miR- 423-5p, hsa-miR-744-5p e hsa-let-7g-5p mostraram correlação com o perfil lipídico dos pacientes que, por sua vez, apresentou influência nos valores de agregação compreendidos no T3 de resposta a ambos os medicamentos. Sendo assim, conclui-se que maiores taxas de agregação plaquetária podem estar indiretamente relacionadas com os padrões de expressão de hsa-miR- 423-3p, hsa-miR-744-5p e hsa-let-7g-5p. Sugere-se que a avaliação do perfil de expressão destes 3 miRNAs no sangue periférico de pacientes com DAC possa predizer resposta terapêutica inadequada ao AAS e ao clopidogrel

Constraining storm-scale forecasts of deep convective initiation with surface weather observations

Thesis (Ph.D.)--University of Washington, 2016-06

Successive nonparametric estimation of conditional distributions

Spatial characterization of non-Gaussian attributes in earth sciences and engineering commonly requires the estimation of their conditional distribution. The indicator and probability kriging approaches of current nonparametric geostatistics provide approximations for estimating conditional distributions. They do not, however, provide results similar to those in the cumbersome implementation of simultaneous cokriging of indicators. This paper presents a new formulation termed successive cokriging of indicators that avoids the classic simultaneous solution and related computational problems, while obtaining equivalent results to the impractical simultaneous solution of cokriging of indicators. A successive minimization of the estimation variance of probability estimates is performed, as additional data are successively included into the estimation process. In addition, the approach leads to an efficient nonparametric simulation algorithm for non-Gaussian random functions based on residual probabilities.

Conditional quantum-state engineering using ancillary squeezed-vacuum states

We investigate an optical scheme to conditionally engineer quantum states using a beam splitter, homodyne detection, and a squeezed vacuum as an ancillar state. This scheme is efficient in producing non-Gaussian quantum states such as squeezed single photons and superpositions of coherent states (SCSs). We show that a SCS with well defined parity and high fidelity can be generated from a Fock state of n

Universal quantum computation with continuous-variable cluster states

We describe a generalization of the cluster-state model of quantum computation to continuous-variable systems, along with a proposal for an optical implementation using squeezed-light sources, linear optics, and homodyne detection. For universal quantum computation, a nonlinear element is required. This can be satisfied by adding to the toolbox any single-mode non-Gaussian measurement, while the initial cluster state itself remains Gaussian. Homodyne detection alone suffices to perform an arbitrary multimode Gaussian transformation via the cluster state. We also propose an experiment to demonstrate cluster-based error reduction when implementing Gaussian operations.

Towards a framework for combining stochastic and deterministic descriptions of nonstationary financial time series

We present in this paper ideas to tackle the problem of analysing and forecasting nonstationary time series within the financial domain. Accepting the stochastic nature of the underlying data generator we assume that the evolution of the generator's parameters is restricted on a deterministic manifold. Therefore we propose methods for determining the characteristics of the time-localised distribution. Starting with the assumption of a static normal distribution we refine this hypothesis according to the empirical results obtained with the methods anc conclude with the indication of a dynamic non-Gaussian behaviour with varying dependency for the time series under consideration.

Regularisation of mixture density networks

Mixture Density Networks are a principled method to model conditional probability density functions which are non-Gaussian. This is achieved by modelling the conditional distribution for each pattern with a Gaussian Mixture Model for which the parameters are generated by a neural network. This thesis presents a novel method to introduce regularisation in this context for the special case where the mean and variance of the spherical Gaussian Kernels in the mixtures are fixed to predetermined values. Guidelines for how these parameters can be initialised are given, and it is shown how to apply the evidence framework to mixture density networks to achieve regularisation. This also provides an objective stopping criteria that can replace the `early stopping' methods that have previously been used. If the neural network used is an RBF network with fixed centres this opens up new opportunities for improved initialisation of the network weights, which are exploited to start training relatively close to the optimum. The new method is demonstrated on two data sets. The first is a simple synthetic data set while the second is a real life data set, namely satellite scatterometer data used to infer the wind speed and wind direction near the ocean surface. For both data sets the regularisation method performs well in comparison with earlier published results. Ideas on how the constraint on the kernels may be relaxed to allow fully adaptable kernels are presented.

Conditional distribution modeling for estimating and exploiting uncertainty in control systems

This paper presents a general methodology for estimating and incorporating uncertainty in the controller and forward models for noisy nonlinear control problems. Conditional distribution modeling in a neural network context is used to estimate uncertainty around the prediction of neural network outputs. The developed methodology circumvents the dynamic programming problem by using the predicted neural network uncertainty to localize the possible control solutions to consider. A nonlinear multivariable system with different delays between the input-output pairs is used to demonstrate the successful application of the developed control algorithm. The proposed method is suitable for redundant control systems and allows us to model strongly non Gaussian distributions of control signal as well as processes with hysteresis.