Bayesian estimation of regression parameters in elliptical measurement error models

The main object of this paper is to discuss the Bayes estimation of the regression coefficients in the elliptically distributed simple regression model with measurement errors. The posterior distribution for the line parameters is obtained in a closed form, considering the following: the ratio of the error variances is known, informative prior distribution for the error variance, and non-informative prior distributions for the regression coefficients and for the incidental parameters. We proved that the posterior distribution of the regression coefficients has at most two real modes. Situations with a single mode are more likely than those with two modes, especially in large samples. The precision of the modal estimators is studied by deriving the Hessian matrix, which although complicated can be computed numerically. The posterior mean is estimated by using the Gibbs sampling algorithm and approximations by normal distributions. The results are applied to a real data set and connections with results in the literature are reported.     

Preliminary test estimation in system regression models in view of asymmetry

In this paper, we consider the system regression model introduced by Arashi and Roozbeh (Comput Stat 30:359–376, 2015) and study the performance of the feasible preliminary test estimator (FPTE) both analytically and computationally, under the assumption that constraints may hold on the vector parameter space. The performance of the FPTE is analysed through a Monte Carlo simulation study under bounded and or asymmetric loss functions. An application of the so-called Cobb–Douglas production function in economic modelling together with the results from the simulation study shows that the bounded linear exponential (BLINEX) loss function outperforms the linear exponential loss function (LINEX) by comparing risk values.     

Efficient MCMC estimation of some elliptical copula regression models through scale mixtures of normals

This paper proposes an efficient estimation method for some elliptical copula regression models by expressing both copula density and marginal density functions as scale mixtures of normals (SMN). Implementing these models using the SMN is novel and allows efficient estimation via Bayesian methods. An innovative algorithm for the case of complex semicontinuous margins is also presented. We utilize the facts that copulas are invariant to the location and scale of the margins; all elliptical distributions have the same correlation structure; and some densities can be represented by the SMN. Two simulation studies, one on continuous margins and the other on semicontinuous margins, highlight the favorable performance of the proposed methods. Two empirical studies, one on the US excess returns and one on the Thai wage earnings, further illustrate the applicability of the proposals.     

Combined regression models

The goal of this paper is two-fold. First, new regression models obtained by combinations of the least squares (LS), minimax (MM), and the least sum of absolute deviations (LSAD) are proposed. Second, measures for assessing the influence of observations on the fitted models are suggested. The paper is interdisciplinary because the theory behind the proposed method draws from results in the operations research area. The methods are illustrated by their application to some examples and graphical illustrations are given.     

Simple homogeneous models

Geometrical stability theory is a powerful set of model-theoretic tools that can lead to structural results on models of a simple first-order theory. Typical results offer a characterization of the groups definable in a model of the theory. The work is carried out in a universal domain of the theory (a saturated model) in which the Stone space topology on ultrafilters of definable relations is compact. Here we operate in the more general setting of homogeneous models, which typically have noncompact Stone topologies. A structure equipped with a class of finitary relations is strongly -homogeneous if orbits under automorphisms of have finite character in the following sense: Given an ordinal and sequences , from , if and have the same orbit, for all and , then for some automorphism of . In this paper strongly -homogeneous models in which the elements of induce a symmetric and transitive notion of independence with bounded character are studied. This notion of independence, defined using a combinatorial condition called ``dividing', agrees with forking independence when is saturated. A concept central to the development of stability theory for saturated structures, namely parallelism, is also shown to be well-behaved in this setting. These results broaden the scope of the methods of geometrical stability theory.

     

Seemingly unrelated regression models

The cross-covariance matrix of observation vectors in two linear statistical models need not be zero matrix. In such a case the problem is to find explicit expressions for the best linear unbiased estimators of both model parameters and estimators of variance components in the simplest structure of the covariance matrix. Univariate and multivariate forms of linear models are dealt with.     

Multivariate spatial regression models

This paper describes the inference procedures required to perform Bayesian inference to some multivariate econometric models. These models have a spatial component built into commonly used multivariate models. In particular, the common component models are addressed and extended to accommodate for spatial dependence. Inference procedures are based on a variety of simulation-based schemes designed to obtain samples from the posterior distribution of model parameters. They are also used to provide a basis to forecast new observations.     

Unitary response regression models

The dependent variable in a regular linear regression is a numerical variable, and in a logistic regression it is a binary or categorical variable. In these models the dependent variable has varying values. However, there are problems yielding an identity output of a constant value which can also be modelled in a linear or logistic regression with this constant in place of a numerical or binary response. In a linear model with a positive response, dividing by its values yields a regression of constant output by the relative shares of individual predictors into the total response. Chemical reaction models use the agents' concentration, summing to a constant 100%. Another example can be found in priority modelling by Thurstone scaling for ranked or paired comparison data. The Thurstone scale can be estimated by probit or logit models with identical output across all the responses. Models with a unitary output can be constructed by software for regular regressions, but they give a different interpretation of results. For instance, the coefficient of multiple determination is not an estimate of the explained variance in the total response variance (which is zero), but a measure of the fitting quality of the constant approximated by an aggregate of predictors.     

Bayesian inference for dependent elliptical measurement error models

In this article we provide a Bayesian analysis for dependent elliptical measurement error models considering nondifferential and differential errors. In both cases we compute posterior distributions for structural parameters by using squared radial prior distributions for the precision parameters. The main result is that the posterior distribution of location parameters, for specific priors, is invariant with respect to changes in the generator function, in agreement with previous results obtained in the literature under different assumptions. Finally, although the results obtained are valid for any elliptical distribution for the error term, we illustrate those results by using the student-t distribution and a real data set.     

Random weightingT-statistics in linear regression models

In this paper, we have constructed a random weighting statistic to approximate the distribution of studentized least square estimator in a linear regression model with ideal accuracyo(n ^–1/2). Thus, we have provided a more practical distribution approximating method.Supported by the Doctoral Program Foundation of the Institute of Higher Education and the National Natural Science Foundation of China.     

CLT in functional linear regression models

We propose in this work to derive a CLT in the functional linear regression model. The main difficulty is due to the fact that estimation of the functional parameter leads to a kind of ill-posed inverse problem. We consider estimators that belong to a large class of regularizing methods and we first show that, contrary to the multivariate case, it is not possible to state a CLT in the topology of the considered functional space. However, we show that we can get a CLT for the weak topology under mild hypotheses and in particular without assuming any strong assumptions on the decay of the eigenvalues of the covariance operator. Rates of convergence depend on the smoothness of the functional coefficient and on the point in which the prediction is made.     

Simple discrete-time malarial models

We investigate dynamics of mosquito population models under two assumptions, respectively, and then formulate simple discrete-time compartmental susceptible-exposed-infective-recovered models for the malaria transmission based on the mosquito population models. We show that the mosquito population models either have robust dynamics or exhibit period-doubling bifurcation depending on the model assumptions. We derive a formula for the reproductive number of infection for the malaria model, which determines the stability of the infection-free fixed point. We then determine the existence of endemic fixed points for the malaria models. Using numerical simulations, we demonstrate that the dynamical characteristics of the mosquito populations, such as the global stability of the endemic fixed point and the appearance of a period-doubling bifurcation, are reflected in the dynamics of the malaria transmission.     

Simple stochastic models of cracks

The stochastic Markov models describe various natural and technological processes. They are often used in rather diverse areas. The Markov models are distinguished with discrete time and relatively few states. In concrete cases, such models allow us to carry out effective calculations. In the article are considered some special models with two, four and eight states. The processes are simulated that are connected with hydraulic fracturing and wave front development.     

二值回归模型中自变量误差的处理

In this paper the binary regression model, in which the response Y is binary, i.e. Y takes the value 0 or 1 only, is considered. The conditional probability P{Y=1|X=x} is assumed to be of the form G(a₀+a₁x), where x is the measurement of the factor, a= (a₀,a₁)′ is the vector of the unknown parameter being estimated and G(a) = (1+e^a)^-1. The measurement x is assu-med to be contaminated when the value of x is out of certain limit. M-estimation is employed to estimate the unknown a of the model. In this paper, the consistent M-estimation is proved to be asymptotically normal and the optimal M-equation is obtained. Finally, a computation method is introduced to solve the M-equation and to get the optimal solution.     

Heteroscedasticity checks for regression models

For checking on heteroscedasticity in regression models, a unified approach is proposed to constructing test statistics in parametric and nonparametric regression models. For nonparametric regression, the test is not affected sensitively by the choice of smoothing parameters which are involved in estimation of the nonparametric regression function. The limiting null distribution of the test statistic remains the same in a wide range of the smoothing parameters. When the covariate is one-dimensional, the tests are, under some conditions, asymptotically distribution-free. In the high-dimensional cases, the validity of bootstrap approximations is investigated. It is shown that a variant of the wild bootstrap is consistent while the classical bootstrap is not in the general case, but is applicable if some extra assumption on conditional variance of the squared error is imposed. A simulation study is performed to provide evidence of how the tests work and compare with tests that have appeared in the literature. The approach may readily be extended to handle partial linear, and linear autoregressive models.     

Semiparametric random coefficient regression models

Linear regression models with random coefficients express the idea that each individual sampled may have a different linear response function. Technically speaking, random coefficient regression encompasses a rich variety of submodels. These include deconvolution or affine-mixture models as well as certain classical linear regression models that have heteroscedastic errors, or errors-in-variables, or random effects. This paper studies minimum distance estimates for the coefficient distributions in a general, semiparametric, random coefficient regression model. The analysis yields goodness-of-fit tests for the semiparametric model, prediction regions for future responses, and confidence regions for the distribution of the random coefficients.This research was supported in part by NSF Grant DMS 9001710.     

Multiple deletion diagnostics in beta regression models

We consider the problem of identifying multiple outliers in a general class of beta regression models proposed by Ferrari and Cribari-Neto (J Appl Stat 31:799–815, 2004). The currently available single-case deletion diagnostic measures, e.g., the standardized weighted residual (SWR), the Cook-like distance (LD), etc., often fail to identify multiple outlying observations, because they suffer from the well-known problems of masking and swamping effects. In this article, we develop group deletion diagnostic measures, such as generalized SWR, generalized LD, generalized DFFITS and generalized DFBETAS, and suggest a simple procedure for identifying multiple outliers using these. The performance of the proposed methods is investigated through simulation studies and two practical examples.     

Testing model assumptions in functional regression models

In the functional regression model where the responses are curves, new tests for the functional form of the regression and the variance function are proposed, which are based on a stochastic process estimating L²-distances. Our approach avoids the explicit estimation of the functional regression and it is shown that normalized versions of the proposed test statistics converge weakly. The finite sample properties of the tests are illustrated by means of a small simulation study. It is also demonstrated that for small samples, bootstrap versions of the tests improve the quality of the approximation of the nominal level.     

Testing in linear composite quantile regression models

Composite quantile regression (CQR) can be more efficient and sometimes arbitrarily more efficient than least squares for non-normal random errors, and almost as efficient for normal random errors. Based on CQR, we propose a test method to deal with the testing problem of the parameter in the linear regression models. The critical values of the test statistic can be obtained by the random weighting method without estimating the nuisance parameters. A distinguished feature of the proposed method is that the approximation is valid even the null hypothesis is not true and power evaluation is possible under the local alternatives. Extensive simulations are reported, showing that the proposed method works well in practical settings. The proposed methods are also applied to a data set from a walking behavior survey.