Modelling heterogeneous space–time occurrences of earthquakes and its residual analysis

Summary. Earthquake intensities are modelled as a function of previous activity whose specific form is based on established empirical laws in seismology, but whose parameter values can vary from place to place. This model is used for characterizing regional features of seismic activities in and around Japan, and also for exploring regions where the actual seismicity rate systematically deviates from that of the modelled rate.     

Residual analysis for spatial point processes (with discussion)

Summary.  We define residuals for point process models fitted to spatial point pattern data, and we propose diagnostic plots based on them. The residuals apply to any point process model that has a conditional intensity; the model may exhibit spatial heterogeneity, interpoint interaction and dependence on spatial covariates. Some existing ad hoc methods for model checking (quadrat counts, scan statistic, kernel smoothed intensity and Berman's diagnostic) are recovered as special cases. Diagnostic tools are developed systematically, by using an analogy between our spatial residuals and the usual residuals for (non-spatial) generalized linear models. The conditional intensity λ plays the role of the mean response. This makes it possible to adapt existing knowledge about model validation for generalized linear models to the spatial point process context, giving recommendations for diagnostic plots. A plot of smoothed residuals against spatial location, or against a spatial covariate, is effective in diagnosing spatial trend or co-variate effects. Q – Q -plots of the residuals are effective in diagnosing interpoint interaction.     

An individual measure of relative survival

Summary.  Relative survival techniques are used to compare survival experience in a study cohort with that expected if background population rates apply. The techniques are especially useful when cause-specific death information is not accurate or not available as they provide a measure of excess mortality in a group of patients with a certain disease. Whereas these methods are based on group comparisons, we present here a transformation approach which instead gives for each individual an outcome measure relative to the appropriate background population. The new outcome measure is easily interpreted and can be analysed by using standard survival analysis techniques. It provides additional information on relative survival and gives new options in regression analysis. For example, one can estimate the proportion of patients who survived longer than a given percentile of the respective general population or compare survival experience of individuals while accounting for the population differences. The regression models for the new outcome measure are different from existing models, thus providing new possibilities in analysing relative survival data. One distinctive feature of our approach is that we adjust for expected survival before modelling. The paper is motivated by a study into the survival of patients after acute myocardial infarction.     

A note on measures of similarity based on centrality

This note presents a measure of similarity between connected nodes in terms of centrality based on Euclidean distances, and compares it to ‘assortative mixing’ [Newman, M.E.J., 2002. Assortative mixing in networks. Physical Review Letters 89, 208701], which is based on Pearson correlation coefficient. This study suggests that the measure based on Euclidean distances may be more appropriate for relatively smaller (N < 500) and denser networks.     

A procedure for testing the equality of k dependent correlation coefficients

A procedure is proposed for testing the equality of k dependent correlation coefficients. The procedure is simulated utilizing Monte Carlo techniques; and, a method for post hoc probing is also suggested.     

Diagnostics of multiple group influential observations for logistic regression models

In this paper, two new multiple influential observation detection methods, GCD.GSPR and mCD*, are introduced for logistic regression. The proposed diagnostic measures are compared with the generalized difference in fits (GDFFITS) and the generalized squared difference in beta (GSDFBETA), which are multiple influential diagnostics. The simulation study is conducted with one, two and five independent variable logistic regression models. The performance of the diagnostic measures is examined for a single contaminated independent variable for each model and in the case where all the independent variables are contaminated with certain contamination rates and intensity. In addition, the performance of the diagnostic measures is compared in terms of the correct identification rate and swamping rate via a frequently referred to data set in the literature.     

Influence diagnostics in gamma ridge regression model

In this article, we proposed some influence diagnostics for the gamma regression model (GRM) and the gamma ridge regression model (GRRM). We assess the impact of influential observations on the GRM and GRRM estimates by extending the work of Pregibon [Logistic regression diagnostics. Ann Stat. 1981;9:705–724] and Walker and Birch [Influence measures in ridge regression. Technometrics. 1988;30:221–227]. Comparison of both models is made and demonstrated with the help of a simulation study and a real data set. We report some momentous results in detecting the influential observations and their effects on the GRM and GRRM estimates.     

基于产品生命周期的多阶段剩余收益项目决策模型

剩余收益模型（Residual Income Model,RIM）是评估公司权益价值的经典模型,特别地,奥尔森系列剩余收益模型由于可以利用历史财务与会计数据而得到广泛运用。本文在一般剩余收益模型的基础上,运用规范研究、逻辑推理等方法,基于产品生命周期理论,即项目处于产品生命周期的不同阶段时,其净资产收益率不同,在增长期的净资产收益率持续上升,衰退期的净资产收益率持续下降,而成熟期的净资产收益率则围绕行业平均水平波动,从理论上就项目投资价值的评估构建了多阶段剩余收益项目决策模型,并就所构建的多阶段剩余收益项目决策模型进行参数赋值分析。研究结果表明该新构建的决策模型具有较好的理论和运用价值,可以作为一种有效的项目投资决策模型。本文的研究一是拓展了剩余收益模型的应用,从主要用于对公司权益价值的估值拓展至对项目投资决策的评估;二是有利于寻求项目投资决策的价值最大化。     

Approximating the risk score for disease diagnosis using MARS

In disease screening and diagnosis, often multiple markers are measured and combined to improve the accuracy of diagnosis. McIntosh and Pepe [Combining several screening tests: optimality of the risk score, Biometrics 58 (2002), pp. 657–664] showed that the risk score, defined as the probability of disease conditional on multiple markers, is the optimal function for classification based on the Neyman–Pearson lemma. They proposed a two-step procedure to approximate the risk score. However, the resulting receiver operating characteristic (ROC) curve is only defined in a subrange (L, h) of false-positive rates in (0,1) and the determination of the lower limit L needs extra prior information. In practice, most diagnostic tests are not perfect, and it is usually rare that a single marker is uniformly better than the other tests. Using simulation, I show that multivariate adaptive regression spline is a useful tool to approximate the risk score when combining multiple markers, especially when ROC curves from multiple tests cross. The resulting ROC is defined in the whole range of (0,1) and is easy to implement and has intuitive interpretation. The sample code of the application is shown in the appendix.     

Computation of finite population parameters and approximate probability values for multi-response randomized block permutation procedures (mrbp)

An algorithm is presented for computing the finite population parameters and the approximate probability values associated with a recently-developed class of statistical inference techniques termed multi-response randomized block permutation procedures (MRBP).