多个正态总体均值相等的VDR检验

本文将随机估计由一维参数扩展至多维参数,基于随机估计的密度函数提出VDR检验.在总体方差已知和未知的两种情形下,本文讨论多个正态总体均值是否相同的VDR检验过程,而且得到精确的检验.单因素方差分析是VDR检验的特例.模拟研究表明,VDR检验是一个普遍适用的方法.     

����~���Ƶ�ָ����ط����Ի��ЧӦģ�͵����Լ���

Homogeneity of variance and correlation coefficients is one of assumptions in the analysis of longitudinal data.However, the assumption can be challenged. In this paper, we mainly propose and analyze nonlinear mixed effects models for longitudinal data with exponential correlation covariance structure, intend to introduce Huber's function in the log likelihood function and get robust estimation (M-estimation) by Fisher scoring method. Score test statistics for homogeneity of variance and correlation coefficient based on M-estimation are then studied. A simulation study is carried to assess the performance of test statistics and the method we proposed in the paper is illustrated by an actual data example.     

具有双线性BL(1,1,1,1)误差的非线性回归模型相关性和方差齐性的检验

在时间序列回归模型分析中,相关性和方差齐性的检验是一个很基本的问题.本文讨论了具有双线性BL(1,1,1,1)误差的非线性回归模型的相关性和方差齐性的检验问题, 用Score检验方法给出了双线性项检验、相关性检验、方差齐性检验、以及相关性和方差齐性同时检验的检验统计量.推广和发展了具有线性序列误差项回归模型的结果.本文还用数值实例说明了检验方法的实用价值.     

随机估计和VDR检验

众所周知统计推断有三种理论:普遍承认的Neyman理论(频率学派),Bayes推断和信仰推断(Fiducial)。Bayes推断基于后验分布,由先验分布和样本分布求得。信仰推断是基于信仰分布(Confidence Distribution,简称CD),直接利用样本求得。两者推断方式一致,都是用分布函数作推断,称为分布推断。从分析传统的参数估计、假设检验特性来看,经典统计推断也可以视为分布推断。通常将置信上限看做置信度的函数。其反函数,即置信度是置信上界的函数,恰是分布函数,该分布恰是近年来引起许多学者兴趣的CD。在本文中,基于随机化估计(其分布是一CD)的概率密度函数,提出VDR检验。常见正态分布期望或方差的检验,多元正态分布期望的Hoteling检验等是其特例。VDR(vertical density representation)检验适合于多元分布参数检验,实现了非正态的多元线性变换分布族的参数检验。VDR构造的参数的置信域有最小Lebesgue测度。     

Empirical likelihood inference for median regression models for censored survival data

Recent advances in median regression model have made it possible to use this model for analyzing a variety of censored survival data. For inference on the model parameter vector, there are now semiparametric procedures based on normal approximation that are valid without strong conditions on the error distribution. However, the accuracy of such procedures can be quite low when the censoring proportion is high. In this paper, we propose an alternative semiparametric procedure based on the empirical likelihood. We define the empirical likelihood ratio for the parameter vector and show that its limiting distribution is a weighted sum of chi-square distributions. Numerical results from a simulation study suggest that the empirical likelihood method is more accurate than the normal approximation based method of Ying et al. (J. Amer. Statist. Assoc. 90 (1995) 178).     

Auxiliary Likelihood-Based Approximate Bayesian Computation in State Space Models

A computationally simple approach to inference in state space models is proposed, using approximate Bayesian computation (ABC). ABC avoids evaluation of an intractable likelihood by matching summary statistics for the observed data with statistics computed from data simulated from the true process, based on parameter draws from the prior. Draws that produce a “match” between observed and simulated summaries are retained, and used to estimate the inaccessible posterior. With no reduction to a low-dimensional set ofsufficient statistics being possible in the state space setting, we define the summaries as the maximum of an auxiliary likelihood function, and thereby exploit the asymptotic sufficiency of this estimator for the auxiliary parameter vector. We derive conditions under which this approach—including a computationally efficient version based on the auxiliary score—achieves Bayesian consistency. To reduce the well-documented inaccuracy of ABC in multiparameter settings, we propose the separate treatment of each parameter dimension using an integrated likelihood technique. Three stochastic volatility models for which exact Bayesian inference is either computationally challenging, or infeasible, are used for illustration. We demonstrate that our approach compares favorably against an extensive set of approximate and exact comparators. An empirical illustration completes the article. Supplementary materials for this article are available online.     

Multimodel inference based on smoothed information criteria

The multimodel inference makes statistical inferences from a set of plausible models rather than from a single model. In this paper, we focus on the multimodel inference based on smoothed information criteria proposed by seminal monographs(see Buckland et al.(1997) and Burnham and Anderson(2003)), which are termed as smoothed Akaike information criterion(SAIC) and smoothed Bayesian information criterion(SBIC)methods. Due to their simplicity and applicability, these methods are very widely used in many fields. By using an illustrative example and deriving limiting properties for the weights in the linear regression, we find that the existing variance estimation for SAIC is not applicable because of a restrictive condition, but for SBIC it is applicable. Especially, we propose a simulation-based inference for SAIC based on the limiting properties. Both the simulation study and the real data example show the promising performance of the proposed simulationbased inference.     

Efficient statistical inference for partially nonlinear errors-in-variables models

In this paper, we consider the partially nonlinear errors-in-variables models when the non- parametric component is measured with additive error. The profile nonlinear least squares estimator of unknown parameter and the estimator of nonparametric component are constructed, and their asymptotic properties are derived under general assumptions. Finite sample performances of the proposed statistical inference procedures are illustrated by Monte Carlo simulation studies.     

Testing statistical hypotheses based on the density power divergence

The family of density power divergences is an useful class which generates robust parameter estimates with high efficiency. None of these divergences require any non-parametric density estimate to carry out the inference procedure. However, these divergences have so far not been used effectively in robust testing of hypotheses. In this paper, we develop tests of hypotheses based on this family of divergences. The asymptotic variances of the estimators are generally different from the inverse of the Fisher information matrix, so that the usual drop-in-divergence type statistics do not lead to standard Chi-square limits. It is shown that the alternative test statistics proposed herein have asymptotic limits which are described by linear combinations of Chi-square statistics. Extensive simulation results are presented to substantiate the theory developed.     

Completeness and unbiased estimation of mean vector in the multivariate group sequential case

We consider estimation after a group sequential test about a multivariate normal mean, such as a χ² test or a sequential version of the Bonferroni procedure. We derive the density function of the sufficient statistics and show that the sample mean remains to be the maximum likelihood estimator but is no longer unbiased. We propose an alternative Rao-Blackwell type unbiased estimator. We show that the family of distributions of the sufficient statistic is not complete, and there exist infinitely many unbiased estimators of the mean vector and none has uniformly minimum variance. However, when restricted to truncation-adaptable statistics, completeness holds and the Rao-Blackwell estimator has uniformly minimum variance.     

Structured probabilistic inference

Probabilistic inference is among the main topics with reasoning in uncertainty in AI. For this purpose, Bayesian Networks (BNs) is one of the most successful and efficient Probabilistic Graphical Model (PGM) so far. Since the mid-90s, a growing number of BNs extensions have been proposed. Object-oriented, entity-relationship and first-order logic are the main representation paradigms used to extend BNs. While entity-relationship and first-order models have been successfully used for machine learning in defining lifted probabilistic inference, object-oriented models have been mostly underused. Structured inference, which exploits the structural knowledge encoded in an object-oriented PGM, is a surprisingly unstudied technique. In this paper we propose a full object-oriented framework for PRM and propose two extensions of the state-of-the-art structured inference algorithm: SPI which removes the major flaws of existing algorithms and SPISBB which largely enhances SPI by using d-separation.     

Heavy tailed capital incomes: Zenga index,statistical inference,and ECHP data analysis

Micro-data of European Union (EU) countries show that capital incomes account for a large part of disparity in populations and follow heavy-tailed distributions in many EU countries. Measuring and comparing the disparity requires incorporating the relative nature of ‘small’ and ‘large,’ and for this reason we employ the newly developed Zenga index of economic inequality. Its non-parametric estimator does not fall into any well known class of statistics. This makes the development of statistical inference a challenge even for light-tailed populations, let alone heavy-tailed ones, as is the case with capital incomes. In this paper we construct a heavy-tailed Zenga estimator, establish its asymptotic distribution, and derive confidence intervals. We explore the performance of the confidence intervals in a simulation study and draw conclusions about capital incomes in EU countries, based on the 2001 wave of the European Community Household Panel (ECHP) survey.     

Bias-corrected statistical inference for partially linear varying coefficient errors-in-variables models with restricted condition

In this paper, we consider the statistical inference for the partially liner varying coefficient model with measurement error in the nonparametric part when some prior information about the parametric part is available. The prior information is expressed in the form of exact linear restrictions. Two types of local bias-corrected restricted profile least squares estimators of the parametric component and nonparametric component are conducted, and their asymptotic properties are also studied under some regularity conditions. Moreover, we compare the efficiency of the two kinds of parameter estimators under the criterion of Lo?ner ordering. Finally, we develop a linear hypothesis test for the parametric component. Some simulation studies are conducted to examine the finite sample performance for the proposed method. A real dataset is analyzed for illustration.     

Causal statistical inference in high dimensions

We present a short selective review of causal inference from observational data, with a particular emphasis on the high-dimensional scenario where the number of measured variables may be much larger than sample size. Despite major identifiability problems, making causal inference from observational data very ill-posed, we outline a methodology providing useful bounds for causal effects. Furthermore, we discuss open problems in optimization, non-linear estimation and for assigning statistical measures of uncertainty, and we illustrate the benefits and limitations of high-dimensional causal inference for biological applications.     

Conceptual issues in understanding the inner logic of statistical inference: Insights from two teaching experiments

We report on a sequence of two classroom teaching experiments that investigated high school students’ understandings as they explored connections among the ideas comprising the inner logic of statistical inference—ideas involving a core image of sampling as a repeatable process, and the organization of its outcomes into a distribution of sample statistics as a basis for making inferences. Students’ responses to post-instruction test questions indicate that despite understanding various individual components of inference—a sample, a population, and a distribution of a sample statistic—their abilities to coordinate and compose these into a coherent and well-connected scheme of ideas were usually tenuous. We argue that the coordination and composition required to assemble these component ideas into a coherent scheme is a major source of difficulty in developing a deep understanding of inference.     

Bayesian Conditional Density Filtering

We propose a conditional density filtering (C-DF) algorithm for efficient online Bayesian inference. C-DF adapts MCMC sampling to the online setting, sampling from approximations to conditional posterior distributions obtained by propagating surrogate conditional sufficient statistics (a function of data and parameter estimates) as new data arrive. These quantities eliminate the need to store or process the entire dataset simultaneously and offer a number of desirable features. Often, these include a reduction in memory requirements and runtime and improved mixing, along with state-of-the-art parameter inference and prediction. These improvements are demonstrated through several illustrative examples including an application to high dimensional compressed regression. In the cases where dimension of the model parameter does not grow with time, we also establish sufficient conditions under which C-DF samples converge to the target posterior distribution asymptotically as sampling proceeds and more data arrive. Supplementary materials of C-DF are available online.     

��ָ��ģ���췽���

??In the additive regression models, the single-index model is considered commonly for high dimensional regression analysis. The specification of this model that it is more flexible compared with a parametric model, and it avoids the curse of dimensionality because the single-index reduces the dimensionality of a standard variable vector (x in the multi-regression) to a univariate index (\beta^\T X in the single-index model). In this paper, we developed a single-index regression model with a functional errors' term that serves in checking the heteroscedasticity. Since the efficient inference of a regression model demands that heteroscedasticity is regarded when it exists, this paper presents the assumptions of testing variance constancy in single-index models. The test statistic is assessing the variance homogeneity stated as a combination of Levene's test and the theories of ANOVA for the infinite factor levels. The test statistic in the simulation studies displays appropriately in all situations compared to a well-known method and applies to a real dataset.     

Testing the Heteroscedasticity in Single-Index Model

In the additive regression models, the single-index model is considered commonly for high dimensional regression analysis. The specification of this model that it is more flexible compared with a parametric model, and it avoids the curse of dimensionality because the single-index reduces the dimensionality of a standard variable vector (x in the multi-regression) to a univariate index (\beta^\T X in the single-index model). In this paper, we developed a single-index regression model with a functional errors' term that serves in checking the heteroscedasticity. Since the efficient inference of a regression model demands that heteroscedasticity is regarded when it exists, this paper presents the assumptions of testing variance constancy in single-index models. The test statistic is assessing the variance homogeneity stated as a combination of Levene's test and the theories of ANOVA for the infinite factor levels. The test statistic in the simulation studies displays appropriately in all situations compared to a well-known method and applies to a real dataset.     

��ɾʧ�����·�λ���ع�Ĺ⻬������Ȼ����

This paper is focused on testing the parameters of the quantile regression models. For complete observation, it is shown in literature that the test statistics, based on empirical likelihood (EL) method and smoothed empirical likelihood (SEL) method, both converge weakly to the standard Chi-square distribution $\chi_M^2$ under the null hypothesis. For right censored data, the statistics in literature, by the EL method, have a weighted Chi-square limiting distribution, but the weights are unknown. In this paper, we show that the statistics based on the EL method and the SEL method also converge weakly to $\chi_M^2$ under the null hypothesis, so there is no need to estimate any weights. As its estimating function is smoothed, the SEL method can be Bartlett corrected. Numerical results show that the SEL method, via Bartlett correction, outperforms some recent methods.     

Detection of outliers in weighted least squares regression

In multiple linear regression model, we have presupposed assumptions (independence, normality, variance homogeneity and so on) on error term. When case weights are given because of variance heterogeneity, we can estimate efficiently regression parameter using weighted least squares estimator. Unfortunately, this estimator is sensitive to outliers like ordinary least squares estimator. Thus, in this paper, we proposed some statistics for detection of outliers in weighted least squares regression.