---

---

---

Robust linear regression is one of the most popular problems in the robust statistics community. The parameters of this method are often estimated via least trimmed squares, which minimizes the sum of the k smallest squared residuals. So, the estimation method in contrast to the common least squares estimation method is very computationally expensive. The main idea of this paper is to propose a new estimation method in partial linear models based on minimizing the sum of the k smallest squared residuals which determines the set of outlier point and provides robust estimators. In this regard, first, difference based method in estimation parameters of partial linear models is introduced. Then the method of obtaining robust difference based estimators in partial linear models is introduced which is based on solving an optimization problem minimizing the sum of the k smallest squared residuals. This method can identify outliers. The simulated example and applied numerical example with real data found the proposed robust difference based estimators in the paper produce highly accurate results in compare to the common difference based estimators in partial linear models.

---

Skew spatial data often are modeled by using skew Gaussian random field. The main problem is that simulations from this random field are very time consuming for some parameter values and large dimensions. Also it is impossible in some cases and requires using of an approximation methods. One a spatial statistics branch often used to determine the natural resources such as oil and gas, is analysis of seismic data by inverse model. Bayesian Gaussian inversion model commonly is used in seismic inversion that the analytical and computational can easily be done for large dimensions. But in practice, we are encountered with the variables that are asymmetric and skewed. They are modeled using skew distributions. In Bayesian Analysis of closed skew Gaussian inversion model, there is an important problem to generate samples from closed skew normal distributions. In this paper, an efficient algorithm for the realization of the Closed Skew Normal Distribution is provided with higher dimensions. Also the Closed Skew T Distribution is offered that include heavy tails in the density function and the simulation algorithm for generating samples from the Closed Skew T Distribution is provided. Finally, the discussion and conclusions are presented.

---

Hierarchical spatio-temporal models are used for modeling space-time responses and temporally and spatially correlations of the data is considered via Gaussian latent random field with Matérn covariance function. The most important interest in these models is estimation of the model parameters and the latent variables, and is predict of the response variables at new locations and times. In this paper, to analyze these models, the Bayesian approach is presented. Because of the complexity of the posterior distributions and the full conditional distributions of these models and the use of Monte Carlo samples in a Bayesian analysis, the computation time is too long. For solving this problem, Gaussian latent random field with Matern covariance function are represented as a Gaussian Markov Random Field (GMRF) through the Stochastic Partial Differential Equations (SPDE) approach. Approximatin Baysian method and Integrated Nested Laplace Approximation (INLA) are used to obtain an approximation of the posterior distributions and to inference about the model. Finally, the presented methods are applied to a case study on rainfall data observed in the weather stations of Semnan in 2013.

---

‎In many fields such as econometrics‎, ‎psychology‎, ‎social sciences‎, ‎medical sciences‎, ‎engineering‎, ‎etc.‎, ‎we face with multicollinearity among the explanatory variables and the existence of outliers in data‎. ‎In such situations‎, ‎the ordinary least-squares estimator leads to an inaccurate estimate‎. ‎The robust methods are used to handle the outliers‎. ‎Also‎, ‎to overcome multicollinearity ridge estimators are suggested‎. ‎On the other hand‎, ‎when the error terms are heteroscedastic or correlated‎, ‎the generalized least squares method is used‎. ‎In this paper‎, ‎a fast algorithm for computation of the feasible generalized least trimmed squares ridge estimator in a semiparametric regression model is proposed and then‎, ‎the performance of the proposed estimators is examined through a Monte Carlo simulation study and a real data set.

---

---

The panel data model is used in many areas, such as economics, social sciences, medicine, and epidemiology. In recent decades, inference on regression coefficients has been developed in panel data models. In this paper, methods are introduced to test the equality models of the panel model among the groups in the data set. First, we present a random quantity that we estimate its distribution by two methods of approximation and parametric bootstrap. We also introduce a pivotal quantity for performing this hypothesis test. In a simulation study, we compare our proposed approaches with an available method based on the type I error and test power. We also apply our method to gasoline panel data as a real data set.

---

The Gaussian random field is commonly used to analyze spatial data. One of the important features of this random field is having essential properties of the normal distribution family, such as closure under linear transformations, marginalization and conditioning, which makes the marginal consistency condition of the Kolmogorov extension theorem. Similarly, the skew-Gaussian random field is used to model skewed spatial data. Although the skew-normal distribution has many of the properties of the normal distribution, in some definitions of the skew-Gaussian random field, the marginal consistency property is not satisfied. This paper introduces a stationery skew-Gaussian random field, and its marginal consistency property is investigated. Then, the spatial correlation model of this skew random field is analyzed using an empirical variogram. Also, the likelihood analysis of the introduced random field parameters is expressed with a simulation study, and at the end, a discussion and conclusion are presented.

---

‎The high-dimensional data analysis using classical regression approaches is not applicable, and the consequences may need to be more accurate.
This study tried to analyze such data by introducing new and powerful approaches such as support vector regression, functional regression, LASSO and ridge regression. On this subject, by investigating two high-dimensional data sets  (riboflavin and simulated data sets) using the suggested approaches, it is progressed to derive the most efficient model based on three criteria (correlation squared, mean squared error and mean absolute error percentage deviation) according to the type of data.

---

Gaussian random field is usually used to model Gaussian spatial data. In practice, we may encounter non-Gaussian data that are skewed. One solution to model skew spatial data is to use a skew random field. Recently, many skew random fields have been proposed to model this type of data, some of which have problems such as complexity, non-identifiability, and non-stationarity. In this article, a flexible class of closed skew-normal distribution is introduced to construct valid stationary random fields, and some important properties of this class such as identifiability and closedness under marginalization and conditioning are examined. The reasons for developing valid spatial models based on these skew random fields are also explained. Additionally, the identifiability of the spatial correlation model based on empirical variogram is investigated in a simulation study with the stationary skew random field as a competing model. Furthermore, spatial predictions using a likelihood approach are presented on these skew random fields and a simulation study is performed to evaluate the likelihood estimation of their parameters. 

---

The spatial generalized linear mixed models are often used, where the latent variables representing spatial correlations are modeled through a Gaussian random field to model the categorical spatial data. The violation of the Gaussian assumption affects the accuracy of predictions and parameter estimates in these models. In this paper, the spatial generalized linear mixed models are fitted and analyzed by utilizing a stationary skew Gaussian random field and employing an approximate Bayesian approach. The performance of the model and the approximate Bayesian approach is examined through a simulation example, and implementation on an actual data set is presented.

---

Spatial regression models are used to analyze quantitative spatial responses based on linear and non-linear relationships with explanatory variables. Usually, the spatial correlation of responses is modeled with a Gaussian random field based on a multivariate normal distribution. However, in practice, we encounter skewed responses, which are analyzed using skew-normal distributions. Closed skew-normal distribution is one of the extended families of skew-normal distributions, which has similar properties to normal distributions. This article presents a hierarchical Bayesian analysis based on a flexible subclass of closed skew-normal distributions. Given the time-consuming nature of Monte Carlo methods in hierarchical Bayes analysis, we have opted to use the variational Bayes approach to approximate the posterior distribution. This decision was made to expedite the analysis process without compromising the accuracy of our results. Then, the proposed model is implemented and analyzed based on the real earthquake data of Iran.