An Efficient Regression Type Estimator in Survey Sampling

For the estimation of population mean in simple random sampling, an efficient regression-type estimator is proposed which is more efficient than the conventional regression estimator and hence than mean per unit estimator, ratio and product estimators and many other estimators proposed by various authors. Some numerical examples are included for illustration.     

A physiological analogy of the niche for projecting the potential distribution of plants

Aim  To develop a physiologically based model of the plant niche for use in species distribution modelling. Location  Europe. Methods  We link the Thornley transport resistance (TTR) model with functions which describe how the TTR’s model parameters are influenced by abiotic environmental factors. The TTR model considers how carbon and nutrient uptake, and the allocation of these assimilates, influence growth. We use indirect statistical methods to estimate the model parameters from a high resolution data set on tree distribution for 22 European tree species. Results  We infer, from distribution data and abiotic forcing data, the physiological niche dimensions of 22 European tree species. We found that the model fits were reasonable (AUC: 0.79–0.964). The projected distributions were characterized by a false positive rate of 0.19 and a false negative rate 0.12. The fitted models are used to generate projections of the environmental factors that limit the range boundaries of the study species. Main conclusions  We show that physiological models can be used to derive physiological niche dimensions from species distribution data. Future work should focus on including prior information on physiological rates into the parameter estimation process. Application of the TTR model to species distribution modelling suggests new avenues for establishing explicit links between distribution and physiology, and for generating hypotheses about how ecophysiological processes influence the distribution of plants.     

On computational efficiency of the hybrid Monte Carlo method applied to the multicanonical ensemble

Abstract

In the present paper, computational efficiency of the hybrid Monte Carlo (HMC) method applied to the multicanonical ensemble is studied; the HMC is an equation of motion guided Monte Carlo method. As in the standard HMC for the canonical ensemble, the multicanonical HMC calculations with high acceptance ratio show better efficiency; about 60% acceptance yields the best performance for the system examined.     

Correction of a bootstrap approach to testing for evolution along lines of least resistance

Testing for an association between the leading vectors of multivariate trait (co)variation within populations (the ‘line of least resistance’) and among populations is an important tool for exploring variational bias in evolution. In a recent study of stickleback fish populations, a bootstrap‐based test was introduced that takes into account estimation error in both vectors and hence improves the previously available bootstrap method. Because this test was implemented incorrectly, however, I here describe the correct test protocol and provide a reanalysis of the original data set. The application of this new test protocol should improve future investigations of evolution along lines of least resistance and other vector comparisons.     

Conditional Maximum Likelihood Estimation Following a Group Sequential Test

We consider estimation after a group sequential test. An estimator that is unbiased or has small bias may have substantial conditional bias (Troendle and Yu, 1999, Coburger and Wassmer, 2001). In this paper we derive the conditional maximum likelihood estimators of both the primary parameter and a secondary parameter, and investigate their properties within a conditional inference framework. The method applies to both the usual and adaptive group sequential test designs. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Detection of native-like models for amino acid sequences of unknown three-dimensional structure in a data base of known protein conformations.

We present an approach which can be used to identify native-like folds in a data base of protein conformations in the absence of any sequence homology to proteins in the data base. The method is based on a knowledge-based force field derived from a set of known protein conformations. A given sequence is mounted on all conformations in the data base and the associated energies are calculated. Using several conformations and sequences from the globin family we show that the native conformation is identified correctly. In fact the resolution of the force field is high enough to discriminate between a native fold and several closely related conformations. We then apply the procedure to several globins of known sequence but unknown three dimensional structure. The homology of these sequences to globins of known structures in the data base ranges from 49 to 17%. With one exception we find that for all globin sequences one of the known globin folds is identified as the most favorable conformation. These results are obtained using a force field derived from a data base devoid of globins of known structure. We briefly discuss useful applications in protein structural research and future development of our approach.     

Interaction of silent and replacement changes in eukaryotic coding sequences

Summary We examined the codon usages in wellconserved and less-well-conserved regions of vertebrate protein genes and found them to be similar. Despite this similarity, there is a statistically significant decrease in codon bias in the less-well-conserved regions. Our analysis suggests that although those codon changes initially fixed under amino acid replacements tend to follow the overall codon usage pattern, they also reduce the bias in codon usage. This decrease in codon bias leads one to predict that the rate of change of synonymous codons should be greater in those regions that are less well conserved at the amino acid level than in the better-conserved regions. Our analysis supports this prediction. Furthermore, we demonstrate a significantly elevated rate of change of synonymous codons among the adjacent codons 5 to amino acid replacement positions. This provides further support for the idea that there are contextual constraints on the choice of synonymous codons in eukaryotes.     

Sampling bias in estimating Design II variance components with S1 families

Summary The use of several S₁ individuals to represent an S₀ individual permits the use of a Design II mating scheme for plants with only one pistillate flower per plant. Estimates of additive (V _A ) and dominance (V _D ) variance from this mating scheme will be biased upwards, when a small number (10) of individuals of each S₁ line are used. This bias can be computed, and the additive and dominance estimates can be corrected. Of particular interest is the observation that the additive genetic variance contributes to bias in estimates of V _D . When S₀ plants are non inbred and their selfedprogeny (S₁ lines) are used to represent them in developing families for use in the Design II, where m₁ is the number of individuals used to represent an S₁ line in developing half sib-families and m₂ is the number of individuals used to represent the S₁ line in making up full sib-families. For example, in a 3×3 Design II, with about 10 individuals used to represent each S₁ line in each cross, m₂ = 10 and m₁ = 30. When m₁ = m₂ = 1, and Joint contribution from Department of Agronomy, University of Nebraska 68583, and the S. S. Cameron Laboratory, Werribee, Victoria 3030, Australia. Published as paper No. 7395, Journal Series     

Lipid vertical motion and related steric effects in bilayer membranes

We present a theoretical model of a lipid bilayer in its gel state which explicity couples the vertical displacements of the lipid chains to their conformational state. In this model the chains are free to move longitudinally under a potential due to the neighbouring chains. The potential is due to a restoring force with a force constant k and thus acts to keep them in the local plane as defined by their nearest neighbours. It is demonstrated that the force constant k is directly related to the internal bilayer pressure, , and that if a value =33 dynes/cm is assumed then k=17.3 dynes/cm. Steric effects are explicity included by allowing chains to twist into free volume created by the vertical displacement of neighbouring chains. The Hamiltonian is expressed in terms of the projection operators, P _ij , describing the displacement of chain i relative to a neighbour j, and G _ij describing the direction of a twist in chain i. The model is solved both analytically and via Monte Carlo simulations for a one-dimensional system. The possibility of phase-transitions in two-dimensions and the relevance to the bilayer pre-transition is discussed.This work was first presented in poster form at the Canadian Biochemical Society Conference held in Banff, Alberta, Canada, April 29–May 4, 1984Work supported in part by the Natural Sciences and Engineering Research Council of Canada, Le Fonds Formation des Chercheurs et Action a la Recherche du Quebec, and the Advisory Research Council of Queen's University