POOR PERFORMANCE OF BOOTSTRAP CONFIDENCE INTERVALS FOR THE LOCATION OF A QUANTITATIVE TRAIT LOUCS

The aim of many genetic studies is to locate the genomic regions (called quantitative trait loci, QTLs) that contribute to variation in a quantitative trait (such as body weight). Confidence intervals for the locations of QTLs are particularly important for the design of further experiments to identify the gene or genes responsible for the effect. Likelihood support intervals are the most widely used method to obtain confidence intervals for QTL location, but the non-parametric bootstrap has also been recommended. Through extensive computer simulation, we show that bootstrap confidence intervals are poorly behaved and so should not be used in this context. The profile likelihood (or LOD curve) for QTL location has a tendency to peak at genetic markers, and so the distribution of the maximum likelihood estimate (MLE) of QTL location has the unusual feature of point masses at genetic markers; this contributes to the poor behavior of the bootstrap. Likelihood support intervals and approximate Bayes credible intervals, on the other hand, are shown to behave appropriately.

Statistical evaluation of time-dependent metabolite concentrations: estimation of post-mortem intervals based on in situ 1H-MRS of the brain

Knowledge of the time interval from death (post-mortem interval, PMI) has an enormous legal, criminological and psychological impact. Aiming to find an objective method for the determination of PMIs in forensic medicine, 1H-MR spectroscopy (1H-MRS) was used in a sheep head model to follow changes in brain metabolite concentrations after death. Following the characterization of newly observed metabolites (Ith et al., Magn. Reson. Med. 2002; 5: 915-920), the full set of acquired spectra was analyzed statistically to provide a quantitative estimation of PMIs with their respective confidence limits. In a first step, analytical mathematical functions are proposed to describe the time courses of 10 metabolites in the decomposing brain up to 3 weeks post-mortem. Subsequently, the inverted functions are used to predict PMIs based on the measured metabolite concentrations. Individual PMIs calculated from five different metabolites are then pooled, being weighted by their inverse variances. The predicted PMIs from all individual examinations in the sheep model are compared with known true times. In addition, four human cases with forensically estimated PMIs are compared with predictions based on single in situ MRS measurements. Interpretation of the individual sheep examinations gave a good correlation up to 250 h post-mortem, demonstrating that the predicted PMIs are consistent with the data used to generate the model. Comparison of the estimated PMIs with the forensically determined PMIs in the four human cases shows an adequate correlation. Current PMI estimations based on forensic methods typically suffer from uncertainties in the order of days to weeks without mathematically defined confidence information. In turn, a single 1H-MRS measurement of brain tissue in situ results in PMIs with defined and favorable confidence intervals in the range of hours, thus offering a quantitative and objective method for the determination of PMIs.

Quantifying the periodicty of Inter Movement Intervals of Leg Movements

ims: Periodic leg movements in sleep (PLMS) are a frequent finding in polysomnography. Most patients with restless legs syndrome (RLS) display PLMS. However, since PLMS are also often recorded in healthy elderly subjects, the clinical significance of PLMS is still discussed controversially. Leg movements are seen concurrently with arousals in obstructive sleep apnoea (OSA) may also appear periodically. Quantitative assessment of the periodicity of LM/PLM as measured by inter movement intervals (IMI) is difficult. This is mainly due to influencing factors like sleep architecture and sleep stage, medication, inter and intra patient variability, the arbitrary amplitude and sequence criteria which tend to broaden the IMI distributions or make them even multi-modal. Methods: Here a statistical method is presented that enables eliminating such effects from the raw data before analysing the statistics of IMI. Rather than studying the absolute size of IMI (measured in seconds) we focus on the shape of their distribution (suitably normalized IMI). To this end we employ methods developed in Random Matrix Theory (RMT). Patients: The periodicity of leg movements (LM) of four patient groups (10 to 15 each) showing LM without PLMS (group 1), OSA without PLMS (group 2), PLMS and OSA (group 3) as well as PLMS without OSA (group 4) are compared. Results: The IMI of patients without PLMS (groups 1 and 2) and with PLMS (groups 3 and 4) are statistically different. In patients without PLMS the distribution of normalized IMI resembles closely the one of random events. In contrary IMI of PLMS patients show features of periodic systems (e.g. a pendulum) when studied in normalized manner. Conclusions: For quantifying PLMS periodicity proper normalization of the IMI is crucial. Without this procedure important features are hidden when grouping LM/PLM over whole nights or across patients. The clinical significance of PLMS might be eluded when properly separating random LM from LM that show features of periodic systems.

Reference intervals for common carotid intima-media thickness measured with echotracking: relation with risk factors

AIMS Common carotid artery intima-media thickness (CCIMT) is widely used as a surrogate marker of atherosclerosis, given its predictive association with cardiovascular disease (CVD). The interpretation of CCIMT values has been hampered by the absence of reference values, however. We therefore aimed to establish reference intervals of CCIMT, obtained using the probably most accurate method at present (i.e. echotracking), to help interpretation of these measures. METHODS AND RESULTS We combined CCIMT data obtained by echotracking on 24 871 individuals (53% men; age range 15-101 years) from 24 research centres worldwide. Individuals without CVD, cardiovascular risk factors (CV-RFs), and BP-, lipid-, and/or glucose-lowering medication constituted a healthy sub-population (n = 4234) used to establish sex-specific equations for percentiles of CCIMT across age. With these equations, we generated CCIMT Z-scores in different reference sub-populations, thereby allowing for a standardized comparison between observed and predicted ('normal') values from individuals of the same age and sex. In the sub-population without CVD and treatment (n = 14 609), and in men and women, respectively, CCIMT Z-scores were independently associated with systolic blood pressure [standardized βs 0.19 (95% CI: 0.16-0.22) and 0.18 (0.15-0.21)], smoking [0.25 (0.19-0.31) and 0.11 (0.04-0.18)], diabetes [0.19 (0.05-0.33) and 0.19 (0.02-0.36)], total-to-HDL cholesterol ratio [0.07 (0.04-0.10) and 0.05 (0.02-0.09)], and body mass index [0.14 (0.12-0.17) and 0.07 (0.04-0.10)]. CONCLUSION We estimated age- and sex-specific percentiles of CCIMT in a healthy population and assessed the association of CV-RFs with CCIMT Z-scores, which enables comparison of IMT values for (patient) groups with different cardiovascular risk profiles, helping interpretation of such measures obtained both in research and clinical settings.

Changes in the PQRST intervals and heart rate variability associated with rewarming in two newborns undergoing hypothermia therapy.

BACKGROUND: Little is known about the effects of hypothermia therapy and subsequent rewarming on the PQRST intervals and heart rate variability (HRV) in term newborns with hypoxic-ischemic encephalopathy (HIE). OBJECTIVES: This study describes the changes in the PQRST intervals and HRV during rewarming to normal core body temperature of 2 newborns with HIE after hypothermia therapy. METHODS: Within 6 h after birth, 2 newborns with HIE were cooled to a core body temperature of 33.5 degrees C for 72 h using a cooling blanket, followed by gradual rewarming (0.5 degrees C per hour) until the body temperature reached 36.5 degrees C. Custom instrumentation recorded the electrocardiogram from the leads used for clinical monitoring of vital signs. Generalized linear mixed models were calculated to estimate temperature-related changes in PQRST intervals and HRV. Results: For every 1 degrees C increase in body temperature, the heart rate increased by 9.2 bpm (95% CI 6.8-11.6), the QTc interval decreased by 21.6 ms (95% CI 17.3-25.9), and low and high frequency HRV decreased by 0.480 dB (95% CI 0.052-0.907) and 0.938 dB (95% CI 0.460-1.416), respectively. CONCLUSIONS: Hypothermia-induced changes in the electrocardiogram should be monitored carefully in future studies.

Reporting of statistical results in prosthodontic and implantology journals: p values or confidence intervals?

PURPOSE Confidence intervals (CIs) are integral to the interpretation of the precision and clinical relevance of research findings. The aim of this study was to ascertain the frequency of reporting of CIs in leading prosthodontic and dental implantology journals and to explore possible factors associated with improved reporting. MATERIALS AND METHODS Thirty issues of nine journals in prosthodontics and implant dentistry were accessed, covering the years 2005 to 2012: The Journal of Prosthetic Dentistry, Journal of Oral Rehabilitation, The International Journal of Prosthodontics, The International Journal of Periodontics & Restorative Dentistry, Clinical Oral Implants Research, Clinical Implant Dentistry and Related Research, The International Journal of Oral & Maxillofacial Implants, Implant Dentistry, and Journal of Dentistry. Articles were screened and the reporting of CIs and P values recorded. Other information including study design, region of authorship, involvement of methodologists, and ethical approval was also obtained. Univariable and multivariable logistic regression was used to identify characteristics associated with reporting of CIs. RESULTS Interrater agreement for the data extraction performed was excellent (kappa = 0.88; 95% CI: 0.87 to 0.89). CI reporting was limited, with mean reporting across journals of 14%. CI reporting was associated with journal type, study design, and involvement of a methodologist or statistician. CONCLUSIONS Reporting of CI in implant dentistry and prosthodontic journals requires improvement. Improved reporting will aid appraisal of the clinical relevance of research findings by providing a range of values within which the effect size lies, thus giving the end user the opportunity to interpret the results in relation to clinical practice.

The effects of type of interval, sensory modality, base duration, and psychophysical task on the discrimination of brief time intervals

The present study was designed to investigate the influences of type of psychophysical task (two-alternative forced-choice [2AFC] and reminder tasks), type of interval (filled vs. empty), sensory modality (auditory vs. visual), and base duration (ranging from 100 through 1,000 ms) on performance on duration discrimination. All of these factors were systematically varied in an experiment comprising 192 participants. This approach allowed for obtaining information not only on the general (main) effect of each factor alone, but also on the functional interplay and mutual interactions of some or all of these factors combined. Temporal sensitivity was markedly higher for auditory than for visual intervals, as well as for the reminder relative to the 2AFC task. With regard to base duration, discrimination performance deteriorated with decreasing base durations for intervals below 400 ms, whereas longer intervals were not affected. No indication emerged that overall performance on duration discrimination was influenced by the type of interval, and only two significant interactions were apparent: Base Duration × Type of Interval and Base Duration × Sensory Modality. With filled intervals, the deteriorating effect of base duration was limited to very brief base durations, not exceeding 100 ms, whereas with empty intervals, temporal discriminability was also affected for the 200-ms base duration. Similarly, the performance decrement observed with visual relative to auditory intervals increased with decreasing base durations. These findings suggest that type of task, sensory modality, and base duration represent largely independent sources of variance for performance on duration discrimination that can be accounted for by distinct nontemporal mechanisms.

Random Sampling with Interspike-Intervals of the Exponential Integrate and Fire Neuron: A Computational Interpretation of UP-States

Oscillations between high and low values of the membrane potential (UP and DOWN states respectively) are an ubiquitous feature of cortical neurons during slow wave sleep and anesthesia. Nevertheless, a surprisingly small number of quantitative studies have been conducted only that deal with this phenomenon’s implications for computation. Here we present a novel theory that explains on a detailed mathematical level the computational benefits of UP states. The theory is based on random sampling by means of interspike intervals (ISIs) of the exponential integrate and fire (EIF) model neuron, such that each spike is considered a sample, whose analog value corresponds to the spike’s preceding ISI. As we show, the EIF’s exponential sodium current, that kicks in when balancing a noisy membrane potential around values close to the firing threshold, leads to a particularly simple, approximative relationship between the neuron’s ISI distribution and input current. Approximation quality depends on the frequency spectrum of the current and is improved upon increasing the voltage baseline towards threshold. Thus, the conceptually simpler leaky integrate and fire neuron that is missing such an additional current boost performs consistently worse than the EIF and does not improve when voltage baseline is increased. For the EIF in contrast, the presented mechanism is particularly effective in the high-conductance regime, which is a hallmark feature of UP-states. Our theoretical results are confirmed by accompanying simulations, which were conducted for input currents of varying spectral composition. Moreover, we provide analytical estimations of the range of ISI distributions the EIF neuron can sample from at a given approximation level. Such samples may be considered by any algorithmic procedure that is based on random sampling, such as Markov Chain Monte Carlo or message-passing methods. Finally, we explain how spike-based random sampling relates to existing computational theories about UP states during slow wave sleep and present possible extensions of the model in the context of spike-frequency adaptation.

LC-MS/MS based assay and reference intervals in children and adolescents for oxysterols elevated in Niemann-Pick diseases

BACKGROUND Niemann-Pick type C (NP-C) is a rare progressive neurodegenerative lipid storage disorder with heterogeneous clinical presentation and challenging diagnostic procedures. Recently oxysterols have been reported to be specific biomarkers for NP-C but knowledge on the intra-individual variation and on reference intervals in children and adolescents are lacking. METHODS We established a LC-MS/MS assay to measure Cholestane-3β, 5α, 6β-triol (C-triol) and 7-Ketocholesterol (7-KC) following Steglich esterification. To assess reference intervals and intra-individual variation we determined oxysterols in 148 children and adolescents from 0 to 18 years and repeat measurements in 19 of them. RESULTS The reported method is linear (r>0.99), sensitive (detection limit of 0.03 ng/mL [0.07 nM] for C-triol, and 0.54 ng/mL [1.35 nM] for 7-KC) and precise, with an intra-day imprecision of 4.8% and 4.1%, and an inter-day imprecision of 7.0% and 11.0% for C-triol (28 ng/ml, 67 nM) and 7-KC (32 ng/ml, 80 nM), respectively. Recoveries for 7-KC and C-triol range between 93% and 107%. The upper reference limit obtained for C-triol is 40.4 ng/mL (95% CI: 26.4-61.7 ng/mL, 96.0 nM, 95% CI: 62.8-146.7 nM) and 75.0 ng/mL for 7-KC (95% CI: 55.5-102.5 ng/mL, 187.2 nM, 95% CI: 138.53-255.8 nM), with no age or gender dependency. Both oxysterols have a broad intra-individual variation of 46%±23% for C-triol and 52%±29% for 7-KC. Nevertheless, all Niemann-Pick patients showed increased C-triol levels including Niemann-Pick type A and B patients. CONCLUSIONS The LC-MS/MS assay is a robust assay to quantify C-triol and 7-KC in plasma with well documented reference intervals in children and adolescents to screen for NP-C in the pediatric population. In addition our results suggest that especially the C-triol is a biomarker for all three Niemann-Pick diseases.