Applying fast initial response features on GWMA control charts for monitoring autocorrelation data

ABSTRACT

A generally weighted moving average (GWMA) control chart with fast initial response (FIR) features is addressed to monitor an autoregressive process mean shift. Numerical simulations based on average run length (ARL) show that the GWMA control chart with additional FIR feature requires less time to detect small or moderate shifts than GWMA control chart at low level of autocorrelation; whereas these two control charts perform similarly at high level of autocorrelation. Regardless of any level of autocorrelation, GWMA control charts provided with additional FIR feature have a good performance in detecting large shifts during the initial stage.     

Minimum Mean Squared Error Estimation of the Noise in Unobserved Component Models

In model-based estimation of unobserved components, the minimum mean squared error estimator of the noise component is different from white noise. In this article, some of the differences are analyzed. It is seen how the variance of the component is always underestimated, and the smaller the noise variance, the larger the underestimation. Estimators of small-variance noise components will also have large autocorrelations. Finally, in the context of an application, the sample autocorrelation function of the estimated noise is seen to perform well as a diagnostic tool, even when the variance is small and the series is of relatively short length.     

Testing the Existence of Change-Point in NHPP Software Reliability Models

In the software testing process, the nature of the failure data is affected by many factors, such as the testing environment, testing strategy, and resource allocation. These factors are unlikely to all be kept stable during the entire process of software testing. As a result, the statistical structure of the failure data is likely to experience major changes. Recently, some useful non homogeneous Poisson process (NHPP) models with change-point are proposed. However, in many realistic situations, whether a change-point exists is unknown. Furthermore, some real data seem to have two or more change-points. In this article we propose test statistics to test the existence of change-point(s). The experimental results of real data show that our tests perform well.     

Pitfalls in the Estimation of a Differenced Model

This article assesses the potential magnitude of the loss of estimation efficiency caused by the adoption of a differenced model when the disturbances of the original (levels) linear regression model follow either a stable (autoregressive) AR(1) process or a fixed start-up random-walk process (hence no filtering is necessary from the standpoint of estimation). The magnitude of the loss, which can be quite large, is found to be affected by both the form of the original model (homogeneous or nonhomogeneous) and the sign and magnitude of the autocorrelation coefficient of the AR(1) disturbance, as well as by the nature of the exogenous variable (smoothly trended or not).     

中国区域经济差异空间分异特征及其机理研究

运用全局空间自相关指数、局部空间自相关指数和空间回归模型,从空间依赖性和异质性的角度分析中国区域经济差异空间分异的过程。研究发现:中国经济发达地区和欠发达地区空间集聚格局日趋显著,且东部和中西部区域分别演变成为经济发达地区和欠发达地区的相对集聚区。造成这一分异现象的成因是各种内生因子和宏观经济环境因子的区域间差异及其循环累积作用、空间自相关性的空间近邻效应及极化——涓滴效应。     

Comparison of spline estimator at various levels of autocorrelation in smoothing spline non parametric regression for longitudinal data

The purpose of this research are: (1) to obtain spline function estimation in non parametric regression for longitudinal data with and without considering the autocorrelation between data of observation within subject, (2) to develop the algorithm that generates simulation data with certain autocorrelation level based on size of sample (N) and error variance (EV), and (3) to establish shape of spline estimator in non parametric regression for longitudinal data to simulation with various level of autocorrelation, as well as compare DM and TM approaches in predicting spline estimator in the data simulation with different of autocorrelation observational data on within subject. The results of the application are as follows: (a) implementation of smoothing spline with penalized weighted least square (PWLS) approach with or without consideration of autocorrelation in general (in all sizes and all error variances levels) provides significantly different spline estimator when the autocorrelation level >0.8; (b) based on size comparison, spline estimator in non parametric regression smoothing spline with PLS approach with (DM), or without (DM) consideration of autocorrelation showed significantly different result in level of autocorrelation > 0.8 (in overall size, moderate and large sample size), and > 0.7 (in small sample size); (c) based on level of variance, spline estimator in non parametric regression smoothing spline with PLS approach with (DM), or without (DM) consideration of autocorrelation showed significantly different result in level of autocorrelation > 0.8 (in overall level of variance, moderate and large variance), and > 0.7 (in small variance).     

A Simulation Study of Ridge Regression Estimators with Autocorrelated Errors

In this paper we run a large number of simulations to study the effects of collinearity and autocorrelated disturbances in the performance of several Ridge Regression estimators. The results suggest that with a fair amount of multicollinearity and of autocorrelation the Ridge Regression estimators which take the autocorrelation into account can perform better than the other methods. Also if the error term is only moderately autocorrelated; then the performance of the Ridge Regression estimators built upon ignoring the autocorrelation can outperform the other estimators.     

Tests for a difference in timing of physiological response between two brain regions measured by using functional magnetic resonance imaging

Summary.  Functional magnetic resonance imaging (FMRI) measures the physiological response of the human brain to experimentally controlled stimulation. In a periodically designed experiment it is of interest to test for a difference in the timing (phase shift) of the response between two anatomically distinct brain regions. We suggest two tests for an interregional difference in phase shift: one based on asymptotic theory and one based on bootstrapping. Whilst the two procedures differ in some of their assumptions, both tests rely on employing the large number of voxels (three-dimensional pixels) in non-activated brain regions to take account of spatial autocorrelation between voxelwise phase shift observations within the activated regions of interest. As an example we apply both tests, and their counterparts assuming spatial independence, to FMRI phase shift data that were acquired from a normal young woman during performance of a periodically designed covert verbal fluency task. We conclude that it is necessary to take account of spatial autocovariance between voxelwise FMRI time series parameter estimates such as the phase shift, and that the most promising way of achieving this is by modelling the spatial autocorrelation structure from a suitably defined base region of the image slice.     

Estimation for the autoregressive moving average process observed with noise

SUMMARY The autoregressive moving average process ARMA (p,q) observed with noise has another ARMA (p,k) representation, where k = max (p,q). Parameters for the ARMA (p,k) representation satisfy some non-linear restrictions. We develop restricted Newton-Raphson estimators of the ARMA (p,k) process which takes advantage of the information given in the non-linear restrictions. The asymptotic relative efficiency of the estimators indicates that the proposed restricted Newton-Raphson estimator is more efficient than the unrestricted Newton-Raphson estimator. In a Monte Carlo experiment, the proposed estimator is shown to perform better than the unrestricted estimator of the ARMA (p,k) process.     

A Spatio-Temporal Model for Functional Magnetic Resonance Imaging Data – with a View to Resting State Networks

Abstract.  Functional magnetic resonance imaging (fMRI) is a technique for studying the active human brain. During the fMRI experiment, a sequence of MR images is obtained, where the brain is represented as a set of voxels. The data obtained are a realization of a complex spatio-temporal process with many sources of variation, both biological and technical. We present a spatio-temporal point process model approach for fMRI data where the temporal and spatial activation are modelled simultaneously. It is possible to analyse other characteristics of the data than just the locations of active brain regions, such as the interaction between the active regions. We discuss both classical statistical inference and Bayesian inference in the model. We analyse simulated data without repeated stimuli both for location of the activated regions and for interactions between the activated regions. An example of analysis of fMRI data, using this approach, is presented.     

A New Test for Short Memory in Long Memory Time Series

This article considers short memory characteristics in a long memory process. We derive new asymptotic results for the sample autocorrelation difference ratios. We used these results to develop a new portmanteau test that determines if short memory parameters are statistically significant. In simulations, the new test can detect short memory components more often than the Ljung-Box test when these short memory components are in fact within a long memory process. Interestingly, our test finds short memory autocorrelations in U.S. inflation rate data, whereas the Ljung-Box test fails to find these autocorrelations. Modeling these short memory autocorrelations of the inflation rate data leads to improved model accuracy and more precise prediction.     

Comparison of data analysis procedures for real-time nanoparticle sampling data using classical regression and ARIMA models

Real-time monitoring is necessary for nanoparticle exposure assessment to characterize the exposure profile, but the data produced are autocorrelated. This study was conducted to compare three statistical methods used to analyze data, which constitute autocorrelated time series, and to investigate the effect of averaging time on the reduction of the autocorrelation using field data. First-order autoregressive (AR(1)) and autoregressive-integrated moving average (ARIMA) models are alternative methods that remove autocorrelation. The classical regression method was compared with AR(1) and ARIMA. Three data sets were used. Scanning mobility particle sizer data were used. We compared the results of regression, AR(1), and ARIMA with averaging times of 1, 5, and 10?min. AR(1) and ARIMA models had similar capacities to adjust autocorrelation of real-time data. Because of the non-stationary of real-time monitoring data, the ARIMA was more appropriate. When using the AR(1), transformation into stationary data was necessary. There was no difference with a longer averaging time. This study suggests that the ARIMA model could be used to process real-time monitoring data especially for non-stationary data, and averaging time setting is flexible depending on the data interval required to capture the effects of processes for occupational and environmental nano measurements.     

The effect of Phase I sample size on the run length performance of control charts for autocorrelated data

Traditional control charts assume independence of observations obtained from the monitored process. However, if the observations are autocorrelated, these charts often do not perform as intended by the design requirements. Recently, several control charts have been proposed to deal with autocorrelated observations. The residual chart, modified Shewhart chart, EWMAST chart, and ARMA chart are such charts widely used for monitoring the occurrence of assignable causes in a process when the process exhibits inherent autocorrelation. Besides autocorrelation, one other issue is the unknown values of true process parameters to be used in the control chart design, which are often estimated from a reference sample of in-control observations. Performances of the above-mentioned control charts for autocorrelated processes are significantly affected by the sample size used in a Phase I study to estimate the control chart parameters. In this study, we investigate the effect of Phase I sample size on the run length performance of these four charts for monitoring the changes in the mean of an autocorrelated process, namely an AR(1) process. A discussion of the practical implications of the results and suggestions on the sample size requirements for effective process monitoring are provided.     

Tracking multiple moving objects in images using Markov Chain Monte Carlo

A new Bayesian state and parameter learning algorithm for multiple target tracking models with image observations are proposed. Specifically, a Markov chain Monte Carlo algorithm is designed to sample from the posterior distribution of the unknown time-varying number of targets, their birth, death times and states as well as the model parameters, which constitutes the complete solution to the specific tracking problem we consider. The conventional approach is to pre-process the images to extract point observations and then perform tracking, i.e. infer the target trajectories. We model the image generation process directly to avoid any potential loss of information when extracting point observations using a pre-processing step that is decoupled from the inference algorithm. Numerical examples show that our algorithm has improved tracking performance over commonly used techniques, for both synthetic examples and real florescent microscopy data, especially in the case of dim targets with overlapping illuminated regions.     

The variability of rainfall acidity

Two years of rainfall acidity data for the eastern United States were analyzed. The data consist of rainfall-event pH measurements from a nine station monitoring network. A spatio-temporal stochastic model, including deterministic components for seasonal variation and rainfall washout, and stochastic components for spatial, temporal, and measurement variation, was fitted to the data. The fitted autocorrelation structure from this model was used, in the process known as Kriging, to obtain BLUE contour maps of seasonal and rainfall adjusted yearly average pH over the monitoring region.     

On a mixture autoregressive model

We generalize the Gaussian mixture transition distribution (GMTD) model introduced by Le and co-workers to the mixture autoregressive (MAR) model for the modelling of non-linear time series. The models consist of a mixture of K stationary or non-stationary AR components. The advantages of the MAR model over the GMTD model include a more full range of shape changing predictive distributions and the ability to handle cycles and conditional heteroscedasticity in the time series. The stationarity conditions and autocorrelation function are derived. The estimation is easily done via a simple EM algorithm and the model selection problem is addressed. The shape changing feature of the conditional distributions makes these models capable of modelling time series with multimodal conditional distributions and with heteroscedasticity. The models are applied to two real data sets and compared with other competing models. The MAR models appear to capture features of the data better than other competing models do.     

Place of living and health inequality: a study for elderly Italians

The aim of this study is to explore if the context matters in explaining socioeconomic inequality in the self-rated health of Italian elderly people. Our hypothesis is that health status perception is associated with existing huge imbalances among Italian areas. A multilevel approach is applied to account for the natural hierarchical structure, as individuals nested in geographical regions. Multilevel logistic regression models are performed including both individual and contextual variables, using data from 2005 Italian Health survey. We prove that individual factors (compositional effect), even representing the most important correlates of health, do not completely explain intra-regional heterogeneity, confirming the existence of an autonomous contextual effect. These territorial differences are present among both Regions and large areas, two geographical aggregations relevant in the domain of health. Moreover, for some Regions, the account for contextual factors explains variations in perceived health, leading to an overthrow of the initial situation: these Regions perform better than expected in the field of health. For other Regions, the contextual elements introduced do not catch the milieu heterogeneity. In this regard, we expect, and solicit, a major effort toward data availability, qualitatively and quantitatively, that might help in explaining residual territorial heterogeneity in health perception, a fundamental starting point for targeting specific policy interventions.     

High dimensional variable selection with clustered data: an application of random multivariate survival forests for detection of outlier medical device components

In many medical studies patients are nested or clustered within doctor. With many explanatory variables, variable selection with clustered data can be challenging. We propose a method for variable selection based on random forest that addresses clustered data through stratified binary splits. Our motivating example involves the detection orthopedic device components from a large pool of candidates, where each patient belongs to a surgeon. Simulations compare the performance of survival forests grown using the stratified logrank statistic to conventional and robust logrank statistics, as well as a method to select variables using a threshold value based on a variable's empirical null distribution. The stratified logrank test performs superior to conventional and robust methods when data are generated to have cluster-specific effects, and when cluster sizes are sufficiently large, perform comparably to the splitting alternatives in the absence of cluster-specific effects. Thresholding was effective at distinguishing between important and unimportant variables.     

Clustering microarray data using model-based double K-means

The microarray technology allows the measurement of expression levels of thousands of genes simultaneously. The dimension and complexity of gene expression data obtained by microarrays create challenging data analysis and management problems ranging from the analysis of images produced by microarray experiments to biological interpretation of results. Therefore, statistical and computational approaches are beginning to assume a substantial position within the molecular biology area. We consider the problem of simultaneously clustering genes and tissue samples (in general conditions) of a microarray data set. This can be useful for revealing groups of genes involved in the same molecular process as well as groups of conditions where this process takes place. The need of finding a subset of genes and tissue samples defining a homogeneous block had led to the application of double clustering techniques on gene expression data. Here, we focus on an extension of standard K-means to simultaneously cluster observations and features of a data matrix, namely double K-means introduced by Vichi (2000). We introduce this model in a probabilistic framework and discuss the advantages of using this approach. We also develop a coordinate ascent algorithm and test its performance via simulation studies and real data set. Finally, we validate the results obtained on the real data set by building resampling confidence intervals for block centroids.