Adaptive Gaussian Markov random fields with applications in human brain mapping

Summary.  Functional magnetic resonance imaging has become a standard technology in human brain mapping. Analyses of the massive spatiotemporal functional magnetic resonance imaging data sets often focus on parametric or non-parametric modelling of the temporal component, whereas spatial smoothing is based on Gaussian kernels or random fields. A weakness of Gaussian spatial smoothing is underestimation of activation peaks or blurring of high curvature transitions between activated and non-activated regions of the brain. To improve spatial adaptivity, we introduce a class of inhomogeneous Markov random fields with stochastic interaction weights in a space-varying coefficient model. For given weights, the random field is conditionally Gaussian, but marginally it is non-Gaussian. Fully Bayesian inference, including estimation of weights and variance parameters, can be carried out through efficient Markov chain Monte Carlo simulation. Although motivated by the analysis of functional magnetic resonance imaging data, the methodological development is general and can also be used for spatial smoothing and regression analysis of areal data on irregular lattices. An application to stylized artificial data and to real functional magnetic resonance imaging data from a visual stimulation experiment demonstrates the performance of our approach in comparison with Gaussian and robustified non-Gaussian Markov random-field models.     

Venezuelan Rainfall Data Analysed by Using a Bayesian Space–time Model

We consider a set of data from 80 stations in the Venezuelan state of Guárico consisting of accumulated monthly rainfall in a time span of 16 years. The problem of modelling rainfall accumulated over fixed periods of time and recorded at meteorological stations at different sites is studied by using a model based on the assumption that the data follow a truncated and transformed multivariate normal distribution. The spatial correlation is modelled by using an exponentially decreasing correlation function and an interpolating surface for the means. Missing data and dry periods are handled within a Markov chain Monte Carlo framework using latent variables. We estimate the amount of rainfall as well as the probability of a dry period by using the predictive density of the data. We considered a model based on a full second-degree polynomial over the spatial co-ordinates as well as the first two Fourier harmonics to describe the variability during the year. Predictive inferences on the data show very realistic results, capturing the typical rainfall variability in time and space for that region. Important extensions of the model are also discussed.     

Zero-inflated sum of Conway-Maxwell-Poissons (ZISCMP) regression

While excess zeros are often thought to cause data over-dispersion (i.e. when the variance exceeds the mean), this implication is not absolute. One should instead consider a flexible class of distributions that can address data dispersion along with excess zeros. This work develops a zero-inflated sum-of-Conway-Maxwell-Poissons (ZISCMP) regression as a flexible analysis tool to model count data that express significant data dispersion and contain excess zeros. This class of models contains several special case zero-inflated regressions, including zero-inflated Poisson (ZIP), zero-inflated negative binomial (ZINB), zero-inflated binomial (ZIB), and the zero-inflated Conway-Maxwell-Poisson (ZICMP). Through simulated and real data examples, we demonstrate class flexibility and usefulness. We further utilize it to analyze shark species data from Australia's Great Barrier Reef to assess the environmental impact of human action on the number of various species of sharks.     

Room for manoeuvre: User participation in water resources management in Burkina Faso

Participation in natural resources management is widely promoted in sub‐Saharan Africa, but faces mounting criticism from social science scholars who rarely engage with it in practice. We use the notion of room for manoeuvre to reflect on a multi‐level participatory approach designed to support the Burkinabè Integrated Water Resources Management policy and propose ways of engaging constructively with local users and policy‐makers. Within an “invited space” of participation, water users’ room for manoeuvre was enhanced through the acquisition of new knowledge on the legal Integrated Water Resources Management (IWRM) framework. This allowed them to interact with decision‐makers and to extend their networks. Power imbalances were discussed, but this did not lead to modifying existing power structures that hinge on broader societal dynamics.     

Faking it: Personality and individual difference predictors of willingness to buy counterfeit goods

Counterfeiting is now widely regarded as a serious social, economic, and political issue. This study examined demographic, personality, and individual difference predictors of willingness to buy counterfeit goods (WBCG) in a community sample of British adults. Two-hundred and thirty-seven participants completed a questionnaire that assessed their willingness to buy 19 types of counterfeit goods, attitudes towards counterfeiting, material values, Big Five personality traits, and demographics. Structural equation modelling showed that attitudes towards counterfeiting were the strongest predictors of WBCG. In addition, material values predicted both WBCG directly and indirectly through attitudes towards counterfeiting. Older participants showed lower WBCG, although this effect was moderated by participants’ conscientiousness, material values, and attitudes towards counterfeiting. There were no sex differences in WBCG once participants’ age and conscientiousness had been taken into account. These results are discussed in relation to the extant literature on the profiles of consumers who buy counterfeit products.     

Estimation and inference of the joint conditional distribution for multivariate longitudinal data using nonparametric copulas

In this paper we study estimating the joint conditional distributions of multivariate longitudinal outcomes using regression models and copulas. For the estimation of marginal models, we consider a class of time-varying transformation models and combine the two marginal models using nonparametric empirical copulas. Our models and estimation method can be applied in many situations where the conditional mean-based models are not good enough. Empirical copulas combined with time-varying transformation models may allow quite flexible modelling for the joint conditional distributions for multivariate longitudinal data. We derive the asymptotic properties for the copula-based estimators of the joint conditional distribution functions. For illustration we apply our estimation method to an epidemiological study of childhood growth and blood pressure.     

On generalised estimating equations for vector regression

Generalised estimating equations (GEE) for regression problems with vector‐valued responses are examined. When the response vectors are of mixed type (e.g. continuous–binary response pairs), the GEE approach is a semiparametric alternative to full‐likelihood copula methods, and is closely related to Prentice & Zhao's mean‐covariance estimation equations approach. When the response vectors are of the same type (e.g. measurements on left and right eyes), the GEE approach can be viewed as a ‘plug‐in’ to existing methods, such as the vglm function from the state‐of‐the‐art VGAM package in R. In either scenario, the GEE approach offers asymptotically correct inferences on model parameters regardless of whether the working variance–covariance model is correctly or incorrectly specified. The finite‐sample performance of the method is assessed using simulation studies based on a burn injury dataset and a sorbinil eye trial dataset. The method is applied to data analysis examples using the same two datasets, as well as to a trivariate binary dataset on three plant species in the Hunua ranges of Auckland.     

Spatio‐temporal Ornstein–Uhlenbeck Processes: Theory,Simulation and Statistical Inference

Spatio‐temporal modelling is an increasingly popular topic in Statistics. Our paper contributes to this line of research by developing the theory, simulation and inference for a spatio‐temporal Ornstein–Uhlenbeck process. We conduct detailed simulation studies and demonstrate the practical relevance of these processes in an empirical study of radiation anomaly data. Finally, we describe how predictions can be carried out in the Gaussian setting.     

Modelling the travel time of transit vehicles in real‐time through a GTFS‐based road network using GPS vehicle locations

Predicting the arrival time of a transit vehicle involves not only knowledge of its current position and schedule adherence, but also traffic conditions along the remainder of the route. Road networks are dynamic and can quickly change from free‐flowing to highly congested, which impacts the arrival time of transit vehicles, particularly buses which often share the road with other vehicles, so reliable predictions need to account for real‐time and future traffic conditions. The first step in this process is to construct a framework with which road state (traffic conditions) can be estimated using real‐time transit vehicle position data. Our proposed framework implements a vehicle model using a particle filter to estimate road travel times, which are used in a second model to estimate real‐time traffic conditions. Although development and testing took place in Auckland, New Zealand, we generalised each component to make the framework compatible with other public transport systems around the world. We demonstrate the real‐time feasibility and performance of our approach in real‐time, where a combination of R and C++ was used to obtain the necessary performance results. Future work will use these estimated traffic conditions in combination with historical data to obtain reliable arrival time predictions of transit vehicles.     

A dynamic business modelling approach to design and experiment new business venture strategies

Business Modelling has evolved as a key activity to reflect new business venture strategy by framing the way a firm will operate and how it will function in achieving its goals (e.g., profitability, growth, innovation, social impact). However, scholars and practitioners have criticized the adoption of a too static perspective in the design and use of conventional Business Model representations. Such a static perspective prevents nascent entrepreneurs experimenting with their Business Models and, as a result, identifying the most effective strategies, especially in terms of business sustainability and profitability. In this paper, we argue that combining conventional Business Model schemas with System Dynamics modelling results in a strategy design tool that may overcome several limitations related to a static view of Business Model representation. Mapping the different key elements underlying value creation processes into a system of causal interdependencies – through the use of simulation – allows strategy analysts and entrepreneurs to experiment and learn how the business reacts to strategic and organizational changes in terms of performance, innovation and value creation. As such, Dynamic Business Models provide useful insights to strategy formulation and business venturing by capturing how critical Business Model elements interact to produce enduring competitive advantages over time.