Bayesian analysis of threshold autoregressions

A nonasymptotic Bayesian approach is developed for analysis of data from threshold autoregressive processes with two regimes. Using the conditional likelihood function, the marginal posterior distribution for each of the parameters is derived along with posterior means and variances. A test for linear functions of the autoregressive coefficients is presented. The approach presented uses a posterior p-value averaged over the values of the threshold. The one-step ahead predictive distribution is derived along with the predictive mean and variance. In addition, equivalent results are derived conditional upon a value of the threshold. A numerical example is presented to illustrate the approach.     

On trends and constants in periodic autoregressions

Periodic autoregressions are characterised by autoregressive structures that vary with the season. If a time series is periodically integrated, one needs a seasonally varying differencing filter to remove the stochastic trend. When the periodic regression model contains constants and trends with unrestricted parameters, the data can show diverging seasonal deterministic trends. In this paper we derive explicit expressions for parameter restrictions that result in common deterministic trends under periodic trend stationarity and periodic integration.     

On trends and constants in periodic autoregressions

Periodic autoregressions are characterised by autoregressive structures that vary with the season. If a time series is periodically integrated, one needs a seasonally varying differencing filter to remove the stochastic trend. When the periodic regression model contains constants and trends with unrestricted parameters, the data can show diverging seasonal deterministic trends. In this paper we derive explicit expressions for parameter restrictions that result in common deterministic trends under periodic trend stationarity and periodic integration.     

Bootstrap confidence regions computed from autoregressions of arbitrary order

Given a linear time series, e.g. an autoregression of infinite order, we may construct a finite order approximation and use that as the basis for confidence regions. The sieve or autoregressive bootstrap, as this method is often called, is generally seen as a competitor with the better-understood block bootstrap approach. However, in the present paper we argue that, for linear time series, the sieve bootstrap has significantly better performance than blocking methods and offers a wider range of opportunities. In particular, since it does not corrupt second-order properties then it may be used in a double-bootstrap form, with the second bootstrap application being employed to calibrate a basic percentile method confidence interval. This approach confers second-order accuracy without the need to estimate variance. That offers substantial benefits, since variances of statistics based on time series can be difficult to estimate reliably, and—partly because of the relatively small amount of information contained in a dependent process—are notorious for causing problems when used to Studentize. Other advantages of the sieve bootstrap include considerably greater robustness against variations in the choice of the tuning parameter, here equal to the autoregressive order, and the fact that, in contradistinction to the case of the block bootstrap, the percentile t version of the sieve bootstrap may be based on the 'raw' estimator of standard error. In the process of establishing these properties we show that the sieve bootstrap is second order correct.     

Estimating multivariate autoregressive moving average models by fitting long autoregressions

Some simple methods for the estimation of mixed multivariate autoregressive moving average time series models are introduced. The methods require the fitting of a long autoregression to the data and the computation of consistent initial estimates for the parameters of the model. After these preliminaries the estimators of the paper are obtained by applying weighted least squares to a multivariate auxiliary regression model. Two types of weight matrices are considered. Both of them yield estimators which are strongly consistent and asymptotically normally distributed. The first estimators are also asymptotically efficient while the second ones are not fully efficient but computationally simple. A simulation study is performed to illustrate the behaviour of the estimators in finite samples.     

Modelling daily multivariate pollutant data at multiple sites

Summary. This paper considers the spatiotemporal modelling of four pollutants measured daily at eight monitoring sites in London over a 4-year period. Such multiple-pollutant data sets measured over time at multiple sites within a region of interest are typical. Here, the modelling was carried out to provide the exposure for a study investigating the health effects of air pollution. Alternative objectives include the design problem of the positioning of a new monitoring site, or for regulatory purposes to determine whether environmental standards are being met. In general, analyses are hampered by missing data due, for example, to a particular pollutant not being measured at a site, a monitor being inactive by design (e.g. a 6-day monitoring schedule) or because of an unreliable or faulty monitor. Data of this type are modelled here within a dynamic linear modelling framework, in which the dependences across time, space and pollutants are exploited. Throughout the approach is Bayesian, with implementation via Markov chain Monte Carlo sampling.     

Demographic policy in the light of the resolutions of the Twenty-sixth Congress of the Communist Party

The USSR is moving from extensive to intensive type of population reproduction, not because of any governmental measures but as a reaction to objective circumstances. At present there is a reduction in the population cohort of working age. This has an adverse effect on production. And the number of the elderly is increasing; we must involve them more in the production process, in monitoring the service sphere, for instance. The infant mortality rate has declined sharply in the USSR, especially since World War II. The birthrate has also dropped, but it is still higher than in the US, France, or West Germany. Some authorities think that this (the ZPG) is a good idea, but this author thinks it is a bad idea, meaning a drop in society's productive forces. We need an active demographic policy, meaning one which ensures optimal reproduction of the population with an improvement in its health, culture, and living standard. The policy must be differentiated, as conditions vary from one part of the country to another. An example is the 1981 decree on assistance to families with children which provides for payments of 35-50 rubles/month to mothers who must stay home to look after children. The number of children per family to ensure reproduction must be 2.6 on the average. In the Central Asian republics family sizes of 4-5 children are anticipated, but in the RSFSR, Ukraine, Belorussia, and the Baltic republics the anticipated size of families is rather low. Measures must be taken to alter this by providing assistance to families with children.     

流动人口社会融合指标体系内在关系研究

利用2014年流动人口社会融合专题调查数据,分别以经济融合、社区融合、文化接纳、自我认同中任何一个变量为因变量,剩余三个变量为自变量,构建由经济模型、社会模型、文化模型和自我认同模型组成的联立方程组,并采用多元线性回归和logistic回归为主要方法,引入流动性别、年龄、受教育程度、户口性质、流动区域、民族等六个变量为控制变量,对流动人口社会融合的四个维度的内在关系进行了分析。研究发现,经济收入与社区融合统计关系不显著,与文化接纳和身份认同统计关系显著,表明收入越高文化接纳越好,身份认同却越差;经济收入对心理文化及社区参与的影响要略微强于心理文化及社区参与对经济收入的影响,经济融合有一定的独立性特征,而社区参与、文化接纳与身份认同存在更高的一致性关系。     

Detecting multiple confounders

This paper proposes an approach for detecting multiple confounders which combines the advantages of two causal models, the potential outcome model and the causal diagram. The approach need not use a complete causal diagram as long as it is known that a known covariate set Z$Z$ contains the parent set of the exposure E  . On the other hand, whether a covariate is or not a confounder may depend on its categorization. We introduce uniform non-confounding which implies non-confounding in any subpopulation defined by the interval of a covariate (or any pooled level for a discrete covariate). We show that the conditions in Miettinen and Cook's criteria for non-confounding also imply uniform non-confounding. Further we present an algorithm for deleting non-confounders from the potential confounder set Z$Z$, which extends Greenland et al.'s [1999a. Causal diagrams for epidemiologic research. Epidemiology 10, 37–48] approach by splitting Z$Z$ into a series of potential confounder subsets. We also discuss conditions for non-confounding bias in the subpopulations in which we are interested, where the subpopulations may be defined by non-confounders.     

Testing multiple slippages

A class of invariant Bayes rules is derived for testing homogeneity of k (≥2) different populations against (^k_t) slippage alternatives that some (unknown) subset of size t of the given populations has parameter larger than the remaining k-t, where t is a given integer between 1 and k-1. For a similar problem in nonparametric situations, locally best tests based on ranks are derived.     

MIMCA: multiple imputation for categorical variables with multiple correspondence analysis

We propose a multiple imputation method to deal with incomplete categorical data. This method imputes the missing entries using the principal component method dedicated to categorical data: multiple correspondence analysis (MCA). The uncertainty concerning the parameters of the imputation model is reflected using a non-parametric bootstrap. Multiple imputation using MCA (MIMCA) requires estimating a small number of parameters due to the dimensionality reduction property of MCA. It allows the user to impute a large range of data sets. In particular, a high number of categories per variable, a high number of variables or a small number of individuals are not an issue for MIMCA. Through a simulation study based on real data sets, the method is assessed and compared to the reference methods (multiple imputation using the loglinear model, multiple imputation by logistic regressions) as well to the latest works on the topic (multiple imputation by random forests or by the Dirichlet process mixture of products of multinomial distributions model). The proposed method provides a good point estimate of the parameters of the analysis model considered, such as the coefficients of a main effects logistic regression model, and a reliable estimate of the variability of the estimators. In addition, MIMCA has the great advantage that it is substantially less time consuming on data sets of high dimensions than the other multiple imputation methods.     

Nonparametric multiple testing procedures

Assume that we have n_i independent observations from each of k independent populations. Each population has the same distribution except for a translation parameter. We are interested in specific pairwise differences of the parameters in various settings, such as treatment vs. control, change point or all pairwise differences. We propose new multiple testing procedures for the pairwise differences. The new procedures are based on ranks and they have desirable practical properties not shared by existing procedures. These include tests that satisfy the interval property. Furthermore, the test method provides an interval that serves as an estimate of the difference in the parameters of interest.     

Rank-based multiple change-point detection

Abstract

A nonparametric procedure is proposed to estimate multiple change-points of location changes in a univariate data sequence by using ranks instead of the raw data. While existing rank-based multiple change-point detection methods are mostly based on sequential tests, we treat it as a model selection problem. We derive the corresponding Schwarz’s information criterion for rank-statistics, theoretically prove the consistency of the change-point estimator and use a pruned dynamic programing algorithm to achieve the change-point estimator. Simulation studies show our method’s robustness, effectiveness and efficiency in detecting mean-changes. We also apply the method to a gene dataset as an illustration.     

Bootstrap multiple test procedure

A method for multiple comparisons is presented. It is based on the bootstrap and takes care of the logical structure among the hypotheses to be tested. The test procedure is very general according to number of hypotheses, sample distribution and number of observations.     

Inferring multiple graphical structures

Gaussian Graphical Models provide a convenient framework for representing dependencies between variables. Recently, this tool has received a high interest for the discovery of biological networks. The literature focuses on the case where a single network is inferred from a set of measurements. But, as wetlab data is typically scarce, several assays, where the experimental conditions affect interactions, are usually merged to infer a single network. In this paper, we propose two approaches for estimating multiple related graphs, by rendering the closeness assumption into an empirical prior or group penalties. We provide quantitative results demonstrating the benefits of the proposed approaches. The methods presented in this paper are embeded in the R package simone from version 1.0-0 and later.     

A revision of the tukey multiple comparisons procedure to control the probability of committing at most one type i error

Halperin et al. (1988) suggested an approach which allows for k Type I errors while using Scheffe's method of multiple comparisons for linear combinations of p means. In this paper we apply the same type of error control to Tukey's method of multiple pairwise comparisons. In fact, the variant of the Tukey (1953) approach discussed here defines the error control objective as assuring with a specified probability that at most one out of the p(p-l)/2 comparisons between all pairs of the treatment means is significant in two-sided tests when an overall null hypothesis (all p means are equal) is true or, from a confidence interval point of view, that at most one of a set of simultaneous confidence intervals for all of the pairwise differences of the treatment means is incorrect. The formulae which yield the critical values needed to carry out this new procedure are derived and the critical values are tabulated. A Monte Carlo study was conducted and several tables are presented to demonstrate the experimentwise Type I error rates and the gains in power furnished by the proposed procedure     

Intrinsic Bayesian inference on a Poisson rate and on the ratio of two Poisson rates

Our purpose is to explore the intrinsic Bayesian inference on the rate of a Poisson distribution and on the ratio of the rates of two independent Poisson distributions, with the natural conjugate family of priors in the first case and the semi-conjugate family of priors defined by Laurent and Legrand (2011) in the second case. Intrinsic Bayesian inference is derived from the Bayesian decision theory framework based on the intrinsic discrepancy loss function. We cover in particular the case of some objective Bayesian procedures suggested by Bernardo when considering reference priors.     

Joint impact of multiple observations on a subset of variables in multiple linear regression

In multiple linear regression analysis, each observation affects the fitted regression equation differently and has varying influences on the regression coefficients of the different variables. Chatterjee & Hadi (1988) have proposed some measures such as DSSE_ij (Impact on Residual Sum of Squares of simultaneously omitting the ith observation and the jth variable), F_j (Partial F-test for the jth variable) and F_j(i) (Partial F-test for the jth variable omitting the ith observation) to show the joint impact and the interrelationship that exists among a variable and an observation. In this paper we have proposed more extended form of those measures DSSE_IJ, F_J and F_J(I) to deal with the interrelationships that exist among the multiple observations and a subset of variables by monitoring the effects of the simultaneous omission of multiple variables and multiple observations.     

Fitting a multiple regression function

Consider the p-dimensional unit cube [0,1]^p, p≥1. Partition [0, 1]^p into n regions, R_1,n,…,R_n,n such that the volume Δ(R_j,n) is of order n^?1,j=1,…,n. Select and fix a point in each of these regions so that we have x⁽ⁿ⁾₁,…,x⁽ⁿ⁾_n. Suppose that associated with the j-th predictor vector x⁽ⁿ⁾_j there is an observable variable Y⁽ⁿ⁾_j, j=1,…,n, satisfying the multiple regression model $\text{Y}{}^{\mathrm{(n)}}{}_{\mathrm{j}}\text{=g(}\text{x}{}^{\mathrm{(n)}}{}_{\mathrm{j}}\text{)+e}{}^{\mathrm{(n)}}{}_{\mathrm{j}}$, where g is an unknown function defined on [0, 1]^pand {e⁽ⁿ⁾_j} are independent identically distributed random variables with Ee⁽ⁿ⁾₁=0 and Var e⁽ⁿ⁾₁=σ²<∞. This paper proposes as an estimator of g(x), where k(u) is a known p-dimensional bounded density and {a_n} is a sequence of reals converging to 0 asn→∞. Weak and strong consistency of g_n(x) and rates of convergence are obtained. Asymptoticnormality of the estimator is established. Also proposed is $\text{σ}{}^2{}_{\mathrm{n}}\text{=n}{}^{\mathrm{?1}}\text{Σ}{}^{\mathrm{n}}{}_{\mathrm{j=1}}\text{(Y}{}^{\mathrm{(n)}}{}_{\mathrm{j}}\text{?g}{}_{\mathrm{n}}\text{(}\text{x}{}^{\mathrm{(n)}}{}_{\mathrm{j}}\text{))}{}^2$ as a consistent estimate of σ².