Tests for independence in a bivariate negative binomial model

The score test and LR test statistic for testing independence are proposed in a bivariate negative binomial regression model. We also propose an adjusted score test in order to enhance the efficiency of the score test. This study is an extension of the work in a univariate model by Dean and Lawless [Dean, C., Lawless, F. (1989). Tests for detecting overdispersion in Poisson regression models. Journal of the American Statistical Association, 84, 467–472]. The adjusted score test proposed in this study is more efficient than the complicated LR test.     

Interval estimation of the negative binomial dispersion parameter

Inference concerning the negative binomial dispersion parameter, denoted by c, is important in many biological and biomedical investigations. Properties of the maximum-likelihood estimator of c and its bias-corrected version have been studied extensively, mainly, in terms of bias and efficiency [W.W. Piegorsch, Maximum likelihood estimation for the negative binomial dispersion parameter, Biometrics 46 (1990), pp. 863–867; S.J. Clark and J.N. Perry, Estimation of the negative binomial parameter κ by maximum quasi-likelihood, Biometrics 45 (1989), pp. 309–316; K.K. Saha and S.R. Paul, Bias corrected maximum likelihood estimator of the negative binomial dispersion parameter, Biometrics 61 (2005), pp. 179–185]. However, not much work has been done on the construction of confidence intervals (C.I.s) for c. The purpose of this paper is to study the behaviour of some C.I. procedures for c. We study, by simulations, three Wald type C.I. procedures based on the asymptotic distribution of the method of moments estimate (mme), the maximum-likelihood estimate (mle) and the bias-corrected mle (bcmle) [K.K. Saha and S.R. Paul, Bias corrected maximum likelihood estimator of the negative binomial dispersion parameter, Biometrics 61 (2005), pp. 179–185] of c. All three methods show serious under-coverage. We further study parametric bootstrap procedures based on these estimates of c, which significantly improve the coverage probabilities. The bootstrap C.I.s based on the mle (Boot-MLE method) and the bcmle (Boot-BCM method) have coverages that are significantly better (empirical coverage close to the nominal coverage) than the corresponding bootstrap C.I. based on the mme, especially for small sample size and highly over-dispersed data. However, simulation results on lengths of the C.I.s show evidence that all three bootstrap procedures have larger average coverage lengths. Therefore, for practical data analysis, the bootstrap C.I. Boot-MLE or Boot-BCM should be used, although Boot-MLE method seems to be preferable over the Boot-BCM method in terms of both coverage and length. Furthermore, Boot-MLE needs less computation than Boot-BCM.     

Parallelizing computation of expected values in recombinant binomial trees

Recombinant binomial trees are binary trees where each non-leaf node has two child nodes, but adjacent parents share a common child node. Such trees arise in option pricing in finance. For example, an option can be valued by evaluating the expected payoffs with respect to random paths in the tree. The cost to exactly compute expected values over random paths grows exponentially in the depth of the tree, rendering a serial computation of one branch at a time impractical. We propose a parallelization method that transforms the calculation of the expected value into an embarrassingly parallel problem by mapping the branches of the binomial tree to the processes in a multiprocessor computing environment. We also discuss a parallel Monte Carlo method and verify the convergence and the variance reduction behavior by simulation study. Performance results from R and Julia implementations are compared on a distributed computing cluster.     

A transient analysis of Markov fluid models with jumps

In this paper, we consider Markov fluid models with jumps which are useful for e.g. insurance risk modeling and the performance analysis of high-speed data networks. Recently, Ahn and Ramaswami [Ahn, S. & Ramaswami, V. (2004). Transient analysis of fluid flow models via stochastic coupling to a queue. Stochastic Models, 20 (1) 71–101] provided a transient analysis of the Markov modulated fluid flow model using stochastic coupling to a queueing model. Here we extend their results and provide a transient analysis of Markov fluid models with jumps. We also present some numerical examples.     

Inference in segmented line regression: a simulation study

A segmented line regression model has been used to describe changes in cancer incidence and mortality trends [Kim, H.-J., Fay, M.P., Feuer, E.J. and Midthune, D.N., 2000, Permutation tests for joinpoint regression with applications to cancer rates. Statistics in Medicine, 19, 335–351. Kim, H.-J., Fay, M.P., Yu, B., Barrett., M.J. and Feuer, E.J., 2004, Comparability of segmented line regression models. Biometrics, 60, 1005–1014.]. The least squares fit can be obtained by using either the grid search method proposed by Lerman [Lerman, P.M., 1980, Fitting segmented regression models by grid search. Applied Statistics, 29, 77–84.] which is implemented in Joinpoint 3.0 available at http://srab.cancer.gov/joinpoint/index.html, or by using the continuous fitting algorithm proposed by Hudson [Hudson, D.J., 1966, Fitting segmented curves whose join points have to be estimated. Journal of the American Statistical Association, 61, 1097–1129.] which will be implemented in the next version of Joinpoint software. Following the least squares fitting of the model, inference on the parameters can be pursued by using the asymptotic results of Hinkley [Hinkley, D.V., 1971, Inference in two-phase regression. Journal of the American Statistical Association, 66, 736–743.] and Feder [Feder, P.I., 1975a, On asymptotic distribution theory in segmented regression Problems-Identified Case. The Annals of Statistics, 3, 49–83.] Feder [Feder, P.I., 1975b, The log likelihood ratio in segmented regression. The Annals of Statistics, 3, 84–97.] Via simulations, this paper empirically examines small sample behavior of these asymptotic results, studies how the two fitting methods, the grid search and the Hudson's algorithm affect these inferential procedures, and also assesses the robustness of the asymptotic inferential procedures.     

A Strong Limit Theorem for Weighted Sums of Negatively Dependent Random Variables

In this article, we establish a new complete convergence theorem for weighted sums of negatively dependent random variables. As corollaries, many results on the almost sure convergence and complete convergence for weighted sums of negatively dependent random variables are obtained. In particular, the results of Jing and Liang (2008 Jing, B.Y., Liang, H.Y. (2008). Strong limit theorems for weighted sums of negatively associated random variables. J. Theor. Probab. 21:890–909.[Crossref], [Web of Science ®] , [Google Scholar]), Sung (2012 Sung, S.H. (2012). Complete convergence for weighted sums of negatively dependent random variables. Stat. Pap. 53:73–82.[Crossref], [Web of Science ®] , [Google Scholar]), and Wu (2010) can be obtained.     

Simulation of extremes of diffusions

The authors consider the problem of simulating the times of events such as extremes and barrier crossings in diffusion processes. They develop a rejection sampler based on Shepp [Shepp, Journal of Applied Probability 1979; 16:423–427] for simulating an extreme of a Brownian motion and use it in a general recursive scheme for more complex simulations, including simultaneous simulation of the minimum and maximum and application to more general diffusions. They price exotic options that are difficult to price analytically: a knock‐out barrier option with a modified payoff function, a lookback option that includes discounting at the risk‐free interest rate, and a chooser option where the choice is made at the time of a barrier crossing. The Canadian Journal of Statistics 38: 738–755; 2010 © 2010 Statistical Society of Canada     

Book reviews

S.M.Kendall:Multivariate Analysis.Charles Griffin & Co. Ltd., London and High Wycombe 1975, 210 pp

C.T.Leondes (ed.):Control and Dynamic Systems, Advances in Theory and Applications. Vol. 11, Academic Press, New York and London 1974, 516 pp., $ 24.50.

CH. R.Nelson:Applied Time Series Analysis for Managerial Forecasting. Holden Day, Inc., San Francisco 1973, 231 pp., $ 14.95.

B.DE Finetti:Theory of Probability. Vol. 1, 2, John Wiley & Sons, New York, London, Sydney, Toronto. Vol. 1, 1974, 300 pp., £ 7.50 - Vol. 2 1975, 375 pp., £ 10.50.

P.Erdördos, J. Spencer:Probabilistic Methods in Combinatorics. Akadémic Press, New York and London; Akadémiai Kiadó, Budapest 1974, 106 pp., $ 11.75.

J.S.R.Ustagi:Variational Methods in Statistics.Academic Press, New York and London 1975.

J.S.Rustagi:Optimizing Methods in Statistics.Academic Press, New York and London 1971,488 pp., $ 17.00.

Karl V. Bury:Statistical Models in Applied Science. John Wiley & Sons, New York-London-Sydney-Toronto 1975,.625 pp., £ 15.60; $ 28.00.

Michael R. Anderberg:Cluster Analysis for Applications. Academic Press, New York-San Francisco-London 1973, 359 pp., $ 27,–.

J.L.Fleiss:Statistical Methods for Rates and Proportions. John Wiley & Sons, New York-London-Sydney-Toronto 1973. 223 pp., £ 6.50.

J.Tanur et al. (Ed.):Statistics:A Guide to the Unknown. Holden Day, Inc., San Francisco 1972, 430 pp.

H.VÁliaho, T. Pekkonen:A Procedure for Stepwise Regression Analysis.Akademie-Verlag, Berlin 1976, 90 pp., 18,– M.

M.Reinfeldt, U. TrÁnkle:Signifikanztabellen statistischer Testvertellungen. R. Oldenbourg Verlag, Mnchen, Wien 1976, 151 S., DM 44,–.

R.E.Barlow, D.J.Batholomew, J.M.Bremner, H.D.Brunk:Statistical Inference Under Order Restrictions.(The Theory and Applications of Isotonic Regression.) John Wiley & Sons, New York 1972, 388 pp., £ 7.50.

H.J.Larson:Introduction to Probability Theory and Statistical Inference. Wiley, New York 1974, 430 pp., £ 6.85.

R.A.Carlson:Statistics. Holden Day, Inc., San Francisco 1973, 393 pp.

E. Page:Queueing Theory in OR. Butterworths, London 1972, 187 pp., £ 3.60.

H.Krampe, J.Kubat, W.Runge:Bedienungsmodelle. Ein Leitfaden für die praktische Anwendung, Verlag Die Wirtachaft, Berlin 1973, 512 S., 79,– M.

G.S.Fishman:Concepts and Methods in Discrete Events Digital Simulation. John Wiley & Sons, New York 1973, 385 pp., £ 8,75.     

Improved likelihood inference for the shape parameter in Weibull regression

We obtain adjustments to the profile likelihood function in Weibull regression models with and without censoring. Specifically, we consider two different modified profile likelihoods: (i) the one proposed by Cox and Reid [Cox, D.R. and Reid, N., 1987, Parameter orthogonality and approximate conditional inference. Journal of the Royal Statistical Society B, 49, 1–39.], and (ii) an approximation to the one proposed by Barndorff–Nielsen [Barndorff–Nielsen, O.E., 1983, On a formula for the distribution of the maximum likelihood estimator. Biometrika, 70, 343–365.], the approximation having been obtained using the results by Fraser and Reid [Fraser, D.A.S. and Reid, N., 1995, Ancillaries and third-order significance. Utilitas Mathematica, 47, 33–53.] and by Fraser et al. [Fraser, D.A.S., Reid, N. and Wu, J., 1999, A simple formula for tail probabilities for frequentist and Bayesian inference. Biometrika, 86, 655–661.]. We focus on point estimation and likelihood ratio tests on the shape parameter in the class of Weibull regression models. We derive some distributional properties of the different maximum likelihood estimators and likelihood ratio tests. The numerical evidence presented in the paper favors the approximation to Barndorff–Nielsen's adjustment.     

Assessing the performance of variational methods for mixed logistic regression models

We present a variational estimation method for the mixed logistic regression model. The method is based on a lower bound approximation of the logistic function [Jaakkola, J.S. and Jordan, M.I., 2000, Bayesian parameter estimation via variational methods. Statistics & Computing, 10, 25–37.]. Based on the approximation, an EM algorithm can be derived that results in a considerable simplification of the maximization problem in that it does not require the numerical evaluation of integrals over the random effects. We assess the performance of the variational method for the mixed logistic regression model in a simulation study and an empirical data example, and compare it to Laplace's method. The results indicate that the variational method is a viable choice for estimating the fixed effects of the mixed logistic regression model under the condition that the number of outcomes within each cluster is sufficiently high.     

Book reviews

E. Kofler, G. Menges: Entscheidungen bei unvollständiger Information. Lecture Notes in Economics and Mathematical Systems 136. Springer-Verlag, Berlin-Heidelberg-NewYork 1976, 357 S., DM 31.

B. S. Everitt: The Analysis ot Contingency Tables. Chapman and Han, London 1977, 128 S., £ 3.75.

O. Barndorff-Nielsen: Information and Exponential Families in Statistical Theory. J. Wiley & Sons, Chichester-New York-Brisbane-Toronto 1978, 247 S., £ 13.50; $ 28.50.

W. Gilchrist: Statisticai Forecasting. J. Wiley & Sons, London-New York-Sydney-Toronto 1976, 321 S., £ 9.50; $ 19.00.

J. M. Chambers: Computational Methods for Data Analysis. J. Wiley & Sons, New York Chichester-Brisbane-Toronto 1977, 279 pp., £ 11.30.

A. Hughes, D. Grawoig: Statistics: A Foundation for Analysis. Addison-Wesley Publishing Company, Reading-Menlo-Park London-Don Mills 1971, 525 S., $ 11.50.

K.Krickeberg, H.Ziezold: Stochastiche Methoden. Springer Verlag, Berlin-Heidelberg-New York 1977, 201 S., 13 Abb., DM 28.

Ch. Schneeweiss: Inventory-Production Theory. A Linear Policy Approach. Lecture Notes in Economics and Mathematical Systems 151. Springer Verlag, Berlin-Heidelberg-NewYork 1977, 122 S., 13 Abb., 9 Tab., DM 18.

S. Dworatschek: Grundlagen der Datenverarbeltung, 6., völlig neu bearb. u. erw. Aufl., Walter de Gruyter, Berlin-New York 1977, 538 S., 200 Abb., 212 Üb.Aufg., 59 Fotos, DM 38.

H. E. Steinhagen, S. Fuchs:Objekterkennung. Einführung in die mathematischen Methoden der Zeichenerkennung. VEB Verlag Technik, Berlin 1976, 436 S., 165 Abb., 32 Tab., 49,–.M.

G. Tinhofer: Mathematik für Studlenanfänger, Carl Hanser Verlag, München 1977, 464 S., 191 Abb., DM 38.     

A comparison of two varying coefficient meta-analysis methods for an average risk difference

Two interval estimation methods for a general linear function of binomial proportions have been proposed. One method [Zou GY, Huang W, Zhang X. A note on confidence interval estimation for a linear function of binomial proportions. Comput Statist Data Anal. 2009;53:1080–1085] combines Wilson interval estimates of individual proportions, and the other method [Price RM, Bonett DG. An improved confidence interval for a linear function of binomial proportions. Comput Statist Data Anal. 2004;45:449–456] uses an adjusted Wald interval. Both methods are appropriate in varying coefficient meta-analysis models where the risk differences are allowed to vary across studies. The two methods were compared in a simulation study under realistic meta-analysis conditions and the adjusted Wald method was found to have the best performance characteristics.     

Equivalent Conditions of Complete Moment Convergence of Weighted Sums for ϕ-Mixing Sequence of Random Variables

In this article, the complete moment convergence of weighted sums for ?-mixing sequence of random variables is investigated. By applying moment inequality and truncation methods, the equivalent conditions of complete moment convergence of weighted sums for ?-mixing sequence of random variables are established. These results promote and improve the corresponding results obtained by Li et al. (1995 Li, D.L., Rao, M.B., Jiang, T.F., Wang, X.C. (1995). Complete convergence and almost sure convergence of weighted sums of random variables. J. Theoret. Probab. 8:49–76.[Crossref], [Web of Science ®] , [Google Scholar]) and Gut (1993 Gut, A. (1993). Complete convergence and Cesàro summation for i.i.d. random variables. Probab. Theory Related Fields 97:169–178.[Crossref], [Web of Science ®] , [Google Scholar]) from i.i.d. to ?-mixing setting. Moreover, we obtain the complete moment convergence of moving average processes based on ?-mixing random variables, which extends the result of Kim et al. (2008 Kim, T.S., Ko, M.H. (2008). Complete moment convergence of moving average processes under dependence assumptions. Statist. Probab. Lett. 78:839–846.[Crossref], [Web of Science ®] , [Google Scholar]) in the sense that it does not require a specific mixing rate.     

On Complete Convergence in Marcinkiewicz-Zygmund Type SLLN for END Random Variables and Its Applications

Abstract

In this paper, the complete convergence for maximal weighted sums of extended negatively dependent (END, for short) random variables is investigated. Some sufficient conditions for the complete convergence and some applications to a nonparametric model are provided. The results obtained in the paper generalize and improve the corresponding ones of Wang et al. (2014 Wang, X. J., X. Deng, L. L. Zheng, and S. H. Hu. 2014. Complete convergence for arrays of rowwise negatively superadditive-dependent random variables and its applications. A Journal of Theoretical and Applied Statistics 48(4):834–50. [Google Scholar]b) and Shen, Xue, and Wang (2017 Shen, A., M. Xue, and W. Wang. 2017. Complete convergence for weighted sums of extended negatively dependent random variables. Communications in Statistics – Theory and Methods 46(3):1433–44.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]).     

Fair Valuation of Life Insurance Contracts Under a Two-Sided Jump Diffusion Model

In this article, we study the fair valuation of participating life insurance contract, which is one of the most common life insurance products, under the two-sided jump diffusion model with the consideration of default risk. The participating life insurance contracts considered here can be expressed as portfolios of options as shown by Grosen and Jøgrgensen (1997 Grosen , A. , Jøgrgensen , P. (1997). Valuation of early exercisable interest rate guarantees. J. Risk Ins. 64:481–503.[Crossref], [Web of Science ®] , [Google Scholar]). We can give the Laplace transforms for these options under the two-sided jump diffusion model, and then price these options by inverting Laplace transforms.     

Geometric characterization of planar regression

The geometric characterization of linear regression in terms of the ‘concentration ellipse’ by Galton [Galton, F., 1886, Family likeness in stature (with Appendix by Dickson, J.D.H.). Proceedings of the Royal Society of London, 40, 42–73.] and Pearson [Pearson, K., 1901, On lines and planes of closest fit to systems of points in space. Philosophical Magazine, 2, 559–572.] was extended to the case of unequal variances of the presumably uncorrelated errors in the experimental data [McCartin, B.J., 2003, A geometric characterization of linear regression. Statistics, 37(2), 101–117.]. In this paper, this geometric characterization is further extended to planar (and also linear) regression in three dimensions where a beautiful interpretation in terms of the concentration ellipsoid is developed.     

A Unified Approach to Ordering Comparison of GPS Distributions with their Mixtures

Recently, several authors have been concerned with ordering comparison of known distributions of the family of generalized power series (GPS) distributions with their mixtures in various senses. In this article, we shall employ a unified approach and obtain similar results, more generally, for all members of the class of the GPS distributions. Some of the previous findings of Misra et al. (2003 Misra , N. , Singh , H. , Harner , E. J. ( 2003 ). Stochastic comparisons of Poisson and binomial random variables with their mixtures . Statist. Probab. Lett. 65 : 279 – 290 .[Crossref], [Web of Science ®] , [Google Scholar]), Alamatsaz and Abbasi (2008 Alamatsaz , M. H. , Abbasi , S. ( 2008 ). Ordering comparison of negative binomial random variables with their mixtures . Statist. Probab. Lett. 78 : 2234 – 2239 . [Google Scholar]), and Aghababaei Jazi and Alamatsaz (2010 Aghababaei Jazi , M. , Alamatsaz , M. H. ( 2010 ). Ordering comparison of logarithmic series random variables with their mixtures . Commun. Statist. Theor. Meth. 39 : 3252 – 3263 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) in this connection, then, follow as corollaries. Further, we have derived some more ordering comparison results.     

Revisiting Beal's Confidence Intervals for the Difference of Two Binomial Proportions

Confidence interval construction for the difference of two independent binomial proportions is a well-known problem with a full panoply of proposed solutions. In this paper, we focus largely on the family of intervals proposed by Beal (1987 Beal , S. ( 1987 ). Asymptotic confidence intervals for the difference between two binomial parameters for use with small samples . Biometrics 43 : 941 – 950 . [CSA] [CROSSREF] [Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). This family, which includes the Haldane and Jeffreys–Perks intervals as special cases, assumes a symmetric prior distribution for the population proportions p ₁ and p ₂. We propose new methods that allow the currently observed data to set the prior distribution by taking a parametric empirical-Bayes approach; in addition, we also provide an investigation of the new interval' behaviors in small-sample situations. Unlike other solutions, our intervals can be used adaptively for experiments conducted in multiple stages over time. We illustrate this notion using data from an Argentinean study involving the Mal Rio Cuarto virus and its transmission to susceptible maize crops.     

Dependent Lindeberg central limit theorem for the fidis of empirical processes of cluster functionals

Drees H. and Rootzén H. [Limit theorems for empirical processes of cluster functionals (EPCF). Ann Stat. 2010;38(4):2145–2186] have proven central limit theorems (CLTs) for EPCF built from β-mixing processes. However, this family of β-mixing processes is quite restrictive. We expand some of those results, for the finite-dimensional marginal distributions (fidis), to a more general dependent processes family, known as weakly dependent processes in the sense of Doukhan P. and Louhichi S. [A new weak dependence condition and applications to moment inequalities. Stoch. Proc. Appl. 1999;84:313–342]. In this context, the CLT for the fidis of EPCF is sufficient in some applications. For instance, we prove the convergence without mixing conditions of the extremogram estimator, including a small example with simulation of the extremogram of a weakly dependent random process but nonmixing, in order to confirm the efficacy of our result.