Novel statistical approaches for non-normal censored immunological data: analysis of cytokine and gene expression data

Background: For several immune-mediated diseases, immunological analysis will become more complex in the future with datasets in which cytokine and gene expression data play a major role. These data have certain characteristics that require sophisticated statistical analysis such as strategies for non-normal distribution and censoring. Additionally, complex and multiple immunological relationships need to be adjusted for potential confounding and interaction effects. Objective: We aimed to introduce and apply different methods for statistical analysis of non-normal censored cytokine and gene expression data. Furthermore, we assessed the performance and accuracy of a novel regression approach in order to allow adjusting for covariates and potential confounding. Methods: For non-normally distributed censored data traditional means such as the Kaplan-Meier method or the generalized Wilcoxon test are described. In order to adjust for covariates the novel approach named Tobit regression on ranks was introduced. Its performance and accuracy for analysis of non-normal censored cytokine/gene expression data was evaluated by a simulation study and a statistical experiment applying permutation and bootstrapping. Results: If adjustment for covariates is not necessary traditional statistical methods are adequate for non-normal censored data. Comparable with these and appropriate if additional adjustment is required, Tobit regression on ranks is a valid method. Its power, type-I error rate and accuracy were comparable to the classical Tobit regression. Conclusion: Non-normally distributed censored immunological data require appropriate statistical methods. Tobit regression on ranks meets these requirements and can be used for adjustment for covariates and potential confounding in large and complex immunological datasets

Endoscopic internal drainage for the management of leak, fistula, and collection after sleeve gastrectomy: our experience in 617 consecutive patients

Background: Endoscopy plays a pivotal role in the management of adverse events (AE) following bariatric surgery. Leaks, fistulae, and post-operative collection after sleeve gastrectomy (SG) may occur in up to 10% of cases. Objectives: To evaluate the efficacy and safety of endoscopic internal drainage (EID) for the management of leak, fistula, and collection following SG. Setting: Retrospective, observational, single center study on patients referred from several bariatric surgery departments to an endoscopic referral center. Methods: EID was used as first-line treatment for the management of leaks, fistulae, and collections. Leaks and fistulae were treated with double pigtail stent (DPS) deployment in order to guarantee internal drainage and second intention cavity obliteration. Collections were treated with endoscropic ultrasound (EUS)–guided deployment of DPS or lumen apposing metal stents. Results: A total of 617 patients (83.3% female; mean age, 43.1 yr) were enrolled in the study for leak (n = 300, 48.6%), fistula (n = 285, 46.2%), and collection (n = 32, 5.2%). Median follow-up was 19.5 months. Overall clinical success was 84.7% whereas 15.3% of cases required revisional surgery after EID failure. Clinical success according to type of AE was 89.5%, 78.5%, and 90% for leak, fistula, and collection, respectively. A total of 10 of 547 (1.8%) presented a recurrence during follow-up. A total of 28 (4.5%) AE related to the endoscopic treatment occurred. At univariate logistic regression predictors of failure were: fistula (OR 2.012), combined endoscopic approach (OR 2.319), need for emergency surgery (OR 1.755), and previous endoscopic treatment (OR 4.818). Conclusion: Early EID for the management of leak, fistula, and post-operative collection after SG seems a safe and effective first-line approach with good long-term results

Search for R-parity Violating Supersymmetry in Dimuon and Four-Jets Channel

We present results of a search for R-parity-violating decay of the neutralino chi_1^0, taken to be the Lightest Supersymmetric Particle. It is assumed that this decay proceeds through one of the lepton-number violating couplings lambda-prime_2jk (j=1,2; k=1,2,3). This search is based on 77.5 pb-1 of data, collected by the D0 experiment at the Fermilab Tevatron in ppbar collisions at a center of mass energy of 1.8 TeV in 1992-1995.Comment: 10 pages, 3 figure

A measurement of the W boson mass using large rapidity electrons

We present a measurement of the W boson mass using data collected by the D0 experiment at the Fermilab Tevatron during 1994--1995. We identify W bosons by their decays to e-nu final states where the electron is detected in a forward calorimeter. We extract the W boson mass, Mw, by fitting the transverse mass and transverse electron and neutrino momentum spectra from a sample of 11,089 W -> e nu decay candidates. We use a sample of 1,687 dielectron events, mostly due to Z -> ee decays, to constrain our model of the detector response. Using the forward calorimeter data, we measure Mw = 80.691 +- 0.227 GeV. Combining the forward calorimeter measurements with our previously published central calorimeter results, we obtain Mw = 80.482 +- 0.091 GeV

Spin Correlation in tt-bar Production from pp-bar Collisions at sqrt(s)=1.8 TeV

The D0 collaboration has performed a study of spin correlation in tt-bar production for the process tt-bar to bb-bar W^+W^-, where the W bosons decay to e-nu or mu-nu. A sample of six events was collected during an exposure of the D0 detector to an integrated luminosity of approximately 125 pb^-1 of sqrt{s}=1.8 TeV pp-bar collisions. The standard model (SM) predicts that the short lifetime of the top quark ensures the transmission of any spin information at production to the tt-bar decay products. The degree of spin correlation is characterized by a correlation coefficient k. We find that k>-0.25 at the 68% confidence level, in agreement with the SM prediction of k=0.88.Comment: Submitted to PRL, Added references, minor changes to tex

Search for the Higgs boson in events with missing transverse energy and b quark jets produced in proton-antiproton collisions at s**(1/2)=1.96 TeV

We search for the standard model Higgs boson produced in association with an electroweak vector boson in events with no identified charged leptons, large imbalance in transverse momentum, and two jets where at least one contains a secondary vertex consistent with the decay of b hadrons. We use ~1 fb-1 integrated luminosity of proton-antiproton collisions at s**(1/2)=1.96 TeV recorded by the CDF II experiment at the Tevatron. We find 268 (16) single (double) b-tagged candidate events, where 248 +/- 43 (14.4 +/- 2.7) are expected from standard model background processes. We place 95% confidence level upper limits on the Higgs boson production cross section for several Higgs boson masses ranging from 110 GeV/c2 to 140 GeV/c2. For a mass of 115 GeV/c2 the observed (expected) limit is 20.4 (14.2) times the standard model prediction.Comment: 8 pages, 2 figures, submitted to Phys. Rev. Let

Search for right-handed W bosons in top quark decay

We present a measurement of the fraction f+ of right-handed W bosons produced in top quark decays, based on a candidate sample of $t\bar{t}$ events in the lepton+jets decay mode. These data correspond to an integrated luminosity of 230pb^-1, collected by the DO detector at the Fermilab Tevatron $p\bar{p}$ Collider at sqrt(s)=1.96 TeV. We use a constrained fit to reconstruct the kinematics of the $t\bar{t}$ and decay products, which allows for the measurement of the leptonic decay angle $\theta^*$ for each event. By comparing the $\cos\theta^*$ distribution from the data with those for the expected background and signal for various values of f+, we find f+=0.00+-0.13(stat)+-0.07(syst). This measurement is consistent with the standard model prediction of f+=3.6x10^-4.Comment: Submitted to Physical Review D Rapid Communications 7 pages, 3 figure

Measurement of Semileptonic Branching Fractions of B Mesons to Narrow D** States

Using the data accumulated in 2002-2004 with the DO detector in proton-antiproton collisions at the Fermilab Tevatron collider with centre-of-mass energy 1.96 TeV, the branching fractions of the decays B -> \bar{D}_1^0(2420) \mu^+ \nu_\mu X and B -> \bar{D}_2^{*0}(2460) \mu^+ \nu_\mu X and their ratio have been measured: BR(\bar{b}->B) \cdot BR(B-> \bar{D}_1^0 \mu^+ \nu_\mu X) \cdot BR(\bar{D}_1^0 -> D*- pi+) = (0.087+-0.007(stat)+-0.014(syst))%; BR(\bar{b}->B)\cdot BR(B->D_2^{*0} \mu^+ \nu_\mu X) \cdot BR(\bar{D}_2^{*0} -> D*- \pi^+) = (0.035+-0.007(stat)+-0.008(syst))%; and (BR(B -> \bar{D}_2^{*0} \mu^+ \nu_\mu X)BR(D2*0->D*- pi+)) / (BR(B -> \bar{D}_1^{0} \mu^+ \nu_\mu X)\cdot BR(\bar{D}_1^{0}->D*- \pi^+)) = 0.39+-0.09(stat)+-0.12(syst), where the charge conjugated states are always implied.Comment: submitted to Phys. Rev. Let

Improved W boson mass measurement with the D0 detector

We have measured the W boson mass using the D0 detector and a data sample of 82 pb^-1 from the Tevatron collider. This measurement used W -> e nu decays, where the electron is close to a boundary of a central electromagnetic calorimeter module. Such 'edge' electrons have not been used in any previous D0 analysis, and represent a 14% increase in the W boson sample size. For these electrons, new response and resolution parameters are determined, and revised backgrounds and underlying event energy flow measurements are made. When the current measurement is combined with previous D0 W boson mass measurements, we obtain M_W = 80.483 +/- 0.084 GeV. The 8% improvement from the previous D0 measurement is primarily due to the improved determination of the response parameters for non-edge electrons using the sample of Z bosons with non-edge and edge electrons.Comment: submitted to Phys. Rev. D; 20 pages, 18 figures, 9 table

Measurement of the Lifetime Difference in the B_s^0 System

We present a study of the decay B_s^0 -> J/psi phi We obtain the CP-odd fraction in the final state at time zero, R_perp = 0.16 +/- 0.10 (stat) +/- 0.02 (syst), the average lifetime of the (B_s, B_sbar) system, tau (B_s^0) =1.39^{+0.13}_{-0.16} (stat) ^{+0.01}_{-0.02} (syst) ps, and the relative width difference between the heavy and light mass eigenstates, Delta Gamma/Gamma = (Gamma_L - Gamma_H)/Gamma =0.24^{+0.28}_{-0.38} (stat) ^{+0.03}_{-0.04} (syst). With the additional constraint from the world average of the B_s^0$lifetime measurements using semileptonic decays, we find tau (B_s^0)= 1.39 +/- 0.06 ~ps and Delta Gamma/\Gamma = 0.25^{+0.14}_{-0.15}. For the ratio of the B_s^0 and B^0 lifetimes we obtain tau(B_s^0)/tau(B^0)} = 0.91 +/- 0.09 (stat) +/- 0.003 (syst).Comment: submitted to Phys. Rev. Lett. FERMILAB-PUB-05-324-