Survival analysis statistics applied to threshold data obtained from the ascending forced-choice method of limits

Sensory thresholds are often collected through ascending forced-choice methods. Group thresholds are important for comparing stimuli or populations; yet, the method has two problems. An individual may correctly guess the correct answer at any concentration step and might detect correctly at low concentrations but become adapted or fatigued at higher concentrations. The survival analysis method deals with both issues. Individual sequences of incorrect and correct answers are adjusted, taking into account the group performance at each concentration. The technique reduces the chance probability where there are consecutive correct answers. Adjusted sequences are submitted to survival analysis to determine group thresholds. The technique was applied to an aroma threshold and a taste threshold study. It resulted in group thresholds similar to ASTM or logarithmic regression procedures. Significant differences in taste thresholds between younger and older adults were determined. The approach provides a more robust technique over previous estimation methods.

Tests of sunspot number sequences: 3. Effects of regression procedures on the calibration of historic sunspot data

We use sunspot group observations from the Royal Greenwich Observatory (RGO) to investigate the effects of intercalibrating data from observers with different visual acuities. The tests are made by counting the number of groups RB above a variable cut-off threshold of observed total whole-spot area (uncorrected for foreshortening) to simulate what a lower acuity observer would have seen. The synthesised annual means of RB are then re-scaled to the full observed RGO group number RA using a variety of regression techniques. It is found that a very high correlation between RA and RB (rAB > 0.98) does not prevent large errors in the intercalibration (for example sunspot maximum values can be over 30 % too large even for such levels of rAB). In generating the backbone sunspot number (RBB), Svalgaard and Schatten (2015, this issue) force regression fits to pass through the scatter plot origin which generates unreliable fits (the residuals do not form a normal distribution) and causes sunspot cycle amplitudes to be exaggerated in the intercalibrated data. It is demonstrated that the use of Quantile-Quantile (“Q  Q”) plots to test for a normal distribution is a useful indicator of erroneous and misleading regression fits. Ordinary least squares linear fits, not forced to pass through the origin, are sometimes reliable (although the optimum method used is shown to be different when matching peak and average sunspot group numbers). However, other fits are only reliable if non-linear regression is used. From these results it is entirely possible that the inflation of solar cycle amplitudes in the backbone group sunspot number as one goes back in time, relative to related solar-terrestrial parameters, is entirely caused by the use of inappropriate and non-robust regression techniques to calibrate the sunspot data.

EMG amplitude and frequency parameters of muscular activity: Effect of resistance training based on electromyographic fatigue threshold

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Anaerobic threshold estimation by statistical modelling.

Anaerobic threshold (AT) is usually estimated as a change point problem by visual analysis of the cardiorespiratory response to incremental dynamic exercise. In this study, two phase linear (TPL) models of the linear-linear and linear-quadratic type were used for the estimation of AT. The correlation coefficient between the classical and statistical approaches was 0.88, and 0.89 after outlier exclusion. The TPL models provide a simple method for estimating AT that can be easily implemented using a digital computer for the automatic pattern recognition of AT.

A general threshold stress hybrid hazard model for lifetime data

In this paper we propose a hybrid hazard regression model with threshold stress which includes the proportional hazards and the accelerated failure time models as particular cases. To express the behavior of lifetimes the generalized-gamma distribution is assumed and an inverse power law model with a threshold stress is considered. For parameter estimation we develop a sampling-based posterior inference procedure based on Markov Chain Monte Carlo techniques. We assume proper but vague priors for the parameters of interest. A simulation study investigates the frequentist properties of the proposed estimators obtained under the assumption of vague priors. Further, some discussions on model selection criteria are given. The methodology is illustrated on simulated and real lifetime data set.

Blood glucose and prognosis in children with presumed severe malaria: is there a threshold for 'hypoglycaemia'?

OBJECTIVES: Hypoglycaemia (glucose <2.2 mmol/l) is a defining feature of severe malaria, but the significance of other levels of blood glucose has not previously been studied in children with severe malaria. METHODS: A prospective study of 437 consecutive children with presumed severe malaria was conducted in Mali. We defined hypoglycaemia as <2.2 mmol/l, low glycaemia as 2.2-4.4 mmol/l and hyperglycaemia as >8.3 mmol/l. Associations between glycaemia and case fatality were analysed for 418 children using logistic regression models and a receiver operator curve (ROC). RESULTS: There was a significant difference between blood glucose levels in children who died (median 4.6 mmol/l) and survivors (median 7.6 mmol/l, P < 0.001). Case fatality declined from 61.5% of the hypoglycaemic children to 46.2% of those with low glycaemia, 13.4% of those with normal glycaemia and 7.6% of those with hyperglycaemia (P < 0.001). Logistic regression showed an adjusted odds ratio (AOR) of 0.75 (0.64-0.88) for case fatality per 1 mmol/l increase in baseline blood glucose. Compared to a normal blood glucose, hypoglycaemia and low glycaemia both significantly increased the odds of death (AOR 11.87, 2.10-67.00; and 5.21, 1.86-14.63, respectively), whereas hyperglycaemia reduced the odds of death (AOR 0.34, 0.13-0.91). The ROC [area under the curve at 0.753 (95% CI 0.684-0.820)] indicated that glycaemia had a moderate predictive value for death and identified an optimal threshold at glycaemia <6.1 mmol/l, (sensitivity 64.5% and specificity 75.1%). CONCLUSIONS: If there is a threshold of blood glucose which defines a worse prognosis, it is at a higher level than the current definition of 2.2 mmol/l.

Migration and population expansion of Abies, Fagus, Picea, and Quercus since 15000 years in and across the Alps, based on pollen-percentage threshold values

Aims: The aim of this study is to explore the migration (colonization of new areas) and subsequent population expansion (within an area) since 15 ka cal BP of Abies, Fagus, Picea, and Quercus into and through the Alps solely on the basis of high-quality pollen data. Methods: Chronologies of 101 pollen sequences are improved or created. Data from the area delimited by 45.5–48.1°N and 6–14°E are summarized in three ways: (1) in a selection of pollen-percentage threshold maps (thresholds 0.5%, 1%, 2%, 4%, 8%, 16%, and 32% of land pollen); (2) in graphic summaries of 250-year time slices and geographic segments (lengthwise and transverse in relation to the main axis of the Alps) as pollen-percentage curves, pollen-percentage difference curves, and pollen-percentage threshold ages cal BP graphed against both the length and the transverse Alpine axes; and (3) in tables showing statistical relationships of either pollen-percentage threshold ages cal BP or pollen expansion durations (=time lapse between different pollen-percentage threshold ages cal BP) with latitude, longitude, and elevation; to establish these relationships we used both simple linear regression and multiple linear regression after stepwise-forward selection. Results: The statistical results indicate that (a) the use of pollen-percentage thresholds between 0.5% and 8% yield mostly similar directions of tree migration, so the method is fairly robust, (b) Abies migrated northward, Fagus southward, Picea westward, and Quercus northward; more detail does not emerge due to an extreme scarcity of high-quality data especially along the southern foothills of the Alps and in the eastern Alps. This scarcity allows the reconstruction of one immigration route only of Abies into the southern Alps. The speed of population expansion (following arrival) of Abies increased and of Picea decreased during the Holocene, of Fagus it decreased especially during the later Holocene, and of Quercus it increased especially at the start of the Holocene.

Quantile Inferences for Inflation and Its Variability: Does a Threshold Inflation Rate Exist?

Using quantile regressions and cross-sectional data from 152 countries, we examine the relationship between inflation and its variability. We consider two measures of inflation - the mean and median - and three different measures of inflation variability - the standard deviation, coefficient of variation, and median deviation. Using the mean and standard deviation or the median and the median deviation, the results support both the hypothesis that higher inflation creates more inflation variability and that inflation variability raises inflation across quantiles. Moreover, higher quantiles in both cases lead to larger marginal effects of inflation (inflation variability) on inflation variability (inflation). Using the mean and the coefficient of variation, however, the findings largely support no correlation between inflation and its variability. Finally, we also consider whether thresholds for inflation rate or inflation variability exist before finding such positive correlations. We find evidence of thresholds for inflation rates below 3 percent, but mixed results for thresholds for inflation variability.

Clinically-relevant threshold of preformed donor-specific anti-HLA antibodies in kidney transplantation.

BACKGROUND: Pretransplant anti-HLA donor-specific antibodies (DSA) are recognized as a risk factor for acute antibody-mediated rejection (AMR) in kidney transplantation. The predictive value of C4d-fixing capability by DSA or of IgG DSA subclasses for acute AMR in the pretransplant setting has been recently studied. In addition DSA strength assessed by mean fluorescence intensity (MFI) may improve risk stratification. We aimed to analyze the relevance of preformed DSA and of DSA MFI values. METHODS: 280 consecutive patients with negative complement-dependent cytotoxicity crossmatches received a kidney transplant between 01/2008 and 03/2014. Sera were screened for the presence of DSA with a solid-phase assays on a Luminex flow analyzer, and the results were correlated with biopsy-proven acute AMR in the first year and survival. RESULTS: Pretransplant anti-HLA antibodies were present in 72 patients (25.7%) and 24 (8.6%) had DSA. There were 46 (16.4%) acute rejection episodes, 32 (11.4%) being cellular and 14 (5.0%) AMR. The incidence of acute AMR was higher in patients with pretransplant DSA (41.7%) than in those without (1.6%) (p<0.001). The median cumulative MFI (cMFI) of the group DSA+/AMR+ was 5680 vs 2208 in DSA+/AMR- (p=0.058). With univariate logistic regression a threshold value of 5280 cMFI was predictive for acute AMR. DSA cMFI's ability to predict AMR was also explored by ROC analysis. AUC was 0.728 and the best threshold was a cMFI of 4340. Importantly pretransplant DSA>5280 cMFI had a detrimental effect on 5-year graft survival. CONCLUSIONS: Preformed DSA cMFI values were clinically-relevant for the prediction of acute AMR and graft survival in kidney transplantation. A threshold of 4300-5300 cMFI was a significant outcome predictor.

Threshold values of sagittal abdominal diameter for the detection of cardio-metabolic risk factors in Northeastern Mexico: a cross-sectional study

Background: The use of sagittal abdominal diameter (SAD) has been proposed for screening cardio-metabolic risk factors; however, its accuracy can be influenced by the choice of thresholds values. Aim: To determine the SAD threshold values for cardio-metabolic risk factors in Mexican adults; to assess whether parallel and serial SAD testing can improve waist circumference (WC) sensitivity and specificity; and to analyze the effect of considering SAD along with WC and body mass index (BMI) in detecting cardio-metabolic risk. Methods: This cross-sectional study was conducted during 2012-2014 in Northeast Mexico (n = 269). Data on anthropometric, clinical, and biochemical measurements were collected. Sex-adjusted receiver-operating characteristic curves (ROC) were obtained using hypertension, dysglycemia, dyslipidemia and insulin resistance as individual outcomes and metabolic syndrome as a composite outcome. Age-adjusted odds ratios and 95% confidence intervals (CI) were estimated using logistic regression. Results: The threshold value for SAD with acceptable combination of sensitivity and specificity was 24.6 cm in men and 22.5 cm in women. Parallel SAD testing improved WC sensitivity and serial testing improved WC specificity. The co-occurrence of high WC/high SAD increased the risk for insulin resistance by 2.4-fold (95% CI: 1.1-5.3), high BMI/high SAD by 4.3-fold (95% CI: 1.7-11.9) and SAD alone by 2.2-fold (95% CI: 1.2.-4.2). Conclusions: The use of SAD together with traditional obesity indices such as WC and BMI has advantages over using either of these indices alone. SAD may be a powerful screening tool for interventions for high-risk individuals.

Structural Connectivity Of The Default Mode Network And Cognition In Alzheimer׳s Disease.

Disconnectivity between the Default Mode Network (DMN) nodes can cause clinical symptoms and cognitive deficits in Alzheimer׳s disease (AD). We aimed to examine the structural connectivity between DMN nodes, to verify the extent in which white matter disconnection affects cognitive performance. MRI data of 76 subjects (25 mild AD, 21 amnestic Mild Cognitive Impairment subjects and 30 controls) were acquired on a 3.0T scanner. ExploreDTI software (fractional Anisotropy threshold=0.25 and the angular threshold=60°) calculated axial, radial, and mean diffusivities, fractional anisotropy and streamline count. AD patients showed lower fractional anisotropy (P=0.01) and streamline count (P=0.029), and higher radial diffusivity (P=0.014) than controls in the cingulum. After correction for white matter atrophy, only fractional anisotropy and radial diffusivity remained significantly lower in AD compared to controls (P=0.003 and P=0.05). In the parahippocampal bundle, AD patients had lower mean and radial diffusivities (P=0.048 and P=0.013) compared to controls, from which only radial diffusivity survived for white matter adjustment (P=0.05). Regression models revealed that cognitive performance is also accounted for by white matter microstructural values. Structural connectivity within the DMN is important to the execution of high-complexity tasks, probably due to its relevant role in the integration of the network.

Censored Linear Regression Models For Irregularly Observed Longitudinal Data Using The Multivariate-t Distribution.

In acquired immunodeficiency syndrome (AIDS) studies it is quite common to observe viral load measurements collected irregularly over time. Moreover, these measurements can be subjected to some upper and/or lower detection limits depending on the quantification assays. A complication arises when these continuous repeated measures have a heavy-tailed behavior. For such data structures, we propose a robust structure for a censored linear model based on the multivariate Student's t-distribution. To compensate for the autocorrelation existing among irregularly observed measures, a damped exponential correlation structure is employed. An efficient expectation maximization type algorithm is developed for computing the maximum likelihood estimates, obtaining as a by-product the standard errors of the fixed effects and the log-likelihood function. The proposed algorithm uses closed-form expressions at the E-step that rely on formulas for the mean and variance of a truncated multivariate Student's t-distribution. The methodology is illustrated through an application to an Human Immunodeficiency Virus-AIDS (HIV-AIDS) study and several simulation studies.

Pain Threshold And Sleep Quality In Women With Endometriosis.

Pain is a common complaint in women with endometriosis and can be influenced by many variables, including sleep disorders; however, no data are available on the sleep quality of women with endometriosis or on the correlation between sleep quality and pain. The 510 volunteers included in this study were divided into two groups: 257 women with a laparoscopic and histopathological diagnosis of endometriosis and 253 women with no history of endometriosis and no endometriosis-related symptoms. The volunteers answered two questionnaires: the Post-Sleep Inventory to evaluate sleep quality and the International Physical Activity Questionnaire to assess their level of physical activity. Pain was evaluated using a visual analogue scale (VAS) and women were also submitted to a physical examination, during which their pain threshold was assessed at 20 different body sites. Sleep quality was significantly poorer in women with endometriosis compared to women without the disease. The pain threshold was significantly lower in the greater trochanter and abdomen in women with endometriosis when compared to women without the disease; however, there was no difference in VAS pain score between the groups. The higher the VAS pain score, the lower the Post-Sleep Inventory score. Additionally, there was a significant positive correlation between the pain threshold at some body sites and sleep quality. Sleep quality was poorer and the pain threshold at certain body sites was lower in the group of women with endometriosis.

Determination Of Detergent And Dispensant Additives In Gasoline By Ring-oven And Near Infrared Hypespectral Imaging.

A method using the ring-oven technique for pre-concentration in filter paper discs and near infrared hyperspectral imaging is proposed to identify four detergent and dispersant additives, and to determine their concentration in gasoline. Different approaches were used to select the best image data processing in order to gather the relevant spectral information. This was attained by selecting the pixels of the region of interest (ROI), using a pre-calculated threshold value of the PCA scores arranged as histograms, to select the spectra set; summing up the selected spectra to achieve representativeness; and compensating for the superimposed filter paper spectral information, also supported by scores histograms for each individual sample. The best classification model was achieved using linear discriminant analysis and genetic algorithm (LDA/GA), whose correct classification rate in the external validation set was 92%. Previous classification of the type of additive present in the gasoline is necessary to define the PLS model required for its quantitative determination. Considering that two of the additives studied present high spectral similarity, a PLS regression model was constructed to predict their content in gasoline, while two additional models were used for the remaining additives. The results for the external validation of these regression models showed a mean percentage error of prediction varying from 5 to 15%.

Comparing Near-infrared Conventional Diffuse Reflectance Spectroscopy And Hyperspectral Imaging For Determination Of The Bulk Properties Of Solid Samples By Multivariate Regression: Determination Of Mooney Viscosity And Plasticity Indices Of Natural Rubber.

Conventional reflectance spectroscopy (NIRS) and hyperspectral imaging (HI) in the near-infrared region (1000-2500 nm) are evaluated and compared, using, as the case study, the determination of relevant properties related to the quality of natural rubber. Mooney viscosity (MV) and plasticity indices (PI) (PI0 - original plasticity, PI30 - plasticity after accelerated aging, and PRI - the plasticity retention index after accelerated aging) of rubber were determined using multivariate regression models. Two hundred and eighty six samples of rubber were measured using conventional and hyperspectral near-infrared imaging reflectance instruments in the range of 1000-2500 nm. The sample set was split into regression (n = 191) and external validation (n = 95) sub-sets. Three instruments were employed for data acquisition: a line scanning hyperspectral camera and two conventional FT-NIR spectrometers. Sample heterogeneity was evaluated using hyperspectral images obtained with a resolution of 150 × 150 μm and principal component analysis. The probed sample area (5 cm(2); 24,000 pixels) to achieve representativeness was found to be equivalent to the average of 6 spectra for a 1 cm diameter probing circular window of one FT-NIR instrument. The other spectrophotometer can probe the whole sample in only one measurement. The results show that the rubber properties can be determined with very similar accuracy and precision by Partial Least Square (PLS) regression models regardless of whether HI-NIR or conventional FT-NIR produce the spectral datasets. The best Root Mean Square Errors of Prediction (RMSEPs) of external validation for MV, PI0, PI30, and PRI were 4.3, 1.8, 3.4, and 5.3%, respectively. Though the quantitative results provided by the three instruments can be considered equivalent, the hyperspectral imaging instrument presents a number of advantages, being about 6 times faster than conventional bulk spectrometers, producing robust spectral data by ensuring sample representativeness, and minimizing the effect of the presence of contaminants.