Accounting for baseline differences and measurement error in the analysis of change over time

If change over time is compared in several groups, it is important to take into account baseline values so that the comparison is carried out under the same preconditions. As the observed baseline measurements are distorted by measurement error, it may not be sufficient to include them as covariate. By fitting a longitudinal mixed-effects model to all data including the baseline observations and subsequently calculating the expected change conditional on the underlying baseline value, a solution to this problem has been provided recently so that groups with the same baseline characteristics can be compared. In this article, we present an extended approach where a broader set of models can be used. Specifically, it is possible to include any desired set of interactions between the time variable and the other covariates, and also, time-dependent covariates can be included. Additionally, we extend the method to adjust for baseline measurement error of other time-varying covariates. We apply the methodology to data from the Swiss HIV Cohort Study to address the question if a joint infection with HIV-1 and hepatitis C virus leads to a slower increase of CD4 lymphocyte counts over time after the start of antiretroviral therapy.

The von Mises-Fisher distribution of the first exit point from the hypersphere of the drifted Brownian motion and the density of the first exit time

A characterization is provided for the von Mises–Fisher random variable, in terms of first exit point from the unit hypersphere of the drifted Wiener process. Laplace transform formulae for the first exit time from the unit hypersphere of the drifted Wiener process are provided. Post representations in terms of Bell polynomials are provided for the densities of the first exit times from the circle and from the sphere.

A sequential Cox approach for estimating the causal effect of treatment in the presence of time-dependent confounding applied to data from the Swiss HIV Cohort Study

When estimating the effect of treatment on HIV using data from observational studies, standard methods may produce biased estimates due to the presence of time-dependent confounders. Such confounding can be present when a covariate, affected by past exposure, is both a predictor of the future exposure and the outcome. One example is the CD4 cell count, being a marker for disease progression for HIV patients, but also a marker for treatment initiation and influenced by treatment. Fitting a marginal structural model (MSM) using inverse probability weights is one way to give appropriate adjustment for this type of confounding. In this paper we study a simple and intuitive approach to estimate similar treatment effects, using observational data to mimic several randomized controlled trials. Each 'trial' is constructed based on individuals starting treatment in a certain time interval. An overall effect estimate for all such trials is found using composite likelihood inference. The method offers an alternative to the use of inverse probability of treatment weights, which is unstable in certain situations. The estimated parameter is not identical to the one of an MSM, it is conditioned on covariate values at the start of each mimicked trial. This allows the study of questions that are not that easily addressed fitting an MSM. The analysis can be performed as a stratified weighted Cox analysis on the joint data set of all the constructed trials, where each trial is one stratum. The model is applied to data from the Swiss HIV cohort study.

Ragu: a free tool for the analysis of EEG and MEG event-related scalp field data using global randomization statistics

We present a program (Ragu; Randomization Graphical User interface) for statistical analyses of multichannel event-related EEG and MEG experiments. Based on measures of scalp field differences including all sensors, and using powerful, assumption-free randomization statistics, the program yields robust, physiologically meaningful conclusions based on the entire, untransformed, and unbiased set of measurements. Ragu accommodates up to two within-subject factors and one between-subject factor with multiple levels each. Significance is computed as function of time and can be controlled for type II errors with overall analyses. Results are displayed in an intuitive visual interface that allows further exploration of the findings. A sample analysis of an ERP experiment illustrates the different possibilities offered by Ragu. The aim of Ragu is to maximize statistical power while minimizing the need for a-priori choices of models and parameters (like inverse models or sensors of interest) that interact with and bias statistics.

The risk of early mortality of polytrauma patients associated to ISS, NISS, APACHE II values and prothrombin time

BACKGROUND: The early hemodynamic normalization of polytrauma patients may lead to better survival outcomes. The aim of this study was to assess the diagnostic quality of trauma and physiological scores from widely used scoring systems in polytrauma patients. METHODS: In total, 770 patients with ISS > 16 who were admitted to a trauma center within the first 24 hours after injury were included in this retrospective study. The patients were subdivided into three groups: those who died on the day of admission, those who died within the first three days, and those who survived for longer than three days. ISS, NISS, APACHE II score, and prothrombin time were recorded at admission. RESULTS: The descriptive statistics for early death in polytrauma patients who died on the day of admission, 1--3 days after admission, and > 3 days after admission were: ISS of 41.0, 34.0, and 29.0, respectively; NISS of 50.0, 50.0, and 41.0, respectively; APACHE II score of 30.0, 25.0, and 15.0, respectively; and prothrombin time of 37.0%, 56.0%, and 84%, respectively. These data indicate that prothrombin time (AUC: 0.89) and APACHE II (AUC: 0.88) have the greatest prognostic utility for early death. CONCLUSION: The estimated densities of the scores may suggest a direction for resuscitative procedures in polytrauma patients.Trial registration: "Retrospektive Analysen in der Chirurgischen Intensivmedizin" StV01-2008.http://www.kek.zh.ch/internet/gesundheitsdirektion/kek/de/home.html.

Non-exponential return time distributions for vorticity extremes explained by fractional Poisson processes

Serial correlation of extreme midlatitude cyclones observed at the storm track exits is explained by deviations from a Poisson process. To model these deviations, we apply fractional Poisson processes (FPPs) to extreme midlatitude cyclones, which are defined by the 850 hPa relative vorticity of the ERA interim reanalysis during boreal winter (DJF) and summer (JJA) seasons. Extremes are defined by a 99% quantile threshold in the grid-point time series. In general, FPPs are based on long-term memory and lead to non-exponential return time distributions. The return times are described by a Weibull distribution to approximate the Mittag–Leffler function in the FPPs. The Weibull shape parameter yields a dispersion parameter that agrees with results found for midlatitude cyclones. The memory of the FPP, which is determined by detrended fluctuation analysis, provides an independent estimate for the shape parameter. Thus, the analysis exhibits a concise framework of the deviation from Poisson statistics (by a dispersion parameter), non-exponential return times and memory (correlation) on the basis of a single parameter. The results have potential implications for the predictability of extreme cyclones.

MR Angiography at 3 Tesla to Assess Proximal Internal Carotid Artery Stenoses: Contrast-Enhanced or 3D Time-of-Flight MR Angiography?

PURPOSE The aim of this study was to compare the diagnostic accuracy of 3D time-of-flight (TOF-MRA) and contrast-enhanced (CE-MRA) magnetic resonance angiography at 3 T for detection and quantification of proximal high-grade stenosis using multidetector computed tomography angiography (MDCTA) as reference standard. METHODS The institutional ethics committee approved this prospective study. A total of 41 patients suspected of having internal carotid artery (ICA) stenosis underwent both MDCTA and MRA. CE-MRA and TOF-MRA were performed using a 3.0-T imager with a dedicated eight-element cervical coil. ICA stenoses were measured according to the North American Symptomatic Carotid Endarterectomy Trial criteria and categorized as 0-25 % (minimal), 25-50 % (mild), 50-69 % (moderate), 70-99 % (high grade), and 100 % (occlusion). Sensitivity and specificity for the detection of high-grade ICA stenoses (70-99 %) and ICA occlusions were determined. In addition, intermodality agreement was assessed with κ-statistics for detection of high-grade ICA stenoses (70-99 %) and ICA occlusions. RESULTS A total of 80 carotid arteries of 41 patients were reviewed. Two previously stented ICAs were excluded from analysis. On MDCTA, 7 ICAs were occluded, 12 ICAs presented with and 63 without a high-grade ICA stenosis (70-99 %). For detecting 70-99 % stenosis, both 3D TOF-MRA and CE-MRA were 91.7 % sensitive and 98.5 % specific, respectively. Both MRA techniques were highly sensitive (100 %), and specific (CE-MRA, 100 %; TOF-MRA, 98.7 %) for the detection of ICA occlusion. However, TOF-MRA misclassified one high-grade stenosis as occlusion. Intermodality agreement for detection of 70-99 % ICA stenoses was excellent between TOF-MRA and CE-MRA [κ = 0.902, 95 % confidence interval (CI) = 0.769-1.000], TOF-MRA and MDCTA (κ = 0.902, 95 % CI = 0.769-1.000), and CE-MRA and MDCTA (κ = 0.902, 95 % CI = 0.769-1.000). CONCLUSION Both 3D TOF-MRA and CE-MRA at 3 T are reliable tools for detecting high-grade proximal ICA stenoses (70-99 %). 3D TOF-MRA might misclassify pseudo-occlusions as complete occlusions. If there are no contraindications for CE-MRA, CE-MRA is recommended as primary MR imaging modality.

Comment on “Low time resolution analysis of polar ice cores cannot detect impulsive nitrate events” by D.F. Smart et al.

Smart et al. (2014) suggested that the detection of nitrate spikes in polar ice cores from solar energetic particle (SEP) events could be achieved if an analytical system with sufficiently high resolution was used. Here we show that the spikes they associate with SEP events are not reliably recorded in cores from the same location, even when the resolution is clearly adequate. We explain the processes that limit the effective resolution of ice cores. Liquid conductivity data suggest that the observed spikes are associated with sodium or another nonacidic cation, making it likely that they result from deposition of sea salt or similar aerosol that has scavenged nitrate, rather than from a primary input of nitrate in the troposphere. We consider that there is no evidence at present to support the identification of any spikes in nitrate as representing SEP events. Although such events undoubtedly create nitrate in the atmosphere, we see no plausible route to using nitrate spikes to document the statistics of such events.