---

In this paper we propose a new two-parameters distribution, which is an extension of the Lindley distribution with increasing and bathtub-shaped failure rate, called as the Lindley-logarithmic (LL) distribution. The new distribution is obtained by compounding Lindley (L) and Logarithmic distributions. We obtain several properties of the new distribution such as its probability density function, its failure rate functions, quantiles and moments. The maximum likelihood estimation procedure via a EM-algorithm is presented in this paper. At the end, in order to show the flexibility and potentiality of this new class, some series of real data is used to fit.

---

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.

---

The optimal criteria are used to find the optimal design in the studied model. These kinds of models are included the paired comparison models. In these models, the optimal criteria (D-optimality) determine the optimal paired comparison. In this paper, in addition to introducing the quadratic regression model with random effects, the paired comparison models were presented and the optimal design has been calculated for them.

---

In weighted sampling as a generalization of random sampling, every observation, y, is recorded with probably proportional to a non-negative function of y. In this paper, the normal regression model is investigated under the weighted sampling for a common weight function. Parameters of the model are estimated for known and unknown weight parameters. Using simulation, efficiency of estimators is studied when they have not closed forms. As an application, the data of number of visited  patients by specialist doctors in Social Security Organization of Ahvaz in Iran (SSOAI) are analyzed.

---

Statistical analysis of fractional Brownian motion process is one of the most important issues in the field of stochastic processes. The most important issue in the study of this process is statistical inference about the Hurst parametersof the fractional Brownian motion. One of the methods for estimation of aforementioned parameter is maximum likelihood approach. Due to the computational complexity of this approach to give a closed estimate, it is attempting to derive the parameter estimated through the numerical method approach. Also, the theoretical result of the paper is evaluated in a simulation study for different scenarios.

---

Value-at-Risk and Average Value-at-Risk are tow important risk measures based on statistical methoeds that used to measure the market's risk with quantity structure. Recently, linear regression models such as least squares and quantile methods are introduced to estimate these risk measures. In this paper, these two risk measures are estimated by using omposite quantile regression. To evaluate the performance of the proposed model with the other models, a simulation study was conducted and at the end, applications to real data set from Iran's stock market are illustarted.

---

In this paper, an approximate tolerance interval is presented for the discrete size-biased Poisson-Lindley distribution. This approximate tolerance interval, is constructed based on large sample Wald confidence interval for the parameter of the size-biased Poisson-Lindley distribution. Then, coverage probabilities and expected widths of the proposed tolerance interval is considered. The results show that the coverage probabilities have a better performance for the small values of the parameter and are close to the nominal confidence level, and are conservative for the large values of the parameter. Finally, an applicable example is provided for illustrating approximate tolerance interval.

---

Among all statistical distributions, standard normal distribution has been the most important and practical distribution in which calculation of area under probability density function and cumulative distribution function are required. Unfortunately, the cumulative distribution function of this is, in general, expressed as a definite integral with no closed form or analytical solution. Consequently, it has to be approximated. In this paper, attempts have been made for Winitzki's approximation to be proved by a new approach. Then, the approximation is improved with some modifications and shown that the maximum error resulted from this is less than 0.0000584. Finally, an inverse function for computation of normal distribution quantiles has been derived.

---

In some applied problems we need to choose a population from the given populations and estimate the parameter of the selected population. Suppose k random samples are chosen from k populations with proportional hazard rate model or proportional reversed hazard rate model. According to a specified selection rule, it is desired to estimate the parameter of the best (worst) selected population. In this paper, under the entropy loss function we obtain the  uniformly minimum risk unbiased (UMRU) estimator of  the parameters of the selected population, and derived sufficient conditions for minimaxity of a given estimator. Then we find the class of admissible and inadmissible linear estimators of the parameters of the selected population and determine the class of dominators of a given estimator. We show that the UMRU estimator is inadmissible and compare the obtained estimators by plotting their risk functions.

---

The presence of endogenous variables in the statistical models leads to inconsistent and bias estimators for the parameters. In this case, several approaches have been proposed which are able to tackle the biase and inconsistency problems only in large sample situations. One of these methods is biased on instrumental variables which causes removing endogenous variables. The method of two-stage least squares is another approach in this case that it has more accurate than ordinary least squares. This paper aims to enhance the accuracy of three methods of estimation based upon least square methodology called, two-stage iterative least squares, two-stage Jackknife least squares and also two-stage calibration least squares. In order to evaluate the performance of each method, a simulation study is conducted. Also, using data collected in 1390 related to the cost and revenue in Iran, those methods to estimate parameters are compared.

---

Consider a repairable system which starts operating at t=0. Once the system fails, it is immediately replaced by another one of the same type or it is repaired and back to its working functions. In this paper, the system's activity is studied from t>0 for a fixed period of time w. Different replacement policies are considered. In each cases, for a fixed period of time w, the probability model and likelihood function of repair process, say window censored, are obtained. The obtained results depend on the lifetime distribution of the original system, so, expression for the maximum likelihood estimator and Fisher information are derived, by assuming the lifetime follows an exponential distribution.

---

The collective risk model of insurance company with constant initial capital when process of claims number have the poisson distribution with constant rate is considered. For computing the infinite time ruin probability the stochastic processes and differential equations are used. Also a formula is obtained to compute the Lundberg approximation in finding the approximate of infinite time ruin probability based on the distribution function of claims number. The numerical examples to illustrate these results are given and showed that for any value of initial capital the approximate of our infinite time ruin probability is closer to its real value rather than the ruin probability computed by other authors and has less error.

---

Maximum of the Renyi entropy and the Tsallis entropy are generalization of the maximum entropy for a larger class of Shannon entropy. In this paper we introduce the maximum Renyi entropy and some of the attributes of distributions which have maximum Renyi entropy investigated. The form of distributions with maximum Renyi entropy is power so we state some properties of these distributions and we have a new form of the Renyi entropy. After pointing the topics of minimum Renyi divergence, some other points in this relation have been discussed. An another form of Renyi divergence have also obtained. Therefore we discussed some of the economic applications of the maximum entropy. Meanwhile, the review of the Csiszar information measure, the general form of distributions with minimum Renyi divergence have obtained.

---

We consider the purely sequential procedure for estimating the scale parameter of an exponential distribution, when the risk function is bounded by the known preassigned number. In this paper, we provide explicit formulas for the expectation of the total sample size. Also, we propose how to adjust the stopping variable so that the risk is uniformly bounded by a known preassigned number. In the end, the performances of the proposed methodology are investigated with the help of simulations.

---

In this paper the issue of variable selection with new approach in finite mixture of semi-parametric regression models is studying, although it is supposed that data have Poisson distribution. When we use Poisson distribution, two problems such as overdispersion and excess zeros will happen that can affect on variable selection and parameter estimation. Actually parameter estimation in parametric component of the semi-parametric regression model is done by penalized likelihood approach. However, in nonparametric component after local approximation using Teylor series, the estimation of nonparametric coefficients along with estimated parametric coefficients will be calculated. Using new approach leads to a properly variable selection results. In addition to representing related theories, overdispersion and excess zeros are considered in data simulation section and using EM algorithm in parameter estimation leads to increase the accuracy of end results.

---

The (n-k+1)-out-of-n systems are important types of coherent systems and have many applications in various areas of engineering. In this paper, the general inactivity time of failed components of (n-k+1)-out-of-n system is studied when the system fails at time t>0. First we consider a parallel system including two exchangeable components and then using Farlie-Gumbel-Morgenstern copula, investigate the behavior of mean inactivity time of failed components of the system. In the next part, (n-k+1)-out-of-n systems with exchangeable components are considered and then, some stochastic ordering properties of the general inactivity time of the systems are presented based on one sample or two samples.

---

Penalized estimators for estimating regression parameters have been considered by many authors for many decades. Penalized regression with rectangular norm is one of the mainly used since it does variable selection and estimating parameters, simultaneously. In this paper, we propose some new estimators by employing uncertain prior information on parameters. Superiority of the proposed shrinkage estimators over the least absoluate and shrinkage operator (LASSO) estimator is demonstrated via a Monte Carlo study. The prediction rate of the proposed estimators compared to the LASSO estimator is also studied in the US State Facts and Figures dataset.

---

Usually, we estimate the unknown parameter by observing a random sample and using the usual methods of estimation such as maximum likelihood method. In some situations, we have information about the real parameter in the form of a guess. In these cases, one may shrink the maximum likelihood or other estimators towards a guess value and construct a shrinkage estimator. In this paper, we study the behavior of a Bayes shrinkage estimator for the scale parameter of exponential distribution based on censored samples under an asymmetric and scale invariant loss function. To do this, we propose a Bayes shrinkage estimator and compute the relative efficiency between this estimator and the best linear estimator within a subclass with respect to sample size, hyperparameters of the prior distribution and the vicinity of the guess and real parameter. Also, the obtained results are extended to Weibull and Rayleigh lifetime distributions.

---

In this paper, we introduce a family of bivariate generalized Gompertz-power series distributions. This new class of bivariate distributions contains several models such as: bivariate generalized Gompertz -geometric, -Poisson, - binomial, -logarithmic, -negative binomial and bivariate generalized exponental-power series distributions as special cases. We express the method of construction and derive different properties of the proposed class of distributions. The method of maximum likelihood and EM algorithm are used for estimating the model parameters. Finally, we illustrate the usefulness of the new distributions by means of application to real data sets.