Numerical Methods of Statistics: Contents

This book explains how computer software is designed to perform the tasks required for sophisticated statistical analysis. For statisticians, it examines the nitty-gritty computational problems behind statistical methods. For mathematicians and computer scientists, it looks at the application of mathematical tools to statistical problems. The first half of the book offers a basic background in numerical analysis that emphasizes issues important to statisticians. The next several chapters cover a broad array of statistical tools, such as maximum likelihood and nonlinear regression. The author also treats the application of numerical tools; numerical integration and random number generation are explained in a unified manner reflecting complementary views of Monte Carlo methods. Each chapter contains exercises that range from simple questions to research problems. Most of the examples are accompanied by demonstration and source code available in from the author's Web site. New in this second edition are demonstrations coded in R, as well as new sections on linear programming and the Nelder-Mead search algorithm.

[1]  Forman S. Acton,et al.  Numerical methods that work , 1970 .

[2]  Takuji Nishimura,et al.  Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator , 1998, TOMC.

[3]  D. Boos.  On Generalized Score Tests , 1992 .

[4]  M. Dowell,et al.  A modified regula falsi method for computing the root of an equation , 1971 .

[5]  Nancy E. Bosten,et al.  Incomplete Beta Ratio (Remark on Algorithm 179) , 1974, Commun. ACM.

[6]  Marcello Pagano,et al.  When is an altoregressive scheme stationary , 1973 .

[7]  T. W. Anderson,et al.  Statistical analysis of time series , 1972 .

[8]  William J. Stewart,et al.  A Simultaneous Iteration Algorithm for Real Matrices , 1981, TOMS.

[9]  Brian D. Ripley,et al.  Stochastic Simulation , 2005 .

[10]  Changbao Wu,et al.  Asymptotic Theory of Nonlinear Least Squares Estimation , 1981 .

[11]  J. Geweke,et al.  Bayesian Inference in Econometric Models Using Monte Carlo Integration , 1989 .

[12]  B. Ripley.  Computer generation of random variables: a tutorial , 1983 .

[13]  Maxwell L. King,et al.  IMPROVING THE NUMERICAL TECHNIQUE FOR COMPUTING THE ACCUMULATED DISTRIBUTION OF A QUADRATIC FORM IN NORMAL VARIABLES , 2002 .

[14]  Ted G. Lewis,et al.  Generalized Feedback Shift Register Pseudorandom Number Algorithm , 1973, JACM.

[15]  Dennis D. Boos,et al.  Using extreme value theory to estimate large percentiles , 1984 .

[16]  John F. Monahan.  Extensions of von Neumann’s method for generating random variables , 1979 .

[17]  Thomas J. Santner,et al.  A note on A. Albert and J. A. Anderson's conditions for the existence of maximum likelihood estimates in logistic regression models , 1986 .

[18]  W. Weston Meyer,et al.  Optimal error bounds for cubic spline interpolation , 1976 .

[19]  Ross Ihaka,et al.  Lexical Scope and Statistical Computing , 2000 .

[20]  J. G. Ramage,et al.  Computer methods for sampling from student's t distribution , 1977 .

[21]  J. Roscoe,et al.  An Investigation of the Restraints with Respect to Sample Size Commonly Imposed on the Use of the Chi-Square Statistic , 1971 .

[22]  J. Monahan.  Fully Bayesian analysis of ARMA time series models , 1983 .

[23]  R. Singleton.  An algorithm for computing the mixed radix fast Fourier transform , 1969 .

[24]  E S Eppright,et al.  Nutrition of infants and preschool children in the North Central Region of the United States of America. , 1972, World review of nutrition and dietetics.

[25]  R. Lew,et al.  An Approximation to the Cumulative Normal Distribution with Simple Coefficients , 1981 .

[26]  G. S. Fishman,et al.  A Statistical Evaluation of Multiplicative Congruential Random Number Generators with Modulus 231 — 1 , 1982 .

[27]  B. Efron,et al.  Assessing the accuracy of the maximum likelihood estimator: Observed versus expected Fisher information , 1978 .

[28]  Chris Chatfield,et al.  Introduction to Statistical Time Series. , 1976 .

[29]  L. Devroye.  Non-Uniform Random Variate Generation , 1986 .

[30]  S. E. Hills,et al.  Illustration of Bayesian Inference in Normal Data Models Using Gibbs Sampling , 1990 .

[31]  R. H. Myers.  Classical and modern regression with applications , 1986 .

[32]  S. Tezuka.  Uniform Random Numbers: Theory and Practice , 1995 .

[33]  A. Stroud.  Approximate calculation of multiple integrals , 1973 .

[34]  Irving John Good,et al.  Exact Distributions for χ2 and for the Likelihood-Ratio Statistic for the Equiprobable Multinomial Distribution , 1970 .

[35]  A. Stroud,et al.  Approximate Calculation of Integrals , 1962 .

[36]  John F. Monahan.  An algorithm for generating chi random variables , 1987, TOMS.

[37]  Josef Leydold,et al.  Automatic sampling with the ratio-of-uniforms method , 2000, TOMS.

[38]  Daniel A. Relles.  Variance Reduction Techniques for Monte Carlo Sampling from Student Distributions , 1970 .

[39]  T. A. Ryan,et al.  Algorithm 516: An Algorithm for Obtaining Confidence Intervals and Point Estimates Based on Ranks in the Two-Sample Location Problem [G1] , 1977, TOMS.

[40]  M. Powell.  A New Algorithm for Unconstrained Optimization , 1970 .

[41]  Henry R. Neave,et al.  On Using the Box‐Müller Transformation with Multiplicative Congruential Pseudo‐Random Number Generators , 1973 .

[42]  M. Bartlett.  XX.—On the Theory of Statistical Regression. , 1934 .

[43]  Å. Björck.  Solving linear least squares problems by Gram-Schmidt orthogonalization , 1967 .

[44]  R. R. Hocking,et al.  Algorithm AS 53: Wishart Variate Generator , 1972 .

[45]  R. Davies.  Numerical inversion of a characteristic function , 1973 .

[46]  M. R. B. Clarke,et al.  A Givens Algorithm for Moving from One Linear Model to Another Without Going Back to the Data , 1981 .

[47]  Walter R. Gilks,et al.  Introduction to general state-space Markov chain theory , 1995 .

[48]  I. P. Mysovskikh.  THE APPROXIMATION OF MULTIPLE INTEGRALS BY USING INTERPOLATORY CUBATURE FORMULAE , 1980 .

[49]  Peter Craven,et al.  Smoothing noisy data with spline functions , 1978 .

[50]  G. Marsaglia.  Generating Exponential Random Variables , 1961 .

[51]  David G. Luenberger,et al.  Linear and nonlinear programming , 1984 .

[52]  R. Kronmal,et al.  On the Alias Method for Generating Random Variables From a Discrete Distribution , 1979 .

[53]  S. Chib,et al.  Understanding the Metropolis-Hastings Algorithm , 1995 .

[54]  Roy E. Welsch,et al.  Efficient Computing of Regression Diagnostics , 1981 .

[55]  R. Davies.  The distribution of a linear combination of 2 random variables , 1980 .

[56]  R. Tibshirani.  Regression Shrinkage and Selection via the Lasso , 1996 .

[57]  Philip Rabinowitz,et al.  Methods of Numerical Integration , 1985 .

[58]  Adrian F. M. Smith,et al.  Bayesian Inference for Generalized Linear and Proportional Hazards Models Via Gibbs Sampling , 1993 .

[59]  P. Hall.  On Some Simple Estimates of an Exponent of Regular Variation , 1982 .

[60]  E. Frome.  The analysis of rates using Poisson regression models. , 1983, Biometrics.

[61]  Peter D. Welch,et al.  Historical notes on the fast Fourier transform , 1967 .

[62]  Bradley P. Carlin,et al.  Markov Chain Monte Carlo in Practice: A Roundtable Discussion , 1998 .

[63]  C. Kelley.  Iterative Methods for Linear and Nonlinear Equations , 1987 .

[64]  Gary Simon.  Computer Simulation Swindles, with Applications to Estimates of Location and Dispersion , 1976 .

[65]  A. Albert,et al.  On the existence of maximum likelihood estimates in logistic regression models , 1984 .

[66]  G. Forsythe.  Von Neumann''s comparison method for random sampling from the normal and other distributions. , 1972 .

[67]  M. Stephens.  Use of the Kolmogorov-Smirnov, Cramer-Von Mises and Related Statistics without Extensive Tables , 1970 .

[68]  Armido R. Didonato,et al.  Algorithm 708: Significant digit computation of the incomplete beta function ratios , 1988, TOMS.

[69]  C. Chatfield,et al.  Fourier Analysis of Time Series: An Introduction , 1977, IEEE Transactions on Systems, Man, and Cybernetics.

[70]  Nitin R. Patel,et al.  ALGORITHM 643: FEXACT: a FORTRAN subroutine for Fisher's exact test on unordered r×c contingency tables , 1986, TOMS.

[71]  H. Keng,et al.  Applications of number theory to numerical analysis , 1981 .

[72]  M. E. Muller,et al.  A Note on the Generation of Random Normal Deviates , 1958 .

[73]  John D. Vedder,et al.  An invertible approximation to the normal distribution function , 1993 .

[74]  Alfred V. Aho,et al.  The Design and Analysis of Computer Algorithms , 1974 .

[75]  E. Haeusler,et al.  On Asymptotic Normality of Hill's Estimator for the Exponent of Regular Variation , 1985 .

[76]  Mark Priestley.  Practical Object-Oriented Design , 1996 .

[77]  Alan Jennings,et al.  Matrix Computation for Engineers and Scientists , 1977 .

[78]  C. Geyer,et al.  Discussion: Markov Chains for Exploring Posterior Distributions , 1994 .

[79]  T. A. Bray,et al.  A Convenient Method for Generating Normal Variables , 1964 .

[80]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[81]  S. Fienberg,et al.  Efficient Calculation of All Possible Regressions , 1968 .

[82]  William J. Stewart,et al.  A comparison of numerical techniques in Markov modeling , 1978, CACM.

[83]  Mark B. Wells.  Generation of permutations by transposition , 1961 .

[84]  D. Tritchler.  An Algorithm for Exact Logistic Regression , 1984 .

[85]  F. R. Gantmakher.  The Theory of Matrices , 1984 .

[86]  G. Alefeld,et al.  Introduction to Interval Computation , 1983 .

[87]  M GayDavid,et al.  Algorithm 611: Subroutines for Unconstrained Minimization Using a Model/Trust-Region Approach , 1983 .

[88]  D. G. Watts,et al.  Relative Curvature Measures of Nonlinearity , 1980 .

[89]  W. E. Gentleman.  Least Squares Computations by Givens Transformations Without Square Roots , 1973 .

[90]  G. Casella,et al.  Explaining the Gibbs Sampler , 1992 .

[91]  J. Neyman,et al.  Consistent Estimates Based on Partially Consistent Observations , 1948 .

[92]  M. R. Osborne,et al.  On the LASSO and its Dual , 2000 .

[93]  Andrew F. Siegel,et al.  Unbiased Monte Carlo Integration Methods with Exactness for Low Order Polynomials , 1985 .

[94]  Carl Tim Kelley,et al.  Iterative methods for optimization , 1999, Frontiers in applied mathematics.

[95]  C. A. R. Hoare,et al.  Algorithm 64: Quicksort , 1961, Commun. ACM.

[96]  R. E. Kalman,et al.  Optimum Seeking Methods. , 1964 .

[97]  G. M. Furnival.  All Possible Regressions with Less Computation , 1971 .

[98]  Jorge J. Moré,et al.  Testing Unconstrained Optimization Software , 1981, TOMS.

[99]  Pierre L'Ecuyer,et al.  Random numbers for simulation , 1990, CACM.

[100]  Lauwerens Kuipers,et al.  Uniform distribution of sequences , 1974 .

[101]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[102]  Walter Gautschi,et al.  A Computational Procedure for Incomplete Gamma Functions , 1979, TOMS.

[103]  C. Ansley.  An algorithm for the exact likelihood of a mixed autoregressive-moving average process , 1979 .

[104]  John E. Dennis,et al.  An Adaptive Nonlinear Least-Squares Algorithm , 1977, TOMS.

[105]  R. Koenker,et al.  The Gaussian hare and the Laplacian tortoise: computability of squared-error versus absolute-error estimators , 1997 .

[106]  Harald Niederreiter,et al.  A statistical analysis of generalized feedback shift register pseudorandom number generators , 1987 .

[107]  George Marsaglia,et al.  Uniform Random Number Generators , 1965, JACM.

[108]  J. Monahan,et al.  Spherical-Radial Integration Rules for Bayesian Computation , 1997 .

[109]  L. L. Cam,et al.  On the Asymptotic Theory of Estimation and Testing Hypotheses , 1956 .

[110]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[111]  Shu Tezuka,et al.  The k-distribution of generalized feedback shift register pseudorandom numbers , 1983, CACM.

[112]  M. Pagano,et al.  On Obtaining Permutation Distributions in Polynomial Time , 1983 .

[113]  Gene H. Golub,et al.  Rules and software for detecting rank degeneracy , 1980 .

[114]  B. Silverman,et al.  Some Aspects of the Spline Smoothing Approach to Non‐Parametric Regression Curve Fitting , 1985 .

[115]  Robert F. Ling,et al.  A Study of the Accuracy of Some Approximations for t, χ2, and F Tail Probabilities , 1978 .

[116]  New York Dover,et al.  ON THE CONVERGENCE PROPERTIES OF THE EM ALGORITHM , 1983 .

[117]  B. Silverman.  Density estimation for statistics and data analysis , 1986 .

[118]  D. Goldfarb.  Factorized variable metric methods for unconstrained optimization , 1976 .

[119]  L. Rabiner,et al.  The chirp z-transform algorithm and its application , 1969 .

[120]  L. Mark Berliner,et al.  Subsampling the Gibbs Sampler , 1994 .

[121]  James W. Longley.  An Appraisal of Least Squares Programs for the Electronic Computer from the Point of View of the User , 1967 .

[122]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[123]  R. Tausworthe.  Random Numbers Generated by Linear Recurrence Modulo Two , 1965 .

[124]  Tony F. Chan,et al.  An Improved Algorithm for Computing the Singular Value Decomposition , 1982, TOMS.

[125]  John F. Monahan,et al.  Algorithm 616: fast computation of the Hodges-Lehmann location estimator , 1984, ACM Trans. Math. Softw..

[126]  W. Gilks,et al.  Adaptive Rejection Metropolis Sampling Within Gibbs Sampling , 1995 .

[127]  L. R. Moore,et al.  An Exhaustive Analysis of Multiplicative Congruential Random Number Generators with Modulus $2^{31} - 1$ , 1986 .

[128]  Richard M. Heiberger,et al.  Generation of Random Orthogonal Matrices , 1978 .

[129]  Philip Heidelberger,et al.  Simulation Run Length Control in the Presence of an Initial Transient , 1983, Oper. Res..

[130]  G. C. Tiao,et al.  Bayesian inference in statistical analysis , 1973 .

[131]  G. Wahba,et al.  Gcvpack – routines for generalized cross validation , 1987 .

[132]  T. Hesterberg,et al.  Weighted Average Importance Sampling and Defensive Mixture Distributions , 1995 .

[133]  Alastair J. Walker,et al.  An Efficient Method for Generating Discrete Random Variables with General Distributions , 1977, TOMS.

[134]  A. Raftery,et al.  How Many Iterations in the Gibbs Sampler , 1991 .

[135]  P. A. W. Lewis,et al.  A Pseudo-Random Number Generator for the System/360 , 1969, IBM Syst. J..

[136]  G. Casella,et al.  The Effect of Improper Priors on Gibbs Sampling in Hierarchical Linear Mixed Models , 1996 .

[137]  G. P. Bhattacharjee,et al.  Algorithm AS 63: The Incomplete Beta Integral , 1973 .

[138]  Gene H. Golub,et al.  The simplex method of linear programming using LU decomposition , 1969, Commun. ACM.

[139]  A. F. Smith,et al.  Conjugate likelihood distributions , 1993 .

[140]  G. Stewart.  The Efficient Generation of Random Orthogonal Matrices with an Application to Condition Estimators , 1980 .

[141]  A. Raftery,et al.  Space-time modeling with long-memory dependence: assessing Ireland's wind-power resource. Technical report , 1987 .

[142]  Linus Schrage,et al.  A More Portable Fortran Random Number Generator , 1979, TOMS.

[143]  Jack Dongarra,et al.  LINPACK Users' Guide , 1987 .

[144]  P. D. Coddington,et al.  Analysis of random number generators using Monte Carlo simulation , 1993, cond-mat/9309017.

[145]  John F. Monahan,et al.  Accuracy in random number generation , 1985 .

[146]  M. C. Jones,et al.  Spline Smoothing and Nonparametric Regression. , 1989 .

[147]  G. S. Watson.  Statistics on Spheres , 1983 .

[148]  Michael N. Vrahatis,et al.  Solving systems of nonlinear equations using the nonzero value of the topological degree , 1988, TOMS.

[149]  I. D. Hill,et al.  Remark on Algorithm 179: Incomplete Beta ratio , 1967, CACM.

[150]  James Hardy Wilkinson,et al.  Symmetric decomposition of positive definite band matrices , 1965 .

[151]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[152]  R. Jennrich,et al.  Application of Stepwise Regression to Non-Linear Estimation , 1968 .

[153]  R. Cranley,et al.  Randomization of Number Theoretic Methods for Multiple Integration , 1976 .

[154]  A. W. van der Vaart,et al.  On Profile Likelihood , 2000 .

[155]  S. K. Park,et al.  Random number generators: good ones are hard to find , 1988, CACM.

[156]  Åke Björck,et al.  Numerical Methods , 2021, Markov Renewal and Piecewise Deterministic Processes.

[157]  D. F. Andrews,et al.  Robust Estimates of Location , 1972 .

[158]  Philip E. Gill,et al.  Numerically stable methods for quadratic programming , 1978, Math. Program..

[159]  Ernst Stadlober,et al.  Sampling from Poisson, binomial and hypergeometric distributions: ratio of uniforms as a simple and fast alternative , 1989 .

[160]  John F. Monahan.  An alternative method for computing overflow probabilities , 1987 .

[161]  Douglas M. Bates,et al.  Nonlinear Regression Analysis and Its Applications , 1988 .

[162]  George M. Furnival,et al.  Regressions by leaps and bounds , 2000 .

[163]  Carter Bays,et al.  Improving a Poor Random Number Generator , 1976, TOMS.

[164]  K. Fang,et al.  Number-theoretic methods in statistics , 1993 .

[165]  Petko D. Proinov.  Discrepancy and integration of continuous functions , 1988 .

[166]  S. R. Searle.  Linear Models , 1971 .

[167]  P. J. Huber.  Robust Estimation of a Location Parameter , 1964 .

[168]  T. N. Sriram.  Asymptotics in Statistics–Some Basic Concepts , 2002 .

[169]  Herbert S. Bright,et al.  Quasi-Random Number Sequences from a Long-Period TLP Generator with Remarks on Application to Cryptography , 1979, CSUR.

[170]  John F. Monahan.  Approximating the Log of the Normal Cumulative , 1981 .

[171]  E. R. Golder.  The Spectral Test for the Evaluation of Congruential Pseudo‐Random Generators , 1976 .

[172]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[173]  Roy H. Wampler.  Test procedures and test problems for least squares algorithms , 1980 .

[174]  D. Poirier.  Piecewise Regression Using Cubic Splines , 1973 .

[175]  H. Ruben,et al.  Probability Content of Regions Under Spherical Normal Distributions, III: The Bivariate Normal Integral , 1961 .

[176]  J. Hammersley,et al.  Monte Carlo Methods , 1965 .

[177]  W. Gilks,et al.  Adaptive Rejection Sampling for Gibbs Sampling , 1992 .

[178]  M. F. Hutchinson,et al.  Algorithm 642: A fast procedure for calculating minimum cross-validation cubic smoothing splines , 1986, TOMS.

[179]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .

[180]  J.E. Dennis,et al.  A user's guide to nonlinear optimization algorithms , 1984, Proceedings of the IEEE.

[181]  E. Stadlober,et al.  The ratio of uniforms approach for generating discrete random variates , 1990 .

[182]  Gary D. Knott.  A numbering system for combinations , 1974, CACM.

[183]  L. Tierney,et al.  Accurate Approximations for Posterior Moments and Marginal Densities , 1986 .

[184]  Alan M. Ferrenberg,et al.  Monte Carlo simulations: Hidden errors from "good" random number generators. , 1992, Physical review letters.

[185]  L. Tierney.  Markov Chains for Exploring Posterior Distributions , 1994 .

[186]  B. M. Hill,et al.  A Simple General Approach to Inference About the Tail of a Distribution , 1975 .

[187]  William J. Stewart,et al.  Numerical Solution of Markov Chains , 1993 .

[188]  John F. Hart,et al.  Computer Approximations , 1978 .

[189]  J. L. Hodges,et al.  Estimates of Location Based on Rank Tests , 1963 .

[190]  C. T. Fike.  Computer Evaluation of Mathematical Functions , 1970 .

[191]  Dominique Maisonneuve.  Recherche et Utilisation des “Bons Treillis.” Programmation et Résultats Numériques , 1972 .

[192]  C. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[193]  D. J. Finney,et al.  The estimation from individual records of the relationship between dose and quantal response. , 1947, Biometrika.

[194]  Charles J. Geyer,et al.  Practical Markov Chain Monte Carlo , 1992 .

[195]  Donald B. Rubin,et al.  The Acceptability of Regression Solutions: Another Look at Computational Accuracy , 1976 .

[196]  Gene H. Golub,et al.  Matrix computations , 1983 .

[197]  P. Holland,et al.  Robust regression using iteratively reweighted least-squares , 1977 .

[198]  W. Wong,et al.  The calculation of posterior distributions by data augmentation , 1987 .

[199]  Cecil Hastings,et al.  Approximations for digital computers , 1955 .

[200]  Jorge J. Moré,et al.  The Levenberg-Marquardt algo-rithm: Implementation and theory , 1977 .

[201]  G. Kedem,et al.  A table of good lattice points in three dimensions , 1974 .

[202]  Albert J. Kinderman,et al.  Computer Generation of Random Variables Using the Ratio of Uniform Deviates , 1977, TOMS.

[203]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[204]  Joseph B. Kruskal,et al.  Extremely portable random number generator , 1969, CACM.

[205]  J. Imhof.  Computing the distribution of quadratic forms in normal variables , 1961 .

[206]  Alan Agresti,et al.  Exact conditional tests for cross-classifications: Approximation of attained significance levels , 1979 .

[207]  T. Brubaker,et al.  Nonlinear Parameter Estimation , 1979 .

[208]  H. H. Rosenbrock,et al.  An Automatic Method for Finding the Greatest or Least Value of a Function , 1960, Comput. J..

[209]  B. J. R. Bailey.  Alternatives to Hastings' Approximation to the Inverse of the Normal Cumulative Distribution Function , 1981 .

[210]  J. L. Walsh,et al.  The theory of splines and their applications , 1969 .

[211]  H. C. Williams,et al.  Advanced Linear Programming , 1983, The Mathematical Gazette.

[212]  William M. Waite,et al.  Software manual for the elementary functions , 1980 .

[213]  P. McCullagh,et al.  Generalized Linear Models , 1992 .

[214]  R. C. Elston,et al.  Estimation of Time-Response Curves and Their Confidence Bands , 1962 .

[215]  William G. Poole,et al.  A geometric theory for the QR, LU and power iterations. , 1973 .

[216]  G. W. Snedecor.  Statistical Methods , 1964 .

[217]  G. Marsaglia.  Random numbers fall mainly in the planes. , 1968, Proceedings of the National Academy of Sciences of the United States of America.

[218]  Marcello Pagano,et al.  WHEN IS AN AUTOREGRESSIVE SCHEME STATIONARY , 1973 .

[219]  James Hardy Wilkinson,et al.  Linear algebra , 1971, Handbook for automatic computation.

[220]  Wolfgang Hörmann,et al.  A portable random number generator well suited for the rejection method , 1993, TOMS.

[221]  M. Degroot.  Optimal Statistical Decisions , 1970 .

[222]  J. H. Wilkinson.  The algebraic eigenvalue problem , 1966 .

[223]  G. Marsaglia.  The Structure of Linear Congruential Sequences , 1972 .

[224]  G. Stewart.  The economical storage of plane rotations , 1976 .

[225]  James Durbin,et al.  The fitting of time series models , 1960 .

[226]  L. N. Balaam,et al.  Statistical Theory and Methodology in Science and Engineering , 1966 .

[227]  G. R. Hext,et al.  Sequential Application of Simplex Designs in Optimisation and Evolutionary Operation , 1962 .

[228]  A. Wald.  Note on the Consistency of the Maximum Likelihood Estimate , 1949 .

[229]  David Tritchler,et al.  On Inverting Permutation Tests , 1984 .

[230]  John F. Monahan,et al.  Monte Carlo Comparison of ANOVA, MIVQUE, REML, and ML Estimators of Variance Components , 1984 .

[231]  A. Wald,et al.  On the Choice of the Number of Class Intervals in the Application of the Chi Square Test , 1942 .

[232]  Pierre L'Ecuyer,et al.  Efficient and portable combined random number generators , 1988, CACM.

[233]  Philip E. Gill,et al.  Practical optimization , 1981 .

[234]  M. V. Ramakrishna.  Computing the probability of hash table/urn overflow , 1987 .

[235]  Calyampudi R. Rao,et al.  Linear Statistical Inference and Its Applications. , 1975 .

[236]  G. Stewart,et al.  Matrix Perturbation Theory , 1990 .

[237]  D. R. Divgi.  Calculation of Univariate and Bivariate Normal Probability Functions , 1979 .

[238]  Joseph L. Leva.  A fast normal random number generator , 1992, TOMS.

[239]  R. W. Wedderburn,et al.  On the existence and uniqueness of the maximum likelihood estimates for certain generalized linear models , 1976 .

[240]  G. Wahba.  Interpolating Spline Methods for Density Estimation I. Equi-Spaced Knots , 1975 .

[241]  W. J. Cody,et al.  Rational Chebyshev approximation using linear equations , 1968 .

[242]  I. Weissman.  Estimation of Parameters and Large Quantiles Based on the k Largest Observations , 1978 .

[243]  D. Pollard.  A User's Guide to Measure Theoretic Probability by David Pollard , 2001 .

[244]  Anthony C. Davison,et al.  Bootstrap methods and their application: Tests , 1997 .

[245]  Harald Niederreiter,et al.  Random number generation and Quasi-Monte Carlo methods , 1992, CBMS-NSF regional conference series in applied mathematics.

[246]  Masaaki Sibuya,et al.  On exponential and other random variable generators , 1961 .

[247]  Donald Ervin Knuth,et al.  The Art of Computer Programming , 1968 .

[248]  Gwilym M. Jenkins,et al.  Time series analysis, forecasting and control , 1972 .

[249]  Donald B. Johnson,et al.  Lower Bounds for Selection in X + Y and Other Multisets , 1978, JACM.

[250]  R. Plemmons,et al.  Comparison of Some Direct Methods for Computing Stationary Distributions of Markov Chains , 1984 .

[251]  E. C. Fieller.  SOME PROBLEMS IN INTERVAL ESTIMATION , 1954 .

[252]  John Geweke,et al.  Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments , 1991 .

[253]  J. Monahan.  A note on enforcing stationarity in autoregressive-moving average models , 1984 .

[254]  Gene H. Golub,et al.  Algorithms for Computing the Sample Variance: Analysis and Recommendations , 1983 .

[255]  L. Lecam.  On the Assumptions Used to Prove Asymptotic Normality of Maximum Likelihood Estimates , 1970 .

[256]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[257]  W. F. Trench.  An Algorithm for the Inversion of Finite Toeplitz Matrices , 1964 .

[258]  G. O. Wesolowsky,et al.  On the computation of the bivariate normal integral , 1990 .

[259]  S. Weisberg,et al.  Residuals and Influence in Regression , 1982 .

[260]  Andreas Griewank,et al.  Evaluating derivatives - principles and techniques of algorithmic differentiation, Second Edition , 2000, Frontiers in applied mathematics.

[261]  Alston S. Householder,et al.  The Theory of Matrices in Numerical Analysis , 1964 .

[262]  R. E. Odeh,et al.  Algorithm AS 70: The Percentage Points of the Normal Distribution , 1974 .

[263]  Peter Alfeld,et al.  The exact analysis of sparse rectangular linear systems , 1991, TOMS.

[264]  J. Wolfowitz.  On Wald's Proof of the Consistency of the Maximum Likelihood Estimate , 1949 .

[265]  W. W. Muir,et al.  Regression Diagnostics: Identifying Influential Data and Sources of Collinearity , 1980 .

[266]  A. Gallant,et al.  Nonlinear Statistical Models , 1988 .

[267]  Ronald A. Thisted,et al.  A Remark on as 127. Generation of Random Orthogonal Matrices , 1982 .

[268]  Richard L. Tweedie,et al.  Markov Chains and Stochastic Stability , 1993, Communications and Control Engineering Series.

[269]  Zvi Drezner,et al.  Computation of the bivariate normal integral , 1978 .

[270]  H. Niederreiter.  Quasi-Monte Carlo methods and pseudo-random numbers , 1978 .

[271]  Paul Bratley,et al.  A guide to simulation , 1983 .

[272]  W. G. Cochran.  Some Methods for Strengthening the Common χ 2 Tests , 1954 .

[273]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[274]  Ed Anderson,et al.  LAPACK Users' Guide , 1995 .

[275]  L. Bluestein.  A linear filtering approach to the computation of discrete Fourier transform , 1970 .

[276]  B. Silverman,et al.  Kernel Density Estimation Using the Fast Fourier Transform , 1982 .

[277]  Marcello Pagano,et al.  Efficient Calculation of the Permutation Distribution of Trimmed Means , 1991 .

[278]  Nicholas G. Polson,et al.  On the Geometric Convergence of the Gibbs Sampler , 1994 .

[279]  R. Tapia,et al.  Nonparametric Probability Density Estimation , 1978 .

[280]  Sylvia Richardson,et al.  Markov Chain Monte Carlo in Practice , 1997 .

[281]  C. Reinsch.  Smoothing by spline functions , 1967 .

[282]  F. N. David,et al.  Principles and procedures of statistics. , 1961 .

[283]  D. G. Watts,et al.  Calculation of intrinsic and parameter-effects curvatures for nonlinear regression models , 1983 .

[284]  Jeffrey Scott Vitter,et al.  An efficient algorithm for sequential random sampling , 1987, TOMS.

[285]  I. D. Hill,et al.  An Efficient and Portable Pseudo‐Random Number Generator , 1982 .

[286]  David R. Cox.  The analysis of binary data , 1970 .

[287]  G. Stewart.  Introduction to matrix computations , 1973 .

[288]  M. Pagano,et al.  Methods for Exact Goodness-of-Fit Tests , 1992 .

[289]  G. Wahba.  Smoothing noisy data with spline functions , 1975 .

[290]  H. Bohman.  Numerical inversions of characteristic functions , 1975 .

[291]  George Cybenko,et al.  The Numerical Stability of the Levinson-Durbin Algorithm for Toeplitz Systems of Equations , 1980 .

[292]  J. H. Wilkinson,et al.  AN ESTIMATE FOR THE CONDITION NUMBER OF A MATRIX , 1979 .

[293]  A. Mcintosh.  Fitting Linear Models: An Application of Conjugate Gradient Algorithms , 1982 .

[294]  T. C. Hu,et al.  Combinatorial algorithms , 1982 .

[295]  Nancy Flournoy,et al.  Statistical Multiple Integration , 1991 .

[296]  L. Tierney,et al.  Fully Exponential Laplace Approximations to Expectations and Variances of Nonpositive Functions , 1989 .

[297]  John Monahan,et al.  Stochastic Integration Rules for Infinite Regions , 1998, SIAM J. Sci. Comput..

[298]  A. M. Gross.  A Monte Carlo Swindle for Estimators of Location , 1973 .

[299]  J. Goodnight.  A Tutorial on the SWEEP Operator , 1979 .

[300]  Bradley P. Carlin,et al.  Markov Chain Monte Carlo conver-gence diagnostics: a comparative review , 1996 .

[301]  Donald B. Rubin,et al.  More on Computational Accuracy in Regression: Comment , 1977 .

[302]  A. Stroud,et al.  Gaussian quadrature formulas , 1966 .