A likelihood ratio appropach to family-based association studies with covariates

We introduce a procedure for association based analysis of nuclear families that allows for dichotomous and more general measurements of phenotype and inclusion of covariate information. Standard generalized linear models are used to relate phenotype and its predictors. Our test procedure, based on the likelihood ratio, unifies the estimation of all parameters through the likelihood itself and yields maximum likelihood estimates of the genetic relative risk and interaction parameters. Our method has advantages in modelling the covariate and gene-covariate interaction terms over recently proposed conditional score tests that include covariate information via a two-stage modelling approach. We apply our method in a study of human systemic lupus erythematosus and the C-reactive protein that includes sex as a covariate.

Host niches and defensive extended phenotypes structure parasitoid wasp communities

Oak galls are spectacular extended phenotypes of gallwasp genes in host oak tissues and have evolved complex morphologies that serve, in part, to exclude parasitoid natural enemies. Parasitoids and their insect herbivore hosts have coevolved to produce diverse communities comprising about a third of all animal species. The factors structuring these communities, however, remain poorly understood. An emerging theme in community ecology is the need to consider the effects of host traits, shaped by both natural selection and phylogenetic history, on associated communities of natural enemies. Here we examine the impact of host traits and phylogenetic relatedness on 48 ecologically closed and species-rich communities of parasitoids attacking gall-inducing wasps on oaks. Gallwasps induce the development of spectacular and structurally complex galls whose species- and generation-specific morphologies are the extended phenotypes of gallwasp genes. All the associated natural enemies attack their concealed hosts through gall tissues, and several structural gall traits have been shown to enhance defence against parasitoid attack. Here we explore the significance of these and other host traits in predicting variation in parasitoid community structure across gallwasp species. In particular, we test the "Enemy Hypothesis,'' which predicts that galls with similar morphology will exclude similar sets of parasitoids and therefore have similar parasitoid communities. Having controlled for phylogenetic patterning in host traits and communities, we found significant correlations between parasitoid community structure and several gall structural traits (toughness, hairiness, stickiness), supporting the Enemy Hypothesis. Parasitoid community structure was also consistently predicted by components of the hosts' spatiotemporal niche, particularly host oak taxonomy and gall location (e.g., leaf versus bud versus seed). The combined explanatory power of structural and spatiotemporal traits on community structure can be high, reaching 62% in one analysis. The observed patterns derive mainly from partial niche specialisation of highly generalist parasitoids with broad host ranges (>20 hosts), rather than strict separation of enemies with narrower host ranges, and so may contribute to maintenance of the richness of generalist parasitoids in gallwasp communities. Though evolutionary escape from parasitoids might most effectively be achieved via changes in host oak taxon, extreme conservatism in this trait for gallwasps suggests that selection is more likely to have acted on gall morphology and location. Any escape from parasitoids associated with evolutionary shifts in these traits has probably only been transient, however, due to subsequent recruitment of parasitoid species already attacking other host galls with similar trait combinations.

A novel parallel approach to the likelihood-based estimation of admixture in population genetics

Inferring population admixture from genetic data and quantifying it is a difficult but crucial task in evolutionary and conservation biology. Unfortunately state-of-the-art probabilistic approaches are computationally demanding. Effectively exploiting the computational power of modern multiprocessor systems can thus have a positive impact to Monte Carlo-based simulation of admixture modeling. A novel parallel approach is briefly described and promising results on its message passing interface (MPI)-based C implementation are reported.

A comparison of the active surveillance of scrapie in the European Union

The abattoir and the fallen stock surveys constitute the active surveillance component aimed at improving the detection of scrapie across the European Union. Previous studies have suggested the occurrence of significant differences in the operation of the surveys across the EU. In the present study we assessed the standardisation of the surveys throughout time across the EU and identified clusters of countries with similar underlying characteristics allowing comparisons between them. In the absence of sufficient covariate information to explain the observed variability across countries, we modelled the unobserved heterogeneity by means of non-parametric distributions on the risk ratios of the fallen stock over the abattoir survey. More specifically, we used the profile likelihood method on 2003, 2004 and 2005 active surveillance data for 18 European countries on classical scrapie, and on 2004 and 2005 data for atypical scrapie separately. We extended our analyses to include the limited covariate information available, more specifically, the proportion of the adult sheep population sampled by the fallen stock survey every year. Our results show that the between-country heterogeneity dropped in 2004 and 2005 relative to that of 2003 for classical scrapie. As a consequence, the number of clusters in the last two years was also reduced indicating the gradual standardisation of the surveillance efforts across the EU. The crude analyses of the atypical data grouped all the countries in one cluster and showed non-significant gain in the detection of this type of scrapie by any of the two sources. The proportion of the population sampled by the fallen stock appeared significantly associated with our risk ratio for both types of scrapie, although in opposite directions: negative for classical and positive for atypical. The initial justification for the fallen stock, targeting a high-risk population to increase the likelihood of case finding, appears compromised for both types of scrapie in some countries.

A note on the accuracy of PAC-likelihood inference with microsatellite data

Stephens and Donnelly have introduced a simple yet powerful importance sampling scheme for computing the likelihood in population genetic models. Fundamental to the method is an approximation to the conditional probability of the allelic type of an additional gene, given those currently in the sample. As noted by Li and Stephens, the product of these conditional probabilities for a sequence of draws that gives the frequency of allelic types in a sample is an approximation to the likelihood, and can be used directly in inference. The aim of this note is to demonstrate the high level of accuracy of "product of approximate conditionals" (PAC) likelihood when used with microsatellite data. Results obtained on simulated microsatellite data show that this strategy leads to a negligible bias over a wide range of the scaled mutation parameter theta. Furthermore, the sampling variance of likelihood estimates as well as the computation time are lower than that obtained with importance sampling on the whole range of theta. It follows that this approach represents an efficient substitute to IS algorithms in computer intensive (e.g. MCMC) inference methods in population genetics. (c) 2006 Elsevier Inc. All rights reserved.

Likelihood-based estimation of microsatellite mutation rates

Microsatellites are widely used in genetic analyses, many of which require reliable estimates of microsatellite mutation rates, yet the factors determining mutation rates are uncertain. The most straightforward and conclusive method by which to study mutation is direct observation of allele transmissions in parent-child pairs, and studies of this type suggest a positive, possibly exponential, relationship between mutation rate and allele size, together with a bias toward length increase. Except for microsatellites on the Y chromosome, however, previous analyses have not made full use of available data and may have introduced bias: mutations have been identified only where child genotypes could not be generated by transmission from parents' genotypes, so that the probability that a mutation is detected depends on the distribution of allele lengths and varies with allele length. We introduce a likelihood-based approach that has two key advantages over existing methods. First, we can make formal comparisons between competing models of microsatellite evolution; second, we obtain asymptotically unbiased and efficient parameter estimates. Application to data composed of 118,866 parent-offspring transmissions of AC microsatellites supports the hypothesis that mutation rate increases exponentially with microsatellite length, with a suggestion that contractions become more likely than expansions as length increases. This would lead to a stationary distribution for allele length maintained by mutational balance. There is no evidence that contractions and expansions differ in their step size distributions.

Extended-chain, multinuclear transition metal complexes bridged by cyanodiazenido(1-), [N=N-C=N](-), ligands

Extended-chain complexes containing multiple transition metal centres linked by conjugated mu-cyanodiazenido(1-) ligands [N= N-C N]-have been obtained by reaction of trans-[BrW(dppe)(2)(N2CN)], 1, [dppe = 1,2-bis(diphenylphosphino) ethane] with dirhodium(II) tetra-acetate, bis(benzonitrile) palladium(II) dichloride, and bis(aqua) M(II) bis(hexa. uoroacetylacetonate) (M = Mn, Ni, Cu, Zn): stronger Lewis acids such as tetrakis(acetonitrile) palladium(II) tetra. uoroborate and boron trifl. uoride promote hydrolysis of complex 1, leading to the isolation of a novel carbamoylhydrazido(2-) complex, trans-[BrW(dppe) 2(N2HC=ONH2)](+)[BF4](-).

Implicit memory: a prime example for brand consideration and choice

Three experiments investigated the influence of implicit memory for familiar brand names on consumer choice. Priming was measured using modified preference judgment tasks that comprised both brand consideration and choice components. Experiment 1 used a 'complex choice task' where consideration and choice stages were characterized as acting in sequence. Experiment 2 explored a different formulation whereby consideration and choice were assumed to act in parallel, Both experiments demonstrated that priming had an influence on brand consideration but not on final or preferred choice. Finally, Experiment 3 replicated and extended these findings under more realistic conditions where participants actually received some of the products that they selected. Overall, the experiments suggested that for many decisions involving the consideration of familiar brands prior to choice, previous exposure to brand names can increase the likelihood that they will enter the consumers' consideration set. However, the advantage does not appear to extend to choice itself. Copyright (C) 2004 John Wiley Sons, Ltd.

Cooperative transport of extended payloads

In this paper we present experimental results for the dual robot transport of an extended payload. Two robotic rovers that were designed specifically for the extended payload transport task are described. Each rover incorporates a 4-DOF robot arm incorporating three active joints (one of which is a gripper), a passive wrist, and a mobile base which employs a rocker-bogie design. A set of behaviours has been developed to support the performance of the task, integrating simple sensing with controls. We describe the behaviours and their integration within the overall task structure. The experimental results presented focus on the manipulation elements of the task, but incorporate a complete cycle of pick-up, traversal, and putdown.

A fast identification algorithm for Box-Cox transformation based radial basis function neural network

In this letter, a Box-Cox transformation-based radial basis function (RBF) neural network is introduced using the RBF neural network to represent the transformed system output. Initially a fixed and moderate sized RBF model base is derived based on a rank revealing orthogonal matrix triangularization (QR decomposition). Then a new fast identification algorithm is introduced using Gauss-Newton algorithm to derive the required Box-Cox transformation, based on a maximum likelihood estimator. The main contribution of this letter is to explore the special structure of the proposed RBF neural network for computational efficiency by utilizing the inverse of matrix block decomposition lemma. Finally, the Box-Cox transformation-based RBF neural network, with good generalization and sparsity, is identified based on the derived optimal Box-Cox transformation and a D-optimality-based orthogonal forward regression algorithm. The proposed algorithm and its efficacy are demonstrated with an illustrative example in comparison with support vector machine regression.

A neurofuzzy network knowledge extraction and extended gram-Schmidt algorithm for model subspace decomposition

This paper introduces a new neurofuzzy model construction and parameter estimation algorithm from observed finite data sets, based on a Takagi and Sugeno (T-S) inference mechanism and a new extended Gram-Schmidt orthogonal decomposition algorithm, for the modeling of a priori unknown dynamical systems in the form of a set of fuzzy rules. The first contribution of the paper is the introduction of a one to one mapping between a fuzzy rule-base and a model matrix feature subspace using the T-S inference mechanism. This link enables the numerical properties associated with a rule-based matrix subspace, the relationships amongst these matrix subspaces, and the correlation between the output vector and a rule-base matrix subspace, to be investigated and extracted as rule-based knowledge to enhance model transparency. The matrix subspace spanned by a fuzzy rule is initially derived as the input regression matrix multiplied by a weighting matrix that consists of the corresponding fuzzy membership functions over the training data set. Model transparency is explored by the derivation of an equivalence between an A-optimality experimental design criterion of the weighting matrix and the average model output sensitivity to the fuzzy rule, so that rule-bases can be effectively measured by their identifiability via the A-optimality experimental design criterion. The A-optimality experimental design criterion of the weighting matrices of fuzzy rules is used to construct an initial model rule-base. An extended Gram-Schmidt algorithm is then developed to estimate the parameter vector for each rule. This new algorithm decomposes the model rule-bases via an orthogonal subspace decomposition approach, so as to enhance model transparency with the capability of interpreting the derived rule-base energy level. This new approach is computationally simpler than the conventional Gram-Schmidt algorithm for resolving high dimensional regression problems, whereby it is computationally desirable to decompose complex models into a few submodels rather than a single model with large number of input variables and the associated curse of dimensionality problem. Numerical examples are included to demonstrate the effectiveness of the proposed new algorithm.