Modeling developmental data using U-shaped threshold dose-response curves

This paper develops threshold models for developmental toxicity data. The distinguishing feature of these threshold models is their flexibility in modeling data below threshold with a U-shaped function if the data warrants. The method is applied to actual data from a developmental study which exhibits U-shaped behavior in early dose groups. Results from a simulation study demonstrate the flexibility of the threshold model to pick up on U-shaped trends in the data. In addition, the simulation study reveals important considerations in design of developmental studies.     

Threshold selection for regional peaks-over-threshold data

A hurdle in the peaks-over-threshold approach for analyzing extreme values is the selection of the threshold. A method is developed to reduce this obstacle in the presence of multiple, similar data samples. This is for instance the case in many environmental applications. The idea is to combine threshold selection methods into a regional method. Regionalized versions of the threshold stability and the mean excess plot are presented as graphical tools for threshold selection. Moreover, quantitative approaches based on the bootstrap distribution of the spatially averaged Kolmogorov–Smirnov and Anderson–Darling test statistics are introduced. It is demonstrated that the proposed regional method leads to an increased sensitivity for too low thresholds, compared to methods that do not take into account the regional information. The approach can be used for a wide range of univariate threshold selection methods. We test the methods using simulated data and present an application to rainfall data from the Dutch water board Vallei en Veluwe.     

Analysis of the Threshold Values of Semantic Structure Analysis in Identifying Causal Relationships

This study uses a semantic structure analysis (SSA) method to construct the causal relationships among the criteria from survey data. The literatures provide a predetermined threshold value when the SSA is applied without explanation, but we use a Monte Carlo simulation based on the raw data to determine the threshold values with the significant levels of 0.05 and 0.10 for constructing the causal relationships. The results show that the causal relationships among the criteria using the suggested threshold value are too complicated, while the causal relationships by the simulated threshold values are relatively easy to be understood and used practically.     

Robust wavelet shrinkage using robust selection of thresholds

This paper considers the problem of selecting a robust threshold of wavelet shrinkage. Previous approaches reported in literature to handle the presence of outliers mainly focus on developing a robust procedure for a given threshold; this is related to solving a nontrivial optimization problem. The drawback of this approach is that the selection of a robust threshold, which is crucial for the resulting fit is ignored. This paper points out that the best fit can be achieved by a robust wavelet shrinkage with a robust threshold. We propose data-driven selection methods for a robust threshold. These approaches are based on a coupling of classical wavelet thresholding rules with pseudo data. The concept of pseudo data has influenced the implementation of the proposed methods, and provides a fast and efficient algorithm. Results from a simulation study and a real example demonstrate the promising empirical properties of the proposed approaches.     

Common threshold in quantile regressions with an application to pricing for reputation

The paper develops a systematic estimation and inference procedure for quantile regression models where there may exist a common threshold effect across different quantile indices. We first propose a sup-Wald test for the existence of a threshold effect, and then study the asymptotic properties of the estimators in a threshold quantile regression model under the shrinking threshold effect framework. We consider several tests for the presence of a common threshold value across different quantile indices and obtain their limiting distributions. We apply our methodology to study the pricing strategy for reputation through the use of a data set from Taobao.com. In our economic model, an online seller maximizes the sum of the profit from current sales and the possible future gain from a targeted higher reputation level. We show that the model can predict a jump in optimal pricing behavior, which is considered as “reputation effect” in this paper. The use of threshold quantile regression model allows us to identify and explore the reputation effect and its heterogeneity in data. We find both reputation effects and common thresholds for a range of quantile indices in seller’s pricing strategy in our application.     

Threshold‐free estimation of functional antibody titers of a group B streptococcus opsonophagocytic killing assay

Opsonophagocytic killing assays (OPKA) are routinely used for the quantification of bactericidal antibodies in blood serum samples. Quantification of the OPKA readout, the titer, provides the basis for the statistical analysis of vaccine clinical trials having functional immune response endpoints. Traditional OPKA titers are defined as the maximum serum dilution yielding a predefined bacterial killing threshold value, and they are estimated by fitting a dose‐response model to the dilution‐killing curve. This paper illustrates a novel definition of titer, the threshold‐free titer, which preserves biological interpretability while not depending on any killing threshold or on a postulated shape of the dose‐response curve. These titers are shown to be more precise than the traditional threshold‐based titers when using simulated and experimental group B streptococcus OPKA experimental data. Also, titer linearity is shown to be not measurable when using threshold‐based titers, whereas it becomes measurable using threshold‐free titers. The biological interpretability and operational characteristics demonstrated here indicate that threshold‐free titers are an appropriate tool for the routine analysis of OPKA data. Copyright © 2015 John Wiley & Sons, Ltd.     

Threshold autoregression analysis for finite-range time series of counts with an application on measles data

This article studies the threshold autoregression analysis for the self-exciting threshold binomial autoregressive processes. Parameters' point estimation and interval estimation problems are considered via the empirical likelihood method. A new algorithm to estimate the threshold value of the threshold model is also given. Simulation study is conducted for the evaluation of the developed approach. An application on measles data is provided to show the applicability of the method.     

A Threshold Dose-Response Model with Random Effects in Teratological Experiments

Abstract

Teratological experiments are controlled dose-response studies in which impregnated animals are randomly assigned to various exposure levels of a toxic substance. Subsequently, both continuous and discrete responses are recorded on the litters of fetuses that these animals produce. Discrete responses are usually binary in nature, such as the presence or absence of some fetal anomaly. This clustered binary data usually exhibits over-dispersion (or under-dispersion), which can be interpreted as either variation between litter response probabilities or intralitter correlation. To model the correlation and/or variation, the beta-binomial distribution has been assumed for the number of positive fetal responses within a litter. Although the mean of the beta-binomial model has been linked to dose-response functions, in terms of measuring over-dispersion, it may be a restrictive method in modeling data from teratological studies. Also for certain toxins, a threshold effect has been observed in the dose-response pattern of the data. We propose to incorporate a random effect into a general threshold dose-response model to account for the variation in responses, while at the same time estimating the threshold effect. We fit this model to a well-known data set in the field of teratology. Simulation studies are performed to assess the validity of the random effects threshold model in these types of studies.     

Interlitter Response Variability in a Threshold Dose-Response Model

In this article, we extend a previously formulated threshold dose-response model with random litter effects that was applied to a data set from a developmental toxicity study. The dose-response pattern of the data indicates that a threshold dose level may exist. Additionally, there is noticeable variation between the responses across the dose levels. With threshold estimation being critical, the assumed variability structure should adequately model the variation while not taking away from the estimation of the threshold as well as the other parameters directly involved in the dose-response relationship. In the prior formulation, the random effect was modeled assuming identical variation in the interlitter response probabilities across all dose levels, that is, the model had a single parameter to account for the interlitter variability. In this new model, the random effect is modeled as having different response variability across dose levels, that is, multiple interlitter variability parameters. We performed the likelihood ratio test (LRT) to compare our extended model to the previous model. We conducted a simulation study to compare the bias of each model when fit to data generated with the underlying parametric structure of the opposing model. The extended threshold dose-response model with multiple response variation was less biased.     

Classification rules for identifying individuals at high risk of developing myocardial infarction based on ApoB,ApoA1 and the ratio were determined using a Bayesian approach

We have developed a new approach to determine the threshold of a biomarker that maximizes the classification accuracy of a disease. We consider a Bayesian estimation procedure for this purpose and illustrate the method using a real data set. In particular, we determine the threshold for Apolipoprotein B (ApoB), Apolipoprotein A1 (ApoA1) and the ratio for the classification of myocardial infarction (MI). We first conduct a literature review and construct prior distributions. We then develop classification rules based on the posterior distribution of the location and scale parameters for these biomarkers. We identify the threshold for ApoB and ApoA1, and the ratio as 0.908 (gram/liter), 1.138 (gram/liter) and 0.808, respectively. We also observe that the threshold for disease classification varies substantially across different age and ethnic groups. Next, we identify the most informative predictor for MI among the three biomarkers. Based on this analysis, ApoA1 appeared to be a stronger predictor than ApoB for MI classification. Given that we have used this data set for illustration only, the results will require further investigation for use in clinical applications. However, the approach developed in this article can be used to determine the threshold of any continuous biomarker for a binary disease classification.     

Making predictions from complex longitudinal data, with application to planning monitoring intervals in a national screening programme

When biological or physiological variables change over time, we are often interested in making predictions either of future measurements or of the time taken to reach some threshold value. On the basis of longitudinal data for multiple individuals, we develop Bayesian hierarchical models for making these predictions together with their associated uncertainty. Particular aspects addressed, which include some novel components, are handling curvature in individuals' trends over time, making predictions for both underlying and measured levels, making predictions from a single baseline measurement, making predictions from a series of measurements, allowing flexibility in the error and random-effects distributions, and including covariates. In the context of data on the expansion of abdominal aortic aneurysms over time, where reaching a certain threshold leads to referral for surgery, we discuss the practical application of these models to the planning of monitoring intervals in a national screening programme. Prediction of the time to reach a threshold was too imprecise to be practically useful, and we focus instead on limiting the probability of exceeding the threshold after given time intervals. Although more complex models can be shown to fit the data better, we find that relatively simple models seem to be adequate for planning monitoring intervals.     

The asymptotic distribution of MLE of treatment lag threshold

Semiparametric models with an unknown threshold parameter for hazard rates are formulated for failure time data. Maximum partial likelihood estimates for threshold and regular parameters are shown to be asymptotically consistent. The asymptotic distributions of the estimators are obtained.     

Extended generalised Pareto models for tail estimation

The most popular approach in extreme value statistics is the modelling of threshold exceedances using the asymptotically motivated generalised Pareto distribution. This approach involves the selection of a high threshold above which the model fits the data well. Sometimes, few observations of a measurement process might be recorded in applications and so selecting a high quantile of the sample as the threshold leads to almost no exceedances. In this paper we propose extensions of the generalised Pareto distribution that incorporate an additional shape parameter while keeping the tail behaviour unaffected. The inclusion of this parameter offers additional structure for the main body of the distribution, improves the stability of the modified scale, tail index and return level estimates to threshold choice and allows a lower threshold to be selected. We illustrate the benefits of the proposed models with a simulation study and two case studies.     

我国通货膨胀率的最优目标区间几何？

 内容提要：本文首先提出了面板数据动态门限回归模型的二阶段合并最小二乘（2SPOLS）估计方法;其次,基于中国29省市自治区1978-2008年的面板数据,对中国通货膨胀和经济增长之间关系的实证分析发现,在一定程度上,我国通货膨胀率对经济增长率的作用存在两个门限值的“双门限效应”,其门限值分别为3.2%和15.7%。所以,通货膨胀率位于（0%,3.2%]时,温和通货膨胀对经济增长率存在“托宾效应”。通货膨胀率超过3.2%时,通货膨胀率对经济增长率存在阻碍经济增长的“反托宾效应”,尤其,通货膨胀率高于15.7%后,恶性通货膨胀严重阻碍经济“软扩张”。因此,我国通货膨胀率的最优目标区间是（0%,3.2%]。     

Wavelet Threshold Estimators for Data with Correlated Noise

Wavelet threshold estimators for data with stationary correlated noise are constructed by applying a level-dependent soft threshold to the coefficients in the wavelet transform. A variety of threshold choices is proposed, including one based on an unbiased estimate of mean-squared error. The practical performance of the method is demonstrated on examples, including data from a neurophysiological context. The theoretical properties of the estimators are investigated by comparing them with an ideal but unattainable `bench-mark', that can be considered in the wavelet context as the risk obtained by ideal spatial adaptivity, and more generally is obtained by the use of an `oracle' that provides information that is not actually available in the data. It is shown that the level-dependent threshold estimator performs well relative to the bench-mark risk, and that its minimax behaviour cannot be improved on in order of magnitude by any other estimator. The wavelet domain structure of both short- and long-range dependent noise is considered, and in both cases it is shown that the estimators have near optimal behaviour simultaneously in a wide range of function classes, adapting automatically to the regularity properties of the underlying model. The proofs of the main results are obtained by considering a more general multivariate normal decision theoretic problem.     

Single-Variable Threshold Effects in Ordered Response Models With an Application to Estimating the Income-Happiness Gradient

This short article extends well-known threshold models to the ordered response setting. We consider the case where the sample is endogenously split to estimate regime-dependent coefficients for one variable of interest, while keeping the other coefficients and auxiliary parameters constant across the threshold. We use Monte Carlo methods to examine the behavior of the model. In addition, we derive the formulae for the partial effects associated with the model. We apply our threshold model to the relationship between income and self-reported happiness using data drawn from the U.S. General Social Survey. While the findings suggest the presence of a threshold in the income-happiness gradient at approximately U.S. $76,000, no evidence is found in support of a satiation point. Supplementary materials for this article are available online.     

Adaptive stochastic gradient boosting tree with composite criterion

ABSTRACT

In this paper, we propose an adaptive stochastic gradient boosting tree for classification studies with imbalanced data. The adjustment of cost-sensitivity and the predictive threshold are integrated together with a composite criterion into the original stochastic gradient boosting tree to deal with the issues of the imbalanced data structure. Numerical study shows that the proposed method can significantly enhance the classification accuracy for the minority class with only a small loss in the true negative rate for the majority class. We discuss the relation of the cost-sensitivity to the threshold manipulation using simulations. An illustrative example of the analysis of suboptimal health-state data in traditional Chinese medicine is discussed.     

Bayesian analysis of threshold autoregressions

A nonasymptotic Bayesian approach is developed for analysis of data from threshold autoregressive processes with two regimes. Using the conditional likelihood function, the marginal posterior distribution for each of the parameters is derived along with posterior means and variances. A test for linear functions of the autoregressive coefficients is presented. The approach presented uses a posterior p-value averaged over the values of the threshold. The one-step ahead predictive distribution is derived along with the predictive mean and variance. In addition, equivalent results are derived conditional upon a value of the threshold. A numerical example is presented to illustrate the approach.     

A sickle transition-rate model with starting threshold

Billari (2001) introduced a new type of single-spell parametric transition-rate model: transition-rate models with a starting threshold. In such models, the transition-rate function is composed of two additive terms. The first term is a constant that holds for any given duration; the second is a ‘traditional’ transition-rate function with the threshold as its time origin, and it is added after a certain threshold point. The possibility of allowing for the presence of long-term survivors in the social process has not yet been dealt with, and it is of specific interest in several domains of application. In this paper, we develop the specific case of the sickle model. We discuss its features, its implementation as a starting threshold model, and the estimation of its parameters. The sickle model with starting threshold is then applied to the union formation of Italian men and women, using the Fertility and Family Survey data.     

Continuous threshold models with two-way interactions in survival analysis

Proportional hazards model with the biomarker–treatment interaction plays an important role in the survival analysis of the subset treatment effect. A threshold parameter for a continuous biomarker variable defines the subset of patients who can benefit or lose from a certain new treatment. In this article, we focus on a continuous threshold effect using the rectified linear unit and propose a gradient descent method to obtain the maximum likelihood estimation of the regression coefficients and the threshold parameter simultaneously. Under certain regularity conditions, we prove the consistency, asymptotic normality and provide a robust estimate of the covariance matrix when the model is misspecified. To illustrate the finite sample properties of the proposed methods, we simulate data to evaluate the empirical biases, the standard errors and the coverage probabilities for both the correctly specified models and misspecified models. The proposed continuous threshold model is applied to a prostate cancer data with serum prostatic acid phosphatase as a biomarker.