Modeling continuous self-report measures of perceived emotion using generalized additive mixed models

Emotion research has long been dominated by the “standard method” of displaying posed or acted static images of facial expressions of emotion. While this method has been useful it is unable to investigate the dynamic nature of emotion expression. Although continuous self-report traces have enabled the measurement of dynamic expressions of emotion, a consensus has not been reached on the correct statistical techniques that permit inferences to be made with such measures. We propose Generalized Additive Models and Generalized Additive Mixed Models as techniques that can account for the dynamic nature of such continuous measures. These models allow us to hold constant shared components of responses that are due to perceived emotion across time, while enabling inference concerning linear differences between groups. The mixed model GAMM approach is preferred as it can account for autocorrelation in time series data and allows emotion decoding participants to be modelled as random effects. To increase confidence in linear differences we assess the methods that address interactions between categorical variables and dynamic changes over time. In addition we provide comments on the use of Generalized Additive Models to assess the effect size of shared perceived emotion and discuss sample sizes. Finally we address additional uses, the inference of feature detection, continuous variable interactions, and measurement of ambiguity.

Balancing the (carbon) budget: Using linear inverse models to estimate carbon flows and mass-balance 13C:15N labelling experiments in low oxygen sediments.

Over 1 million km² of seafloor experience permanent low-oxygen conditions within oxygen minimum zones (OMZs). OMZs are predicted to grow as a consequence of climate change, potentially affecting oceanic biogeochemical cycles. The Arabian Sea OMZ impinges upon the western Indian continental margin at bathyal depths (150 - 1500 m) producing a strong depth dependent oxygen gradient at the sea floor. The influence of the OMZ upon the short term processing of organic matter by sediment ecosystems was investigated using in situ stable isotope pulse chase experiments. These deployed doses of ¹³C:¹⁵N labeled organic matter onto the sediment surface at four stations from across the OMZ (water depth 540 - 1100 m; [O₂] = 0.35 - 15 μM). In order to prevent experimentally anoxia, the mesocosms were not sealed. ¹³C and ¹⁵N labels were traced into sediment, bacteria, fauna and ¹³C into sediment porewater DIC and DOC. However, the DIC and DOC flux to the water column could not be measured, limiting our capacity to obtain mass-balance for C in each experimental mesocosm. Linear Inverse Modeling (LIM) provides a method to obtain a mass-balanced model of carbon flow that integrates stable-isotope tracer data with community biomass and biogeochemical flux data from a range of sources. Here we present an adaptation of the LIM methodology used to investigate how ecosystem structure influenced carbon flow across the Indian margin OMZ. We demonstrate how oxygen conditions affect food-web complexity, affecting the linkages between the bacteria, foraminifera and metazoan fauna, and their contributions to benthic respiration. The food-web models demonstrate how changes in ecosystem complexity are associated with oxygen availability across the OMZ and allow us to obtain a complete carbon budget for the stationa where stable-isotope labelling experiments were conducted.

Applying species distribution modelling to identify areas of high conservation value for endangered species: A case study using Margaritifera margaritifera (L.)

The development and implementation of a population supplementation and restoration plan for any endangered species should involve an understanding of the species’ habitat requirements prior to the release of any captive bred individuals. The freshwater pearl mussel, Margaritifera margaritifera, has undergone dramatic declines over the last century and is now globally endangered. In Northern Ireland, the release of captive bred individuals is being used to support wild populations and repatriate the species in areas where it once existed. We employed a combination of maximum entropy modelling (MAXENT) and Generalized Linear Mixed Models (GLMM) to identify ecological parameters necessary to support wild populations using GIS-based landscape scale and ground-truthed habitat scale environmental parameters. The GIS-based landscape scale model suggested that mussel occurrence was associated with altitude and soil characteristics including the carbon, clay, sand, and silt content. Notably, mussels were associated with a relatively narrow band of variance indicating that M. margaritifera has a highly specific landscape niche. The ground-truthed habitat scale model suggested that mussel occurrence was associated with stable consolidated substrates, the extent of bankside trees, presence of indicative macrophyte species and fast flowing water. We propose a three phase conservation strategy for M. margaritifera identifying suitable areas within rivers that (i) have a high conservation value yet needing habitat restoration at a local level, (ii) sites for population supplementation of existing populations and (iii) sites for species reintroduction to rivers where the mussel historically occurred but is now locally extinct. A combined analytical approach including GIS-based landscape scale and ground-truthed habitat scale models provides a robust method by which suitable release sites can be identified for the population supplementation and restoration of an endangered species. Our results will be highly influential in the future management of M. margaritifera in Northern Ireland.

Genetic variation in the interleukin-1 beta-converting enzyme associates with cognitive function. The PROSPER study

Inflammation is thought to play an important role in the development of cognitive decline and dementia in old age. The interleukin-1 signalling pathway may play a prominent role in this process. The gene encoding for interleukin-1 beta-converting enzyme (ICE) is likely to influence IL-1 beta levels. Inhibition of ICE decreases the age-related increase in IL-1 beta levels and may therefore improve memory function. We assessed whether genetic variation in the ICE gene associates with cognitive function in an elderly population. All 5804 participants of the PROspective Study of Pravastatin in the Elderly at Risk (PROSPER) were genotyped for the 10643GC, 9323GA, 8996AG and 5352GA polymorphisms in the ICE gene. Cross-sectional associations between the polymorphisms and cognitive function were assessed with linear regression. Longitudinal associations between polymorphisms, haplotypes and cognitive function were assessed with linear mixed models. All associations were adjusted for sex, age, education, country, treatment with pravastatin and version of test where appropriate. Subjects carrying the variants 10643C and 5352A allele had significantly lower IL-1 beta production levels (P

Identifying appropriate sampling and modelling approaches for analysing distributional patterns of Antarctic terrestrial arthropods along the Victoria Land latitudinal gradient

Biotic communities in Antarctic terrestrial ecosystems are relatively simple and often lack higher trophic levels (e. g. predators); thus, it is often assumed that species' distributions are mainly affected by abiotic factors such as climatic conditions, which change with increasing latitude, altitude and/or distance from the coast. However, it is becoming increasingly apparent that factors other than geographical gradients affect the distribution of organisms with low dispersal capability such as the terrestrial arthropods. In Victoria Land (East Antarctica) the distribution of springtail (Collembola) and mite (Acari) species vary at scales that range from a few square centimetres to regional and continental. Different species show different scales of variation that relate to factors such as local geological and glaciological history, and biotic interactions, but only weakly with latitudinal/altitudinal gradients. Here, we review the relevant literature and outline more appropriate sampling designs as well as suitable modelling techniques (e. g. linear mixed models and eigenvector mapping), that will more adequately address and identify the range of factors responsible for the distribution of terrestrial arthropods in Antarctica.

Molecular pedigree reconstruction and estimation of evolutionary parameters in a wild Atlantic salmon river system with incomplete sampling: a power analysis

Background: Pedigree reconstruction using genetic analysis provides a useful means to estimate fundamental population biology parameters relating to population demography, trait heritability and individual fitness when combined with other sources of data. However, there remain limitations to pedigree reconstruction in wild populations, particularly in systems where parent-offspring relationships cannot be directly observed, there is incomplete sampling of individuals, or molecular parentage inference relies on low quality DNA from archived material. While much can still be inferred from incomplete or sparse pedigrees, it is crucial to evaluate the quality and power of available genetic information a priori to testing specific biological hypotheses. Here, we used microsatellite markers to reconstruct a multi-generation pedigree of wild Atlantic salmon (Salmo salar L.) using archived scale samples collected with a total trapping system within a river over a 10 year period. Using a simulation-based approach, we determined the optimal microsatellite marker number for accurate parentage assignment, and evaluated the power of the resulting partial pedigree to investigate important evolutionary and quantitative genetic characteristics of salmon in the system.

Results: We show that at least 20 microsatellites (ave. 12 alleles/locus) are required to maximise parentage assignment and to improve the power to estimate reproductive success and heritability in this study system. We also show that 1.5 fold differences can be detected between groups simulated to have differing reproductive success, and that it is possible to detect moderate heritability values for continuous traits (h(2) similar to 0.40) with more than 80% power when using 28 moderately to highly polymorphic markers.

Conclusion: The methodologies and work flow described provide a robust approach for evaluating archived samples for pedigree-based research, even where only a proportion of the total population is sampled. The results demonstrate the feasibility of pedigree-based studies to address challenging ecological and evolutionary questions in free-living populations, where genealogies can be traced only using molecular tools, and that significant increases in pedigree assignment power can be achieved by using higher numbers of markers.

Constructing networks of continuous-time velocity-based models

A conventional local model (LM) network consists of a set of affine local models blended together using appropriate weighting functions. Such networks have poor interpretability since the dynamics of the blended network are only weakly related to the underlying local models. In contrast, velocity-based LM networks employ strictly linear local models to provide a transparent framework for nonlinear modelling in which the global dynamics are a simple linear combination of the local model dynamics. A novel approach for constructing continuous-time velocity-based networks from plant data is presented. Key issues including continuous-time parameter estimation, correct realisation of the velocity-based local models and avoidance of the input derivative are all addressed. Application results are reported for the highly nonlinear simulated continuous stirred tank reactor process.

A meta-analysis of latitudinal variations in life-history traits of roach, Rutilus rutilus, over its geographical range: linear or non-linear relationships?

1. We collated information from the literature on life history traits of the roach (a generalist freshwater fish), and analysed variation in absolute fecundity, von Bertalanffy parameters, and reproductive lifespan in relation to latitude, using both linear and non-linear regression models. We hypothesized that because most life history traits are dependent on growth rate, and growth rate is non-linearly related with temperature, it was likely that when analysed over the whole distribution range of roach, variation in key life history traits would show non-linear patterns with latitude.

Performance models of storage contention in cloud environments

We propose simple models to predict the performance degradation of disk requests due to storage device contention in consolidated virtualized environments. Model parameters can be deduced from measurements obtained inside Virtual Machines (VMs) from a system where a single VM accesses a remote storage server. The parameterized model can then be used to predict the effect of storage contention when multiple VMs are consolidated on the same server. We first propose a trace-driven approach that evaluates a queueing network with fair share scheduling using simulation. The model parameters consider Virtual Machine Monitor level disk access optimizations and rely on a calibration technique. We further present a measurement-based approach that allows a distinct characterization of read/write performance attributes. In particular, we define simple linear prediction models for I/O request mean response times, throughputs and read/write mixes, as well as a simulation model for predicting response time distributions. We found our models to be effective in predicting such quantities across a range of synthetic and emulated application workloads. 

Genome-wide association study identifies novel loci associated with resistance to bovine tuberculosis

Tuberculosis (TB) caused by Mycobacterium bovis is a re-emerging disease of livestock that is of major economic importance worldwide, as well as being a zoonotic risk there is significant heritability for host resistance to bovine TB (bTB) in dairy cattle. To identify resistance loci for bTB, we undertook a genome-wide association study in female Holstein-Friesian cattle with 592 cases and 559 age-matched controls from case herds. Cases and controls were categorised into distinct phenotypes: skin test and lesion positive vs skin test negative on multiple occasions, respectively these animals were genotyped with the Illumina BovineHD 700K BeadChip. Genome-wide rapid association using linear and logistic mixed models and regression (GRAMMAR), regional heritability mapping (RHM) and haplotype-sharing analysis identified two novel resistance loci that attained chromosome-wise significance, protein tyrosine phosphatase receptor T (PTPRT; P=4.8 × 10 -7) and myosin IIIB (MYO3B; P=5.4 × 10 -6). We estimated that 21% of the phenotypic variance in TB resistance could be explained by all of the informative single-nucleotide polymorphisms, of which the region encompassing the PTPRT gene accounted for 6.2% of the variance and a further 3.6% was associated with a putative copy number variant in MYO3B the results from this study add to our understanding of variation in host control of infection and suggest that genetic marker-based selection for resistance to bTB has the potential to make a significant contribution to bTB control.

Integrating personality research and animal contest theory:Aggressiveness in the green swordtail xiphophorus helleri

Aggression occurs when individuals compete over limiting resources. While theoretical studies have long placed a strong emphasis on context-specificity of aggression, there is increasing recognition that consistent behavioural differences exist among individuals, and that aggressiveness may be an important component of individual personality. Though empirical studies tend to focus on one aspect or the other, we suggest there is merit in modelling both within- and among-individual variation in agonistic behaviour simultaneously. Here, we demonstrate how this can be achieved using multivariate linear mixed effect models. Using data from repeated mirror trials and dyadic interactions of male green swordtails, Xiphophorus helleri, we show repeatable components of (co)variation in a suite of agonistic behaviour that is broadly consistent with a major axis of variation in aggressiveness. We also show that observed focal behaviour is dependent on opponent effects, which can themselves be repeatable but were more generally found to be context specific. In particular, our models show that within-individual variation in agonistic behaviour is explained, at least in part, by the relative size of a live opponent as predicted by contest theory. Finally, we suggest several additional applications of the multivariate models demonstrated here. These include testing the recently queried functional equivalence of alternative experimental approaches, (e.g., mirror trials, dyadic interaction tests) for assaying individual aggressiveness. © 2011 Wilson et al.

Outlier effects on robust joint modelling of longitudinal and survival data

Robust joint modelling is an emerging field of research. Through the advancements in electronic patient healthcare records, the popularly of joint modelling approaches has grown rapidly in recent years providing simultaneous analysis of longitudinal and survival data. This research advances previous work through the development of a novel robust joint modelling methodology for one of the most common types of standard joint models, that which links a linear mixed model with a Cox proportional hazards model. Through t-distributional assumptions, longitudinal outliers are accommodated with their detrimental impact being down weighed and thus providing more efficient and reliable estimates. The robust joint modelling technique and its major benefits are showcased through the analysis of Northern Irish end stage renal disease patients. With an ageing population and growing prevalence of chronic kidney disease within the United Kingdom, there is a pressing demand to investigate the detrimental relationship between the changing haemoglobin levels of haemodialysis patients and their survival. As outliers within the NI renal data were found to have significantly worse survival, identification of outlying individuals through robust joint modelling may aid nephrologists to improve patient's survival. A simulation study was also undertaken to explore the difference between robust and standard joint models in the presence of increasing proportions and extremity of longitudinal outliers. More efficient and reliable estimates were obtained by robust joint models with increasing contrast between the robust and standard joint models when a greater proportion of more extreme outliers are present. Through illustration of the gains in efficiency and reliability of parameters when outliers exist, the potential of robust joint modelling is evident. The research presented in this thesis highlights the benefits and stresses the need to utilise a more robust approach to joint modelling in the presence of longitudinal outliers.

Is shyness context specific? Relation between shyness and online self-disclosure with and without a live webcam in young adults

The Internet provides a new tool to investigate old questions in experimental social psychology regarding Person x Context interaction. We examined the interaction of self-reported shyness and context on computer-mediated communication measures. Sixty female undergraduates unfamiliar were paired in dyads and engaged in a 10 min free chat conversation on the Internet with and without a live webcam. Free chat conversations were archived, transcripts were objectively coded for communication variables, and a linear mixed model used for data analysis of dyadic interaction was performed on each communication measure. As predicted, increases in self-reported shyness were significantly related to decreases in the number of prompted self-disclosures (after controlling for the number of opportunities to self-disclose) only in the webcam condition. Self-reported shyness was not related to the number of prompted self-disclosures in the no webcam condition, suggesting that shyness was context dependent. The present study appears to be the first to objectively code measures of Internet behaviour in relation to the study of personality in general and shyness in particular. Theoretical and clinical implications for understanding the contextual nature of shyness are discussed. (C) 2006 Elsevier Inc. All rights reserved.

Transient Crystal Size Distribution in DTB and FC Crystallizers

This work deals with the transient analysis of crystal size distribution (CSD) for imperfectly mixed draft tube baffled (DTB) and forced circulation (FC) crystallizers. The DTB and FC crystallizers are described by the Compartmental and Mixed models respectively. Monte Carlo (MC) scheme has been employed for simulation purposes. The simulation results have been compared with the available experimental data of BENNETT and VAN BUREN for continuous urea crystallizers.