Indexes and boundaries for "quantitative significance" in statistical decisions.

Boundaries for delta, representing a "quantitatively significant" or "substantively impressive" distinction, have not been established, analogous to the boundary of alpha, usually set at 0.05, for the stochastic or probabilistic component of "statistical significance". To determine what boundaries are being used for the "quantitative" decisions, we reviewed pertinent articles in three general medical journals. For each contrast of two means, contrast of two rates, or correlation coefficient, we noted the investigators' decisions about stochastic significance, stated in P values or confidence intervals, and about quantitative significance, indicated by interpretive comments. The boundaries between impressive and unimpressive distinctions were best formed by a ratio of greater than or equal to 1.2 for the smaller to the larger mean in 546 comparisons, by a standardized increment of greater than or equal to 0.28 and odds ratio of greater than or equal to 2.2 in 392 comparisons of two rates; and by an r value of greater than or equal to 0.32 in 154 correlation coefficients. Additional boundaries were also identified for "substantially" and "highly" significant quantitative distinctions. Although the proposed boundaries should be kept flexible, indexes and boundaries for decisions about "quantitative significance" are particularly useful when a value of delta must be chosen for calculating sample size before the research is done, and when the "statistical significance" of completed research is appraised for its quantitative as well as stochastic components.

Statistical evaluation of the reproducibility and the influence of paper on the analysis of black gel pen ink using laser desorption ionisation mass spectrometry

Laser desorption ionisation mass spectrometry (LDI-MS) has demonstrated to be an excellent analytical method for the forensic analysis of inks on a questioned document. The ink can be analysed directly on its substrate (paper) and hence offers a fast method of analysis as sample preparation is kept to a minimum and more importantly, damage to the document is minimised. LDI-MS has also previously been reported to provide a high power of discrimination in the statistical comparison of ink samples and has the potential to be introduced as part of routine ink analysis. This paper looks into the methodology further and evaluates statistically the reproducibility and the influence of paper on black gel pen ink LDI-MS spectra; by comparing spectra of three different black gel pen inks on three different paper substrates. Although generally minimal, the influences of sample homogeneity and paper type were found to be sample dependent. This should be taken into account to avoid the risk of false differentiation of black gel pen ink samples. Other statistical approaches such as principal component analysis (PCA) proved to be a good alternative to correlation coefficients for the comparison of whole mass spectra.

Gene finding in the chicken genome.

BACKGROUND: Despite the continuous production of genome sequence for a number of organisms, reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularly true for genomes for which there is not a large collection of known gene sequences, such as the recently published chicken genome. We used the chicken sequence to test comparative and homology-based gene-finding methods followed by experimental validation as an effective genome annotation method. RESULTS: We performed experimental evaluation by RT-PCR of three different computational gene finders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram was computed and each component of it was evaluated. The results showed that de novo comparative methods can identify up to about 700 chicken genes with no previous evidence of expression, and can correctly extend about 40% of homology-based predictions at the 5' end. CONCLUSIONS: De novo comparative gene prediction followed by experimental verification is effective at enhancing the annotation of the newly sequenced genomes provided by standard homology-based methods.

Trabecular bone score improves fracture risk prediction in non-osteoporotic women: the OFELY study.

The use of areal bone mineral density (aBMD) for fracture prediction may be enhanced by considering bone microarchitectural deterioration. Trabecular bone score (TBS) helped in redefining a significant subset of non-osteoporotic women as a higher risk group. INTRODUCTION: TBS is an index of bone microarchitecture. Our goal was to assess the ability of TBS to predict incident fracture. METHODS: TBS was assessed in 560 postmenopausal women from the Os des Femmes de Lyon cohort, who had a lumbar spine (LS) DXA scan (QDR 4500A, Hologic) between years 2000 and 2001. During a mean follow-up of 7.8 ± 1.3 years, 94 women sustained 112 fragility fractures. RESULTS: At the time of baseline DXA scan, women with incident fracture were significantly older (70 ± 9 vs. 65 ± 8 years) and had a lower LS_aBMD and LS_TBS (both -0.4SD, p < 0.001) than women without fracture. The magnitude of fracture prediction was similar for LS_aBMD and LS_TBS (odds ratio [95 % confidence interval] = 1.4 [1.2;1.7] and 1.6 [1.2;2.0]). After adjustment for age and prevalent fracture, LS_TBS remained predictive of an increased risk of fracture. Yet, its addition to age, prevalent fracture, and LS_aBMD did not reach the level of significance to improve the fracture prediction. When using the WHO classification, 39 % of fractures occurred in osteoporotic women, 46 % in osteopenic women, and 15 % in women with T-score > -1. Thirty-seven percent of fractures occurred in the lowest quartile of LS_TBS, regardless of BMD. Moreover, 35 % of fractures that occurred in osteopenic women were classified below this LS_TBS threshold. CONCLUSION: In conclusion, LS_aBMD and LS_TBS predicted fractures equally well. In our cohort, the addition of LS_TBS to age and LS_aBMD added only limited information on fracture risk prediction. However, using the lowest quartile of LS_TBS helped in redefining a significant subset of non-osteoporotic women as a higher risk group which is important for patient management.

Capillaria spp. eggs in Patagonian archaeological sites: statistical analysis of morphometric data

Discriminant analysis was used to identify eggs of Capillaria spp. at specific level found in organic remains from an archaeological site in Patagonia, Argentina, dated of 6,540 ± 110 years before present. In order to distinguish eggshell morphology 149 eggs were measured and grouped into four arbitrary subsets. The analysis used on egg width and length discriminated them into different morphotypes (Wilks' lambda = 0.381, p < 0.05). The correlation analysis suggests that width was the most important variable to discriminate among the Capillaria spp. egg morphotypes (Pearson coefficient = 0.950, p < 0.05). The study of eggshell patterns, the relative frequency in the sample, and the morphometric data allowed us to correlate the four morphotypes with Capillaria species.

Comparison of 2 frailty indexes for prediction of falls, disability, fractures, and death in older women.

BACKGROUND: Frailty, as defined by the index derived from the Cardiovascular Health Study (CHS index), predicts risk of adverse outcomes in older adults. Use of this index, however, is impractical in clinical practice. METHODS: We conducted a prospective cohort study in 6701 women 69 years or older to compare the predictive validity of a simple frailty index with the components of weight loss, inability to rise from a chair 5 times without using arms, and reduced energy level (Study of Osteoporotic Fractures [SOF index]) with that of the CHS index with the components of unintentional weight loss, poor grip strength, reduced energy level, slow walking speed, and low level of physical activity. Women were classified as robust, of intermediate status, or frail using each index. Falls were reported every 4 months for 1 year. Disability (> or =1 new impairment in performing instrumental activities of daily living) was ascertained at 4(1/2) years, and fractures and deaths were ascertained during 9 years of follow-up. Area under the curve (AUC) statistics from receiver operating characteristic curve analysis and -2 log likelihood statistics were compared for models containing the CHS index vs the SOF index. RESULTS: Increasing evidence of frailty as defined by either the CHS index or the SOF index was similarly associated with an increased risk of adverse outcomes. Frail women had a higher age-adjusted risk of recurrent falls (odds ratio, 2.4), disability (odds ratio, 2.2-2.8), nonspine fracture (hazard ratio, 1.4-1.5), hip fracture (hazard ratio, 1.7-1.8), and death (hazard ratio, 2.4-2.7) (P < .001 for all models). The AUC comparisons revealed no differences between models with the CHS index vs the SOF index in discriminating falls (AUC = 0.61 for both models; P = .66), disability (AUC = 0.64; P = .23), nonspine fracture (AUC = 0.55; P = .80), hip fracture (AUC = 0.63; P = .64), or death (AUC = 0.72; P = .10). Results were similar when -2 log likelihood statistics were compared. CONCLUSION: The simple SOF index predicts risk of falls, disability, fracture, and death as well as the more complex CHS index and may provide a useful definition of frailty to identify older women at risk of adverse health outcomes in clinical practice.

Morphological scoring of human pronuclear zygotes for prediction of pregnancy outcome.

BACKGROUND: As embryo selection is not allowed by law in Switzerland, we need a single early scoring system to identify zygotes with high implantation potential and to select zygotes for fresh transfer or cryopreservation. The underlying aim is to maximize the cumulated pregnancy rate while limiting the number of multiple pregnancies. METHODS: In all, 613 fresh and 617 frozen-thawed zygotes were scored for proximity, orientation and centring of the pronuclei, cytoplasmic halo, and number and polarization of the nucleolar precursor bodies. From these individual scores, a cumulated pronuclear score (CPNS) was calculated. Correlation between CPNS and implantation was examined and compared between fresh and frozen-thawed zygotes. The effect of freezing on CPNS was also investigated. RESULTS: CPNS was positively associated with embryo implantation in both fresh and frozen zygotes. With similar CPNS, frozen zygotes presented implantation rates as high as those of fresh zygotes. Nucleolar precursor bodies pattern and cytoplasmic halo appeared as the most important factors predictive of implantation for both types of zygotes, while pronuclei position was specifically relevant for frozen-thawed zygotes. Freezing induced an alteration of most zygote parameters, resulting in a significantly lower CPNS and a lower pregnancy rate. CONCLUSIONS: CPNS may be used as a single prognostic tool for implantation of both fresh and frozen-thawed zygotes. Lower CPNS values of frozen-thawed zygotes may also be indicative of freezing damage to zygotes. Successful implantation of frozen zygotes despite lower CPNS suggests that they may recover after thawing and in vitro culture.

Analysis of Pleistocene paleodrainage evolution in the Po Basin (Italy) by multivariate statistical techniques

In order to obtain a high-resolution Pleistocene stratigraphy, eleven continuouslycored boreholes, 100 to 220m deep were drilled in the northern part of the PoPlain by Regione Lombardia in the last five years. Quantitative provenanceanalysis (QPA, Weltje and von Eynatten, 2004) of Pleistocene sands was carriedout by using multivariate statistical analysis (principal component analysis, PCA,and similarity analysis) on an integrated data set, including high-resolution bulkpetrography and heavy-mineral analyses on Pleistocene sands and of 250 majorand minor modern rivers draining the southern flank of the Alps from West toEast (Garzanti et al, 2004; 2006). Prior to the onset of major Alpine glaciations,metamorphic and quartzofeldspathic detritus from the Western and Central Alpswas carried from the axial belt to the Po basin longitudinally parallel to theSouthAlpine belt by a trunk river (Vezzoli and Garzanti, 2008). This scenariorapidly changed during the marine isotope stage 22 (0.87 Ma), with the onset ofthe first major Pleistocene glaciation in the Alps (Muttoni et al, 2003). PCA andsimilarity analysis from core samples show that the longitudinal trunk river at thistime was shifted southward by the rapid southward and westward progradation oftransverse alluvial river systems fed from the Central and Southern Alps.Sediments were transported southward by braided river systems as well as glacialsediments transported by Alpine valley glaciers invaded the alluvial plain.Kew words: Detrital modes; Modern sands; Provenance; Principal ComponentsAnalysis; Similarity, Canberra Distance; palaeodrainage

Statistical treatment of grain-size curves and empirical distributions: densities as compositions?

The preceding two editions of CoDaWork included talks on the possible considerationof densities as infinite compositions: Egozcue and D´ıaz-Barrero (2003) extended theEuclidean structure of the simplex to a Hilbert space structure of the set of densitieswithin a bounded interval, and van den Boogaart (2005) generalized this to the setof densities bounded by an arbitrary reference density. From the many variations ofthe Hilbert structures available, we work with three cases. For bounded variables, abasis derived from Legendre polynomials is used. For variables with a lower bound, westandardize them with respect to an exponential distribution and express their densitiesas coordinates in a basis derived from Laguerre polynomials. Finally, for unboundedvariables, a normal distribution is used as reference, and coordinates are obtained withrespect to a Hermite-polynomials-based basis.To get the coordinates, several approaches can be considered. A numerical accuracyproblem occurs if one estimates the coordinates directly by using discretized scalarproducts. Thus we propose to use a weighted linear regression approach, where all k-order polynomials are used as predictand variables and weights are proportional to thereference density. Finally, for the case of 2-order Hermite polinomials (normal reference)and 1-order Laguerre polinomials (exponential), one can also derive the coordinatesfrom their relationships to the classical mean and variance.Apart of these theoretical issues, this contribution focuses on the application of thistheory to two main problems in sedimentary geology: the comparison of several grainsize distributions, and the comparison among different rocks of the empirical distribution of a property measured on a batch of individual grains from the same rock orsediment, like their composition

Detection of ingested cocaine-filled packets: comparison of filtered back projection CT with adaptive statistical iterative reconstructed images

Purpose: To evaluate the diagnostic value and image quality of CT with filtered back projection (FBP) compared with adaptive statistical iterative reconstructed images (ASIR) in body stuffers with ingested cocaine-filled packets.Methods and Materials: Twenty-nine body stuffers (mean age 31.9 years, 3 women) suspected for ingestion of cocaine-filled packets underwent routine-dose 64-row multidetector CT with FBP (120kV, pitch 1.375, 100-300 mA and automatic tube current modulation (auto mA), rotation time 0.7sec, collimation 2.5mm), secondarily reconstructed with 30 % and 60 % ASIR. In 13 (44.83%) out of the body stuffers cocaine-filled packets were detected, confirmed by exact analysis of the faecal content including verification of the number (range 1-25). Three radiologists independently and blindly evaluated anonymous CT examinations (29 FBP-CT and 68 ASIR-CT) for the presence and number of cocaine-filled packets indicating observers' confidence, and graded them for diagnostic quality, image noise, and sharpness. Sensitivity, specificity, area under the receiver operating curve (ROC) Az and interobserver agreement between the 3 radiologists for FBP-CT and ASIR-CT were calculated.Results: The increase of the percentage of ASIR significantly diminished the objective image noise (p<0.001). Overall sensitivity and specificity for the detection of the cocaine-filled packets were 87.72% and 76.15%, respectively. The difference of ROC area Az between the different reconstruction techniques was significant (p= 0.0101), that is 0.938 for FBP-CT, 0.916 for 30 % ASIR-CT, and 0.894 for 60 % ASIR-CT.Conclusion: Despite the evident image noise reduction obtained by ASIR, the diagnostic value for detecting cocaine-filled packets decreases, depending on the applied ASIR percentage.

Statistical modeling on (α, +∞)

Observations in daily practice are sometimes registered as positive values larger then a given threshold α. The sample space is in this case the interval (α,+∞), α & 0, which can be structured as a real Euclidean space in different ways. This fact opens the door to alternative statistical models depending not only on the assumed distribution function, but also on the metric which is considered as appropriate, i.e. the way differences are measured, and thus variability

Statistical modeling on coordinates

This paper is a first draft of the principle of statistical modelling on coordinates. Several causes —which would be long to detail—have led to this situation close to the deadline for submitting papers to CODAWORK’03. The main of them is the fast development of the approach along thelast months, which let appear previous drafts as obsolete. The present paper contains the essential parts of the state of the art of this approach from my point of view. I would like to acknowledge many clarifying discussions with the group of people working in this field in Girona, Barcelona, Carrick Castle, Firenze, Berlin, G¨ottingen, and Freiberg. They have given a lot of suggestions and ideas. Nevertheless, there might be still errors or unclear aspects which are exclusively my fault. I hope this contribution serves as a basis for further discussions and new developments

A prediction rule to identify low-risk patients with pulmonary embolism.

BACKGROUND: A simple prognostic model could help identify patients with pulmonary embolism who are at low risk of death and are candidates for outpatient treatment. METHODS: We randomly allocated 15,531 retrospectively identified inpatients who had a discharge diagnosis of pulmonary embolism from 186 Pennsylvania hospitals to derivation (67%) and internal validation (33%) samples. We derived our rule to predict 30-day mortality using classification tree analysis and patient data routinely available at initial examination as potential predictor variables. We used data from a European prospective study to externally validate the rule among 221 inpatients with pulmonary embolism. We determined mortality and nonfatal adverse medical outcomes across derivation and validation samples. RESULTS: Our final model consisted of 10 patient factors (age > or = 70 years; history of cancer, heart failure, chronic lung disease, chronic renal disease, and cerebrovascular disease; and clinical variables of pulse rate > or = 110 beats/min, systolic blood pressure < 100 mm Hg, altered mental status, and arterial oxygen saturation < 90%). Patients with none of these factors were defined as low risk. The 30-day mortality rates for low-risk patients were 0.6%, 1.5%, and 0% in the derivation, internal validation, and external validation samples, respectively. The rates of nonfatal adverse medical outcomes were less than 1% among low-risk patients across all study samples. CONCLUSIONS: This simple prediction rule accurately identifies patients with pulmonary embolism who are at low risk of short-term mortality and other adverse medical outcomes. Prospective validation of this rule is important before its implementation as a decision aid for outpatient treatment.

Improving prediction of recanalization in acute large-vessel occlusive stroke.

BACKGROUND: Recanalization in acute ischemic stroke with large-vessel occlusion is a potent indicator of good clinical outcome. OBJECTIVE: To identify easily available clinical and radiologic variables predicting recanalization at various occlusion sites. METHODS: All consecutive, acute stroke patients from the Acute STroke Registry and Analysis of Lausanne (2003-2011) who had a large-vessel occlusion on computed tomographic angiography (CTA) (< 12 h) were included. Recanalization status was assessed at 24 h (range: 12-48 h) with CTA, magnetic resonance angiography, or ultrasonography. Complete and partial recanalization (corresponding to the modified Treatment in Cerebral Ischemia scale 2-3) were grouped together. Patients were categorized according to occlusion site and treatment modality. RESULTS: Among 439 patients, 51% (224) showed complete or partial recanalization. In multivariate analysis, recanalization of any occlusion site was most strongly associated with endovascular treatment, including bridging therapy (odds ratio [OR] 7.1, 95% confidence interval [CI] 2.2-23.2), and less so with intravenous thrombolysis (OR 1.6, 95% CI 1.0-2.6) and recanalization treatments performed beyond guidelines (OR 2.6, 95% CI 1.2-5.7). Clot location (large vs. intermediate) and tandem pathology (the combination of intracranial occlusion and symptomatic extracranial stenosis) were other variables discriminating between recanalizers and non-recanalizers. For patients with intracranial occlusions, the variables significantly associated with recanalization after 24 h were: baseline National Institutes of Health Stroke Scale (NIHSS) (OR 1.04, 95% CI 1.02-1.1), Alberta Stroke Program Early CT Score (ASPECTS) on initial computed tomography (OR 1.2, 95% CI 1.1-1.3), and an altered level of consciousness (OR 0.2, 95% CI 0.1-0.5). CONCLUSIONS: Acute endovascular treatment is the single most important factor promoting recanalization in acute ischemic stroke. The presence of extracranial vessel stenosis or occlusion decreases recanalization rates. In patients with intracranial occlusions, higher NIHSS score and ASPECTS and normal vigilance facilitate recanalization. Clinical use of these predictors could influence recanalization strategies in individual patients.