A Bootstrap Algorithm for Data from a Periodic Multiplicative Intensity Function

The periodic multiplicative intensity model is considered. A new bootstrap method for non stationary counting processes which intensity function has some periodicity properties is presented. Its main advantage is that it does not destroy the temporal order and the original periodicity of the underlying counting process. The proposed algorithm is used to construct a bootstrap version of the maximum likelihood hazard function estimator. The consistency of the bootstrap method is shown. A possible modification of the proposed bootstrap method is discussed. The bootstrap simultaneous confidence intervals for the hazard function are presented. The telecommunication network traffic real data example is discussed.     

A bayesian analysis of the mixed linear model

The mixed model is defined. The exact posterior distribution for the fixed effect vector is obtained. The exact posterior distribution for the error variance is obtained. The exact posterior mean and variance of a Bayesian estimator for the variances of random effects is also derived. All computations are non-iterative and avoid numerical integrations.     

中国货币政策产业非均衡效应实证研究

 本文采集中国1999—2011年季度数据,构建由不同产业GDP、外汇储备、CPI、不同货币政策中介目标组成的四元SVAR模型系列,测算出不同货币政策冲击对不同产业的具体效应。结果发现：（1）综合渠道冲击的作用力度最大,时滞最短;信贷渠道冲击的作用力度和时滞适中;利率渠道冲击的作用力度最小,时滞最长。（2）在信贷传导渠道乃至综合传导渠道中,第二产业的反应力度最大,第一产业次之,第三产业最小;而在利率传导渠道中,第二产业的反应力度最大,第三产业次之,第一产业最小。（3）不同产业对不同货币政策冲击反应时滞的长短差别较为一致。第三产业最短,第二产业次之,第一产业最长。上述结论对相关政策操作和管理具有重要的启发意义。     

Bayesian hierarchical models for linear networks

The purpose of this study is to highlight dangerous motorways via estimating the intensity of accidents and study its pattern across the UK motorway network. Two methods have been developed to achieve this aim. First, the motorway-specific intensity is estimated by using a homogeneous Poisson process. The heterogeneity across motorways is incorporated using two-level hierarchical models. The data structure is multilevel since each motorway consists of junctions that are joined by grouped segments. In the second method, the segment-specific intensity is estimated. The homogeneous Poisson process is used to model accident data within grouped segments but heterogeneity across grouped segments is incorporated using three-level hierarchical models. A Bayesian method via Markov Chain Monte Carlo is used to estimate the unknown parameters in the models and the sensitivity to the choice of priors is assessed. The performance of the proposed models is evaluated by a simulation study and an application to traffic accidents in 2016 on the UK motorway network. The deviance information criterion (DIC) and the widely applicable information criterion (WAIC) are employed to choose between models.     

Quantile regression in functional linear semiparametric model

This paper proposes nonparametric estimation methods for functional linear semiparametric quantile regression, where the conditional quantile of the scalar responses is modelled by both scalar and functional covariates and an additional unknown nonparametric function term. The slope function is estimated using the functional principal component basis and the nonparametric function is approximated by a piecewise polynomial function. The asymptotic distribution of the estimators of slope parameters is derived and the global convergence rate of the quantile estimator of unknown slope function is established under suitable norm. The asymptotic distribution of the estimator of the unknown nonparametric function is also established. Simulation studies are conducted to investigate the finite-sample performance of the proposed estimators. The proposed methodology is demonstrated by analysing a real data from ADHD-200 sample.     

A practical method to calculate probabilities: illustrative example from the electronic industry business

The real-life environment is made of probabilistic data by nature and the ability to make decisions based on probabilities is crucial in the business world. It is common to have a set of data and the need of calculating the probability of taking a value greater or less than a specific value. It is also common in many companies the unavailability of a statistical software or a specialized professional in statistics. The purpose of this paper is to present a practical and simple method to calculate probabilities from normal or non-normal distributed data set and illustrate it with an application from the electronic industry. The method does not demand statistical knowledge from the user; there is no need of normality assumptions, goodness test or transformations. The proposed method is easy to implement, robust and the experiments have evidenced its quality. The technique is validated with a large variety of instances and compared with the well-known Johnson system of distributions.     

A simple and asymptotically optimal test for the equality of k(≥2) exponential distributions based on type II censored samples

A test for the equality of two or more two-parameter exponential distributions is suggested. It is developed on an intuitive basis and is obtained by combining two independent tests by the Fisher method (1950, pp. 99-101). The test is simple for application and is optimal asymptotically in the sense of Bahadur efficiency (1960). A numerical example is discussed to illustrate its application in a real-world situation. The Monte Carlo simulation is used for calculating its power which is compared with that of the test suggested by Singh and Narayan (1983). The suggested test is found oftener more powerful.     

Estimates of regression coefficients based on the sign covariance matrix

Summary. A new estimator of the regression parameters is introduced in a multivariate multiple-regression model in which both the vector of explanatory variables and the vector of response variables are assumed to be random. The affine equivariant estimate matrix is constructed using the sign covariance matrix (SCM) where the sign concept is based on Oja's criterion function. The influence function and asymptotic theory are developed to consider robustness and limiting efficiencies of the SCM regression estimate. The estimate is shown to be consistent with a limiting multinormal distribution. The influence function, as a function of the length of the contamination vector, is shown to be linear in elliptic cases; for the least squares (LS) estimate it is quadratic. The asymptotic relative efficiencies with respect to the LS estimate are given in the multivariate normal as well as the t -distribution cases. The SCM regression estimate is highly efficient in the multivariate normal case and, for heavy-tailed distributions, it performs better than the LS estimate. Simulations are used to consider finite sample efficiencies with similar results. The theory is illustrated with an example.     

双重广义线性模型在车损险费率厘定中的应用

双重广义线模型是对广义线性模型的扩展,其对反应变量的均值与散度参数同时建立模型,提高了模型运用的灵活性与适应性。将双重广义线性模型应用到车损险费率厘定中,既考虑了费率期望值与费率因子之间的关系,又考虑了变量的分散程度与费率因子之间的关系,并以欧洲一家保险公司的汽车保险损失数据为样本进行实证研究,把无索赔优待等级、地区、车型与年均行驶里程数作为费率因子,建立了费率厘定模型。结果表明,所得到费率结构合理,符合实际。     

Testimator of the scale parameter of the exponential distribution using linex loss function

The present paper investigates the properties of a testimator of scale of an exponential distribution under Linex loss function. The risk function of testimator is derived and compared with that of an admissible estimator relative to Linex loss function. The shrinkage testimator is proposed which is the extension of testimator and its properties have been discussed. The level of significance of testimator is decided on the basis of Akaike information criterion following Hirano (1977, 1978). It is found that the testimator and shrinkage testimator dominates the admissible estimator in terms of risk in certain parametric space.     

A multiple model discrimination procedure

The case of selecting between a set of fixed models is considered. The true model is assumed to be contained in the set of proposed models and errors are taken to be normally distributed. A sequential procedure which yeilds probabilities of incorrect selections is proposed. The procedure is shown to have optimal properties and is extended to the estimated model case by a bootstrap procedure.     

Testing of Suspected Observations in an Exponential Sample With Unknown Origin

The problem of multiple upper outliers in two-parameter exponential sample is considered. A test statistic is proposed to identify the outliers at the upper end of the sample. The null distribution of the test statistic is obtained and the critical values are found. The performance of the test is also compared with the earlier work.     

A nonparametric test for bivariate circular symmetry based on the empirical cdf

A nonparametric test for circular symmetry about 0 in a continuous bivariate distribution is proposed. The test is of the von Mises type, based on the empirical cdf of the sample, expressed in polar co-ordinates. However, the test is independent of the choice of the polar axis. The asymptotic form of the test statistic is obtained by considering the weak convergence of the empirical process to a limiting Gaussian process. The asymptotic distribution of the test statistic is found explicitly, both under the null hypothesis and under simple alternatives. The test is shown to be consistent against all alternatives.     

A new nonparametric bivariate test for two sample location problem

A strictly nonparametric bivariate test for two sample location problem is proposed. The proposed test is easy to apply and does not require the stringent condition of affine-symmetry or elliptical symmetry which is required by some of the major tests available for the same problem. The power function of the proposed test is calculated. The asymptotic distribution of the proposed test statistic is found to be normal. The power of proposed test is compared with some of the well-known tests under various distributions using Monte Carlo simulation technique. The power study shows that the proposed test statistic performs better than most of the test statistics for almost all the distributions considered here. As soon as the underlying population structure deviates from normality, the ability of the proposed test statistic to detect the smallest shift in location increases as compared to its competitors. The application of the test is shown by using a data set.     

Book Reviews

This article uses a Bayesian unit-root test in stochastic volatility models. The time series of interest is the volatility that is unobservable. The unit-root testing is based on the posterior odds ratio, which is approximated by Markov-chain Monte Carlo methods. Simulations show that the testing procedure is efficient for moderate sample size. The unit-root hypothesis is rejected in seven market indexes, and some evidence of nonstationarity is observed in the TWSI of Taiwan.     

Bayesian Analysis of the Prototypal Search Model

Bayesian analysis for a simple but widely applied dynamic programming model is obtained. The setting is the prototypal job-search model. The general case of wage and duration data, with potential censoring, is studied. The optimality condition implied by the dynamic programming setup is fully imposed. The posterior distribution reveals a “ridge” reflecting the characteristic nonstandard nature of the inference problem. Marginal distributions and moments are obtained in a canonical parameterization after a suitable approximation. The adequacy of the approximation is easily assessed. Simulation is applied to study alternative parameterizations and prior robustness and to facilitate prior elicitations. Finally, we illustrate the applicability of our methods by giving posterior distributions for the elasticities of unemployment durations and reemployment wages with respect to unemployment income. Our analysis is easy to implement and all computations are simple to perform.     

Attribute-variable Inspection Policy for Lots Using Resampling Based on EWMA

A new mixed sampling plan which is a combination of the attribute single sampling plan and variables resampling scheme based on EWMA statistic is proposed in this paper. The operating characteristic function of the proposed plan is derived and the plan parameters are determined such that the probability of acceptance of good lot is larger than the specified producer's confidence level and the bad lot acceptance probability is smaller than the consumer's confidence level. The efficiency and the advantages of the proposed plan are discussed over the existing attribute sampling plan. The extensive tables are provided for industrial applications. The use of tables is discussed with the help of a real-time industrial example.     

A multiple warm standby δ-shock system with a repairman having multiple vacations

In this article, a warm standby n-unit system is studied. The system is operational as long as there is one unit normal. The unit online, which has a lifetime distribution governed by a phase-type distribution, is also attacked by a shock from some external causes. Assume that shocks arrive according to a Poisson process. Whenever an interarrival time of shock is less than a threshold, the unit online fails. The lifetimes of the units in warm standby is exponentially distributed. A repairman who can take multiple vacations repairs the failed units based on the “first-in-first-out” rule. The repair times and the vacation times of repairman are governed by different phase-type distributions. For this system, the Markov process governing the system is constructed. The system is studied in a transient and stationary regime; the availability, the reliability, the rates of occurrence of the different types of failures, and the working probability of the repairman are calculated. A numerical application is performed to illustrate the calculations.     

Extreme values identification in regression using a peaks-over-threshold approach

The problem of heavy tail in regression models is studied. It is proposed that regression models are estimated by a standard procedure and a statistical check for heavy tail using residuals is conducted as a tool for regression diagnostic. Using the peaks-over-threshold approach, the generalized Pareto distribution quantifies the degree of heavy tail by the extreme value index. The number of excesses is determined by means of an innovative threshold model which partitions the random sample into extreme values and ordinary values. The overall decision on a significant heavy tail is justified by both a statistical test and a quantile–quantile plot. The usefulness of the approach includes justification of goodness of fit of the estimated regression model and quantification of the occurrence of extremal events. The proposed methodology is supplemented by surface ozone level in the city center of Leeds.