On the Problem of Identification in Multiplicative Intensity-Rate Models with Multiple Interactions

In this paper we examine a multiplicative intensity model in which a covariate interacts with two other covariates in the same model. We demonstrate, analytically, that in such situations a log-linear parameterization based on two pairs of baseline levels cannot be transformed, uniquely, to the, otherwise equivalent, multiplicative parameterization. We show that the problem lies in an oversight of the conditional independence between the two covariates interacting with a common third covariate. As a solution, therefore, we propose an approach that takes due account of such dependence. Our proposed approach uses a common baseline level for the three covariates involved in interaction while estimating the corresponding relative intensities. The issues addressed are illustrated with a demographic data set involving the estimation of rates of transition to parenthood.     

Investigating design for survival models

The aim of this paper is to derive methodology for designing ‘time to event’ type experiments. In comparison to estimation, design aspects of ‘time to event’ experiments have received relatively little attention. We show that gains in efficiency of estimators of parameters and use of experimental material can be made using optimal design theory. The types of models considered include classical failure data and accelerated testing situations, and frailty models, each involving covariates which influence the outcome. The objective is to construct an optimal design based of the values of the covariates and associated model or indeed a candidate set of models. We consider D-optimality and create compound optimality criteria to derive optimal designs for multi-objective situations which, for example, focus on the number of failures as well as the estimation of parameters. The approach is motivated and demonstrated using common failure/survival models, for example, the Weibull distribution, product assessment and frailty models.     

On the Practice of Rescaling Covariates

Whether doing parametric or nonparametric regression with shrinkage, thresholding, penalized likelihood, Bayesian posterior estimators (e.g., ridge regression, lasso, principal component regression, waveshrink or Markov random field ), it is common practice to rescale covariates by dividing by their respective standard errors ρ. The stated goal of this operation is to provide unitless covariates to compare like with like, especially when penalized likelihood or prior distributions are used. We contend that this vision is too simplistic. Instead, we propose to take into account a more essential component of the structure of the regression matrix by rescaling the covariates based on the diagonal elements of the covariance matrix Σ of the maximum-likelihood estimator. We illustrate the differences between the standard ρ- and proposed Σ-rescalings with various estimators and data sets.     

Generalized smooth finite mixtures

We propose a general class of models and a unified Bayesian inference methodology for flexibly estimating the density of a response variable conditional on a possibly high-dimensional set of covariates. Our model is a finite mixture of component models with covariate-dependent mixing weights. The component densities can belong to any parametric family, with each model parameter being a deterministic function of covariates through a link function. Our MCMC methodology allows for Bayesian variable selection among the covariates in the mixture components and in the mixing weights. The model’s parameterization and variable selection prior are chosen to prevent overfitting. We use simulated and real data sets to illustrate the methodology.     

Methods for multicountry studies of corporate governance: Evidence from the BRIKT countries

We discuss empirical challenges in multicountry studies of the effects of firm-level corporate governance on firm value, focusing on emerging markets. We assess the severe data, “construct validity”, and endogeneity issues in these studies, propose methods to respond to those issues, and apply those methods to a study of five major emerging markets—Brazil, India, Korea, Russia, and Turkey. We develop unique time-series datasets on governance in each country. We address construct validity by building country-specific indices which reflect local norms and institutions. These similar-but-not-identical indices predict firm market value in each country, and when pooled across countries, in firm fixed-effects (FE) and random-effects (RE) regressions. In contrast, a “common index”, which uses the same elements in each country, has no predictive power in FE regressions. For the country-specific and pooled indices, FE and RE coefficients on governance are generally lower than in pooled OLS regressions, and coefficients with extensive covariates are generally lower than with limited covariates. These results confirm the value of using FE or RE with extensive covariates to reduce omitted variable bias. We develop lower bounds on our estimates which reflect potential remaining omitted variable bias.     

Isolating the Roles of Individual Covariates in Reweighting Estimation

A host of recent research has used reweighting methods to analyze the extent to which observable characteristics predict between‐group differences in the distribution of an outcome. Less attention has been paid to using reweighting methods to isolate the roles of individual covariates. We analyze two approaches that have been used in previous studies, and we propose a new approach that examines the role of one covariate while holding the marginal distribution of the other covariates constant. We illustrate the differences between the methods with a numerical example and an empirical analysis of black–white wage differentials among males. Copyright © 2015 John Wiley & Sons, Ltd.     

Distributional Effects of a Continuous Treatment with an Application on Intergenerational Mobility

This paper considers the effect of a continuous treatment on the entire distribution of outcomes after adjusting for differences in the distribution of covariates across different levels of the treatment. Our methodology encompasses dose-response functions, counterfactual distributions, and ‘distributional policy effects’ depending on the assumptions invoked by the researcher. We propose a three-step estimator that consists of (i) estimating the distribution of the outcome conditional on the treatment and other covariates using quantile regression; (ii) for each value of the treatment, averaging over a counterfactual distribution of the covariates holding the treatment fixed; (iii) converting the resulting counterfactual distribution into parameters of interest that are easy to interpret. We show that our estimators converge uniformly to Gaussian processes and that the empirical bootstrap can be used to conduct uniformly valid inference across a range of values of the treatment. We use our method to study intergenerational income mobility where we consider effects of parents’ income on features of their child's income distribution such as (i) the fraction of children with income below the poverty line; (ii) the variance of child's income; and (iii) the inter-quantile range of child's income–all as a function of parents’ income.     

Semiparametric estimation of logistic regression model with missing covariates and outcome

We consider a semiparametric method to estimate logistic regression models with missing both covariates and an outcome variable, and propose two new estimators. The first, which is based solely on the validation set, is an extension of the validation likelihood estimator of Breslow and Cain (Biometrika 75:11–20, 1988). The second is a joint conditional likelihood estimator based on the validation and non-validation data sets. Both estimators are semiparametric as they do not require any model assumptions regarding the missing data mechanism nor the specification of the conditional distribution of the missing covariates given the observed covariates. The asymptotic distribution theory is developed under the assumption that all covariate variables are categorical. The finite-sample properties of the proposed estimators are investigated through simulation studies showing that the joint conditional likelihood estimator is the most efficient. A cable TV survey data set from Taiwan is used to illustrate the practical use of the proposed methodology.     

Weighted least squares estimators in possibly misspecified nonlinear regression

The behavior of estimators for misspecified parametric models has been well studied. We consider estimators for misspecified nonlinear regression models, with error and covariates possibly dependent. These models are described by specifying a parametric model for the conditional expectation of the response given the covariates. This is a parametric family of conditional constraints, which makes the model itself close to nonparametric. We study the behavior of weighted least squares estimators both when the regression function is correctly specified, and when it is misspecified and also involves possible additional covariates.     

What are the determinants of financial well-being? A Bayesian LASSO approach

The financial well-being (FWB) of individuals is a topic that is becoming increasingly important across a multitude of disciplines. In this study, we use the 2016 National Financial Well-Being Survey administered by the Consumer Financial Protection Bureau to assess the determinants of an individual's FWB. We identify 144 potential covariates that could explain variation in the FWB score of individuals. The statistical methodology of choice is the Bayesian LASSO, which is a covariate selection algorithm that also allows for the importance ranking of covariates. Out of the 144 potential covariates, we find that 26 have 95% credible intervals that do not contain zero. Broadly speaking, the results show that objective measures of financial competency and psychological and sociological factors contribute the bulk of the explanatory power that help explain an individual's FWB score.     

Doubly robust uniform confidence band for the conditional average treatment effect function

In this paper, we propose a doubly robust method to estimate the heterogeneity of the average treatment effect with respect to observed covariates of interest. We consider a situation where a large number of covariates are needed for identifying the average treatment effect but the covariates of interest for analyzing heterogeneity are of much lower dimension. Our proposed estimator is doubly robust and avoids the curse of dimensionality. We propose a uniform confidence band that is easy to compute, and we illustrate its usefulness via Monte Carlo experiments and an application to the effects of smoking on birth weights.     

Understanding the economic determinants of the severity of operational losses: A regularized generalized Pareto regression approach

We investigate a novel database of 10,217 extreme operational losses from the Italian bank UniCredit. Our goal is to shed light on the dependence between the severity distribution of these losses and a set of macroeconomic, financial, and firm‐specific factors. To do so, we use generalized Pareto regression techniques, where both the scale and shape parameters are assumed to be functions of these explanatory variables. We perform the selection of the relevant covariates with a state‐of‐the‐art penalized‐likelihood estimation procedure relying on L₁‐penalty terms. A simulation study indicates that this approach efficiently selects covariates of interest and tackles spurious regression issues encountered when dealing with integrated time series. Lastly, we illustrate the impact of different economic scenarios on the requested capital for operational risk. Our results have important implications in terms of risk management and regulatory policy.     

Hawkes process modeling of COVID-19 with mobility leading indicators and spatial covariates

Hawkes processes are used in statistical modeling for event clustering and causal inference, while they also can be viewed as stochastic versions of popular compartmental models used in epidemiology. Here we show how to develop accurate models of COVID-19 transmission using Hawkes processes with spatial-temporal covariates. We model the conditional intensity of new COVID-19 cases and deaths in the U.S. at the county level, estimating the dynamic reproduction number of the virus within an EM algorithm through a regression on Google mobility indices and demographic covariates in the maximization step. We validate the approach on both short-term and long-term forecasting tasks, showing that the Hawkes process outperforms several models currently used to track the pandemic, including an ensemble approach and an SEIR-variant. We also investigate which covariates and mobility indices are most important for building forecasts of COVID-19 in the U.S.     

Identification and estimation of nonlinear dynamic panel data models with unobserved covariates

This paper considers nonparametric identification of nonlinear dynamic models for panel data with unobserved covariates. Including such unobserved covariates may control for both the individual-specific unobserved heterogeneity and the endogeneity of the explanatory variables. Without specifying the distribution of the initial condition with the unobserved variables, we show that the models are nonparametrically identified from two periods of the dependent variable _Yit$Y_{it}$ and three periods of the covariate _Xit$X_{it}$. The main identifying assumptions include high-level injectivity restrictions and require that the evolution of the observed covariates depends on the unobserved covariates but not on the lagged dependent variable. We also propose a sieve maximum likelihood estimator (MLE) and focus on two classes of nonlinear dynamic panel data models, i.e., dynamic discrete choice models and dynamic censored models. We present the asymptotic properties of the sieve MLE and investigate the finite sample properties of these sieve-based estimators through a Monte Carlo study. An intertemporal female labor force participation model is estimated as an empirical illustration using a sample from the Panel Study of Income Dynamics (PSID).     

Three Corrected Score Tests for Generalized Linear Models with Dispersion Covariates

We develop three corrected score tests for generalized linear models with dispersion covariates, thus generalizing the results of Cordeiro , Ferrari and Paula (1993) and Cribari-Neto and Ferrari (1995) . We present, in matrix notation, general formulae for the coefficients which define the corrected statistics. The formulae only require simple operations on matrices and can be used to obtain analytically closed-form corrections for score test statistics in a variety of special generalized linear models with dispersion covariates. They also have advantages for numerical purposes since our formulae are readily computable using a language supporting numerical linear algebra. Two examples, namely, iid sampling without covariates on the mean or dispersion parameter oand one-way classification models, are given. We also present some simulations where the three corrected tests perform better than the usual score test, the likelihood ratio test and its Bartlett corrected version. Finally, we present a numerical example for a data set discussed by Simonoff and Tsai (1994) .     

Fisher information in record data

We present some general results on Fisher information (FI) contained in upper (or lower) record values and associated record times generated from a sequence of i.i.d. continuous variables. For the record data obtained from a random sample of fixed size, we establish an interesting relationship between its FI content and the FI in the data consisting of sequential maxima. We also consider the record data from an inverse sampling plan (Samaniego and Whitaker, 1986). We apply the general results to evaluate the FI in upper as well as lower records data from the exponential distribution for both sampling plans. Further, we discuss the implication of our results to statistical inference from these record data. Received: December 2001 Acknowledgements. This research was supported by Fondo Nacional de Desarrollo Cientifico y Tecnologico (FONDECYT) grants 7990089 and 1010222 of Chile. We would like to thank the Department of Statistics at the University of Concepción for its hospitality during the stay of H. N. Nagaraja in Chile in March of 2000, when the initial work was done. We are grateful to the two referees for various comments that let to improvements in the paper.     

Counterfactual decomposition of changes in wage distributions using quantile regression

We propose a method to decompose the changes in the wage distribution over a period of time in several factors contributing to those changes. The method is based on the estimation of marginal wage distributions consistent with a conditional distribution estimated by quantile regression as well as with any hypothesized distribution for the covariates. Comparing the marginal distributions implied by different distributions for the covariates, one is then able to perform counterfactual exercises. The proposed methodology enables the identification of the sources of the increased wage inequality observed in most countries. Specifically, it decomposes the changes in the wage distribution over a period of time into several factors contributing to those changes, namely by discriminating between changes in the characteristics of the working population and changes in the returns to these characteristics. We apply this methodology to Portuguese data for the period 1986–1995, and find that the observed increase in educational levels contributed decisively towards greater wage inequality. Copyright © 2005 John Wiley & Sons, Ltd.     

Analysis of a mixed model for ordinal data by iterative re-weighted REML

An estimation procedure will be presented for a class of threshold models for ordinal data. These models may include both fixed and random effects with associated components of variance on an underlying scale. The residual error distribution on the underlying scale may be rendered greater flexibility by introducing additional shape parameters, e.g. a kurtosis parameter or parameters to model heterogeneous residual variances as a function of factors and covariates. The estimation procedure is an extension of an iterative re-weighted restricted maximum likelihood procedure, originally developed for generalized linear mixed models. This procedure will be illustrated with a practical problem involving damage to potato tubers and with data from animal breeding and medical research from the literature.     

The two‐sample linear regression model with interval‐censored covariates

There are surveys that gather precise information on an outcome of interest, but measure continuous covariates by a discrete number of intervals, in which case the covariates are interval censored. For applications with a second independent dataset precisely measuring the covariates, but not the outcome, this paper introduces a semiparametrically efficient estimator for the coefficients in a linear regression model. The second sample serves to establish point identification. An empirical application investigating the relationship between income and body mass index illustrates the use of the estimator.     

Effect of factors with time-dependent covariates on the survival status of the elderly in Taiwan

The elderly population has been increasing rapidly in recent years because of improvements in medical care and the progress in economics in Taiwan. Now, the population of old people age above 65 has approached 10%, and problems of the elderly have become a major concern for public health. In this study, we use five waves of Survey of Health and Living Status of the Elderly in Taiwan held from 1989 to 2003 to explore 15 variables related to demographic characteristic, health status, health behavior and home condition etc. affect the survival status of the elderly by employing Cox proportional hazard model. The results show that there are nine variables, e.g. age, gender, ethnic group, ADL, self-rated health, physical function, smoking, chewing betel nuts and marital status, were strongly related to the survival status of the elderly. In addition, Aalen’s nonparametric additive model is not only used as an alternative to a proportional hazards to describe the effects of covariates on survival time but assist in detecting and describing the nature of time-dependent effects of covariates used for Cox proportional model. We find there are some time-dependent covariates (e.g. self-rate health, ADL function, and physical function), and Aalen’s nonparametric additive model provides a flexible and nonparametric method for investigating the time-dependent variables through plots of the estimated cumulative regression coefficients, with confidence intervals.