Possible solution of some essential zero problems in compositional data analysis

One of the tantalising remaining problems in compositional data analysis lies in how to deal with data sets in which there are components which are essential zeros. By an essential zero we mean a component which is truly zero, not something recorded as zero simply because the experimental design or the measuring instrument has not been sufficiently sensitive to detect a trace of the part. Such essential zeros occur in many compositional situations, such as household budget patterns, time budgets, palaeontological zonation studies, ecological abundance studies. Devices such as nonzero replacement and amalgamation are almost invariably ad hoc and unsuccessful in such situations. From consideration of such examples it seems sensible to build up a model in two stages, the first determining where the zeros will occur and the second how the unit available is distributed among the non-zero parts. In this paper we suggest two such models, an independent binomial conditional logistic normal model and a hierarchical dependent binomial conditional logistic normal model. The compositional data in such modelling consist of an incidence matrix and a conditional compositional matrix. Interesting statistical problems arise, such as the question of estimability of parameters, the nature of the computational process for the estimation of both the incidence and compositional parameters caused by the complexity of the subcompositional structure, the formation of meaningful hypotheses, and the devising of suitable testing methodology within a lattice of such essential zero-compositional hypotheses. The methodology is illustrated by application to both simulated and real compositional data

System Design: UML Activity and Sequence Diagrams

In this session we look at how to model flow of control and interactions between components using UML Activity and Sequence Diagrams. This is an introductory session and so for Activity Diagrams we only cover branching, forks and joins and swim lanes, and for Sequence we cover lifelines, messages and returns, and alt, par and opt frames.

Identification of core and variable components of the Salmonella enterica subspecies I genome by microarray

We have performed microarray hybridization studies on 40 clinical isolates from 12 common serovars within Salmonella enterica subspecies I to identify the conserved chromosomal gene pool. We were able to separate the core invariant portion of the genome by a novel mathematical approach using a decision tree based on genes ranked by increasing variance. All genes within the core component were confirmed using available sequence and microarray information for S. enterica subspecies I strains. The majority of genes within the core component had conserved homologues in Escherichia coli K-12 strain MG1655. However, many genes present in the conserved set which were absent or highly divergent in K-12 had close homologues in pathogenic bacteria such as Shigella flexneri and Pseudomonas aeruginosa. Genes within previously established virulence determinants such as SPI1 to SPI5 were conserved. In addition several genes within SPI6, all of SPI9, and three fimbrial operons (fim, bcf, and stb) were conserved within all S. enterica strains included in this study. Although many phage and insertion sequence elements were missing from the core component, approximately half the pseudogenes present in S. enterica serovar Typhi were conserved. Furthermore, approximately half the genes conserved in the core set encoded hypothetical proteins. Separation of the core and variant gene sets within S. enterica subspecies I has offered fundamental biological insight into the genetic basis of phenotypic similarity and diversity across S. enterica subspecies I and shown how the core genome of these pathogens differs from the closely related E. coli K-12 laboratory strain.

Variability of dayside high latitude convection associated with a sequence of auroral transients

10 second resolution ionospheric convection data covering the invariant latitude range from 71° to 76°, obtained by using the EISCAT UHF and VHF radars, are combined with optical data from Ny Ålesund during a sequence of auroral transients in the post-noon sector (∼ 15 MLT). Satellite observations of polar cap convection patterns suggest negative BZ and BY components of the interplanetary magnetic field. Burst-like enhancements of westward (sunward) post-noon convection were accompanied by eastward moving auroral forms at higher latitudes, above the convection reversal boundary. In this case the background convection was weak, whereas the integrated potential drop across the radar field-of-view associated with the westward flow bursts was typically ∼ 20-35 kV. The auroral phenomenon consists of a series of similar events with a mean repetition period of 8 min. A close correlation between the auroral activity and convection enhancements in the cleft ionosphere is demonstrated.

The effect of clamping sequence on dimensional variability in sheet metal assembly

Finite element (FE) modelling techniques have become a popular tool for exploring welding and clamping sequence dependence in sheet metal assemblies. In the current paper, the dimensional variability associated with different assembly clamping sequences is investigated with a FE contact modelling approach implemented in the commercial code Abaqus. A simplified channel section assembly consisting of a top hat and bottom plate is the case study investigated. Expected variation modes of bow and twist were used to simulate key variability sources in the main structural component under investigation; the top hat of the channel section. It was found that final assembly variability can change considerably depending on clamp sequence selection. It was also found that different clamp sequences can control particular modes of variation better than others, and that there is not one particular clamping sequence that is the best for containing all variation modes. An adaptable assembly process is therefore suggested, where given the shape of input components the best available clamping sequence is selected. Comparison of the performance of the proposed adaptable clamping sequence to traditional fixed clamping sequences shows improvements for the dimensional control of variability in non-rigid components. While introducing such a method in production would require inspection of each component being assembled and investigation of the alternative clamping sequences, given access to fast and detailed dimensional inspection technology such as optical coordinate measuring machines (OCMM's), the approach shows promise for future application.

A unified learning algorithm to extract principal and minor components

Recently, many unified learning algorithms have been developed to solve the task of principal component analysis (PCA) and minor component analysis (MCA). These unified algorithms can be used to extract principal component and if altered simply by the sign, it can also serve as a minor component extractor. This is of practical significance in the implementations of algorithms. Convergence of the existing unified algorithms is guaranteed only under the condition that the learning rates of algorithms approach zero, which is impractical in many practical applications. In this paper, we propose a unified PCA & MCA algorithm with a constant learning rate, and derive the sufficient conditions to guarantee convergence via analyzing the discrete-time dynamics of the proposed algorithm. The achieved theoretical results lay a solid foundation for the applications of our proposed algorithm.

Multi-envelope HIV-1 vaccine devoid of SIV components controls disease in macaques challenged with heterologous pathogenic SHIV

A central obstacle to the design of a global HIV-1 vaccine is virus diversity. Pathogen diversity is not unique to HIV-1, and has been successfully conquered in other fields by the creation of vaccine cocktails. Here we describe the testing of an HIV-1 envelope cocktail vaccine. Six macaques received the vaccine, delivered by successive immunizations with recombinant DNA, recombinant vaccinia virus and recombinant envelope proteins. Following vaccination, animals developed a diversity of anti-envelope antibody binding and neutralizing activities toward proteins and viruses that were not represented by sequence in the vaccine. T-cells were also elicited, as measured by gamma-interferon production assays with envelope-derived peptide pools. Vaccinated and control animals were then challenged with the heterologous pathogenic SHIV, 89.6P. Vaccinated monkeys experienced significantly lower virus titers and better maintenance of CD4+ T-cells than unvaccinated controls. The B- and T-cell immune responses were far superior post-challenge in the vaccinated group. Four of six vaccinated animals and only one of six control animals survived a 44-week observation period post-challenge. The present report is the first to describe pathogenic SHIV disease control mediated by a heterologous HIV-1 vaccine, devoid of 89.6 or SIV derivatives.

Rapid detection of SMARCB1 sequence variation using high resolution melting

Background : Rhabdoid tumors are rare cancers of early childhood arising in the kidney, central nervous system and other organs. The majority are caused by somatic inactivating mutations or deletions affecting the tumor suppressor locus SMARCB1 [OMIM 601607]. Germ-line SMARCB1 inactivation has been reported in association with rhabdoid tumor, epitheloid sarcoma and familial schwannomatosis, underscoring the importance of accurate mutation screening to ascertain recurrence and transmission risks. We describe a rapid and sensitive diagnostic screening method, using high resolution melting (HRM), for detecting sequence variations in SMARCB1. Methods : Amplicons, encompassing the nine coding exons of SMARCB1, flanking splice site sequences and the 5' and 3' UTR, were screened by both HRM and direct DNA sequencing to establish the reliability of HRM as a primary mutation screening tool. Reaction conditions were optimized with commercially available HRM mixes. Results : The false negative rate for detecting sequence variants by HRM in our sample series was zero. Nine amplicons out of a total of 140 (6.4%) showed variant melt profiles that were subsequently shown to be false positive. Overall nine distinct pathogenic SMARCB1 mutations were identified in a total of 19 possible rhabdoid tumors. Two tumors had two distinct mutations and two harbored SMARCB1 deletion. Other mutations were nonsense or frame-shifts. The detection sensitivity of the HRM screening method was influenced by both sequence context and specific nucleotide change and varied from 1: 4 to 1:1000 (variant to wild-type DNA). A novel method involving digital HRM, followed by re-sequencing, was used to confirm mutations in tumor specimens containing associated normal tissue. Conclusions : This is the first report describing SMARCB1 mutation screening using HRM. HRM is a rapid, sensitive and inexpensive screening technology that is likely to be widely adopted in diagnostic laboratories to facilitate whole gene mutation screening.

A sequence based dynamic SOM model for text clustering

Text clustering can be considered as a four step process consisting of feature extraction, text representation, document clustering and cluster interpretation. Most text clustering models consider text as an unordered collection of words. However the semantics of text would be better captured if word sequences are taken into account.

In this paper we propose a sequence based text clustering model where four novel sequence based components are introduced in each of the four steps in the text clustering process.

Experiments conducted on the Reuters dataset and Sydney Morning Herald (SMH) news archives demonstrate the advantage of the proposed sequence based model, in terms of capturing context with semantics, accuracy and speed, compared to clustering of documents based on single words and n-gram based models.

The metabolic syndrome and cancer: is the metabolic syndrome useful for predicting cancer risk above and beyond its individual components?

AIMS: The metabolic syndrome (MetS) is a risk factor for cancer. However, it is not known if the MetS confers a greater cancer risk than the sum of its individual components, which components drive the association, or if the MetS predicts future cancer risk. MATERIALS AND METHODS: We linked 20,648 participants from the Australian and New Zealand Diabetes and Cancer Collaboration with complete data on the MetS to national cancer registries and used Cox proportional hazards models to estimate associations of the MetS, the number of positive MetS components, and each of the five MetS components separately with the risk for overall, colorectal, prostate and breast cancer. Hazard ratios (HR) and 95% confidence intervals (95%CI) are reported. We assessed predictive ability of the MetS using Harrell's c-statistic. RESULTS: The MetS was inversely associated with prostate cancer (HR 0.85; 95% CI 0.72-0.99). We found no evidence of an association between the MetS overall, colorectal and breast cancers. For those with five positive MetS components the HR was 1.12 (1.02-1.48) and 2.07 (1.26-3.39) for overall, and colorectal cancer, respectively, compared with those with zero positive MetS components. Greater waist circumference (WC) (1.38; 1.13-1.70) and elevated blood pressure (1.29; 1.01-1.64) were associated with colorectal cancer. Elevated WC and triglycerides were (inversely) associated with prostate cancer. MetS models were only poor to moderate discriminators for all cancer outcomes. CONCLUSIONS: We show that the MetS is (inversely) associated with prostate cancer, but is not associated with overall, colorectal or breast cancer. Although, persons with five positive components of the MetS are at a 1.2 and 2.1 increased risk for overall and colorectal cancer, respectively, and these associations appear to be driven, largely, by elevated WC and BP. We also demonstrate that the MetS is only a moderate discriminator of cancer risk.

The contribution of 700,000 ORF sequence tags to the definition of the human transcriptome

open reading frame expressed sequences tags (ORESTES) differ from conventional ESTs by providing sequence data from the central protein coding portion of transcripts. We generated a total of 696,745 ORESTES sequences from 24 human tissues and used a subset of the data that correspond to a set of 15,095 full-length mRNAs as a means of assessing the efficiency of the strategy and its potential contribution to the definition of the human transcriptome. We estimate that ORESTES sampled over 80% of all highly and moderately expressed, and between 40% and 50% of rarely expressed, human genes. In our most thoroughly sequenced tissue, the breast, the 130,000 ORESTES generated are derived from transcripts from an estimated 70% of all genes expressed in that tissue, with an equally efficient representation of both highly and poorly expressed genes. In this respect, we find that the capacity of the ORESTES strategy both for gene discovery and shotgun transcript sequence generation significantly exceeds that of conventional ESTs. The distribution of ORESTES is such that many human transcripts are now represented by a scaffold of partial sequences distributed along the length of each gene product. The experimental joining of the scaffold components, by reverse transcription-PCR, represents a direct route to transcript finishing that may represent a useful alternative to full-length cDNA cloning.

Podzolization as a deferralitization process: dynamics and chemistry of ground and surface waters in an Acrisol-Podzol sequence of the upper Amazon Basin

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Isolation and sequence determination of peptides in the venom of the spider wasp (Cyphononyx dorsalis) guided by matrix-assisted laser desorption/ionization time of flight (MALDI-TOF) mass spectrometry

Micro-scale (sub-pmol) isolation and sequence determination of three peptides from the venom of the solitary spider wasp Cyphononyx dorsalis is described. We isolated two novel peptides Cd-125 and Cd-146 and a known peptide Thr(6)-bradykinin from only two venom sacs of solitary spider wasp Cyphononyx dorsalis without bioassay-guided fractionation. but instead guided by MALDI-TOF MS. The MALDI-TOF MS analysis of each fraction showed the purity and molecular weight of the components, which led to the isolation of the peptides virtually without loss of sample amount. The sequences of the novel peptides Cd-125 (Asp-Thr-Ala-Arg-Leu-Lys-Trp-His) and Cd-146 (Ser-Glu-Thr-Gly-Asn-Thr-Val-Thr-Val-Lys-Gly-Phe-Ser-Pro-Leu-Arg) were determined by Edman degradation together with mass spectrometry. and finally corroborated by solid-phase synthesis. The known peptide Thr(6)-bradykinin (Arg-Pro-Pro-Gly-Phe-Thr-Pro-Phe-Arg) was identified by comparison with the synthetic authentic specimen. This is the first example for any kinins to be found in Pompilidae wasp venoms. The procedure reported here can be applicable to studies on many other components of solitary wasp venoms with limited sample availability. (C) 2001 Elsevier B.V. Ltd. All rights reserved.

Evaluation of the components of production of the ricinus according of rates of calcareous rock and match

The soil of the Paraiba, in generality, are acid and with low levels of available match, seriously limiting the income of the plants. A research in vegetation house was lead, with the objective to evaluate doses of calcareous soil and match in the culture of ricinus. The treatments had been composites for five levels of match: 0.0; 80; 160; 240 and 360 kg ha(-1) of P(2)O(5) and four dosis of calcareous rock: 0.0; 2.5; 3.5 and 4.5 t ha(-1) in experimental delineation of blocks casualized with 4 factorial project x 3 (four levels of match and three doses of calcareous soil) more two treatments you add zero of match (0.0 P(2)O(5): 2,5 t ha(-1) of calcareous soil) and calcareous soil zero (240 P(2)O(5): 0.0 t ha(-1) of calcareous soil), with three repetitions. The calcareous soil reduced the acidity of the ground and effect of the application of the match (>= 80 kg ha(-1) of P(2)O(5)) in the presence of the collagen (>= 2.5 t ha(-1) of CaCO(3)) with of the calcareous soil (>= 2.5 t ha(-1) of CaCO(3)) in the presence of the match was not verified well (>= 80 kg ha(-1) of P(2)O(5)). The match applied in the absence of the calcareous soil was limited to the income of the culture of what the calcareous soil in the absence of the match.

Zero mode effect generalization for the electromagnetic current in the light front

We consider in this work the electromagnetic current for a system composed of two charged bosons and show that it has a structure of many bodies even in the impulse approximation, when described in the light-front time x+. In terms of the two-body component for the bound state, the current contains two-body operators. We consider the photon interacting with two bosons and the process of pair creation connected to this interaction, interpreting it as a zero mode contribution to the current and discuss the consequences of this pair creation to the components of currents in the light front. © 2013 American Physical Society.