Equivalence of truncated count mixture distributions and mixtures of truncated count distributions

This article is about modeling count data with zero truncation. A parametric count density family is considered. The truncated mixture of densities from this family is different from the mixture of truncated densities from the same family. Whereas the former model is more natural to formulate and to interpret, the latter model is theoretically easier to treat. It is shown that for any mixing distribution leading to a truncated mixture, a (usually different) mixing distribution can be found so. that the associated mixture of truncated densities equals the truncated mixture, and vice versa. This implies that the likelihood surfaces for both situations agree, and in this sense both models are equivalent. Zero-truncated count data models are used frequently in the capture-recapture setting to estimate population size, and it can be shown that the two Horvitz-Thompson estimators, associated with the two models, agree. In particular, it is possible to achieve strong results for mixtures of truncated Poisson densities, including reliable, global construction of the unique NPMLE (nonparametric maximum likelihood estimator) of the mixing distribution, implying a unique estimator for the population size. The benefit of these results lies in the fact that it is valid to work with the mixture of truncated count densities, which is less appealing for the practitioner but theoretically easier. Mixtures of truncated count densities form a convex linear model, for which a developed theory exists, including global maximum likelihood theory as well as algorithmic approaches. Once the problem has been solved in this class, it might readily be transformed back to the original problem by means of an explicitly given mapping. Applications of these ideas are given, particularly in the case of the truncated Poisson family.

The amphioxus FoxQ1 gene is expressed in the developing endostyle

The FoxQ1 genes form a distinct group within the Fox (also known as forkhead) gene family. We have isolated a gene from the amphioxus Branchiostoma floridae that encodes a forkhead domain with high identity to FoxQ1 genes in other chordates. Molecular phylogenetic analysis places AmphiFoxQ1 in a robust grouping with vertebrate FoxQ1 genes and with Ciona intestinalis Ci-FoxQ1. This group is separate from that containing AmphiFoxQ2, which instead groups with other invertebrate Fox genes. The expression of AmphiFoxQ1 was analysed by whole mount in situ hybridisation. The results show that AmphiFoxQ1 expression is confined to the developing endoderm, and specifically marks the endostyle and associated peripharyngeal bands of amphioxus larvae. Ci-FoxQ1 is also expressed in the endostyle, highlighting this as a conserved site of FoxQ1 gene expression in basal chordates. (C) 2004 Elsevier B.V. All rights reserved.

Oestrogenic activity of p-hydroxybenzoic acid (common metabolite of paraben esters) and methylparaben in human breast cancer cell lines

This paper addresses the question of whether p-hydroxybenzoic acid, the common metabolite of parabens, possesses oestrogenic activity in human breast cancer cell lines. The alkyl esters of p-hydroxybenzoic acid (parabens) are used widely as preservatives in consumer products to which the human population is exposed and have been shown previously to possess oestrogenic activity and to be present in human breast tumour tissue, which is an oestrogen-responsive tissue. Recent work has shown p-hydroxybenzoic acid to give an oestrogenic response in the rodent uterotrophic assay. We report here that p-hydroxybenzoic acid possesses oestrogenic activity in a panel of assays in human breast cancer cell lines. p-Hydroxybenzoic acid was able to displace [H-3]oestradiol from cytosolic oestrogen receptor of MCF7 human breast cancer cells by 54% at 5 x 10(6)-fold molar excess and by 99% at 10(7)-fold molar excess. It was able to increase the expression of a stably integrated oestrogen responsive reporter gene (ERE-CAT) at a concentration of 5 x 10(-4) M in MCF7 cells after 24 h and 7 days, which could be inhibited by the anti-oestrogen ICI 182 780 (Faslodex, fulvestrant). Proliferation of two human breast cancer cell lines (MCF7, ZR-75-1) could be increased by 10(-5) M p-hydroxybenzoic acid. Following on from previous studies showing a decrease in oestrogenic activity of parabens with shortening of the linear alkyl chain length, this study has compared the oestrogenic activity of p-hydroxybenzoic acid where the alkyl grouping is no longer present with methylparaben, which has the shortest alkyl group. Intrinsic oestrogenic activity of p-hydroxybenzoic acid was similar to that of methylparaben in terms of relative binding to the oestrogen receptor but its oestrogenic activity on gene expression and cell proliferation was lower than that of methylparaben. It can be concluded that removal of the ester group from parabens does not abrogate its oestrogenic activity and that p-hydroxybenzoic acid can give oestrogenic responses in human breast cancer cells. Copyright (C) 2005 John Wiley & Sons, Ltd.

Comparative evaluation of a new effective population size estimator based on approximate Bayesian computation

We describe and evaluate a new estimator of the effective population size (N-e), a critical parameter in evolutionary and conservation biology. This new "SummStat" N-e. estimator is based upon the use of summary statistics in an approximate Bayesian computation framework to infer N-e. Simulations of a Wright-Fisher population with known N-e show that the SummStat estimator is useful across a realistic range of individuals and loci sampled, generations between samples, and N-e values. We also address the paucity of information about the relative performance of N-e estimators by comparing the SUMMStat estimator to two recently developed likelihood-based estimators and a traditional moment-based estimator. The SummStat estimator is the least biased of the four estimators compared. In 32 of 36 parameter combinations investigated rising initial allele frequencies drawn from a Dirichlet distribution, it has the lowest bias. The relative mean square error (RMSE) of the SummStat estimator was generally intermediate to the others. All of the estimators had RMSE > 1 when small samples (n = 20, five loci) were collected a generation apart. In contrast, when samples were separated by three or more generations and Ne less than or equal to 50, the SummStat and likelihood-based estimators all had greatly reduced RMSE. Under the conditions simulated, SummStat confidence intervals were more conservative than the likelihood-based estimators and more likely to include true N-e. The greatest strength of the SummStat estimator is its flexible structure. This flexibility allows it to incorporate any, potentially informative summary statistic from Population genetic data.

Phylogenetic relationships within the moss family Bryaceae based on chloroplast DNA evidence

The Bryaceae are a large cosmopolitan family of mosses containing genera of considerable taxonomic difficulty. Phylogenetic relationships within the family were inferred using data from chloroplast DNA sequences (rps4 and trnL-trnF region). Parsimony and maximum likelihood optimality criteria, and Bayesian phylogenetic inference procedures were employed to reconstruct relationships. The genera Bryum and Brachymenium are not monophyletic groups. A clade comprising Plagiobryum, Acidodontium, Mielichhoferia macrocarpa, Bryum sects. Bryum, Apalodictyon, Limbata, Leucodontium, Caespiticia, Capillaria (in part: sect. Capillaria), and Brachymenium sect. Dicranobryum, is well supported in all analyses and represents a major lineage within the family. Section Dicranobryum of Brachymenium is more closely related to section Bryum than to the other sections of Brachymenium, as are Mielichhoferia macrocarpa and M. himalayana. Species of Acidodontium form a clade with Anomobryum julaceum. The grouping of species with a rosulate gametophytic growth form suggests the presence of a 'rosulate' clade similar in circumscription to the genus Rosulabryum. Mielichhoferia macrocarpa and M. himalayana are transferred to Bryum as B. porsildii and B. caucasicum, respectively.

Oestrogenic activity of benzylparaben

Previous work has demonstrated that the alkyl esters of p-hydroxybenzoic acid (parabens) possess oestrogenic activity, which increases with length of alkyl chain from methylparaben to n-butylparaben and with branching in the alkyl chain from n-butylparaben to isobutylparaben. This study reports on the oestrogenic activity of benzylparaben in a variety of assays in vitro and in vivo. Benzylparaben was able to displace [H-3]oestradiol from cytosolic oestrogen receptor (ER) of MCF7 human breast cancer cells by 22% at 1000-fold molar excess, by 40% at 10000-fold molar excess, by 57% at 100000-fold molar excess and by 100% at 1000000-fold molar excess. It was able to increase expression of a stably transfected oestrogen responsive reporter gene (ERE-CAT) in MCF7 cells after 24 h at 10(-5)M/10(-4)M and after 7 days at 10(-6)M/10(-5)M/10(-4)M. Proliferation of MCF7 cells could be increased by 10(-6)M/10(-5)M benzylparaben and this could be inhibited by 10(-7)M pure anti-oestrogen ICI 182,780, indicating that growth effects were ER mediated. Further evidence for ER-mediation was provided from the ability of benzylparaben to increase the growth of a second oestrogen-dependent human breast cancer cell line ZR-75-1, but not the oestrogen-insensitive NIDA-MB-231 cell line. When tested in the presence of 10(-10)M 17beta-oestradiol, benzylparaben gave no antagonist response on the growth of either MCF7 or ZR-75-1 cells. Finally, benzylparaben could increase uterine weight in the immature mouse following topical application of three daily doses of 33 mg to dorsal skin. These results demonstrate that the oestrogenicity of methylparaben can be increased by the addition of an aryl group as well as by lengthening or branching the alkyl grouping. Copyright (C) 2003 John Wiley Sons, Ltd.

Formate-dependent growth and homoacetogenic fermentation by a bacterium from human feces: Description of Bryantella formatexigens gen. nov., sp nov.

Formate stimulates growth of a new bacterium from human feces. With high formate, it ferments glucose to acetate via the Wood-Ljungdahl pathway. The original isolate fermented vegetable cellulose and carboxymethylcellulose, but it lost this ability after storage at -76degreesC. 16S rRNA gene sequencing identifies it as a distinct line within the Clostridium coccoides supra-generic rRNA grouping. We propose naming it Bryantella formatexigens gen. nov., sp. nov.

Hespellia stercorisuis gen. nov., sp nov and Hespellia porcina sp nov., isolated from swine manure storage pits

Four Gram-positive-staining, strictly anaerobic, non-spore-forming, rod-shaped organisms were isolated from a pig manure storage pit. Comparative 16S rRNA gene sequence analysis revealed that the isolates belonged to two related but distinct groups. Sequence analysis showed that the two groups of isolates were highly related to each other (approx. 97% 16S rRNA gene sequence similarity), forming a distinct cluster within the Clostridium coccoides suprageneric rDNA grouping. Biochemical and physiological studies confirmed the division of the isolates into two related, albeit distinct, groups. Based on both phenotypic and phylogenetic evidence, it is proposed that the unidentified rod-shaped isolates from pig manure should be classified in a novel genus, Hespellia gen. nov., as Hespellia stercorisuis sp. nov. and Hespellia porcina sp. nov. The type species of the novel genus is H. stercorisuis (type strain, PC18(T) = NRRL B-23456(T) = CCUG 46279(T) = ATCC BAA-677(T)) and the type strain of H. porcina is PC80(T) (= NRRL B-23458(T) = ATCC BAA-674(T)).

Colonic microbiota signatures across five northern European countries

The composition of the colonic microbiota of 91 northern Europeans was characterized by fluorescent in situ hybridization using 18 phylogenetic probes. On average 75% of the bacteria were identified, and large interindividual variations were observed. Clostridium coccoides and Clostridium leptum were the dominant groups (28.0% and 25.2%), followed by the Bacteroides (8.5%). According to principal component analysis, no significant grouping with respect to geographic origin, age, or gender was observed.

Fastidiosipila sanguinis gen. nov., sp nov., a new Gram-positive, coccus-shaped organism from human blood

Phenotypic and phylogenetic studies were performed on two strains of an unidentified Gram-positive, fastidious, non-spore-forming, coccus-shaped bacterium recovered from human blood. The organism was catalase-negative and grew under strictly anaerobic conditions and in the presence of 2 and 6% O-2. Comparative 16S rRNA gene sequencing demonstrated that the unidentified bacterium was, phylogenetically, far removed from peptostreptococci and related Gram-positive coccus-shaped organisms, but exhibited a phylogenetic association with Clostridium rRNA cluster III [as defined by Collins et al, Int J Syst Bacteriol 44 (1994), 812-826]. Sequence divergence values of 15% or more were observed between the unidentified bacterium and all other recognized species within this and related rRINIA clostridial clusters. Treeing analysis showed that the unknown bacterium formed a deep line branching at the periphery of rRNA cluster III and represents a hitherto unknown genus within this supra-generic grouping. On the basis of both phylogenetic and phenotypic evidence, it is proposed that the unknown bacterium from blood be classified in a new genus, Fastidiosipila gen. nov., as Fastidiosipila sanguinis sp, nov. The type strain of Fastidiosipila sanguinis is CCUG 47711(T) (= CIP 108292(T)).

A sparse kernel density estimation algorithm using forward constrained regression

Using the classical Parzen window (PW) estimate as the target function, the sparse kernel density estimator is constructed in a forward constrained regression manner. The leave-one-out (LOO) test score is used for kernel selection. The jackknife parameter estimator subject to positivity constraint check is used for the parameter estimation of a single parameter at each forward step. As such the proposed approach is simple to implement and the associated computational cost is very low. An illustrative example is employed to demonstrate that the proposed approach is effective in constructing sparse kernel density estimators with comparable accuracy to that of the classical Parzen window estimate.

Robust classification of SAR imagery

In this work the G(A)(0) distribution is assumed as the universal model for amplitude Synthetic Aperture (SAR) imagery data under the Multiplicative Model. The observed data, therefore, is assumed to obey a G(A)(0) (alpha; gamma, n) law, where the parameter n is related to the speckle noise, and (alpha, gamma) are related to the ground truth, giving information about the background. Therefore, maps generated by the estimation of (alpha, gamma) in each coordinate can be used as the input for classification methods. Maximum likelihood estimators are derived and used to form estimated parameter maps. This estimation can be hampered by the presence of corner reflectors, man-made objects used to calibrate SAR images that produce large return values. In order to alleviate this contamination, robust (M) estimators are also derived for the universal model. Gaussian Maximum Likelihood classification is used to obtain maps using hard-to-deal-with simulated data, and the superiority of robust estimation is quantitatively assessed.

Example uses of a learning object

A self study course for learning to program using the C programming language has been developed. A Learning Object approach was used in the design of the course. One of the benefits of the Learning Object approach is that the learning material can be reused for different purposes. 'Me course developed is designed so that learners can choose the pedagogical approach most suited to their personal learning requirements. For all learning approaches a set of common Assessment Learning Objects (ALOs or tests) have been created. The design of formative assessments with ALOs can be carried out by the Instructional Designer grouping ALOs to correspond to a specific assessment intention. The course is non-credit earning, so there is no summative assessment, all assessment is formative. In this paper examples of ALOs and their uses is presented together with their uses as decided by the Instructional Designer and learner. Personalisation of the formative assessment of skills can be decided by the Instructional Designer or the learner using a repository of pre-designed ALOs. The process of combining ALOs can be carried out manually or in a semi-automated way using metadata that describes the ALO and the skill it is designed to assess.

M-estimator, and D-optimality model construction using orthogonal forward regression

This correspondence introduces a new orthogonal forward regression (OFR) model identification algorithm using D-optimality for model structure selection and is based on an M-estimators of parameter estimates. M-estimator is a classical robust parameter estimation technique to tackle bad data conditions such as outliers. Computationally, The M-estimator can be derived using an iterative reweighted least squares (IRLS) algorithm. D-optimality is a model structure robustness criterion in experimental design to tackle ill-conditioning in model Structure. The orthogonal forward regression (OFR), often based on the modified Gram-Schmidt procedure, is an efficient method incorporating structure selection and parameter estimation simultaneously. The basic idea of the proposed approach is to incorporate an IRLS inner loop into the modified Gram-Schmidt procedure. In this manner, the OFR algorithm for parsimonious model structure determination is extended to bad data conditions with improved performance via the derivation of parameter M-estimators with inherent robustness to outliers. Numerical examples are included to demonstrate the effectiveness of the proposed algorithm.

Rejoinder to 'Performance of different synchronisation measures in real data: A case study on electroencephalograhic signals

We agree with Duckrow and Albano [Phys. Rev. E 67, 063901 (2003)] and Quian Quiroga et al. [Phys. Rev. E 67, 063902 (2003)] that mutual information (MI) is a useful measure of dependence for electroencephalogram (EEG) data, but we show that the improvement seen in the performance of MI on extracting dependence trends from EEG is more dependent on the type of MI estimator rather than any embedding technique used. In an independent study we conducted in search for an optimal MI estimator, and in particular for EEG applications, we examined the performance of a number of MI estimators on the data set used by Quian Quiroga et al. in their original study, where the performance of different dependence measures on real data was investigated [Phys. Rev. E 65, 041903 (2002)]. We show that for EEG applications the best performance among the investigated estimators is achieved by k-nearest neighbors, which supports the conjecture by Quian Quiroga et al. in Phys. Rev. E 67, 063902 (2003) that the nearest neighbor estimator is the most precise method for estimating MI.