曲线分类建模方法及其在GDP预测中的应用

传统方法解决大规模时序曲线的预测建模问题,需要对每条曲线逐一建模,使得建模工作量相当庞大,在实际应用中缺乏可操作性。文章提出一种解决此问题的新方法——曲线分类建模方法。该方法先减少曲线的模型种类,再进行曲线分类和分类建模,在尽可能保留原始信息的前提下较大程度地降低了建模的工作量。文章阐述了该方法的原理和计算过程,并以应用于多地区GDP曲线的预测案例说明该方法的实用性和有效性。     

Graphical representation of survival curves in the presence of time-dependent categorical covariates with application to liver transplantation

Graphical representation of survival curves is often used to illustrate associations between exposures and time-to-event outcomes. However, when exposures are time-dependent, calculation of survival probabilities is not straightforward. Our aim was to develop a method to estimate time-dependent survival probabilities and represent them graphically. Cox models with time-dependent indicators to represent state changes were fitted, and survival probabilities were plotted using pre-specified times of state changes. Time-varying hazard ratios for the state change were also explored. The method was applied to data from the Adult-to-Adult Living Donor Liver Transplantation Cohort Study (A2ALL). Survival curves showing a ‘split’ at a pre-specified time t allow for the qualitative comparison of survival probabilities between patients with similar baseline covariates who do and do not experience a state change at time t. Time since state change interactions can be visually represented to reflect changing hazard ratios over time. A2ALL study results showed differences in survival probabilities among those who did not receive a transplant, received a living donor transplant, and received a deceased donor transplant. These graphical representations of survival curves with time-dependent indicators improve upon previous methods and allow for clinically meaningful interpretation.     

Inflation projections for monetary policy decision making

Many central banks set inflation targets over the medium term and inflation projections are a key input for monetary policy decision making. In this paper, we present the procedures used by the Spanish Central Bank staff to project consumer price inflation. We also provide some illustrations of their policy uses, such as fan charts, deflation probabilities and the monitoring of inflation targets.     

Analysis of microtubule dynamics using growth curve models

Microtubules are part of the structural network within a cell's cytoplasm, providing structural support as well as taking part in many of the cellular processes. A large body of data provide evidence that dynamics of microtubules in a cell is reponsible for the performance of many critical cellular functions such as cell division. In this article, we study the effect of four different isoforms of a protein tau on microtubule dynamics using growth curve models. The results show that a linear growth curve model is sufficient to explain the data. Moreover, we find that a mutated version of a 3-repeat tau protein has a similar effect as a 4-repeat tau protein on microtubule dynamics. The latter findings conform with the biological understanding of the effect of the protein tau on microtubule dynamics.     

An alternative method for global and partial comparison of two diagnostic systems based on ROC curves

In this paper, an alternative method for the comparison of two diagnostic systems based on receiver operating characteristic (ROC) curves is presented. ROC curve analysis is often used as a statistical tool for the evaluation of diagnostic systems. However, in general, the comparison of ROC curves is not straightforward, in particular, when they cross each other. A similar difficulty is also observed in the multi-objective optimization field where sets of solutions defining fronts must be compared with a multi-dimensional space. Thus, the proposed methodology is based on a procedure used to compare the performance of distinct multi-objective optimization algorithms. In general, methods based on the area under the ROC curves are not sensitive to the existence of crossing points between the curves. The new approach can deal with this situation and also allows the comparison of partial portions of ROC curves according to particular values of sensitivity and specificity of practical interest. Simulations results are presented. For illustration purposes, considering real data from newborns with very low birthweight, the new method was applied in order to discriminate the better index for evaluating the risk of death.     

Effects of the Mechanism of Receptor-Mediated Gene Expression on the Shape of the Dose-Response Curve

A mathematical model of receptor-mediated gene expression that includes receptor binding of natural and xenobiotic ligands, protein synthesis and degradation, and metabolism of the xenobiotic ligand was created to identify the determinants of the shape of the dose-response profile. Values of the model's parameters were varied to reflect alternative mechanisms of expression of the protein. These assumptions had dramatic effects on the computed response to a bolus dose of the xenobiotic ligand. If all processes in the model exhibit hyperbolic kinetics, the dose-response curves can appear sigmoidal but actually be linear with a positive slope at low doses. The slope of the curve only approached zero at low dose, indicative of a threshold for response, if binding of the xenobiotic ligand to the receptor exhibited positive cooperativity (ligand binding at one site increases the affinity for ligand at another binding site on the receptor). Positive cooperativity in the rate-limiting step of protein synthesis produced dose-response curves which were "U-shaped" at low doses, also indicative of a threshold. Positive cooperativity in the metabolism of the xenobiotic ligand produced dose-response curves that increased more rapidly than linearly with increasing dose. The model illustrates the fact that response cannot be predicted from qualitative mechanistic arguments alone; any assessment of risk to health from xenobiotic chemicals must be based on a detailed quantitative examination of the kinetic behavior of each chemical species individually.     

Sequential experimental design and response optimisation

We consider the situation where one wants to maximise a functionf(θ,x) with respect tox, with θ unknown and estimated from observationsy _k . This may correspond to the case of a regression model, where one observesy _k =f(θ,x _k )+ε _k , with ε _k some random error, or to the Bernoulli case wherey _k ∈{0, 1}, with Pr[y _k =1|θ,x _k |=f(θ,x _k ). Special attention is given to sequences given by , with an estimated value of θ obtained from (x₁, y₁),...,(x _k ,y _k ) andd _k (x) a penalty for poor estimation. Approximately optimal rules are suggested in the linear regression case with a finite horizon, where one wants to maximize ∑ _i=1 ^N w _i f(θ, x _i ) with {w _i } a weighting sequence. Various examples are presented, with a comparison with a Polya urn design and an up-and-down method for a binary response problem.     

Growth and Heat Resistance Kinetic Variation Among Various Isolates of Salmonella and its Application to Risk Assessment

The abilities of cells of a particular type of bacteria to leave lag phase and begin the process of dividing or surviving heat treatment can depend on the serotypes or strains of the bacteria. This article reports an investigation of serotype-specific differences in growth and heat resistance kinetics of clinical and food isolates of Salmonella. Growth kinetics at 19 degrees C and 37 degrees C were examined in brain heart infusion broth and heat resistance kinetics for 60 degrees C were examined in beef gravy using a submerged coil heating apparatus. Estimates of the parameters of the growth curves suggests a small between-serotype variance of the growth kinetics. However, for inactivation, the results suggest a significant between-serotype effect on the asymptotic D-values, with an estimated between-serotype CV of about 20%. In microbial risk assessment, predictive microbiology is used to estimate growth and inactivation of pathogens. Often the data used for estimating the growth or inactivation kinetics are based on measurements on a cocktail--a mixture of approximately equal proportions of several serotypes or strains of the pathogen being studied. The expected growth or inactivation rates derived from data using cocktails are biased, reflecting the characteristics of the fastest growing or most heat resistant serotype of the cocktail. In this article, an adjustment to decrease this possible bias in a risk assessment is offered. The article also presents discussion of the effect on estimating growth when stochastic assumptions are incorporated in the model. In particular, equations describing the variation of relative growth are derived, accounting for the stochastic variations of the division of cells. For small numbers of cells, the expected value of the relative growth is not an appropriate "representative" value for actual relative growths that might occur.     

The estimation of derivatives of a nonparametric regression function when the data are correlated

A kernel estimator of a derivative of arbitrary order of a nonparametric average population curve is considered for a correlated-errors model with balanced replicate measurements at each design point. Asymptotic expansions of the mean squared error are derived for two classes of correlation functions in the model. Consistency, choice of smoothing parameter, and rates of convergence are examined for the important special cases of estimating the first and second derivatives.     

Assessing program impact using latent growth modeling: a primer for the evaluator

Latent growth modeling (LGM) has emerged as a flexible analytic technique for modeling change over time because it can describe developmental processes at both the inter- and intra-individual levels. The LGM method can also provide a means for testing the contribution of other variables in order to explain variability in growth trajectories. This paper didactically illustrates the use of LGM as an analytical tool in program evaluation. Specifically, a hypothetical evaluation of a high school drug prevention program was used to demonstrate: (a) how LGM can be used to assess the longitudinal impact of a prevention program by comparing treatment and control populations with respect to individual differences in initial status and in rate of change; and (b) how predictors of initial status (post-intervention) and growth selected on the basis of a particular program theory can be incorporated in the model to explain program impact. Some advantages and limitations of using LGM in program evaluation are highlighted.