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.

A new statistical framework for the quantification of covariate associations with species distributions

Identifying processes that shape species geographical ranges is a prerequisite for understanding environmental change. Currently, species distribution modelling methods do not offer credible statistical tests of the relative influence of climate factors and typically ignore other processes (e.g. biotic interactions and dispersal limitation). We use a hierarchical model fitted with Markov Chain Monte Carlo to combine ecologically plausible niche structures using regression splines to describe unimodal but potentially skewed response terms. We apply spatially explicit error terms that account for (and may help identify) missing variables. Using three example distributions of European bird species, we map model results to show sensitivity to change in each covariate. We show that the overall strength of climatic association differs between species and that each species has considerable spatial variation in both the strength of the climatic association and the sensitivity to climate change. Our methods are widely applicable to many species distribution modelling problems and enable accurate assessment of the statistical importance of biotic and abiotic influences on distributions.

Cross-Cultural Patterns in Dynamic Ratings of Positive and Negative Natural Emotional Behaviour

Background: Studies of cross-cultural variations in the perception of emotion have typically compared rates of recognition of static posed stimulus photographs. That research has provided evidence for universality in the recognition of a range of emotions but also for some systematic cross-cultural variation in the interpretation of emotional expression. However, questions remain about how widely such findings can be generalised to real life emotional situations. The present study provides the first evidence that the previously reported interplay between universal and cultural influences extends to ratings of natural, dynamic emotional stimuli.

Methodology/Principal Findings: Participants from Northern Ireland, Serbia, Guatemala and Peru used a computer based tool to continuously rate the strength of positive and negative emotion being displayed in twelve short video sequences by people from the United Kingdom engaged in emotional conversations. Generalized additive mixed models were developed to assess the differences in perception of emotion between countries and sexes. Our results indicate that the temporal pattern of ratings is similar across cultures for a range of emotions and social contexts. However, there are systematic differences in intensity ratings between the countries, with participants from Northern Ireland making the most extreme ratings in the majority of the clips.

Conclusions/Significance: The results indicate that there is strong agreement across cultures in the valence and patterns of ratings of natural emotional situations but that participants from different cultures show systematic variation in the intensity with which they rate emotion. Results are discussed in terms of both ‘in-group advantage’ and ‘display rules’ approaches. This study indicates that examples of natural spontaneous emotional behaviour can be used to study cross-cultural variations in the perception of emotion.

Regression analysis of spatial data

Many of the most interesting questions ecologists ask lead to analyses of spatial data. Yet, perhaps confused by the large number of statistical models and fitting methods available, many ecologists seem to believe this is best left to specialists. Here, we describe the issues that need consideration when analysing spatial data and illustrate these using simulation studies. Our comparative analysis involves using methods including generalized least squares, spatial filters, wavelet revised models, conditional autoregressive models and generalized additive mixed models to estimate regression coefficients from synthetic but realistic data sets, including some which violate standard regression assumptions. We assess the performance of each method using two measures and using statistical error rates for model selection. Methods that performed well included generalized least squares family of models and a Bayesian implementation of the conditional auto-regressive model. Ordinary least squares also performed adequately in the absence of model selection, but had poorly controlled Type I error rates and so did not show the improvements in performance under model selection when using the above methods. Removing large-scale spatial trends in the response led to poor performance. These are empirical results; hence extrapolation of these findings to other situations should be performed cautiously. Nevertheless, our simulation-based approach provides much stronger evidence for comparative analysis than assessments based on single or small numbers of data sets, and should be considered a necessary foundation for statements of this type in future.

Role of Predalors, Habitat Attributes, and Spatial Autocorrelation on the Distribution of Eggs in the Northern Spectacled Salamander (Salamandrina perspicillata)

We predicted that the probability of egg occurrence of salamander Salamandrina perspicillata depended on stream features and predation by native crayfish Austropotamobius fulcisianus and the introduced trout Salmo trutta. We assessed the presence of S. perspicillata at 54 sites within a natural reserve of southern Tuscany, Italy. Generalized linear models with binomial errors were constructed using egg presence/absence and altitude, stream mean size and slope, electrical conductivity, water pH and temperature, and a predation factor, defined according to the presence/absence of crayfish and trout. Some competing models also included an autocovariate term, which estimated how much the response variable at any one sampling point reflected response values at surrounding points. The resulting models were compared using Akaike's information criterion. Model selection led to a subset of 14 models with Delta AIC(c) <7 (i.e., models ranging from substantial support to considerably less support), and all but one of these included an effect of predation. Models with the autocovariate term had considerably more support than those without the term. According to multimodel inference, the presence of trout and crayfish reduced the probability of egg occurrence from a mean level of 0.90 (SE limits: 0.98-0.55) to 0.12 (SE limits: 0.34-0.04). The presence of crayfish alone had no detectable effects (SE limits: 0.86-0.39). The results suggest that introduced trout have a detrimental effect on the reproductive output of S. perspicillata and confirm the fundamental importance of distinguishing the roles of endogenous and exogenous forces that act on population distribution.

Large-scale spatial patterns in the distribution of Collembola (Hexapoda) species in Antarctic terrestrial ecosystems

We tested whether the distribution of three common springtail species (Gressittacantha terranova, Gomphiocephalus hodgsoni and Friesea grisea) in Victoria Land (Antarctica) could be modelled as a function of latitude, longitude, altitude and distance from the sea.

Victoria Land, Ross Dependency, Antarctica.

Generalized linear models were constructed using species presence/absence data relative to geographical features (latitude, longitude, altitude, distance from sea) across the species' entire ranges. Model results were then integrated with the known phylogeography of each species and hypotheses were generated on the role of climate as a major driver of Antarctic springtail distribution.

Based on model selection using Akaike's information criterion, the species' distributions were: hump-shaped relative to longitude and monotonic with altitude for Gressittacantha terranova; hump-shaped relative to latitude and monotonic with altitude for Gomphiocephalus hodgsoni; and hump-shaped relative to longitude and monotonic with latitude, altitude and distance from the sea for Friesea grisea.

No single distributional pattern was shared by the three species. While distributions were partially a response to climatic spatial clines, the patterns observed strongly suggest that past geological events have influenced the observed distributions. Accordingly, present-day spatial patterns are likely to have arisen from the interaction of historical and environmental drivers. Future studies will need to integrate a range of spatial and temporal scales to further quantify their respective roles.

Spatial variation in adult sex ratio across multiple scales in the invasive golden apple snail, Pomacea canaliculata

Adult sex ratio (ASR) has critical effects on behavior and life history and has implications for population demography, including the invasiveness of introduced species. ASR exhibits immense variation in nature, yet the scale dependence of this variation is rarely analyzed. In this study, using the generalized multilevel models, we investigated the variation in ASR across multiple nested spatial scales and analyzed the underlying causes for an invasive species, the golden apple snail Pomacea canaliculata. We partitioned the variance in ASR to describe the variations at different scales and then included the explanatory variables at the individual and group levels to analyze the potential causes driving the variation in ASR. We firstly determined there is a significant female-biased ASR for this species when accounting for the spatial and temporal autocorrelations of sampling. We found that, counter to nearly equal distributed variation at plot, habitat and region levels, ASR showed little variation at the town level. Temperature and precipitation at the region level were significantly positively associated with ASR, whereas the individual weight, the density characteristic, and sampling time were not significant factors influencing ASR. Our study suggests that offspring sex ratio of this species may shape the general pattern of ASR in the population level while the environmental variables at the region level translate the unbiased offspring sex ratio to the female-biased ASR. Future research should consider the implications of climate warming on the female-biased ASR of this invasive species and thus on invasion pattern.

Generalized social trust across 22 European countries: A multi-level analysis of four competing models

Concern with what can explain variation in generalized social trust has led to an abundance of theoretical models. Defining generalized social trust as a belief in human benevolence, we focus on the emancipation theory and social capital theory as well as the ethnic diversity and economic development models of trust. We then determine which dimensions of individuals’ behavior and attitudes as well as of their national context are the most important predictors. Using data from 20 countries that participated in round one of the European Social Survey, we test these models at their respective level of analysis, individual and/or national. Our analysis revealed that individuals’ own trust in the political system as a moral and competent institution was the most important predictor of generalized social trust at the individual level, while a country’s level of affluence was the most important contextual predictor, indicating that different dimensions are significant at the two levels of analysis. This analysis also raised further questions as to the meaning of social capital at the two levels of analysis and the conceptual equivalence of its civic engagement dimension across cultures.

Physical Models of Wind Instruments : A Generalized Excitation Coupled with a Modular Tube Simulation Platform

To develop real-time simulations of wind instruments, digital waveguides filters can be used as an efficient representation of the air column. Many aerophones are shaped as horns which can be approximated using conical sections. Therefore the derivation of conical waveguide filters is of special interest. When these filters are used in combination with a generalized reed excitation, several classes of wind instruments can be simulated. In this paper we present the methods for transforming a continuous description of conical tube segments to a discrete filter representation. The coupling of the reed model with the conical waveguide and a simplified model of the termination at the open end are described in the same way. It turns out that the complete lossless conical waveguide requires only one type of filter.Furthermore, we developed a digital reed excitation model, which is purely based on numerical integration methods, i.e., without the use of a look-up table.

Generalized compound transport of charged particles in turbulent magnetized plasmas

The transport of charged particles in partially turbulent magnetic systems is investigated from first principles. A generalized compound transport model is proposed, providing an explicit relation between the mean-square deviation of the particle parallel and perpendicular to a magnetic mean field, and the mean-square deviation which characterizes the stochastic field-line topology. The model is applied in various cases of study, and the relation to previous models is discussed.

Simulated predator extinctions: Predator identity affects survival and recruitment of oysters

The rate of species loss is increasing at a global scale, and human-induced extinctions are biased toward predator species. We examined the effects of predator extinctions on a foundation species, the eastern oyster (Crassostrea virginica). We performed a factorial experiment manipulating the presence and abundance of three of the most common predatory crabs, the blue crab (Callinectes sapidus), stone crab (Menippe mercenaria), and mud crab (Panopeus herbstii) in estuaries in the eastern United States. We tested the effects of species richness and identity of predators on juvenile oyster survival, oyster recruitment, and organic matter content of sediment. We also manipulated the density of each of the predators and controlled for the loss of biomass of species by maintaining a constant mass of predators in one set of treatments and simultaneously using an additive design. This design allowed us to test the density dependence of our results and test for functional compensation by other species. The identity of predator species, but not richness, affected oyster populations. The loss of blue crabs, alone or in combination with either of the other species, affected the survival rate of juvenile oysters. Blue crabs and stone crabs both affected oyster recruitment and sediment organic matter negatively. Mud crabs at higher than ambient densities, however, could fulfill some of the functions of blue and stone crabs, suggesting a level of ecological redundancy. Importantly, the strong effects of blue crabs in all processes measured no longer occurred when individuals were present at higher-than-ambient densities. Their role as dominant predator is, therefore, dependent on their density within the system and the density of other species within their guild (e.g., mud crabs). Our findings support the hypothesis that the effects of species loss at higher trophic levels are determined by predator identity and are subject to complex intraguild interactions that are largely density dependent. Understanding the role of biodiversity in ecosystem functioning or addressing practical concerns, such as loss of predators owing to overharvesting, remains complicated because accurate predictions require detailed knowledge of the system and should be drawn from sound experimental evidence, not based on observations or generalized models.

Alternative Factor Models and Heritability of the Short Leyton Obsessional Inventory-Children's Version

An alternative models framework was used to test three confirmatory factor analytic models for the Short Leyton Obsessional Inventory-Children's Version (Short LOI-CV) in a general population sample of 517 young adolescent twins (11-16 years). A one-factor model as implicit in current classification systems of Obsessive-Compulsive Disorder (OCD), a two-factor obsessions and compulsions model, and a multidimensional model corresponding to the three proposed subscales of the Short LOI-CV (labelled Obsessions/Incompleteness, Numbers/Luck and Cleanliness) were considered. The three-factor model was the only model to provide an adequate explanation of the data. Twin analyses suggested significant quantitative sex differences in heritability for both the Obsessions/Incompleteness and Numbers/Luck dimensions with these being significantly heritable in males only (heritability of 60% and 65% respectively). The correlation between the additive genetic effects for these two dimensions in males was 0.95 suggesting they largely share the same genetic risk factors.

Algorithms for Hidden Markov Models With Imprecisely Specified Parameters

Hidden Markov models (HMMs) are widely used models for sequential data. As with other probabilistic graphical models, they require the specification of precise probability values, which can be too restrictive for some domains, especially when data are scarce or costly to acquire. We present a generalized version of HMMs, whose quantification can be done by sets of, instead of single, probability distributions. Our models have the ability to suspend judgment when there is not enough statistical evidence, and can serve as a sensitivity analysis tool for standard non-stationary HMMs. Efficient inference algorithms are developed to address standard HMM usage such as the computation of likelihoods and most probable explanations. Experiments with real data show that the use of imprecise probabilities leads to more reliable inferences without compromising efficiency.

Generalized Loopy 2U: A New Algorithm for Approximate Inference in Credal Networks

Credal nets generalize Bayesian nets by relaxing the requirement of precision of probabilities. Credal nets are considerably more expressive than Bayesian nets, but this makes belief updating NP-hard even on polytrees. We develop a new efficient algorithm for approximate belief updating in credal nets. The algorithm is based on an important representation result we prove for general credal nets: that any credal net can be equivalently reformulated as a credal net with binary variables; moreover, the transformation, which is considerably more complex than in the Bayesian case, can be implemented in polynomial time. The equivalent binary credal net is updated by L2U, a loopy approximate algorithm for binary credal nets. Thus, we generalize L2U to non-binary credal nets, obtaining an accurate and scalable algorithm for the general case, which is approximate only because of its loopy nature. The accuracy of the inferences is evaluated by empirical tests.