The Bayesian conditional independence model for measurement error: applications in ecology

The measurement error model is a well established statistical method for regression problems in medical sciences, although rarely used in ecological studies. While the situations in which it is appropriate may be less common in ecology, there are instances in which there may be benefits in its use for prediction and estimation of parameters of interest. We have chosen to explore this topic using a conditional independence model in a Bayesian framework using a Gibbs sampler, as this gives a great deal of flexibility, allowing us to analyse a number of different models without losing generality. Using simulations and two examples, we show how the conditional independence model can be used in ecology, and when it is appropriate.

An infinite-horizon robust filter for uncertain hidden Markov Models with conditional relative entropy constraints

We consider a robust filtering problem for uncertain discrete-time, homogeneous, first-order, finite-state hidden Markov models (HMMs). The class of uncertain HMMs considered is described by a conditional relative entropy constraint on measures perturbed from a nominal regular conditional probability distribution given the previous posterior state distribution and the latest measurement. Under this class of perturbations, a robust infinite horizon filtering problem is first formulated as a constrained optimization problem before being transformed via variational results into an unconstrained optimization problem; the latter can be elegantly solved using a risk-sensitive information-state based filtering.

A conditional retransmission enabled transport protocol for real-time networked control systems

Real-time networked control systems (NCSs) over data networks are being increasingly implemented on a massive scale in industrial applications. Along with this trend, wireless network technologies have been promoted for modern wireless NCSs (WNCSs). However, popular wireless network standards such as IEEE 802.11/15/16 are not designed for real-time communications. Key issues in real-time applications include limited transmission reliability and poor transmission delay performance. Considering the unique features of real-time control systems, this paper develops a conditional retransmission enabled transport protocol (CRETP) to improve the delay performance of the transmission control protocol (TCP) and also the reliability performance of the user datagram protocol (UDP) and its variants. Key features of the CRETP include a connectionless mechanism with acknowledgement (ACK), conditional retransmission and detection of ineffective data packets on the receiver side.

Using Bayesian methods for the estimation of uncertainty in complex statistical models

The research objectives of this thesis were to contribute to Bayesian statistical methodology by contributing to risk assessment statistical methodology, and to spatial and spatio-temporal methodology, by modelling error structures using complex hierarchical models. Specifically, I hoped to consider two applied areas, and use these applications as a springboard for developing new statistical methods as well as undertaking analyses which might give answers to particular applied questions. Thus, this thesis considers a series of models, firstly in the context of risk assessments for recycled water, and secondly in the context of water usage by crops. The research objective was to model error structures using hierarchical models in two problems, namely risk assessment analyses for wastewater, and secondly, in a four dimensional dataset, assessing differences between cropping systems over time and over three spatial dimensions. The aim was to use the simplicity and insight afforded by Bayesian networks to develop appropriate models for risk scenarios, and again to use Bayesian hierarchical models to explore the necessarily complex modelling of four dimensional agricultural data. The specific objectives of the research were to develop a method for the calculation of credible intervals for the point estimates of Bayesian networks; to develop a model structure to incorporate all the experimental uncertainty associated with various constants thereby allowing the calculation of more credible credible intervals for a risk assessment; to model a single day’s data from the agricultural dataset which satisfactorily captured the complexities of the data; to build a model for several days’ data, in order to consider how the full data might be modelled; and finally to build a model for the full four dimensional dataset and to consider the timevarying nature of the contrast of interest, having satisfactorily accounted for possible spatial and temporal autocorrelations. This work forms five papers, two of which have been published, with two submitted, and the final paper still in draft. The first two objectives were met by recasting the risk assessments as directed, acyclic graphs (DAGs). In the first case, we elicited uncertainty for the conditional probabilities needed by the Bayesian net, incorporated these into a corresponding DAG, and used Markov chain Monte Carlo (MCMC) to find credible intervals, for all the scenarios and outcomes of interest. In the second case, we incorporated the experimental data underlying the risk assessment constants into the DAG, and also treated some of that data as needing to be modelled as an ‘errors-invariables’ problem [Fuller, 1987]. This illustrated a simple method for the incorporation of experimental error into risk assessments. In considering one day of the three-dimensional agricultural data, it became clear that geostatistical models or conditional autoregressive (CAR) models over the three dimensions were not the best way to approach the data. Instead CAR models are used with neighbours only in the same depth layer. This gave flexibility to the model, allowing both the spatially structured and non-structured variances to differ at all depths. We call this model the CAR layered model. Given the experimental design, the fixed part of the model could have been modelled as a set of means by treatment and by depth, but doing so allows little insight into how the treatment effects vary with depth. Hence, a number of essentially non-parametric approaches were taken to see the effects of depth on treatment, with the model of choice incorporating an errors-in-variables approach for depth in addition to a non-parametric smooth. The statistical contribution here was the introduction of the CAR layered model, the applied contribution the analysis of moisture over depth and estimation of the contrast of interest together with its credible intervals. These models were fitted using WinBUGS [Lunn et al., 2000]. The work in the fifth paper deals with the fact that with large datasets, the use of WinBUGS becomes more problematic because of its highly correlated term by term updating. In this work, we introduce a Gibbs sampler with block updating for the CAR layered model. The Gibbs sampler was implemented by Chris Strickland using pyMCMC [Strickland, 2010]. This framework is then used to consider five days data, and we show that moisture in the soil for all the various treatments reaches levels particular to each treatment at a depth of 200 cm and thereafter stays constant, albeit with increasing variances with depth. In an analysis across three spatial dimensions and across time, there are many interactions of time and the spatial dimensions to be considered. Hence, we chose to use a daily model and to repeat the analysis at all time points, effectively creating an interaction model of time by the daily model. Such an approach allows great flexibility. However, this approach does not allow insight into the way in which the parameter of interest varies over time. Hence, a two-stage approach was also used, with estimates from the first-stage being analysed as a set of time series. We see this spatio-temporal interaction model as being a useful approach to data measured across three spatial dimensions and time, since it does not assume additivity of the random spatial or temporal effects.

Expert elicitation into Bayesian networks, Part 4 of 4: eliciting a conditional probability table

An introduction to eliciting a conditional probability table in a Bayesian Network model, highlighting three efficient methods for populating a CPT.

Comparison of three-dimensional profiles over time

In this paper, we describe an analysis for data collected on a three-dimensional spatial lattice with treatments applied at the horizontal lattice points. Spatial correlation is accounted for using a conditional autoregressive model. Observations are defined as neighbours only if they are at the same depth. This allows the corresponding variance components to vary by depth. We use the Markov chain Monte Carlo method with block updating, together with Krylov subspace methods, for efficient estimation of the model. The method is applicable to both regular and irregular horizontal lattices and hence to data collected at any set of horizontal sites for a set of depths or heights, for example, water column or soil profile data. The model for the three-dimensional data is applied to agricultural trial data for five separate days taken roughly six months apart in order to determine possible relationships over time. The purpose of the trial is to determine a form of cropping that leads to less moist soils in the root zone and beyond.We estimate moisture for each date, depth and treatment accounting for spatial correlation and determine relationships of these and other parameters over time.

A Bayesian analysis of an agricultural field trial with three spatial dimensions

Modern technology now has the ability to generate large datasets over space and time. Such data typically exhibit high autocorrelations over all dimensions. The field trial data motivating the methods of this paper were collected to examine the behaviour of traditional cropping and to determine a cropping system which could maximise water use for grain production while minimising leakage below the crop root zone. They consist of moisture measurements made at 15 depths across 3 rows and 18 columns, in the lattice framework of an agricultural field. Bayesian conditional autoregressive (CAR) models are used to account for local site correlations. Conditional autoregressive models have not been widely used in analyses of agricultural data. This paper serves to illustrate the usefulness of these models in this field, along with the ease of implementation in WinBUGS, a freely available software package. The innovation is the fitting of separate conditional autoregressive models for each depth layer, the ‘layered CAR model’, while simultaneously estimating depth profile functions for each site treatment. Modelling interest also lay in how best to model the treatment effect depth profiles, and in the choice of neighbourhood structure for the spatial autocorrelation model. The favoured model fitted the treatment effects as splines over depth, and treated depth, the basis for the regression model, as measured with error, while fitting CAR neighbourhood models by depth layer. It is hierarchical, with separate onditional autoregressive spatial variance components at each depth, and the fixed terms which involve an errors-in-measurement model treat depth errors as interval-censored measurement error. The Bayesian framework permits transparent specification and easy comparison of the various complex models compared.

Cellular settings mediating Src Substrate switching between Focal Adhesion Kinase Tyrosine 861 and CUB-domain-containing protein 1 (CDCP1) Tyrosine 734

Reciprocal interactions between Src family kinases (SFKs) and focal adhesion kinase (FAK) are critical during changes in cell attachment. Recently it has been recognized that another SFK substrate, CUB-domain-containing protein 1 (CDCP1), is differentially phosphorylated during these events. However, the molecular processes underlying SFK-mediated phosphorylation of CDCP1 are poorly understood. Here we identify a novel mechanism in which FAK tyrosine 861 and CDCP1-Tyr-734 compete as SFK substrates and demonstrate cellular settings in which SFKs switch between these sites. Our results show that stable CDCP1 expression induces robust SFK-mediated phosphorylation of CDCP1-Tyr-734 with concomitant loss of p-FAK-Tyr-861 in adherent HeLa cells. SFK substrate switching in these cells is dependent on the level of expression of CDCP1 and is also dependent on CDCP1-Tyr-734 but is independent of CDCP1-Tyr-743 and -Tyr-762. In HeLa CDCP1 cells, engagement of SFKs with CDCP1 is accompanied by an increase in phosphorylation of Src-Tyr-416 and a change in cell morphology to a fibroblastic appearance dependent on CDCP1-Tyr-734. SFK switching between FAK-Tyr-861 and CDCP1-Tyr-734 also occurs during changes in adhesion of colorectal cancer cell lines endogenously expressing these two proteins. Consistently, increased p-FAK-Tyr-861 levels and a more epithelial morphology are seen in colon cancer SW480 cells silenced for CDCP1. Unlike protein kinase Cδ, FAK does not appear to form a trimeric complex with Src and CDCP1. These data demonstrate novel aspects of the dynamics of SFK-mediated cell signaling that may be relevant during cancer progression.

Conditional corruption

In this paper we discuss whether corruption is contagious and whether conditional cooperation matters. We use the notion of “conditional corruption” for these effects. We analyze whether the justifiability to be corrupt is influenced by the perceived activities of others. Moreover, we also explore whether – and to what extent – group dynamics or socialization and past experiences affect corruption. We present evidence using two data sets at the micro level and a large macro level international panel data set. The results indicate that the willingness to engage in corruption is influenced by the perceived activities of peers and other individuals. Moreover, the panel data set at the macro level indicates that the past level of corruption has a strong impact on the current corruption level.

Sufficiently informative measurements for stability of approximate conditional mean estimates

This paper establishes sufficient conditions to bound the error in perturbed conditional mean estimates derived from a perturbed model (only the scalar case is shown in this paper but a similar result is expected to hold for the vector case). The results established here extend recent stability results on approximating information state filter recursions to stability results on the approximate conditional mean estimates. The presented filter stability results provide bounds for a wide variety of model error situations.

Interactions and consequences of inertia and switching costs

Purpose: Service research typically relates switching costs to customer loyalty, and portrays them as effective switching deterrents that engender harmful word-of-mouth (WOM). Rather than to customer loyalty, this paper aims to relate switching costs to consumer inertia, and show that while switching costs may result in customer retention, they can engender positive and negative WOM. This depends on whether the inertia stems from satisfaction or indifference. Design/methodology/approach: A mall-intercept survey investigated 518 customers' perceptions of their mobile phone service providers. Structural equation modelling fitted the data to the conceptual model. Findings: Switching costs deterred switching and engendered negative WOM, but only with low-inertia customers. With high-inertia customers, retention and WOM behaviours depended on whether the inertia stemmed from satisfaction or indifference. Satisfied customers with high switching costs tended to stay, gave more positive and less negative WOM. With indifferent customers, switching costs were unrelated to retention or WOM behaviours. Research limitations/implications: While they may be perceived negatively, switching costs can engender PWOM. Hence, research should not consider switching costs alone without considering the context that produces them. Practical implications: Service providers should segment their customers into low-inertia, high-inertia/satisfied and high-inertia/indifferent, and target each segment differently. By converting customers into the high-inertia/satisfied segment, service providers can make the best use of switching costs – not only in the traditional sense as a barrier to defection, but also as a way of generating positive WOM. Originality/value: This study is the first to consider the role of inertia with switching costs, positive WOM, and negative WOM. The findings suggest that past studies portraying switching costs as negative impediments that evoke only negative WOM might be misleading.