Sampling phylogenetic tree space with the generalized Gibbs sampler

The generalized Gibbs sampler (GGS) is a recently developed Markov chain Monte Carlo (MCMC) technique that enables Gibbs-like sampling of state spaces that lack a convenient representation in terms of a fixed coordinate system. This paper describes a new sampler, called the tree sampler, which uses the GGS to sample from a state space consisting of phylogenetic trees. The tree sampler is useful for a wide range of phylogenetic applications, including Bayesian, maximum likelihood, and maximum parsimony methods. A fast new algorithm to search for a maximum parsimony phylogeny is presented, using the tree sampler in the context of simulated annealing. The mathematics underlying the algorithm is explained and its time complexity is analyzed. The method is tested on two large data sets consisting of 123 sequences and 500 sequences, respectively. The new algorithm is shown to compare very favorably in terms of speed and accuracy to the program DNAPARS from the PHYLIP package.

Segmenting eukaryotic genomes with the generalized Gibbs sampler

Eukaryotic genomes display segmental patterns of variation in various properties, including GC content and degree of evolutionary conservation. DNA segmentation algorithms are aimed at identifying statistically significant boundaries between such segments. Such algorithms may provide a means of discovering new classes of functional elements in eukaryotic genomes. This paper presents a model and an algorithm for Bayesian DNA segmentation and considers the feasibility of using it to segment whole eukaryotic genomes. The algorithm is tested on a range of simulated and real DNA sequences, and the following conclusions are drawn. Firstly, the algorithm correctly identifies non-segmented sequence, and can thus be used to reject the null hypothesis of uniformity in the property of interest. Secondly, estimates of the number and locations of change-points produced by the algorithm are robust to variations in algorithm parameters and initial starting conditions and correspond to real features in the data. Thirdly, the algorithm is successfully used to segment human chromosome 1 according to GC content, thus demonstrating the feasibility of Bayesian segmentation of eukaryotic genomes. The software described in this paper is available from the author's website (www.uq.edu.au/similar to uqjkeith/) or upon request to the author.

A generalized Markov sampler

A recent development of the Markov chain Monte Carlo (MCMC) technique is the emergence of MCMC samplers that allow transitions between different models. Such samplers make possible a range of computational tasks involving models, including model selection, model evaluation, model averaging and hypothesis testing. An example of this type of sampler is the reversible jump MCMC sampler, which is a generalization of the Metropolis-Hastings algorithm. Here, we present a new MCMC sampler of this type. The new sampler is a generalization of the Gibbs sampler, but somewhat surprisingly, it also turns out to encompass as particular cases all of the well-known MCMC samplers, including those of Metropolis, Barker, and Hastings. Moreover, the new sampler generalizes the reversible jump MCMC. It therefore appears to be a very general framework for MCMC sampling. This paper describes the new sampler and illustrates its use in three applications in Computational Biology, specifically determination of consensus sequences, phylogenetic inference and delineation of isochores via multiple change-point analysis.

A Bayesian hierarchical model for categorical longitudinal data from a social survey of immigrants

The paper investigates a Bayesian hierarchical model for the analysis of categorical longitudinal data from a large social survey of immigrants to Australia. Data for each subject are observed on three separate occasions, or waves, of the survey. One of the features of the data set is that observations for some variables are missing for at least one wave. A model for the employment status of immigrants is developed by introducing, at the first stage of a hierarchical model, a multinomial model for the response and then subsequent terms are introduced to explain wave and subject effects. To estimate the model, we use the Gibbs sampler, which allows missing data for both the response and the explanatory variables to be imputed at each iteration of the algorithm, given some appropriate prior distributions. After accounting for significant covariate effects in the model, results show that the relative probability of remaining unemployed diminished with time following arrival in Australia.

A Bayesian approach to imposing curvature on distance functions

The estimated parameters of output distance functions frequently violate the monotonicity, quasi-convexity and convexity constraints implied by economic theory, leading to estimated elasticities and shadow prices that are incorrectly signed, and ultimately to perverse conclusions concerning the effects of input and output changes on productivity growth and relative efficiency levels. We show how a Bayesian approach can be used to impose these constraints on the parameters of a translog output distance function. Implementing the approach involves the use of a Gibbs sampler with data augmentation. A Metropolis-Hastings algorithm is also used within the Gibbs to simulate observations from truncated pdfs. Our methods are developed for the case where panel data is available and technical inefficiency effects are assumed to be time-invariant. Two models-a fixed effects model and a random effects model-are developed and applied to panel data on 17 European railways. We observe significant changes in estimated elasticities and shadow price ratios when regularity restrictions are imposed. (c) 2004 Elsevier B.V. All rights reserved.

Role of generalized and episode specific memories in the word frequency effect in recognition

Frequency, recency, and type of prior exposure to very low-and high-frequency words were manipulated in a 3-phase (i.e., familiarization training, study, and test) design. Increasing the frequency with which a definition for a very low-frequency word was provided during familiarization facilitated the word's recognition in both yes-no (Experiment 1) and forced-choice paradigms (Experiment 2). Recognition of very low-frequency words not accompanied by a definition during familiarization first increased, then decreased as familiarization frequency increased (Experiment I). Reasons for these differences were investigated in Experiment 3 using judgments of recency and frequency. Results suggested that prior familiarization of a very low-frequency word with its definition may allow a more adequate episodic representation of the word to be formed during a subsequent study trial. Theoretical implications of these results for current models of memory are discussed.

Observation of generalized synchronization of chaos in a driven chaotic system

We report on the experimental observation of the generalized synchronization of chaos in a real physical system. We show that under a nonlinear resonant interaction, the chaotic dynamics of a single mode laser can become functionally related to that of a chaotic driving signal and furthermore as the coupling strength is further increased, the chaotic dynamics of the laser approaches that of the driving signal.

Febrile seizures and generalized epilepsy associated with a mutation in the Na+-channel beta 1 subunit gene SCN1B

Febrile seizures affect approximately 3% of all children under six years of age and are by far the most common seizure disorder(1). A small proportion of children with febrile seizures later develop ongoing epilepsy with afebrile seizures(2). Segregation analysis suggests the majority of cases have complex inheritance(3) but rare families show apparent autosomal dominant: inheritance. Two putative loci have been mapped (FEB1 and FEB2), but specific genes have not yet been identified(4,5). We recently described a clinical subset, termed generalized epilepsy with febrile seizures plus (GEFS(+)), in which many family members have seizures with fever that may persist beyond six years of age or be associated with afebrile generalized seizures(6). We now report linkage, in another large GEFS(+) family, to chromosome region 19q13.1 and identification of a mutation in the voltage-gated sodium (Na+)-channel beta 1 subunit gene (SCN1B). The mutation changes a conserved cysteine residue disrupting a putative disulfide bridge which normally maintains an extracellular immunoglobulin-like fold. Go-expression of the mutant pr subunit with a brain Na+-channel alpha subunit in Xenopus laevis oocytes demonstrates that the mutation interferes with the ability of the subunit to modulate channel-gating kinetics consistent with a loss-of-function allele. This observation develops the theme that idiopathic epilepsies are a family of channelopathies and raises the possibility of involvement of other Na+-channel subunit genes in febrile seizures and generalized epilepsies with complex inheritance patterns.

Generalized optimal lattice covering of finite-dimensional Euclidean space

A generalization of the classical problem of optimal lattice covering of R-n is considered. Solutions to this generalized problem are found in two specific classes of lattices. The global optimal solution of the generalization is found for R-2. (C) 1998 Elsevier Science Inc. All rights reserved.

Risk premiums and benefit measures for generalized-expected-utility theories

Standard tools for the analysis of economic problems involving uncertainty, including risk premiums, certainty equivalents and the notions of absolute and relative risk aversion, are developed without making specific assumptions on functional form beyond the basic requirements of monotonicity, transitivity, continuity, and the presumption that individuals prefer certainty to risk. Individuals are not required to display probabilistic sophistication. The approach relies on the distance and benefit functions to characterize preferences relative to a given state-contingent vector of outcomes. The distance and benefit functions are used to derive absolute and relative risk premiums and to characterize preferences exhibiting constant absolute risk aversion (CARA) and constant relative risk aversion (CRRA). A generalization of the notion of Schur-concavity is presented. If preferences are generalized Schur concave, the absolute and relative risk premiums are generalized Schur convex, and the certainty equivalents are generalized Schur concave.

Neuronal sodium-channel alpha 1-subunit mutations in generalized epilepsy with febrile seizures plus

Generalized epilepsy with febrile seizures plus (GEFS+) is a familial epilepsy syndrome characterized by the presence of febrile and afebrile seizures. The first gene, GEFS1, was mapped to chromosome 19q and was identified as the sodium-channel beta1-subunit, SCN1B. A second locus on chromosome 2q, GEFS2, was recently identified as the sodium-channel alpha1-subunit, SCN1A. Single-stranded conformation analysis (SSCA) of SCN1A was performed in 53 unrelated index cases to estimate the frequency of mutations in patients with GEFS+. No mutations were found in 17 isolated cases of GEFS+. Three novel SCN1A mutations-D188V, V1353L, and I1656M-were found in 36 familial cases; of the remaining 33 families, 3 had mutations in SCN1B. On the basis of SSCA, the combined frequency of SCN1A and SCN1B mutations in familial cases of GEFS+ was found to be 17%.

Truncation of the GABA(A)-receptor gamma 2 subunit in a family with generalized epilepsy with febrile seizures plus

Recent findings from studies of two families have shown that mutations in the GABA(A)-receptor gamma2 subunit are associated with generalized epilepsies and febrile seizures. Here we describe a family that has generalized epilepsy with febrile seizures plus (GEFS(+)), including an individual with severe myoclonic epilepsy of infancy, in whom a third GABA(A)-receptor gamma2-subunit mutation was found. This mutation lies in the intracellular loop between the third and fourth transmembrane domains of the GABA(A)-receptor gamma2 subunit and introduces a premature stop codon at Q351 in the mature protein. GABA sensitivity in Xenopus laevis oocytes expressing the mutant gamma2(Q351X) subunit is completely abolished, and fluorescent-microscopy studies have shown that receptors containing GFP-labeled gamma2(Q351X) protein are retained in the lumen of the endoplasmic reticulum. This finding reinforces the involvement of GABA(A) receptors in epilepsy.

Generalized epilepsy with febrile seizures plus: Mutation of the sodium channel subunit SCN1B

Generalized epilepsy with febrile seizures plus (GEFS(+)) is an important childhood genetic epilepsy syndrome with heterogeneous phenotypes, including febrile seizures (FS) and generalized epilepsies of variable severity. Forty unrelated GEFS(+) and FS patients were screened for mutations in the sodium channel beta-subunits SCN1B and SCN2B, and the second GEFS(+) family with an SCN1B mutation is described here. The family had 19 affected individuals: 16 with typical GEFS(+) phenotypes and three with other epilepsy phenotypes. Site-specific mutation within SCN1B remains a rare cause of GEFS(+), and the authors found no evidence to implicate SCN2B in this syndrome.

Production under uncertainty and choice under uncertainty in the emergence of generalized expected utility theory

This paper presents a personal view of the interaction between the analysis of choice under uncertainty and the analysis of production under uncertainty. Interest in the foundations of the theory of choice under uncertainty was stimulated by applications of expected utility theory such as the Sandmo model of production under uncertainty. This interest led to the development of generalized models including rank-dependent expected utility theory. In turn, the development of generalized expected utility models raised the question of whether such models could be used in the analysis of applied problems such as those involving production under uncertainty. Finally, the revival of the state-contingent approach led to the recognition of a fundamental duality between choice problems and production problems.

Background risk in generalized expected utility theory

In this paper, it is shown that, for a wide range of risk-averse generalized expected utility preferences, independent risks are complementary, contrary to the results for expected utility preferences satisfying conditions such as proper and standard risk aversion.