Predicting infectious disease outbreak risk via migratory waterfowl vectors

The spread of an emerging infectious disease is a major public health threat. Given the uncertainties associated with vector-borne diseases, in terms of vector dynamics and disease transmission, it is critical to develop statistical models that address how and when such an infectious disease could spread throughout a region such as the USA. This paper considers a spatio-temporal statistical model for how an infectious disease could be carried into the USA by migratory waterfowl vectors during their seasonal migration and, ultimately, the risk of transmission of such a disease to domestic fowl. Modeling spatio-temporal data of this type is inherently difficult given the uncertainty associated with observations, complexity of the dynamics, high dimensionality of the underlying process, and the presence of excessive zeros. In particular, the spatio-temporal dynamics of the waterfowl migration are developed by way of a two-tiered functional temporal and spatial dimension reduction procedure that captures spatial and seasonal trends, as well as regional dynamics. Furthermore, the model relates the migration to a population of poultry farms that are known to be susceptible to such diseases, and is one of the possible avenues toward transmission to domestic poultry and humans. The result is a predictive distribution of those counties containing poultry farms that are at the greatest risk of having the infectious disease infiltrate their flocks assuming that the migratory population was infected. The model naturally fits into the hierarchical Bayesian framework.     

Correlation is first order independent of transformation

We show that the correlation between the estimates of two parameters is almost unchanged if they are each transformed in an arbitrary way. To be more specific, the correlation of two estimates is invariant (except for a possible sign change) up to a first order approximation, to smooth transformations of the estimates. There is a sign change if exactly one of the transformations is decreasing in a neighborhood of its parameter. In addition, we approximate the variance, covariance and correlation between functions of sample means and moments.     

A class of integral operators generated by random variables

H. Kres Statistisehe Tafeln zur multlvariaten Analysis. Springer-Verlag, Berlin- Heidel-berg-New York 1975, XVIII, 431 S., 26 Tab., DM 48.

D. Rasch: Einführung in die mathematische Statistik - WahrscheinUcllkeitsrechnung und Grundlagen der mathematlsehan Statistlk. VEB Deutscher Verlag delr Wissenschaften, Berlin 1976, 371 S., 37 Abb., 46 'I'ab., 40,– M.

D. Rasch: Einführung in die muthematisehe Statlstik - II .Anweuduugen, VEB Deutscher Verlag der Wissenschaften, Berlin 1976.

Donald L. Snyder: Random Point Processes. -JohnWiley &; Sons, New York 1975,485 S.     

Ridge Autoregression Estimation: LS Method

In an AR (p)-model, least-squares estimation of the parameters is considered when it is suspected that the parameters may belong to a linear subspace and the estimated covariance matrix is ill-conditioned. Accordingly, we define five estimators and study their properties in an asymptotic setup to discover dominance properties based on asymptotic distributional bias (ADB), MSE (ADMSE) matrices, and under quadratic risks (ADQR).     

On the Use of a Log-Rate Model for Survey-Weighted Categorical Data

For the analysis of survey-weighted categorical data, one recommended method of analysis is a log-rate model. For each cell in a contingency table, the survey weights are averaged across subjects and incorporated into an offset for a loglinear model. Supposedly, one can then proceed with the analysis of unweighted observed cell counts. We provide theoretical and simulation-based evidence to show that the log-rate analysis is not an effective statistical analysis method and should not be used in general. The root of the problem is in its failure to properly account for variability in the individual weights within cells of a contingency table. This results in goodness-of-fit tests that have higher-than-nominal error rates and confidence intervals for odds ratios that have lower-than-nominal coverage.     

Combining an Internal Pilot with an Interim Analysis for Single Degree of Freedom Tests

An internal pilot with interim analysis (IPIA) design combines interim power analysis (an internal pilot) with interim data analysis (two-stage group sequential). We provide IPIA methods for single df hypotheses within the Gaussian general linear model, including one and two group t tests. The design allows early stopping for efficacy and futility while also re-estimating sample size based on an interim variance estimate. Study planning in small samples requires the exact and computable forms reported here. The formulation gives fast and accurate calculations of power, Type I error rate, and expected sample size.     

Time series analysis of stratospheric ozone

Time series analyses of monthly average total ozone measured at 37 stations throughout the world were used to estimate the extent to which the average ozone trend correlates with the depletion curve hypothesized as due to chlorofluorocarbons (CFCs). Statistical characteristics of stations in the ensemble were used to help define appropriate model and station selection criteria. The maximum likelihood procedure developed herein estimates the weighted average trend, its. variance, and the intra- and inter-station variance components of the trend. Correlations among trends at different stations are also taken into account. The models were subjected to much checking and criticism. Variations in statistical methodology are used to show that the results are insensitive to details of the model selection criteria. The method does not discriminate well between the hypothesized CFC trend and a linear trend. The trend estimates represent the sum of all long-term global effects. The variance includes all effects that differ from station-to-station. The estimated trend and 2α limits for 14 stations with 20-year records (1958-79) is an ozone increase through 1979 of (1.5+1.0) percent. At the 23 stations with shorter records, the trend is (1.0=1.7) percent. It is concluded that no significant depletion in stratospheric ozone has occurred from any cause through the end of 1979.     

Pseudo maximum likelihood estimation for the dirichlet-multinomial distribution

Pseudo maximum likelihood estimation (PML) for the Dirich-let-multinomial distribution is proposed and examined in this pa-per. The procedure is compared to that based on moments (MM) for its asymptotic relative efficiency (ARE) relative to the maximum likelihood estimate (ML). It is found that PML, requiring much less computational effort than ML and possessing considerably higher ARE than MM, constitutes a good compromise between ML and MM. PML is also found to have very high ARE when an estimate for the scale parameter in the Dirichlet-multinomial distribution is all that is needed.     

Bayesian credible intervals for response surface optima

In response surface methodology, one is usually interested in estimating the optimal conditions based on a small number of experimental runs which are designed to optimally sample the experimental space. Typically, regression models are constructed from the experimental data and interrogated in order to provide a point estimate of the independent variable settings predicted to optimize the response. Unfortunately, these point estimates are rarely accompanied with uncertainty intervals. Though classical frequentist confidence intervals can be constructed for unconstrained quadratic models, higher order, constrained or nonlinear models are often encountered in practice. Existing techniques for constructing uncertainty estimates in such situations have not been implemented widely, due in part to the need to set adjustable parameters or because of limited or difficult applicability to constrained or nonlinear problems. To address these limitations a Bayesian method of determining credible intervals for response surface optima was developed. The approach shows good coverage probabilities on two test problems, is straightforward to implement and is readily applicable to the kind of constrained and/or nonlinear problems that frequently appear in practice.