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.

The population pharmacokinetics of citalopram after deliberate self-poisoning: A Bayesian approach

Defining the pharmacokinetics of drugs in overdose is complicated. Deliberate self-poisoning is generally impulsive and associated with poor accuracy in dose history. In addition, early blood samples are rarely collected to characterize the whole plasma-concentration time profile and the effect of decontamination on the pharmacokinetics is uncertain. The aim of this study was to explore a fully Bayesian methodology for population pharmacokinetic analysis of data that arose from deliberate self-poisoning with citalopram. Prior information on the pharmacokinetic parameters was elicited from 14 published studies on citalopram when taken in therapeutic doses. The data set included concentration-time data from 53 patients studied after 63 citalopram overdose events (dose range: 20-1700 mg). Activated charcoal was administered between 0.5 and 4 h after 17 overdose events. The clinical investigator graded the veracity of the patients' dosing history on a 5-point ordinal scale. Inclusion of informative priors stabilised the pharmacokinetic model and the population mean values could be estimated well. There were no indications of non-linear clearance after excessive doses. The final model included an estimated uncertainty of the dose amount which in a simulation study was shown to not affect the model's ability to characterise the effects of activated charcoal. The effect of activated charcoal on clearance and bioavailability was pronounced and resulted in a 72% increase and 22% decrease, respectively. These findings suggest charcoal administration is potentially beneficial after citalopram overdose. The methodology explored seems promising for exploring the dose-exposure relationship in the toxicological settings.

Endoparasite communities of five fish species (Labridae : Cheilininae) from Lizard Island: how important is the ecology and phylogeny of the hosts?

The parasite community of animals is generally influenced by host physiology, ecology, and phylogeny. Therefore, sympatric and phylogenetically related hosts with similar ecologies should have similar parasite communities. To test this hypothesis we surveyed the endoparasites of 5 closely related cheilinine fishes (Labridae) from the Great Barrier Reef. They were Cheilinus chlorounts, C. trilobatus, C. fasciatils, Epibulus insidiator and OxYcheilinus diagrainnia. VVe examined the relationship between parasitological variables (richness, abundance and diversity) and host characteristics (bodv weight, diet and phuylogeny). The 5 fishes had 31 parasite species with 9-18 parasite species per fish species. Cestode larvae (mostly Tetraphyllidea) were the most abundant and prevalent parasites followed by nematodes and digeneans. Parasites, body size and diet of hosts differed between fish species. In general, body weight, diet and host phylogeny each explained some of the variation in richness and composition of parasites among the fishes. The 2 most closely related species, Cheilinus chlorourus and C. trilobatus, had broadly similar parasites but the Other fish species differed significantly in all variables. However, there was no all -encompassing pattern. This may, be because different lineages of parasites may react differently to ecological variables. We also argue that adult parasites may respond principally to host diet. In contrast, larval parasite composition may respond both to host diet and predator-prey interactions because this is the path by which many, parasites complete their life-cycles. Finally, variation in parasite phylogeny and parasite life-cycles among hosts likely increase the complexity of the system making it difficult to find all-encompassing patterns between host characteristics and parasites, particularly when all the species in rich parasite communities are considered.

Phylogeny, evolution and biogeography of the Quadrifoliovariinae Yamaguti, 1965 (Digenea : Lecithasteridae)

The Quadrifoliovariinae is revised and three new species of Quadrifoliovarium Yamaguit, 1965 from acanthurid fishes of the genus Naso from waters of the Indo-Pacific are described: Q, maceria n. sp. from N. tonganus, N. annulatus, N. fageni and N. brevirostris; Q. simplex n. sp. from N. tonganus and N. quannulatus; and Q. quattuordecim n. sp. from N. tonganus. Amendments are made to the characterisation of the Quadrifoliovariinae, Quadrifoliovarium, Bilacinia Manter, 1969 and Unilacinia Manter, 1969 in light of observations on type and new material. A molecular phylogeny based on ITS2 and 28S regions of the ribosomal DNA is proposed. The phylogeny suggests that U. asymmetrica is the most basal taxon and Q. simplex n. sp. and Q. quattuordecim n. sp. the most derived. Evolution of morphological traits within the Quadrifoliovariinae are discussed in light of the molecular phylogeny. Molecular sequences of the ITS2 rDNA were identical between specimens of Q. pritchardae collected off Exmouth (Indian Ocean), Heron Island and Lizard Island (Western Pacific) and Moorea (far Eastern Indo-Pacific), indicating a broad Indo-Pacific distribution. All members of the subfamily are recorded only from the acanthurid genus Naso, with the exception of B. lobatum (Yamaguti, 1970), which has been recorded from a pomacanthid. The restricted host range of the group is discussed in the light of the phylogeny of the host genus Naso.

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.

Motor unit number estimation - A Bayesian approach

All muscle contractions are dependent on the functioning of motor units. In diseases such as amyotrophic lateral sclerosis (ALS), progressive loss of motor units leads to gradual paralysis. A major difficulty in the search for a treatment for these diseases has been the lack of a reliable measure of disease progression. One possible measure would be an estimate of the number of surviving motor units. Despite over 30 years of motor unit number estimation (MUNE), all proposed methods have been met with practical and theoretical objections. Our aim is to develop a method of MUNE that overcomes these objections. We record the compound muscle action potential (CMAP) from a selected muscle in response to a graded electrical stimulation applied to the nerve. As the stimulus increases, the threshold of each motor unit is exceeded, and the size of the CMAP increases until a maximum response is obtained. However, the threshold potential required to excite an axon is not a precise value but fluctuates over a small range leading to probabilistic activation of motor units in response to a given stimulus. When the threshold ranges of motor units overlap, there may be alternation where the number of motor units that fire in response to the stimulus is variable. This means that increments in the value of the CMAP correspond to the firing of different combinations of motor units. At a fixed stimulus, variability in the CMAP, measured as variance, can be used to conduct MUNE using the "statistical" or the "Poisson" method. However, this method relies on the assumptions that the numbers of motor units that are firing probabilistically have the Poisson distribution and that all single motor unit action potentials (MUAP) have a fixed and identical size. These assumptions are not necessarily correct. We propose to develop a Bayesian statistical methodology to analyze electrophysiological data to provide an estimate of motor unit numbers. Our method of MUNE incorporates the variability of the threshold, the variability between and within single MUAPs, and baseline variability. Our model not only gives the most probable number of motor units but also provides information about both the population of units and individual units. We use Markov chain Monte Carlo to obtain information about the characteristics of individual motor units and about the population of motor units and the Bayesian information criterion for MUNE. We test our method of MUNE on three subjects. Our method provides a reproducible estimate for a patient with stable but severe ALS. In a serial study, we demonstrate a decline in the number of motor unit numbers with a patient with rapidly advancing disease. Finally, with our last patient, we show that our method has the capacity to estimate a larger number of motor units.

A review of placoderm scales, and their significance in placoderm phylogeny

An historical review of the literature relating to placoderm scales preserved in association with articulated dermal plates, or as isolated units in microvertebrate assemblages, is followed by a discussion of their relevance in phylogenetic analyses of the Placodermi. The dentinous tissue forming the tubercles of Early Devonian acanthothoracid scales and dermal bone is similar to that of the dermal bone ornament of some osteostracans, and denticles of the vertebrate Skiichthys from the Ordovician Harding Sandstone. This similarity supports the proposition that the gnathostomes are the sister-group of the Osteostraci, with the Placodermi branching earliest within the gnathostomes, and the Acanthothoraci branching earliest within the Placodermi. The meso-semidentine in acanthothoracid tubercles, rather than semidentine (sensu stricto), is most likely to be synapomorphic for the Placodermi.

A bayesian artificial neural network method to characterise laminar defects using dynamic measurements

This paper reports on the development of an artificial neural network (ANN) method to detect laminar defects following the pattern matching approach utilizing dynamic measurement. Although structural health monitoring (SHM) using ANN has attracted much attention in the last decade, the problem of how to select the optimal class of ANN models has not been investigated in great depth. It turns out that the lack of a rigorous ANN design methodology is one of the main reasons for the delay in the successful application of the promising technique in SHM. In this paper, a Bayesian method is applied in the selection of the optimal class of ANN models for a given set of input/target training data. The ANN design method is demonstrated for the case of the detection and characterisation of laminar defects in carbon fibre-reinforced beams using flexural vibration data for beams with and without non-symmetric delamination damage.

Bayesian population pharmacokinetic analysis of sirolimus

Objective: It is usual that data collected from routine clinical care is sparse and unable to support the more complex pharmacokinetic (PK) models that may have been reported in previous rich data studies. Informative priors may be a pre-requisite for model development. The aim of this study was to estimate the population PK parameters of sirolimus using a fully Bayesian approach with informative priors. Methods: Informative priors including prior mean and precision of the prior mean were elicited from previous published studies using a meta-analytic technique. Precision of between-subject variability was determined by simulations from a Wishart distribution using MATLAB (version 6.5). Concentration-time data of sirolimus retrospectively collected from kidney transplant patients were analysed using WinBUGS (version 1.3). The candidate models were either one- or two-compartment with first order absorption and first order elimination. Model discrimination was based on computation of the posterior odds supporting the model. Results: A total of 315 concentration-time points were obtained from 25 patients. Most data were clustered at trough concentrations with range of 1.6 to 77 hours post-dose. Using informative priors, either a one- or two-compartment model could be used to describe the data. When a one-compartment model was applied, information was gained from the data for the value of apparent clearance (CL/F = 18.5 L/h), and apparent volume of distribution (V/F = 1406 L) but no information was gained about the absorption rate constant (ka). When a two-compartment model was fitted to the data, the data were informative about CL/F, apparent inter-compartmental clearance, and apparent volume of distribution of the peripheral compartment (13.2 L/h, 20.8 L/h, and 579 L, respectively). The posterior distribution of the volume distribution of central compartment and ka were the same as priors. The posterior odds for the two-compartment model was 8.1, indicating the data supported the two-compartment model. Conclusion: The use of informative priors supported the choice of a more complex and informative model that would otherwise have not been supported by the sparse data.