Out-of-Sample Performance of Discrete-Time Spot Interest Rate Models

We provide a comprehensive analysis of the out-of-sample performance of a wide variety of spot rate models in forecasting the probability density of future interest rates. Although the most parsimonious models perform best in forecasting the conditional mean of many financial time series, we find that the spot rate models that incorporate conditional heteroscedasticity and excess kurtosis or heavy tails have better density forecasts. Generalized autoregressive conditional heteroscedasticity significantly improves the modeling of the conditional variance and kurtosis, whereas regime switching and jumps improve the modeling of the marginal density of interest rates. Our analysis shows that the sophisticated spot rate models in the existing literature are important for applications involving density forecasts of interest rates.     

Data and projections of HIV and AIDS in Portugal

Projections of AIDS incidence are critical for assessing future healthcare needs. This paper focuses on the method of back-calculation for obtaining forecasts. The first problem faced was the need to account for delays and underreporting in reporting of cases and to adjust the incidence data. The method used to estimate the reporting delay distribution is based on Poisson regression and involves cross-classifying each reported case by calendar time of diagnosis and reporting delay. The adjusted AIDS incidence data are then used to obtain short-term projections and lower bounds on the size of the AIDS epidemic. The estimation procedure 'back-calculates' from AIDS incidence data using the incubation period distribution to obtain estimates of the numbers previously infected. These numbers are then projected forward. The problem can be shown to reduce to estimating the size of a multinomial population. The expectation-maximization (EM) algorithm is used to obtain maximum-likelihood estimates when the density of infection times is parametrized as a step function. The methodology is applied to AIDS incidence data in Portugal for four different transmission categories: injecting drug users, sexual transmission (homosexual/bisexual and heterosexual contact) and other, mainly haemophilia and blood transfusion related, to obtain short-term projections and an estimate of the minimum size of the epidemic.     

Stochastic population forecasts based on conditional expert opinions

The paper develops and applies an expert-based stochastic population forecasting method, which can also be used to obtain a probabilistic version of scenario-based official forecasts. The full probability distribution of population forecasts is specified by starting from expert opinions on the future development of demographic components. Expert opinions are elicited as conditional on the realization of scenarios, in a two-step (or multiple-step) fashion. The method is applied to develop a stochastic forecast for the Italian population, starting from official scenarios from the Italian National Statistical Office.     

Coherent mortality forecasting by the weighted multilevel functional principal component approach

In human mortality modelling, if a population consists of several subpopulations it can be desirable to model their mortality rates simultaneously while taking into account the heterogeneity among them. The mortality forecasting methods tend to result in divergent forecasts for subpopulations when independence is assumed. However, under closely related social, economic and biological backgrounds, mortality patterns of these subpopulations are expected to be non-divergent in the future. In this article, we propose a new method for coherent modelling and forecasting of mortality rates for multiple subpopulations, in the sense of nondivergent life expectancy among subpopulations. The mortality rates of subpopulations are treated as multilevel functional data and a weighted multilevel functional principal component (wMFPCA) approach is proposed to model and forecast them. The proposed model is applied to sex-specific data for nine developed countries, and the results show that, in terms of overall forecasting accuracy, the model outperforms the independent model and the Product-Ratio model as well as the unweighted multilevel functional principal component approach.     

基于均衡关系的中国人口死亡率预测模型

人口死亡率反映了人口的死亡程度,准确预测死亡率是人口科学及人口经济学研究的重点之一,同时也是长寿风险测量的重要数据基础。基于Lee-Carter模型,探索中国大陆与台湾地区死亡率的相关性,通过协整分析考虑两地死亡率的长期均衡关系,创新性地建立基于相关性的向量误差修正模型(VECM),克服传统自回归移动平均模型(ARIMA)使用有限数据进行预测的局限性;均方预测误差作为检验标准,结果表明:基于VECM模型的预测效果比传统的预测效果更佳;基于中国大陆地区和台湾地区的死亡率长期均衡关系,可以为两地联合长寿债券的定价提供重要参考。     

On future household structure

Summary.  We develop a method for computing probabilistic household forecasts which quantifies uncertainty in the future number of households of various types in a country. A probabilistic household forecast helps policy makers, planners and other forecast users in the fields of housing, energy, social security etc. in taking appropriate decisions, because some household variables are more uncertain than others. Deterministic forecasts traditionally do not quantify uncertainty. We apply the method to data from Norway. We find that predictions of future numbers of married couples, cohabiting couples and one-person households are more certain than those of lone parents and other private households. Our method builds on an existing method for computing probabilistic population forecasts, combining such a forecast with a random breakdown of the population according to household position (single, cohabiting, living with a spouse, living alone etc.). In this application, uncertainty in the total numbers of households of different types derives primarily from random shares, rather than uncertain future population size. A similar method could be applied to obtain probabilistic forecasts for other divisions of the population, such as household size, health or disability status, region of residence and labour market status.     

A Bayesian approach to mixture cure models with spatial frailties for population‐based cancer relative survival data

As the treatments of cancer progress, a certain number of cancers are curable if diagnosed early. In population‐based cancer survival studies, cure is said to occur when mortality rate of the cancer patients returns to the same level as that expected for the general cancer‐free population. The estimates of cure fraction are of interest to both cancer patients and health policy makers. Mixture cure models have been widely used because the model is easy to interpret by separating the patients into two distinct groups. Usually parametric models are assumed for the latent distribution for the uncured patients. The estimation of cure fraction from the mixture cure model may be sensitive to misspecification of latent distribution. We propose a Bayesian approach to mixture cure model for population‐based cancer survival data, which can be extended to county‐level cancer survival data. Instead of modeling the latent distribution by a fixed parametric distribution, we use a finite mixture of the union of the lognormal, loglogistic, and Weibull distributions. The parameters are estimated using the Markov chain Monte Carlo method. Simulation study shows that the Bayesian method using a finite mixture latent distribution provides robust inference of parameter estimates. The proposed Bayesian method is applied to relative survival data for colon cancer patients from the Surveillance, Epidemiology, and End Results (SEER) Program to estimate the cure fractions. The Canadian Journal of Statistics 40: 40–54; 2012 © 2012 Statistical Society of Canada     

Scenarios, Uncertainty and Conditional Forecasts of the World Population

Current official population forecasts differ little from those that Whelpton made 50 years ago either in the cohort–component methodology used or in the arguments used to motivate the assumptions. However, Whelpton produced some of the most erroneous forecasts of this century. This suggests that current forecasters should ensure that they give users an assessment of the uncertainty of their forecasts. We show how simple statistical methods can be combined with expert judgment to arrive at an overall predictive distribution for the future population. We apply the methods to a world population forecast that was made in 1994. Accepting that point forecast, we find that the probability is only about 2% that the world population in the year 2030 will be less than the low scenario of 8317 million. The probability that the world population will exceed the high scenario of 10 736 million is about 13%. Similarly, the probability is only about 51% that the high–low interval of a recent United Nations (UN) forecast will contain the true population in the year 2025. Even if we consider the UN high–low intervals as conditional on the possible future policies of its member states, they appear to have a relatively small probability of encompassing the future population.     

Multilevel modelling of the geographical distributions of diseases

Multilevel modelling is used on problems arising from the analysis of spatially distributed health data. We use three applications to demonstrate the use of multilevel modelling in this area. The first concerns small area all-cause mortality rates from Glasgow where spatial autocorrelation between residuals is examined. The second analysis is of prostate cancer cases in Scottish counties where we use a range of models to examine whether the incidence is higher in more rural areas. The third develops a multiple-cause model in which deaths from cancer and cardiovascular disease in Glasgow are examined simultaneously in a spatial model. We discuss some of the issues surrounding the use of complex spatial models and the potential for future developments.     

Performance of control charts for autoregressive conditional heteroscedastic processes

This paper examines the robustness of control schemes to data conditional heteroscedasticity. Overall, the results show that the control schemes which do not account for heteroscedasticity fail in providing reliable information on the status of the process. Consequently, incorrect conclusions will be drawn by applying these procedures in the presence of data conditional heteroscedasticity. Control charts with time-varying control limits are shown to be useful in that context.     

A Note on Forecaster Discord and Consensus Prediction Error

In using survey data, the cross-sectional dispersion of forecasts is commonly used as a proxy for uncertainty. This article proposes a test for the validity of this practice that is a modification of Amemiya's (1977) procedure for estimating heteroscedasticity functions. To deal with the common practical problem of overlapping forecasts, a moving average specification is employed. The technique is applied to the Livingston forecasts of several variables, and the results suggest a significant positive relationship between cross-forecaster dispersion (a measure of forecaster discord) and the variance of ex post prediction errors.     

Uncertain population forecasting

"Errors in population forecasts arise from errors in the jump-off population and errors in the predictions of future vital rates. The propagation of these errors through the linear (Leslie) growth model is studied, and prediction intervals for future population are developed. For U.S. national forecasts, the prediction intervals are compared with the U.S. Census Bureau's high-low intervals." In order to assess the accuracy of the predictions of vital rates, the authors "derive the predictions from a parametric statistical model and estimate the extent of model misspecification and errors in parameter estimates. Subjective, expert opinion, so important in real forecasting, is incorporated with the technique of mixed estimation. A robust regression model is used to assess the effects of model misspecification."     

Forecasting murder within a population of probationers and parolees: a high stakes application of statistical learning

Summary.  Forecasts of future dangerousness are often used to inform the sentencing decisions of convicted offenders. For individuals who are sentenced to probation or paroled to community supervision, such forecasts affect the conditions under which they are to be supervised. The statistical criterion for these forecasts is commonly called recidivism, which is defined as a charge or conviction for any new offence, no matter how minor. Only rarely do such forecasts make distinctions on the basis of the seriousness of offences. Yet seriousness may be central to public concerns, and judges are increasingly required by law and sentencing guidelines to make assessments of seriousness. At the very least, information about seriousness is essential for allocating scarce resources for community supervision of convicted offenders. The paper focuses only on murderous conduct by individuals on probation or parole. Using data on a population of over 60000 cases from Philadelphia's Adult Probation and Parole Department, we forecast whether each offender will be charged with a homicide or attempted homicide within 2 years of beginning community supervision. We use a statistical learning approach that makes no assumptions about how predictors are related to the outcome. We also build in the costs of false negative and false positive charges and use half of the data to build the forecasting model, and the other half of the data to evaluate the quality of the forecasts. Forecasts that are based on this approach offer the possibility of concentrating rehabilitation, treatment and surveillance resources on a small subset of convicted offenders who may be in greatest need, and who pose the greatest risk to society.     

Error models for official mortality forecasts

"The Office of the Actuary, U.S. Social Security Administration, produces alternative forecasts of mortality to reflect uncertainty about the future.... In this article we identify the components and assumptions of the official forecasts and approximate them by stochastic parametric models. We estimate parameters of the models from past data, derive statistical intervals for the forecasts, and compare them with the official high-low intervals. We use the models to evaluate the forecasts rather than to develop different predictions of the future. Analysis of data from 1972 to 1985 shows that the official intervals for mortality forecasts for males or females aged 45-70 have approximately a 95% chance of including the true mortality rate in any year. For other ages the chances are much less than 95%."     

Estimation and interpretation of models of absolute risk from epidemiologic data,including family-based studies

Absolute risk is the chance that a person with given risk factors and free of the disease of interest at age a will be diagnosed with that disease in the interval (a, a + τ]. Absolute risk is sometimes called cumulative incidence. Absolute risk is a “crude” risk because it is reduced by the chance that the person will die of competing causes of death before developing the disease of interest. Cohort studies admit flexibility in modeling absolute risk, either by allowing covariates to affect the cause-specific relative hazards or to affect the absolute risk itself. An advantage of cause-specific relative risk models is that various data sources can be used to fit the required components. For example, case–control data can be used to estimate relative risk and attributable risk, and these can be combined with registry data on age-specific composite hazard rates for the disease of interest and with national data on competing hazards of mortality to estimate absolute risk. Family-based designs, such as the kin-cohort design and collections of pedigrees with multiple affected individuals can be used to estimate the genotype-specific hazard of disease. Such analyses must be adjusted for ascertainment, and failure to take into account residual familial risk, such as might be induced by unmeasured genetic variants or by unmeasured behavioral or environmental exposures that are correlated within families, can lead to overestimates of mutation-specific absolute risk in the general population.     

Empirical Bayes procedures for stabilizing maps of U.S. cancer mortality rates

"The geographic mapping of age-standardized, cause-specific death rates is a powerful tool for identifying possible etiologic factors, because the spatial distribution of mortality risks can be examined for correlations with the spatial distribution of disease-specific risk factors. This article presents a two-stage empirical Bayes procedure for calculating age-standardized cancer death rates, for use in mapping, which are adjusted for the stochasticity of rates in small area populations. Using the adjusted rates helps isolate and identify spatial patterns in the rates. The model is applied to sex-specific data on U.S. county cancer mortality in the white population for 15 cancer sites for three decades: 1950-1959, 1960-1969, and 1970-1979. Selected results are presented as maps of county death rates for white males."     

Temporal disaggregation and restricted forecasting of multiple population time series

This article presents some applications of time-series procedures to solve two typical problems that arise when analyzing demographic information in developing countries: (1) unavailability of annual time series of population growth rates (PGRs) and their corresponding population time series and (2) inappropriately defined population growth goals in official population programs. These problems are considered as situations that require combining information of population time series. Firstly, we suggest the use of temporal disaggregation techniques to combine census data with vital statistics information in order to estimate annual PGRs. Secondly, we apply multiple restricted forecasting to combine the official targets on future PGRs with the disaggregated series. Then, we propose a mechanism to evaluate the compatibility of the demographic goals with the annual data. We apply the aforementioned procedures to data of the Mexico City Metropolitan Zone divided by concentric rings and conclude that the targets established in the official program are not feasible. Hence, we derive future PGRs that are both in line with the official targets and with the historical demographic behavior. We conclude that growth population programs should be based on this kind of analysis to be supported empirically. So, through specialized multivariate time-series techniques, we propose to obtain first an optimal estimate of a disaggregate vector of population time series and then, produce restricted forecasts in agreement with some data-based population policies here derived.     

A marginalized diffusion model for estimating age at first lower endoscopy use from current status data

Summary.  We propose an approach for estimating the age at first lower endoscopy examination from current status data that were collected via two series of cross-sectional surveys. To model the national probability of ever having a lower endoscopy examination, we incorporate birth cohort effects into a mixed influence diffusion model. We link a state-specific model to the national level diffusion model by using a marginalized modelling approach. In future research, results from our model will be used as microsimulation model inputs to estimate the contribution of endoscopy examinations to observed changes in colorectal cancer incidence and mortality.     

Bayesian Random Segmentation Models to Identify Shared Copy Number Aberrations for Array CGH Data

Array-based comparative genomic hybridization (aCGH) is a high-resolution high-throughput technique for studying the genetic basis of cancer. The resulting data consists of log fluorescence ratios as a function of the genomic DNA location and provides a cytogenetic representation of the relative DNA copy number variation. Analysis of such data typically involves estimation of the underlying copy number state at each location and segmenting regions of DNA with similar copy number states. Most current methods proceed by modeling a single sample/array at a time, and thus fail to borrow strength across multiple samples to infer shared regions of copy number aberrations. We propose a hierarchical Bayesian random segmentation approach for modeling aCGH data that utilizes information across arrays from a common population to yield segments of shared copy number changes. These changes characterize the underlying population and allow us to compare different population aCGH profiles to assess which regions of the genome have differential alterations. Our method, referred to as BDSAcgh (Bayesian Detection of Shared Aberrations in aCGH), is based on a unified Bayesian hierarchical model that allows us to obtain probabilities of alteration states as well as probabilities of differential alteration that correspond to local false discovery rates. We evaluate the operating characteristics of our method via simulations and an application using a lung cancer aCGH data set.