An Easily Implemented Generalised Inverse Gaussian Generator

Using the ‘ratio’ method an easily implemented algorithm is derived for the generalised inverse Gaussian distribution. Computer timings and efficiency calculations show that the procedure is fast over a wide range of distribution parameter values.     

EGARCH Model with Weighted Liquidity

We analyze a variant of the EGARCH model which captures the variation of the intra-day price. We study the asymptotic behavior of the estimators for the parameters of the model. We also illustrate our theoretical results by empirical studies.     

Reliability of global sensitivity indices

Uncertainty and sensitivity analysis is an essential ingredient of model development and applications. For many uncertainty and sensitivity analysis techniques, sensitivity indices are calculated based on a relatively large sample to measure the importance of parameters in their contributions to uncertainties in model outputs. To statistically compare their importance, it is necessary that uncertainty and sensitivity analysis techniques provide standard errors of estimated sensitivity indices. In this paper, a delta method is used to analytically approximate standard errors of estimated sensitivity indices for a popular sensitivity analysis method, the Fourier amplitude sensitivity test (FAST). Standard errors estimated based on the delta method were compared with those estimated based on 20 sample replicates. We found that the delta method can provide a good approximation for the standard errors of both first-order and higher-order sensitivity indices. Finally, based on the standard error approximation, we also proposed a method to determine a minimum sample size to achieve the desired estimation precision for a specified sensitivity index. The standard error estimation method presented in this paper can make the FAST analysis computationally much more efficient for complex models.     

Lattice path combinatorics for multiple product identities

Lattice paths are enumerated as walks on a lattice under the Delannoy criterion of vertical, horizontal and upward diagonal steps. Delannoy recursions are generalized for arbitrary weights in each of these directions, where the row elements of the Delannoy triangles correspond to lattice paths for all walks at a given length.     

Response predictions in regressions on panel data

In this paper we consider the problem of obtaining best linear unbiased estimators of individual response coefficients in a Random Coefficient Linear Regression (RCR} Model, comparing alternative estimators for these response vectors through a simulation study. We also provide an empirical example that illustrates the estimation procedure proposed here.     

Semiparametric Analysis of Network Formation

ABSTRACT

We consider a statistical model for directed network formation that features both node-specific parameters that capture degree heterogeneity and common parameters that reflect homophily among nodes. The goal is to perform statistical inference on the homophily parameters while treating the node-specific parameters as fixed effects. Jointly estimating all parameters leads to incidental-parameter bias and incorrect inference. As an alternative, we develop an approach based on a sufficient statistic that separates inference on the homophily parameters from estimation of the fixed effects. The estimator is easy to compute and can be applied to both dense and sparse networks, and is shown to have desirable asymptotic properties under sequences of growing networks. We illustrate the improvements of this estimator over maximum likelihood and bias-corrected estimation in a series of numerical experiments. The technique is applied to explain the import and export patterns in a dense network of countries and to estimate a more sparse advice network among attorneys in a corporate law firm.     

A Note on Functional Forms of Order Statistics

Functional forms of order statistics, as the solution of a system of equations, are studied. The case of the smaller and the larger of two random variables is discussed in detail. Some applications for normal and binomial distributions are presented.     

Simulation of the p-generalized Gaussian distribution

In this paper, we introduce the p-generalized polar methods for the simulation of the p-generalized Gaussian distribution. On the basis of geometric measure representations, the well-known Box–Muller method and the Marsaglia–Bray rejecting polar method for the simulation of the Gaussian distribution are generalized to simulate the p-generalized Gaussian distribution, which fits much more flexibly to data than the Gaussian distribution and has already been applied in various fields of modern sciences. To prove the correctness of the p-generalized polar methods, we give stochastic representations, and to demonstrate their adequacy, we perform a comparison of six simulation techniques w.r.t. the goodness of fit and the complexity. The competing methods include adapted general methods and another special method. Furthermore, we prove stochastic representations for all the adapted methods.     

The mean of linear signed rank statistics under alternatives

An expression is presented for the mean of a linear signed rank statistic for the one sample location model. Several examples are given to illustrate its application.     

Nonparametric regression with cross‐classified responses

In this article, we develop regression models with cross‐classified responses. Conditional independence structures can be explored/exploited through the selective inclusion/exclusion of terms in a certain functional ANOVA decomposition, and the estimation is done nonparametrically via the penalized likelihood method. A cohort of computational and data analytical tools are presented, which include cross‐validation for smoothing parameter selection, Kullback–Leibler projection for model selection, and Bayesian confidence intervals for odds ratios. Random effects are introduced to model possible correlations such as those found in longitudinal and clustered data. Empirical performances of the methods are explored in simulation studies of limited scales, and a real data example is presented using some eyetracking data from linguistic studies. The techniques are implemented in a suite of R functions, whose usage is briefly described in the appendix. The Canadian Journal of Statistics 39: 591–609; 2011. © 2011 Statistical Society of Canada