Adaptive testing in arch models

Specification tests for conditional heteroskedasticity that are derived under the assumption that the density of the innovation is Gaussian may not be powerful in light of the recent empirical results that the density is not Gaussian. We obtain specification tests for conditional heteroskedasticity under the assumption that the innovation density is a member of a general family of densities. Our test statistics maximize asymptotic local power and weighted average power criteria for the general family of densities. We establish both first-order and second-order theory for our procedures. Simulations indicate that asymptotic power gains are achievable in finite samples.     

A Two Sample Test for Mean Vectors with Unequal Covariance Matrices

In this paper, we consider testing the equality of two mean vectors with unequal covariance matrices. In the case of equal covariance matrices, we can use Hotelling’s T² statistic, which follows the F distribution under the null hypothesis. Meanwhile, in the case of unequal covariance matrices, the T² type test statistic does not follow the F distribution, and it is also difficult to derive the exact distribution. In this paper, we propose an approximate solution to the problem by adjusting the degrees of freedom of the F distribution. Asymptotic expansions up to the term of order N^{? 2} for the first and second moments of the U statistic are given, where N is the total sample size minus two. A new approximate degrees of freedom and its bias correction are obtained. Finally, numerical comparison is presented by a Monte Carlo simulation.     

A Comparison of Spatiotemporal Surveillance Methods for Nonhomogeneous Change Size

Spatiotemporal surveillance, especially in detection of emerging outbreaks is of particular importance. When an outbreak spreads across some areas, the incidence rate at the center of the outbreak area might be expected to be much higher than the rate at its edge. However, to the best of our knowledge, all existing methods assume a uniformly increasing rate across the entire area of the outbreak. The purpose of this study is to compare the performance of the spatiotemporal surveillance methods such as multivariate cumulative sum (MCUSUM) or multivariate exponentially weighted moving average (MEWMA) when the changes in size are nonhomogeneous. Monte Carlo simulations were conducted to examine the properties of these spatiotemporal surveillance methods and compared them in terms of the detection speed and the identification rate under various scenarios. The results showed that when nonhomogeneous change sizes are involved, the MCUSUM method taking into account spatial nonhomogeneity of increase rates yields a better identification than the method ignoring such change size pattern although the detection speeds are similar. Further, a case study for the detection of male thyroid cancer data in New Mexico in the United States was performed to demonstrate the applicability of these methods.     

Patterns of consent: evidence from a general household survey

Summary.  We analyse patterns of consent and consent bias in the context of a large general household survey, the 'Improving survey measurement of income and employment' survey, also addressing issues that arise when there are multiple consent questions. A multivariate probit regression model for four binary outcomes with two incidental truncations is used. We show that there are biases in consent to data linkage with benefit and tax credit administrative records that are held by the Department for Work and Pensions, and with wage and employment data held by employers. There are also biases in respondents' willingness and ability to supply their national insurance number. The biases differ according to the question that is considered. We also show that modelling questions on consent independently rather than jointly may lead to misleading inferences about consent bias. A positive correlation between unobservable individual factors affecting consent to Department for Work and Pensions record linkage and consent to employer record linkage is suggestive of a latent individual consent propensity.     

Case-deletion Influence Measures for the Data from Multivariate t Distributions

For the data from multivariate t distributions, it is very hard to make an influence analysis based on the probability density function since its expression is intractable. In this paper, we present a technique for influence analysis based on the mixture distribution and EM algorithm. In fact, the multivariate t distribution can be considered as a particular Gaussian mixture by introducing the weights from the Gamma distribution. We treat the weights as the missing data and develop the influence analysis for the data from multivariate t distributions based on the conditional expectation of the complete-data log-likelihood function in the EM algorithm. Several case-deletion measures are proposed for detecting influential observations from multivariate t distributions. Two numerical examples are given to illustrate our methodology.     

Forecasting composite indicators with anticipated information: an application to the industrial production index

Summary. Many economic and social phenomena are measured by composite indicators computed as weighted averages of a set of elementary time series. Often data are collected by means of large sample surveys, and processing takes a long time, whereas the values of some elementary component series may be available a considerable time before the others and may be used for forecasting the composite index. This problem is addressed within the framework of prediction theory for stochastic processes. A method is proposed for exploiting anticipated information to minimize the mean-square forecast error, and for selecting the most useful elementary series. An application to the Italian general industrial production index is illustrated, which demonstrates that knowledge of anticipated values of some, or even just one, component series may reduce the forecast error considerably.     

A test for the presence of conditional heteroskedasticity within arch-m framework

This paper is concerned with testing the presence of ARCH within the ARCH-M model as the alternative hypothesis. Standard testing procedures are inapplicable since a nuisance parameter is unidentified under the null hypothesis. Nonetheless, the diagnostic tests for the presence of the conditional variance is very important since any misspecification in the conditional variance equation leads to inconsistent estimates of the conditional mean parameters. BTo resolve the problem of unidentified nuisance parameter, 'Ne apply Davies' approach, and investigate its finite sample performance through a Monte Carlo study.     

A multivariate finite mixture latent trajectory model with application to dementia studies

Dementia patients exhibit considerable heterogeneity in individual trajectories of cognitive decline, with some patients showing rapid decline following diagnoses while others exhibiting slower decline or remaining stable for several years. Dementia studies often collect longitudinal measures of multiple neuropsychological tests aimed to measure patients’ decline across a number of cognitive domains. We propose a multivariate finite mixture latent trajectory model to identify distinct longitudinal patterns of cognitive decline simultaneously in multiple cognitive domains, each of which is measured by multiple neuropsychological tests. EM algorithm is used for parameter estimation and posterior probabilities are used to predict latent class membership. We present results of a simulation study demonstrating adequate performance of our proposed approach and apply our model to the Uniform Data Set from the National Alzheimer's Coordinating Center to identify cognitive decline patterns among dementia patients.     

The Combined Model: A Tool for Simulating Correlated Counts with Overdispersion

The combined model as introduced by Molenberghs et al. (2007 Molenberghs, G., Verbeke, G., Demétrio, C. (2007). An extended random-effects approach to modeling repeated, overdispersed count data. Lifetime Data Analysis 13:513–531.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar], 2010 Molenberghs, G., Verbeke, G., Demétrio, C., Vieira, A. (2010). A family of generalized linear models for repeated measures with normal and conjugate random effects. Statistical Science 25:325–347.[Crossref], [Web of Science ®] , [Google Scholar]) has been shown to be an appealing tool for modeling not only correlated or overdispersed data but also for data that exhibit both these features. Unlike techniques available in the literature prior to the combined model, which use a single random-effects vector to capture correlation and/or overdispersion, the combined model allows for the correlation and overdispersion features to be modeled by two sets of random effects. In the context of count data, for example, the combined model naturally reduces to the Poisson-normal model, an instance of the generalized linear mixed model in the absence of overdispersion and it also reduces to the negative-binomial model in the absence of correlation. Here, a Poisson model is specified as the parent distribution of the data conditional on a normally distributed random effect at the subject or cluster level and/or a gamma distribution at observation level. Importantly, the development of the combined model and surrounding derivations have relevance well beyond mere data analysis. It so happens that the combined model can also be used to simulate correlated data. If a researcher is interested in comparing marginal models via Monte Carlo simulations, a necessity to generate suitable correlated count data arises. One option is to induce correlation via random effects but calculation of such quantities as the bias is then not straightforward. Since overdispersion and correlation are simultaneous features of longitudinal count data, the combined model presents an appealing framework for generating data to evaluate statistical properties, through a pre-specification of the desired marginal mean (possibly in terms of the covariates and marginal parameters) and a marginal variance-covariance structure. By comparing the marginal mean and variance of the combined model to the desired or pre-specified marginal mean and variance, respectively, the implied hierarchical parameters and the variance-covariance matrices of the normal and Gamma random effects are then derived from which correlated Poisson data are generated. We explore data generation when a random intercept or random intercept and slope model is specified to induce correlation. The data generator, however, allows for any dimension of the random effects although an increase in the random-effects dimension increases the sensitivity of the derived random effects variance-covariance matrix to deviations from positive-definiteness. A simulation study is carried out for the random-intercept model and for the random intercept and slope model, with or without the normal and Gamma random effects. We also pay specific attention to the case of serial correlation.