Display of probability densities for data from a continuous distribution

Based on cumulative distribution functions, Fourier series expansion and Kolmogorov tests, we present a simple method to display probability densities for data drawn from a continuous distribution. It is often more efficient than using histograms.Comment: 5 pages, 4 figures, presented at Computer Simulation Studies XXIV, Athens, GA, 201

Prevalence of lipid abnormalities before and after introduction of lipid modifying therapy among Swedish patients with dyslipidemia (PRIMULA)

<p>Abstract</p> <p>Background</p> <p>Data on the prevalence of dyslipidemia and attainment of goal/normal lipid levels in a Swedish population are scarce. The objective of this study is to estimate the prevalence of dyslipidemia and attainment of goal/normal lipid levels in patients treated with lipid modifying therapy (LMT).</p> <p>Methods</p> <p>This longitudinal retrospective observational study covers time periods before and after treatment. Data were collected from 1994-2007 electronic patient records in public primary healthcare centers in Uppsala County, Sweden. Patients were included if they had been treated with LMT and had at least one lipid abnormality indicating dyslipidemia and if complete lipid profile data were available. Thresholds levels for lipids were defined as per Swedish guidelines.</p> <p>Results</p> <p>Among 5,424 patients included, at baseline, the prevalence of dyslipidemia (≥1 lipid abnormality) was by definition 100%, while this figure was 82% at follow-up. At baseline, 60% had elevated low-density lipoprotein (LDL-C) combined with low high-density lipoprotein (HDL-C) and/or elevated triglycerides (TG s), corresponding figure at follow-up was 36%. Low HDL-C and/or elevated TGs at follow-up remained at 69% for patients with type 2 diabetes mellitus (T2DM), 50% among patients with coronary heart disease (CHD) and 66% among patients with 10 year CHD risk >20%. Of the total sample, 40% attained goal levels of LDL-C and 18% attained goal/normal levels on all three lipid parameters.</p> <p>Conclusions</p> <p>Focusing therapy on LDL-C reduction allows 40% of patients to achieve LDL-C goal and helps reducing triglyceride levels. Almost 60% of patients experience persistent HDL-C and/or triglyceride abnormality independently of LDL-C levels and could be candidates for additional treatments.</p

Double-check: validation of diagnostic statistics for PLS-DA models in metabolomics studies

Partial Least Squares-Discriminant Analysis (PLS-DA) is a PLS regression method with a special binary ‘dummy’ y-variable and it is commonly used for classification purposes and biomarker selection in metabolomics studies. Several statistical approaches are currently in use to validate outcomes of PLS-DA analyses e.g. double cross validation procedures or permutation testing. However, there is a great inconsistency in the optimization and the assessment of performance of PLS-DA models due to many different diagnostic statistics currently employed in metabolomics data analyses. In this paper, properties of four diagnostic statistics of PLS-DA, namely the number of misclassifications (NMC), the Area Under the Receiver Operating Characteristic (AUROC), Q2 and Discriminant Q2 (DQ2) are discussed. All four diagnostic statistics are used in the optimization and the performance assessment of PLS-DA models of three different-size metabolomics data sets obtained with two different types of analytical platforms and with different levels of known differences between two groups: control and case groups. Statistical significance of obtained PLS-DA models was evaluated with permutation testing. PLS-DA models obtained with NMC and AUROC are more powerful in detecting very small differences between groups than models obtained with Q2 and Discriminant Q2 (DQ2). Reproducibility of obtained PLS-DA models outcomes, models complexity and permutation test distributions are also investigated to explain this phenomenon. DQ2 and Q2 (in contrary to NMC and AUROC) prefer PLS-DA models with lower complexity and require higher number of permutation tests and submodels to accurately estimate statistical significance of the model performance. NMC and AUROC seem more efficient and more reliable diagnostic statistics and should be recommended in two group discrimination metabolomic studies

Studying the Underlying Event in Drell-Yan and High Transverse Momentum Jet Production at the Tevatron

We study the underlying event in proton-antiproton collisions by examining the behavior of charged particles (transverse momentum pT > 0.5 GeV/c, pseudorapidity |\eta| < 1) produced in association with large transverse momentum jets (~2.2 fb-1) or with Drell-Yan lepton-pairs (~2.7 fb-1) in the Z-boson mass region (70 < M(pair) < 110 GeV/c2) as measured by CDF at 1.96 TeV center-of-mass energy. We use the direction of the lepton-pair (in Drell-Yan production) or the leading jet (in high-pT jet production) in each event to define three regions of \eta-\phi space; toward, away, and transverse, where \phi is the azimuthal scattering angle. For Drell-Yan production (excluding the leptons) both the toward and transverse regions are very sensitive to the underlying event. In high-pT jet production the transverse region is very sensitive to the underlying event and is separated into a MAX and MIN transverse region, which helps separate the hard component (initial and final-state radiation) from the beam-beam remnant and multiple parton interaction components of the scattering. The data are corrected to the particle level to remove detector effects and are then compared with several QCD Monte-Carlo models. The goal of this analysis is to provide data that can be used to test and improve the QCD Monte-Carlo models of the underlying event that are used to simulate hadron-hadron collisions.Comment: Submitted to Phys.Rev.

Observation of the Baryonic Flavor-Changing Neutral Current Decay Lambda_b -> Lambda mu+ mu-

We report the first observation of the baryonic flavor-changing neutral current decay Lambda_b -> Lambda mu+ mu- with 24 signal events and a statistical significance of 5.8 Gaussian standard deviations. This measurement uses ppbar collisions data sample corresponding to 6.8fb-1 at sqrt{s}=1.96TeV collected by the CDF II detector at the Tevatron collider. The total and differential branching ratios for Lambda_b -> Lambda mu+ mu- are measured. We find B(Lambda_b -> Lambda mu+ mu-) = [1.73+-0.42(stat)+-0.55(syst)] x 10^{-6}. We also report the first measurement of the differential branching ratio of B_s -> phi mu+ mu- using 49 signal events. In addition, we report branching ratios for B+ -> K+ mu+ mu-, B0 -> K0 mu+ mu-, and B -> K*(892) mu+ mu- decays.Comment: 8 pages, 2 figures, 4 tables. Submitted to Phys. Rev. Let

Measurement of Lifetime and Decay-Width Difference in B0s -> J/psi phi Decays

We measure the mean lifetime, tau=2/(Gamma_L+Gamma_H), and the width difference, DeltaGamma=Gamma_L-Gamma_H, of the light and heavy mass eigenstates of the B0s meson, B0sL and B0sH, in B0s -> J/psi phi decays using 1.7 fb^-1 of data collected with the CDF II detector at the Fermilab Tevatron ppbar collider. Assuming CP conservation, a good approximation for the B0s system in the Standard Model, we obtain DeltaGamma = 0.076^+0.059_-0.063 (stat.) +- 0.006 (syst.) ps^-1 and tau = 1.52 +- 0.04 (stat.) +- 0.02 (syst.) ps, the most precise measurements to date. Our constraints on the weak phase and DeltaGamma are consistent with CP conservation. Dedicated to the memory of our dear friend and colleague, Michael P. Schmid

Forward-Backward Asymmetry in Top Quark Production in ppbar Collisions at sqrt{s}=1.96 TeV

Reconstructable final state kinematics and charge assignment in the reaction ppbar->ttbar allows tests of discrete strong interaction symmetries at high energy. We define frame dependent forward-backward asymmetries for the outgoing top quark in both the ppbar and ttbar rest frames, correct for experimental distortions, and derive values at the parton-level. Using 1.9/fb of ppbar collisions at sqrt{s}=1.96 TeV recorded with the CDF II detector at the Fermilab Tevatron, we measure forward-backward top quark production asymmetries in the ppbar and ttbar rest frames of A_{FB,pp} = 0.17 +- 0.08 and A_{FB,tt} = 0.24 +- 0.14.Comment: 7 pages, 2 figures, submitted to Phys.Rev.Lett, corrected references and change of tex

Measurement of the W+W−W^+W^- Production Cross Section and Search for Anomalous WWγWW\gamma and WWZWWZ Couplings in ppˉp \bar p Collisions at s=1.96\sqrt{s} = 1.96 TeV

This Letter describes the current most precise measurement of the $W$ boson pair production cross section and most sensitive test of anomalous $WW\gamma$ and $WWZ$ couplings in $p \bar p$ collisions at a center-of-mass energy of 1.96 TeV. The $WW$ candidates are reconstructed from decays containing two charged leptons and two neutrinos, where the charged leptons are either electrons or muons. Using data collected by the CDF II detector from 3.6 fb$^{-1}$ of integrated luminosity, a total of 654 candidate events are observed with an expected background contribution of $320 \pm 47$ events. The measured total cross section is $\sigma (p \bar p \to W^+ W^- + X) = 12.1 \pm 0.9 \textrm{(stat)} ^{+1.6}_{-1.4} \textrm{(syst)}$ pb, which is in good agreement with the standard model prediction. The same data sample is used to place constraints on anomalous $WW\gamma$ and $WWZ$ couplings.Comment: submitted to Phys. Rev. Let

Search for Pair Production of Scalar Top Quarks Decaying to a tau Lepton and a b Quark in ppbar Collisions at sqrt{s}=1.96 TeV

We search for pair production of supersymmetric top quarks (~t_1), followed by R-parity violating decay ~t_1 -> tau b with a branching ratio beta, using 322 pb^-1 of ppbar collisions at sqrt{s}=1.96 TeV collected by the CDF II detector at Fermilab. Two candidate events pass our final selection criteria, consistent with the standard model expectation. We set upper limits on the cross section sigma(~t_1 ~tbar_1)*beta^2 as a function of the stop mass m(~t_1). Assuming beta=1, we set a 95% confidence level limit m(~t_1)>153 GeV/c^2. The limits are also applicable to the case of a third generation scalar leptoquark (LQ_3) decaying LQ_3 -> tau b.Comment: 7 pages, 2 eps figure