Modeling Clustered Ordered Categorical Data: A Survey

This article surveys various strategies for modeling ordered categorical (ordinal) response variables when the data have some type of clustering, extending a similar survey for binary data by Pendergast, Gange, Newton, Lindstrom, Palta & Fisher (1996). An important special case is when repeated measurement occurs at various occasions for each subject, such as in longitudinal studies. A much greater variety of models and fitting methods are available than when a similar survey for repeated ordinal response data was prepared a decade ago (Agresti, 1989). The primary emphasis of the review is on two classes of models, marginal models for which effects are averaged over all clusters at particular levels of predictors, and cluster-specific models for which effects apply at the cluster level. We present the two types of models in the ordinal context, review the literature for each, and discuss connections between them. Then, we summarize some alternative modeling approaches and ways of estimating parameters, including a Bayesian approach. We also discuss applications and areas likely to be popular for future research, such as ways of handling missing data and ways of modeling agreement and evaluating the accuracy of diagnostic tests. Finally, we review the current availability of software for using the methods discussed in this article.     

Semiparametric analysis of clustered survival data under nonparametric frailty

In this article a novel approach to analyze clustered survival data that are subject to extravariation encountered through clustering of survival times is proposed. This is accomplished by extending the Cox proportional hazard model to a frailty model where the cluster-specific shared frailty is modeled nonparametrically. We assume a nonparametric Dirichlet process for the distribution of frailty. In such a semiparametric setup, we propose a hybrid method to draw model-based inferences. In the framework of the proposed hybrid method, the estimation of parameters is performed by implementing Monte Carlo expected conditional maximization algorithm. A simulation study is conducted to study the efficiency of our methodology. The proposed methodology is, thereafter, illustrated by a real-life data on recurrence time to infections in kidney patient.     

Joint and marginal specification tests for conditional mean and variance models

This article proposes a class of joint and marginal spectral diagnostic tests for parametric conditional means and variances of linear and nonlinear time series models. The use of joint and marginal tests is motivated from the fact that marginal tests for the conditional variance may lead to misleading conclusions when the conditional mean is misspecified. The new tests are based on a generalized spectral approach and do not need to choose a lag order depending on the sample size or to smooth the data. Moreover, the proposed tests are robust to higher order dependence of unknown form, in particular to conditional skewness and kurtosis. It turns out that the asymptotic null distributions of the new tests depend on the data generating process. Hence, we implement the tests with the assistance of a wild bootstrap procedure. A simulation study compares the finite sample performance of the proposed and competing tests, and shows that our tests can play a valuable role in time series modeling. Finally, an application to the S&P 500 highlights the merits of our approach.     

Können Unternehmen durch Cluster Wettbewerbsvorteile realisieren? Eine ressourcenbasierte Analyse des Uhrenclusters Glashütte

Single firms as parts of clusters are in the focus of this paper. They are analysed with regard to their contribution to the generation of firm- and cluster-specific competitive advantages. The resource-based view is extended to ensure applicability in cluster research. This approach serves as the basis for the analysis of a case study undertaken in 2008 in the watch cluster in Glashütte, Germany.     

Nonnegativity constraints and intratemporal uncertainty in a multi-good life-cycle model

In the standard multi-good life-cycle consumption model (with intertemporal additive utility) the intratemporal relations between the marginal utilities of the different goods are deterministic. However, these deterministic identities will not usually be satisfied by the data. To avoid these deterministic relationships, we apply an approach which consists of introducing intratemporal uncertainty, and which is, in particular, interesting when additional nonnegativity constraints are present. We estimate some simple versions of the model with this so-called intratemporal uncertainty. The estimation results are, in general, in accordance with the theory, and most versions of the model are not rejected by Hansen and Singleton's misspecification test.     

Agricultural productivity,shadow wages and off-farm labor decisions in Nicaragua

The objective of this study is to analyze the decision to work off-farm by male and female farmers in Nicaragua using a three-year unbalanced panel dataset. Shadow income and shadow wages are also estimated. Moreover, to mitigate biases from unobserved individual and farm time-invariant characteristics as well as from sample selection, we apply a semiparametric approach for panel data. Our main findings suggest that shadow wages play a major role in off-farm labor decisions for both males and females. This implies that less labor is allocated to off-farm activities as the opportunity cost for agricultural work goes up. In addition, as the on-farm marginal productivity of households (i.e., shadow income) rises, both males and females reduce the hours they allocate to off-farm activities.     

Modeling non-normal corporate bond yield spreads by copula

This research focuses on modeling for how corporate bond yield spreads are affected by explanatory variables such as equity volatility, interest rate volatility, r, slope, rating, liquidity, coupon rate, and maturity. The existing literature assumes normality and linearity in the analysis, which is not the case in our sample. Thus, through a powerful and flexible copula approach, we study the dependence at the mean of the joint distribution by using the Gaussian copula marginal regression method and the dependence structure at the tails by using various copula functions. To our knowledge, this is the first application of the copula marginal regression model to bond market data. In addition, we employ several copula functions to test for the tail dependence between yield spreads and other explanatory variables. We find stronger tail dependence in the joint upper tail for the relation between equity volatility and yield spreads, among others. This result indicates the positive effect of equity volatility on yield spreads in the upper tail is greater than that in the low tail. This finding should be useful to practitioners, such as investors. By relying on better-fitting, more meaningful statistical models, this paper contributes to the extant literature on how corporate bond yield spreads are determined.     

On solving bias‐corrected non‐linear estimation equations with an application to the dynamic linear model

In a seminal paper, Mak, Journal of the Royal Statistical Society B, 55, 1993, 945, derived an efficient algorithm for solving non‐linear unbiased estimation equations. In this paper, we show that when Mak's algorithm is applied to biased estimation equations, it results in the estimates that would come from solving a bias‐corrected estimation equation, making it a consistent estimator if regularity conditions hold. In addition, the properties that Mak established for his algorithm also apply in the case of biased estimation equations but for estimates from the bias‐corrected equations. The marginal likelihood estimator is obtained when the approach is applied to both maximum likelihood and least squares estimation of the covariance matrix parameters in the general linear regression model. The new approach results in two new estimators when applied to the profile and marginal likelihood functions for estimating the lagged dependent variable coefficient in the dynamic linear regression model. Monte Carlo simulation results show the new approach leads to a better estimator when applied to the standard profile likelihood. It is therefore recommended for situations in which standard estimators are known to be biased.     

Bayesian estimation and model selection for spatial Durbin error model with finite distributed lags

In this paper we investigate a spatial Durbin error model with finite distributed lags and consider the Bayesian MCMC estimation of the model with a smoothness prior. We study also the corresponding Bayesian model selection procedure for the spatial Durbin error model, the spatial autoregressive model and the matrix exponential spatial specification model. We derive expressions of the marginal likelihood of the three models, which greatly simplify the model selection procedure. Simulation results suggest that the Bayesian estimates of high order spatial distributed lag coefficients are more precise than the maximum likelihood estimates. When the data is generated with a general declining pattern or a unimodal pattern for lag coefficients, the spatial Durbin error model can better capture the pattern than the SAR and the MESS models in most cases. We apply the procedure to study the effect of right to work (RTW) laws on manufacturing employment.     

Stick-breaking autoregressive processes

This paper considers the problem of defining a time-dependent nonparametric prior for use in Bayesian nonparametric modelling of time series. A recursive construction allows the definition of priors whose marginals have a general stick-breaking form. The processes with Poisson-Dirichlet and Dirichlet process marginals are investigated in some detail. We develop a general conditional Markov Chain Monte Carlo (MCMC) method for inference in the wide subclass of these models where the parameters of the marginal stick-breaking process are nondecreasing sequences. We derive a generalised Pólya urn scheme type representation of the Dirichlet process construction, which allows us to develop a marginal MCMC method for this case. We apply the proposed methods to financial data to develop a semi-parametric stochastic volatility model with a time-varying nonparametric returns distribution. Finally, we present two examples concerning the analysis of regional GDP and its growth.     

Testing a linear dynamic panel data model against nonlinear alternatives

The most popular econometric models in the panel data literature are the class of linear panel data models with unobserved individual- and/or time-specific effects. The consistency of parameter estimators and the validity of their economic interpretations as marginal effects depend crucially on the correct functional form specification of the linear panel data model. In this paper, a new class of residual-based tests is proposed for checking the validity of dynamic panel data models with both large cross-sectional units and time series dimensions. The individual and time effects can be fixed or random, and panel data can be balanced or unbalanced. The tests can detect a wide range of model misspecifications in the conditional mean of a dynamic panel data model, including functional form and lag misspecification. They check a large number of lags so that they can capture misspecification at any lag order asymptotically. No common alternative is assumed, thus allowing for heterogeneity in the degrees and directions of functional form misspecification across individuals. Thanks to the use of panel data with large $N$ and $T$, the proposed nonparametric tests have an asymptotic normal distribution under the null hypothesis without requiring the smoothing parameters to grow with the sample sizes. This suggests better nonparametric asymptotic approximation for the panel data than for time series or cross sectional data. This is confirmed in a simulation study. We apply the new tests to test linear specification of cross-country growth equations and found significant nonlinearities in mean for OECD countries’ growth equation for annual and quintannual panel data.     

Retirement Choices in Italy: What an Option Value Model Tells Us*

Using Italian data, we estimate an option value model to quantify the effect of financial incentives on retirement choices. As far as we know, this is the first empirical study to estimate the conditional multiple‐years model put forward by Stock and Wise (1990) . This implies that we account for dynamic self‐selection bias. We also present an extended version of this model in which the marginal value of leisure is random. For the female sample, the model is able to predict almost perfectly the age‐specific hazard rates. For the male sample, we obtain a good fit. Dynamic self‐selection results in a downward bias in the estimate of the marginal utility of leisure. We perform a simulation study to gauge the effects of a dramatic pension reform. Underestimation of the value of leisure translates into sizeable over‐prediction of the impact of reform. Due to lack of data, results for males should be interpreted with caution since we are not able to fully correct for dynamic self‐selection bias.     

A residual-based multivariate constant correlation test

We propose a new multivariate constant correlation test based on residuals. This test takes into account the whole correlation matrix instead of the considering merely marginal correlations between bivariate data series. In financial markets, it is unrealistic to assume that the marginal variances are constant. This motivates us to develop a constant correlation test which allows for non-constant marginal variances in multivariate time series. However, when the assumption of constant marginal variances is relaxed, it can be shown that the residual effect leads to nonstandard limit distributions of the test statistics based on residual terms. The critical values of the test statistics are not directly available and we use a bootstrap approximation to obtain the corresponding critical values for the test. We also derive the limit distribution of the test statistics based on residuals under the null hypothesis. Monte Carlo simulations show that the test has appealing size and power properties in finite samples. We also apply our test to the stock returns in Euro Stoxx 50 and integrate the test into a binary segmentation algorithm to detect multiple break points.     

An alternative quasi likelihood approach,Bayesian analysis and data-based inference for model specification

This paper studies an alternative quasi likelihood approach under possible model misspecification. We derive a filtered likelihood from a given quasi likelihood (QL), called a limited information quasi likelihood (LI-QL), that contains relevant but limited information on the data generation process. Our LI-QL approach, in one hand, extends robustness of the QL approach to inference problems for which the existing approach does not apply. Our study in this paper, on the other hand, builds a bridge between the classical and Bayesian approaches for statistical inference under possible model misspecification. We can establish a large sample correspondence between the classical QL approach and our LI-QL based Bayesian approach. An interesting finding is that the asymptotic distribution of an LI-QL based posterior and that of the corresponding quasi maximum likelihood estimator share the same “sandwich”-type second moment. Based on the LI-QL we can develop inference methods that are useful for practical applications under possible model misspecification. In particular, we can develop the Bayesian counterparts of classical QL methods that carry all the nice features of the latter studied in  White (1982). In addition, we can develop a Bayesian method for analyzing model specification based on an LI-QL.     

Predicting default risk under asymmetric binary link functions

In this article we propose the use of an asymmetric binary link function to extend the proportional hazard model for predicting loan default. The rationale behind this approach is that the symmetry assumption that has been widely used in the literature could be considered as quite restrictive, especially during periods of financial distress. In our approach we allow for a flexible level of asymmetry in the probability of default by the use of the skewed logit distribution. This enable us to estimate the actual level of asymmetry that is associated with the data at hand. We implement our approach to both simulated data and a rich micro dataset of consumer loan accounts. Our results provide clear evidence that ignoring the actual level of asymmetry leads to seriously biased estimates of the slope coefficients, inaccurate marginal effects of the covariates of the model, and overestimation of the probability of default. Regarding the predictive power of the covariates of the model, we have found that loan-specific covariates contain considerably more information about the loan default than macroeconomic covariates, which are often used in practice to carry out macroprudential stress testing.     

Modelling oil price expectations: Evidence from survey data

Using Consensus Forecast survey data on WTI oil price expectations for 3- and 12-month horizons over the period November 1989 to December 2008, we find that the rational expectation hypothesis is rejected and that none of the traditional extrapolative, regressive and adaptive processes fits the data by itself. We suggest a mixed expectation model defined as a linear combination of these traditional processes, which we interpret as the aggregation of individual mixing behavior and of heterogenous groups of agents using these simple processes. This approach is consistent with the economically rational expectations theory. We show that the target oil price included in the regressive component of this model depends on the long-run marginal cost of crude oil production and on short term macroeconomic fundamentals whose effects are subject to structural changes. For the two horizons, estimation results provide evidence for our mixed expectation model incorporating this break-dependent target price.     

List augmentation with model based multiple imputation: a case study using a mixed-outcome factor model

This study concerns list augmentation in direct marketing. List augmentation is a special case of missing data imputation. We review previous work on the mixed outcome factor model and apply it for the purpose of list augmentation. The model deals with both discrete and continuous variables and allows us to augment the data for all subjects in a company's transaction database with soft data collected in a survey among a sample of those subjects. We propose a bootstrap-based imputation approach, which is appealing to use in combination with the factor model, since it allows one to include estimation uncertainty in the imputation procedure in a simple, yet adequate manner. We provide an empirical case study of the performance of the approach to a transaction data base of a bank.     

Estimating household demand for housing attributes in rent-controlled markets

We estimate the household’s marginal willingness to pay for housing attributes in the rent-controlled sector, so where rents are not freely market determined. The application of hedonic price approaches to obtain estimates of the household’s value of housing characteristics is then invalid. We apply an alternative estimation approach based on residential mobility. In our application, we focus on the households’ willingness to pay for number of rooms as well as the willingness to pay to avoid a long commuting distance. Our estimates appear plausible. For example, for households in the rent-controlled sector are willing to pay about 7% of their household income for an additional room. The implied marginal costs of commuting are about €0.17–€0.23 per (one-way) kilometre.     

Modelling Strategies for Repeated Multiple Response Data

This article discusses modelling strategies for repeated measurements of multiple response variables. Such data arise in the context of categorical variables where one can select more than one of the categories as the response. We consider each of the multiple responses as a binary outcome and use a marginal (or population‐averaged) modelling approach to analyse its means. Generalized estimating equations are used to account for different correlation structures, both over time and between items. We also discuss an alternative approach using a generalized linear mixed model with conditional interpretations. We illustrate the methods using data from a panel study in Australia called the Household, Income, and Labour Dynamics Survey.     

Distribution-free tests of mean vectors and covariance matrices for multivariate paired data

We study a permutation procedure to test the equality of mean vectors, homogeneity of covariance matrices, or simultaneous equality of both mean vectors and covariance matrices in multivariate paired data. We propose to use two test statistics for the equality of mean vectors and the homogeneity of covariance matrices, respectively, and combine them to test the simultaneous equality of both mean vectors and covariance matrices. Since the combined test has composite null hypothesis, we control its type I error probability and theoretically prove the asymptotic unbiasedness and consistency of the combined test. The new procedure requires no structural assumption on the covariances. No distributional assumption is imposed on the data, except that the permutation test for mean vector equality assumes symmetric joint distribution of the paired data. We illustrate the good performance of the proposed approach with comparison to competing methods via simulations. We apply the proposed method to testing the symmetry of tooth size in a dental study and to finding differentially expressed gene sets with dependent structures in a microarray study of prostate cancer.