Theoretical perspectives of trade-off analysis using DEA

Because of the piecewise linear nature of the frontier in data envelopment analysis (DEA), estimated marginal rates of substitution are only valid for infinitesimal, or small finite, changes in one or more variables. Analysing the impacts of very small changes is not adequate for many situations, where the consequences of substantial changes in variables are of interest to for instance managers. This paper develops methods for evaluating larger, i.e. non-marginal, tradeoffs between variables in DEA. The methods are capable of handling both scalar and additive changes. Finally, the approaches for evaluating basic pairwise tradeoffs are generalised in order to enable assessment of the impact on one or more throughputs of the change in one or more of the other throughputs.     

An estimating equation for parametric shared frailty models with marginal additive hazards

Summary.  Multivariate failure time data arise when data consist of clusters in which the failure times may be dependent. A popular approach to such data is the marginal proportional hazards model with estimation under the working independence assumption. In some contexts, however, it may be more reasonable to use the marginal additive hazards model. We derive asymptotic properties of the Lin and Ying estimators for the marginal additive hazards model for multivariate failure time data. Furthermore we suggest estimating equations for the regression parameters and association parameters in parametric shared frailty models with marginal additive hazards by using the Lin and Ying estimators. We give the large sample properties of the estimators arising from these estimating equations and investigate their small sample properties by Monte Carlo simulation. A real example is provided for illustration.     

Continuous-time capture-recapture models with time variation and behavioural response

This paper develops a likelihood-based inference procedure for continuous-time capture-recapture models. The first-capture and recapture intensities are assumed to be in constant proportion but may otherwise vary arbitrarily through time. The full likelihood is partitioned into two factors, one of which is analogous to the likelihood in a special type of multiplicative intensity model arising in failure time analysis. The remaining factor is free of the non-parametric nuisance parameter and is easily maximized. This factor provides an estimator of population size and an asymptotic variance under a counting process framework. The resulting estimation procedure is shown to be equivalent to that derived from a martingale-based estimating function approach. Simulation results are presented to examine the performance of the proposed estimators.     

Robust estimation in generalized linear mixed models

Generalized linear mixed models (GLMMs) are widely used to analyse non-normal response data with extra-variation, but non-robust estimators are still routinely used. We propose robust methods for maximum quasi-likelihood and residual maximum quasi-likelihood estimation to limit the influence of outlying observations in GLMMs. The estimation procedure parallels the development of robust estimation methods in linear mixed models, but with adjustments in the dependent variable and the variance component. The methods proposed are applied to three data sets and a comparison is made with the nonparametric maximum likelihood approach. When applied to a set of epileptic seizure data, the methods proposed have the desired effect of limiting the influence of outlying observations on the parameter estimates. Simulation shows that one of the residual maximum quasi-likelihood proposals has a smaller bias than those of the other estimation methods. We further discuss the equivalence of two GLMM formulations when the response variable follows an exponential family. Their extensions to robust GLMMs and their comparative advantages in modelling are described. Some possible modifications of the robust GLMM estimation methods are given to provide further flexibility for applying the method.     

A Lack-of-Fit Test for Heteroscedastic Regression Models via Cosine-Series Smoothers

In this paper, a test is derived to assess the validity of heteroscedastic nonlinear regression models by a non‐parametric cosine regression method. For order selection, the paper proposes a data‐driven method that uses the parametric null model optimal order. This method yields a test that is asymptotically normally distributed under the null hypothesis and is consistent against any fixed alternative. Simulation studies that test the lack of fit of a generalized linear model are conducted to compare the performance of the proposed test with that of an existing non‐parametric kernel test. A dataset of esterase levels is used to demonstrate the proposed method in practice.     

从“边界冲突”到“关系自觉”—论费孝通如何用“场”修正“差序格局”

“差序格局”原本是费孝通概括中国传统社会亲属关系和社会结构形态的理论工具,也是分析中国人“群己之别”、“圈层意识”的“私”道德的功能概念,更蕴含着人际关系紧张、秩序隐患和边界冲突等负面影响。利用“场”的扩散、交叉性和非边界、融合等特性,来补充“差序格局”,既是费孝通晚年对差序格局内在“推己及人”达成“和”的社会秩序的文化反思,也是对消除现实生活中人际关系隔阂,实现关系自觉的实践关怀。中国人际关系结构的“差序”和“场”二重性,行动者“差序”和“边际”的双重人格,构成中国人际关系交往的实践逻辑,也是费孝通用“场”补充“差序格局”的理论启示。推动费孝通借助“场”来补充和修正“差序格局”的文化意义和研究范畴,仍需进一步开掘。     

Asymptotic properties of maximum quasi-likelihood estimators in generalized linear models with adaptive designs

In this paper, we present the asymptotic properties of maximum quasi-likelihood estimators (MQLEs) in generalized linear models with adaptive designs under some mild regular conditions. The existence of MQLEs in quasi-likelihood equation is discussed. The rate of convergence and asymptotic normality of MQLEs are also established. The results are illustrated by Monte-Carlo simulations.     

Longitudinal poisson regression with disturbed random intercept

Consider repeated event-count data from a sequence of exposures, during each of which a subject can experience some number of events, which is reported at ‘visits’ following each exposure. Within-subject heterogeneity not accounted for by visit-varying covariates is called ‘visit-level’ heterogeneity. Using generalized linear mixed models with log link for longitudinal Poisson regression, I model visit-level heterogeneity by cumulatively adding ‘disturbances’ to the random intercept of each subject over visits to create a ‘disturbed-random-intercept$rsquo; model. I also create a ‘disturbed-random-slope’ model, where the slope is over visits, and both intercept and slope are random but only the slope is disturbed. Simulation studies compare fixed-effect estimation for these models in data with 15 visits, large visit-level heterogeneity, and large multiplicative overdispersion. These studies show statistically significant superiority of the disturbed-random-intercept model. Examples with epidemiological data compare results of this model with those from other published models.     

The Trace Restriction: An Alternative Identification Strategy for the Bayesian Multinomial Probit Model

Previous authors have made Bayesian multinomial probit models identifiable by fixing a parameter on the main diagonal of the covariance matrix. The choice of which element one fixes can influence posterior predictions. Thus, we propose restricting the trace of the covariance matrix, which we achieve without computational penalty. This permits a prior that is symmetric to permutations of the nonbase outcome categories. We find in real and simulated consumer choice datasets that the trace-restricted model is less prone to making extreme predictions. Further, the trace restriction can provide stronger identification, yielding marginal posterior distributions that are more easily interpreted.