A model for extreme wind gusts

Estimates of the largest wind gust that will occur at a given location over a specified period are required by civil engineers. Estimation is usually based on models which are derived from the limiting distributions of maxima of stationary time series and which are fitted to data on extreme gusts. In this paper we develop a model for maximum gusts which also incorporates data on hourly mean speeds through a distributional relationship between maxima and means. This joint model is closely linked to the physical processes which generate the most extreme values and thus provides a mechanism by which data on means can augment those on gusts. It is argued that this increases the credibility of extrapolation in estimates of long period return gusts. The model is shown to provide a good fit to data obtained at a location in northern England and is compared with a more traditional modelling approach, which also performs well for this site.     

A hierarchical model for extreme wind speeds

Summary.  A typical extreme value analysis is often carried out on the basis of simplistic inferential procedures, though the data being analysed may be structurally complex. Here we develop a hierarchical model for hourly gust maximum wind speed data, which attempts to identify site and seasonal effects for the marginal densities of hourly maxima, as well as for the serial dependence at each location. A Gaussian model for the random effects exploits the meteorological structure in the data, enabling increased precision for inferences at individual sites and in individual seasons. The Bayesian framework that is adopted is also exploited to obtain predictive return level estimates at each site, which incorporate uncertainty due to model estimation, as well as the randomness that is inherent in the processes that are involved.     

A local moment type estimator for an extreme quantile in regression with random covariates

A conditional extreme quantile estimator is proposed in the presence of random covariates. It is based on an adaptation of the moment estimator introduced by Dekkers et al. (1989 Dekkers, A.L.M., Einmahl, J.H.J., de Haan, L. (1989). A moment estimator for the index of an extreme-value distribution. Ann. Statist. 17:1833–1855.[Crossref], [Web of Science ®] , [Google Scholar]) in the classical univariate setting, and thus it is valid in the domain of attraction of the extreme value distribution, i.e., whatever the sign of the extreme value index is. Asymptotic normality of the estimator is established under suitable assumptions, and its finite sample behavior is evaluated with a small simulation study, where a comparison with an alternative estimator already proposed in the literature is provided. An illustration to a real dataset concerning the world catalogue of earthquake magnitudes is also proposed.     

A generalized linear model approach to seasonal aspects of wind speed modeling

The aim of the article is to identify the intraday seasonality in a wind speed time series. Following the traditional approach, the marginal probability law is Weibull and, consequently, we consider seasonal Weibull law. A new estimation and decision procedure to estimate the seasonal Weibull law intraday scale parameter is presented. We will also give statistical decision-making tools to discard or not the trend parameter and to validate the seasonal model.     

Some remarks on experimental design for estimating the expectation of a stationary random process

For estimation of the mean of a stationary random process the variance-optimal choice of the observation points (the so called experimental design) is considered. For this discrete and continuous designs are introduced, some known results of process statistics are interpreted to experimental design and a proposal for simplification of the minimization problem is offered, moreover it is proved, that for monotone decreasing eovarianee functions a design, for which the points near the ends of the observation interval are more dense than in the middle, is better than the equidistant design.     

A note on the limiting properties of the least squares estimation for the random coefficient autoregressive model

This note is concerned with the limiting properties of the least squares estimation for the random coefficient autoregressive model. In contrast with existing results, ours is applicable to a wide range of models under more general assumptions.     

A power comparison of three tests for design effects in a random effects covariance model

A preliminary testing procedure for design ettecta in a ran-dom effects covariance model is Compared with the usual procedure to see if the power of the latter can be improved. A procedure which ignores the random covariate effects is included for comparison and for study of misspecification effects. Methodology is based on Roebruck's (1982) results for regular linear models.     

A method for estimating age-specific reference intervals ('normal ranges') based on fractional polynomials and exponential transformation

The age-specific reference interval is an important screening tool in medicine. Put crudely, an individual whose value of a variable of interest lies outside certain extreme centiles may be suspected of abnormality. We propose a parametric method for constructing such intervals. It provides smooth centile curves and explicit formulae for the centile estimates and for standard deviation (SD) scores (age-standardized values). Each parameter of an exponential–normal or modulus–exponential–normal density is modelled as a fractional polynomial function of age. Estimation is by maximum likelihood. These three- and four-parameter models involve transformations of the data towards normality which remove non-normal skewness and/or kurtosis. Fractional polynomials provide more flexible curve shapes than do conventional polynomials. The method easily accommodates binary covariates facilitating, for example, parsimonious modelling of age- and sex-specific centile curves. A method of calculating precision profiles for centile estimates is proposed. Goodness of fit is assessed by using Q–Q -plots and Shapiro–Wilk W -tests of the SD scores, and likelihood ratio tests of the parameters of an enlarged model. Four substantial real data sets are used to illustrate the method. Comparisons are made with the semiparametric LMS method of Cole and Green.     

Bayesian inference for generalized additive mixed models based on Markov random field priors

Most regression problems in practice require flexible semiparametric forms of the predictor for modelling the dependence of responses on covariates. Moreover, it is often necessary to add random effects accounting for overdispersion caused by unobserved heterogeneity or for correlation in longitudinal or spatial data. We present a unified approach for Bayesian inference via Markov chain Monte Carlo simulation in generalized additive and semiparametric mixed models. Different types of covariates, such as the usual covariates with fixed effects, metrical covariates with non-linear effects, unstructured random effects, trend and seasonal components in longitudinal data and spatial covariates, are all treated within the same general framework by assigning appropriate Markov random field priors with different forms and degrees of smoothness. We applied the approach in several case-studies and consulting cases, showing that the methods are also computationally feasible in problems with many covariates and large data sets. In this paper, we choose two typical applications.     

A modified nested-error regression model for small area estimation

A nested-error regression model having both fixed and random effects is introduced to estimate linear parameters of small areas. The model is applicable to data having a proportion of domains where the variable of interest cannot be described by a standard linear mixed model. Algorithms and formulas to fit the model, to calculate EBLUP and to estimate mean-squared errors are given. A Monte Carlo simulation experiment is presented to illustrate the gain of precision obtained by using the proposed model and to obtain some practical conclusions. A motivating application to Spanish Labour Force Survey data is also given.     

A new stratified three-stage unrelated randomized response model for estimating a rare sensitive attribute based on the Poisson distribution

This article suggests an efficient method of estimating a rare sensitive attribute which is assumed following Poisson distribution by using three-stage unrelated randomized response model instead of the Land et al. model (2011 Land, M., S. Singh, and S. A. Sedory. 2011. Estimation of a rare sensitive attribute using poisson distribution. Statistics 46 (3):351–60. doi:10.1080/02331888.2010.524300.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) when the population consists of some different sized clusters and clusters selected by probability proportional to size(:pps) sampling. A rare sensitive parameter is estimated by using pps sampling and equal probability two-stage sampling when the parameter of a rare unrelated attribute is assumed to be known and unknown.

We extend this method to the case of stratified population by applying stratified pps sampling and stratified equal probability two-stage sampling. An empirical study is carried out to show the efficiency of the two proposed methods when the parameter of a rare unrelated attribute is assumed to be known and unknown.