Extreme value and cluster analysis of European daily temperature series

Time series of daily mean temperature obtained from the European Climate Assessment data set is analyzed with respect to their extremal properties. A time-series clustering approach which combines Bayesian methodology, extreme value theory and classification techniques is adopted for the analysis of the regional variability of temperature extremes. The daily mean temperature records are clustered on the basis of their corresponding predictive distributions for 25-, 50- and 100-year return values. The results of the cluster analysis show a clear distinction between the highest altitude stations, for which the return values are lowest, and the remaining stations. Furthermore, a clear distinction is also found between the northernmost stations in Scandinavia and the stations in central and southern Europe. This spatial structure of the return period distributions for 25-, 50- and 100-years seems to be consistent with projected changes in the variability of temperature extremes over Europe pointing to a different behavior in central Europe than in northern Europe and the Mediterranean area, possibly related to the effect of soil moisture and land-atmosphere coupling.     

A three-state continuous time Markov chain model for HIV disease burden

Plasma HIV viral load (VL) is the clinical indicator used to evaluate disease burden for HIV-infected patients. We developed a covariate-adjusted, three-state, homogenous continuous time Markov chain model for HIV/AIDS disease burden among subgroups. We defined Detectable and Undetectable HIV VL levels as two transient states and Death as the third absorbing state. We implemented the exact maximum likelihood method to estimate the parameters with related asymptotic distribution to conduct hypothesis testing. We evaluated the proposed model using HIV-infected individuals from South Carolina (SC) HIV surveillance data. Using the developed model, we estimated and compared the transition hazards, transition probabilities, and the state-specific duration for HIV-infected individuals. We examined gender, race/ethnicity, age, CD4 count, place of residence, and antiretroviral treatment regimen prescribed at the beginning of the study period. We found that patients with a higher CD4 count, increased age, heterosexual orientation, white, and single tablet regimen users were associated with reduced risk of transitioning to a Detectable VL from an Undetectable VL, whereas shorter time since diagnosis, being male, and injection drug use increased the risk of the same transition.     

Challenges and Opportunities of Open Access: A Panel Discussion

This year's panel discussion, moderated by Beverly Dowdy, offered three perspectives on the current and future trends in open access (OA) from David Crotty, a senior editor from a prominent university publisher; Beth Bernhardt, an assistant dean for collection development and scholarly communications at an academic library; and Dr. Paige Hall Smith, an associate professor at a large state university. Each discussed different opportunities and challenges that open access presents. A question-and-answer period followed, prompting further discussion of the points brought up during the presentations, as well as related issues.     

ORCID Pop-ups: Starting Small to Disseminate at Large

Abstract

In this lightning talk for the North Carolina Serials Conference, Lynnee Argabright describes the efforts of the ORCID team at the University of North Carolina at Chapel Hill. In addition to discussing the team’s pop-up outreach planning and documentation, Argabright shares advice from the team’s experience on how to run an effective ORCID pop-up event to facilitate faculty ORCID registration.     

Communities of Knowledge: Creating and Connecting Resource Metadata

A growing number of academic libraries and consortia are developing their own e-resource knowledge bases. The structure, format, and ultimate goals of their projects vary. If they are successful, these community-driven knowledge bases have the potential to make a positive impact on e-resource management and discovery. A situation analysis of the current state of affairs considers the problems these new communities of knowledge are attempting to solve, compares the scope and approaches of selected projects, suggests opportunities for connectivity across communities, and outlines challenges to their success. This article is the outgrowth of a 2013 North Carolina Serials Conference discussion.     

A Bayesian model for multiple change point to extremes,with application to environmental and financial data

Abrupt changes often occur for environmental and financial time series. Most often, these changes are due to human intervention. Change point analysis is a statistical tool used to analyze sudden changes in observations along the time series. In this paper, we propose a Bayesian model for extreme values for environmental and economic datasets that present a typical change point behavior. The model proposed in this paper addresses the situation in which more than one change point can occur in a time series. By analyzing maxima, the distribution of each regime is a generalized extreme value distribution. In this model, the change points are unknown and considered parameters to be estimated. Simulations of extremes with two change points showed that the proposed algorithm can recover the true values of the parameters, in addition to detecting the true change points in different configurations. Also, the number of change points was a problem to be considered, and the Bayesian estimation can correctly identify the correct number of change points for each application. Environmental and financial data were analyzed and results showed the importance of considering the change point in the data and revealed that this change of regime brought about an increase in the return levels, increasing the number of floods in cities around the rivers. Stock market levels showed the necessity of a model with three different regimes.     

Analyzing cross-sectional and longitudinal lung-function measurements: The effects of age

The problem of interpreting lung-function measurements in industrial workers is examined. The data under discussion pertain to FEV₁ and FVC measurements in smoking and in nonsmoking groups of grain-elevator workers in British Columbia and of workers in Vancouver City Hall. Initial observations have now been enriched by longitudinal follow up data on the same groups after three and after six years. It is shown that interesting selection phenomena, favouring “fit” individuals, take place over time, with regard both to lung symptoms and lung functions. Thus cross-sectional and longitudinal studies refer to somewhat different populations. It also appears that longitudinal studies are considerably more sensitive to identifying cumulative lung damage than are corresponding cross-sectional studies. The nonlinearity of the effect of age on lung functions is noted in the longitudinal data in a number of cases, lending support to the hypothesis of association between quadratic age effect and cumulative exposure to lung insults.     

Assessing safety at the end of clinical trials using system organ classes: A case and comparative study

Recent approaches to the statistical analysis of adverse event (AE) data in clinical trials have proposed the use of groupings of related AEs, such as by system organ class (SOC). These methods have opened up the possibility of scanning large numbers of AEs while controlling for multiple comparisons, making the comparative performance of the different methods in terms of AE detection and error rates of interest to investigators. We apply two Bayesian models and two procedures for controlling the false discovery rate (FDR), which use groupings of AEs, to real clinical trial safety data. We find that while the Bayesian models are appropriate for the full data set, the error controlling methods only give similar results to the Bayesian methods when low incidence AEs are removed. A simulation study is used to compare the relative performances of the methods. We investigate the differences between the methods over full trial data sets, and over data sets with low incidence AEs and SOCs removed. We find that while the removal of low incidence AEs increases the power of the error controlling procedures, the estimated power of the Bayesian methods remains relatively constant over all data sizes. Automatic removal of low-incidence AEs however does have an effect on the error rates of all the methods, and a clinically guided approach to their removal is needed. Overall we found that the Bayesian approaches are particularly useful for scanning the large amounts of AE data gathered.     

Me and My Shadow: Observation,Documentation, and Analysis of Serials and Electronic Resources Workflow

Librarians in the Metadata & Cataloging Department at the North Carolina State University (NCSU) Libraries embarked on a comprehensive project to document and review the Libraries' serials and electronic resources workflow with the goals of introducing greater efficiency, clarity, and simplicity into processes across the serials unit and of increasing training opportunities for department-wide understanding of the serials lifecycle. This paper examines the method used to observe and collect information about serials and electronic resources work, the process of synthesizing that information into a graphical representation of the serials lifecycle, and the workflow analysis undertaken to introduce improved serials processes. The value of the project is demonstrated through a discussion of efforts for workflow education, analysis, and improvement at the NCSU Libraries directly instigated by review of the process documentation by technical services staff members.     

Estimating the variance of estimated trends in proportions when there is no unique subject identifier

Summary.  Longitudinal population-based surveys are widely used in the health sciences to study patterns of change over time. In many of these data sets unique patient identifiers are not publicly available, making it impossible to link the repeated measures from the same individual directly. This poses a statistical challenge for making inferences about time trends because repeated measures from the same individual are likely to be positively correlated, i.e., although the time trend that is estimated under the naïve assumption of independence is unbiased, an unbiased estimate of the variance cannot be obtained without knowledge of the subject identifiers linking repeated measures over time. We propose a simple method for obtaining a conservative estimate of variability for making inferences about trends in proportions over time, ensuring that the type I error is no greater than the specified level. The method proposed is illustrated by using longitudinal data on diabetes hospitalization proportions in South Carolina.     

Me and My Shadow: Observation,Documentation, and Analysis of Serials and Electronic Resources Workflow

Abstract

Librarians in the Metadata & Cataloging Department at the North Carolina State University (NCSU) Libraries embarked on a comprehensive project to document and review the Libraries' serials and electronic resources workflow with the goals of introducing greater efficiency, clarity, and simplicity into processes across the serials unit and of increasing training opportunities for department-wide understanding of the serials lifecycle. This paper examines the method used to observe and collect information about serials and electronic resources work, the process of synthesizing that information into a graphical representation of the serials lifecycle, and the workflow analysis undertaken to introduce improved serials processes. The value of the project is demonstrated through a discussion of efforts for workflow education, analysis, and improvement at the NCSU Libraries directly instigated by review of the process documentation by technical services staff members.     

Little Magazine Interview Index 2003

Abstract

Reports on the 2005 Potomac Technical Processing Librarians Conference, Annapolis, Maryland, and the 25th Annual Charleston Conference in South Carolina.     

Overdispersed poisson regression models for studies of air pollution and human health

This paper presents results from a simulation study motivated by a recent study of the relationships between ambient levels of air pollution and human health in the community of Prince George, British Columbia. The simulation study was designed to evaluate the performance of methods based on overdispersed Poisson regression models for the analysis of series of count data. Aspects addressed include estimation of the dispersion parameter, estimation of regression coefficients and their standard errors, and the performance of model selection tests. The effects of varying amounts of overdispersion and differing underlying variance structure on this performance were of particular interest. This study is related to work reported by Breslow (1990) although the context is quite different. Preliminary work led to the conclusion that estimation of the dispersion parameter should be based on Pearson's chi-square statistic rather than the Poisson deviance. Regression coefficients are well estimated, even in the présence of substantial overdispersion and when the model for the variance function is incorrectly specified. Despite potential greater variability, the empirical estimator of the covariance matrix is preferred because the model-based estimator is unreliable in general. When the model for the variance function is incorrect, model-based test statistics may perform poorly, in sharp contrast to empirical test statistics, which performed very well in this study.     

Alternating event processes during lifetimes: population dynamics and statistical inference

In the literature studying recurrent event data, a large amount of work has been focused on univariate recurrent event processes where the occurrence of each event is treated as a single point in time. There are many applications, however, in which univariate recurrent events are insufficient to characterize the feature of the process because patients experience nontrivial durations associated with each event. This results in an alternating event process where the disease status of a patient alternates between exacerbations and remissions. In this paper, we consider the dynamics of a chronic disease and its associated exacerbation-remission process over two time scales: calendar time and time-since-onset. In particular, over calendar time, we explore population dynamics and the relationship between incidence, prevalence and duration for such alternating event processes. We provide nonparametric estimation techniques for characteristic quantities of the process. In some settings, exacerbation processes are observed from an onset time until death; to account for the relationship between the survival and alternating event processes, nonparametric approaches are developed for estimating exacerbation process over lifetime. By understanding the population dynamics and within-process structure, the paper provide a new and general way to study alternating event processes.     

IMPROVING POPULATION-SPECIFIC ALLELE FREQUENCY ESTIMATES BY ADAPTING SUPPLEMENTAL DATA: AN EMPIRICAL BAYES APPROACH

Estimation of the allele frequency at genetic markers is a key ingredient in biological and biomedical research, such as studies of human genetic variation or of the genetic etiology of heritable traits. As genetic data becomes increasingly available, investigators face a dilemma: when should data from other studies and population subgroups be pooled with the primary data? Pooling additional samples will generally reduce the variance of the frequency estimates; however, used inappropriately, pooled estimates can be severely biased due to population stratification. Because of this potential bias, most investigators avoid pooling, even for samples with the same ethnic background and residing on the same continent. Here, we propose an empirical Bayes approach for estimating allele frequencies of single nucleotide polymorphisms. This procedure adaptively incorporates genotypes from related samples, so that more similar samples have a greater influence on the estimates. In every example we have considered, our estimator achieves a mean squared error (MSE) that is smaller than either pooling or not, and sometimes substantially improves over both extremes. The bias introduced is small, as is shown by a simulation study that is carefully matched to a real data example. Our method is particularly useful when small groups of individuals are genotyped at a large number of markers, a situation we are likely to encounter in a genome-wide association study.     

Highs and lows: Some properties of the extremes of a diffusion and applications in finance

Observations on security prices, currency exchange rates, interest rates, and other financial time series usually include not only an open and close, but also a high and low price for the period. For Brown‐ian motion and certain diffusion processes, the information on high and low prices is of considerable value, particularly for estimating volatility, correlations between processes, and in the pricing of look‐back and barrier options. For pricing more general derivatives, this information is useful to the extent that change in volatility is an important ingredient in the price. The author gives a simple geometric device for generating the extremes of Brownian motion, and geometric Brownian motion; he then uses these extremes in the estimation of the volatility of the processes and to study survivorship bias.     

Regression Parameter Estimation from Panel Counts

This paper considers a study where each subject may experience multiple occurrences of an event and the rate of the event occurrences is of primary interest. Specifically, we are concerned with the situations where, for each subject, there are only records of the accumulated counts for the event occurrences at a finite number of time points over the study period. Sets of observation times may vary from subject to subject and differ between groups. We model the mean of the event occurrence number over time semiparametrically, and estimate the regression parameter. The proposed estimation procedures are illustrated with data from a bladder cancer study ( Byar, 1980 ). Both asymptotics and simulation studies on the estimators are presented.     

Modeling of rates over a hierarchical health administrative structure

Identifying the distribution of the incidence rate of a disease over a region is a prediction problem where area‐specific random effects are to be estimated. The authors consider the inclusion of such effects at different levels of a hierarchical health administrative structure. They develop inference procedures for this type of multi‐level model and show that the predicted rates are approximately weighted sums of the crude rates obtained by pooling data on each level of the hierarchy. Their techniques are illustrated using infant mortality data from British Columbia.     

Estimating the probability of simultaneous rainfall extremes within a region: a spatial approach

In this paper we investigate the impact of model mis-specification, in terms of the dependence structure in the extremes of a spatial process, on the estimation of key quantities that are of interest to hydrologists and engineers. For example, it is often the case that severe flooding occurs as a result of the observation of rainfall extremes at several locations in a region simultaneously. Thus, practitioners might be interested in estimates of the joint exceedance probability of some high levels across these locations. It is likely that there will be spatial dependence present between the extremes, and this should be properly accounted for when estimating such probabilities. We compare the use of standard models from the geostatistics literature with max-stables models from extreme value theory. We find that, in some situations, using an incorrect spatial model for our extremes results in a significant under-estimation of these probabilities which – in flood defence terms – could lead to substantial under-protection.