An empirical evaluation of imputation accuracy for association statistics reveals increased type-I error rates in genome-wide associations

Background: Genome wide association studies (GWAS) are becoming the approach of choice to identify genetic determinants of complex phenotypes and common diseases. The astonishing amount of generated data and the use of distinct genotyping platforms with variable genomic coverage are still analytical challenges. Imputation algorithms combine directly genotyped markers information with haplotypic structure for the population of interest for the inference of a badly genotyped or missing marker and are considered a near zero cost approach to allow the comparison and combination of data generated in different studies. Several reports stated that imputed markers have an overall acceptable accuracy but no published report has performed a pair wise comparison of imputed and empiric association statistics of a complete set of GWAS markers. Results: In this report we identified a total of 73 imputed markers that yielded a nominally statistically significant association at P < 10(-5) for type 2 Diabetes Mellitus and compared them with results obtained based on empirical allelic frequencies. Interestingly, despite their overall high correlation, association statistics based on imputed frequencies were discordant in 35 of the 73 (47%) associated markers, considerably inflating the type I error rate of imputed markers. We comprehensively tested several quality thresholds, the haplotypic structure underlying imputed markers and the use of flanking markers as predictors of inaccurate association statistics derived from imputed markers. Conclusions: Our results suggest that association statistics from imputed markers showing specific MAF (Minor Allele Frequencies) range, located in weak linkage disequilibrium blocks or strongly deviating from local patterns of association are prone to have inflated false positive association signals. The present study highlights the potential of imputation procedures and proposes simple procedures for selecting the best imputed markers for follow-up genotyping studies.

Measurement of the underground atmospheric muon charge ratio using the MINOS Near Detector

The magnetized MINOS Near Detector, at a depth of 225 mwe, is used to measure the atmospheric muon charge ratio. The ratio of observed positive to negative atmospheric muon rates, using 301 days of data, is measured to be 1.266 +/- 0.001(stat)(-0.014)(+0.015)(syst). This measurement is consistent with previous results from other shallow underground detectors and is 0.108 +/- 0.019(stat + syst) lower than the measurement at the functionally identical MINOS Far Detector at a depth of 2070 mwe. This increase in charge ratio as a function of depth is consistent with an increase in the fraction of muons arising from kaon decay for increasing muon surface energies.

Measurement of neutrino oscillations with the MINOS detectors in the NuMI beam

This Letter reports new results from the MINOS experiment based on a two-year exposure to muon neutrinos from the Fermilab NuMI beam. Our data are consistent with quantum-mechanical oscillations of neutrino flavor with mass splitting vertical bar Delta m(2)vertical bar = (2.43 +/- 0.13) x 10(-3) eV(2) (68% C.L.) and mixing angle sin(2)(2 theta) > 0.90 (90% C.L.). Our data disfavor two alternative explanations for the disappearance of neutrinos in flight: namely, neutrino decays into lighter particles and quantum decoherence of neutrinos, at the 3.7 and 5.7 standard-deviation levels, respectively.

Measurement of the Neutrino Mass Splitting and Flavor Mixing by MINOS

Measurements of neutrino oscillations using the disappearance of muon neutrinos from the Fermilab NuMI neutrino beam as observed by the two MINOS detectors are reported. New analysis methods have been applied to an enlarged data sample from an exposure of 7.25 x 10(20) protons on target. A fit to neutrino oscillations yields values of vertical bar Delta m(2)vertical bar = (2.32(-0.08)(+0.12) x 10(-3) eV(2) for the atmospheric mass splitting and sin(2)(2 theta) > 0.90 (90% C.L.) for the mixing angle. Pure neutrino decay and quantum decoherence hypotheses are excluded at 7 and 9 standard deviations, respectively.

Measurement of neutral mesons in p plus p collisions at root s=200 GeV and scaling properties of hadron production

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the invariant differential cross section for production of K(S)(0), omega, eta', and phi mesons in p + p collisions at root s 200 GeV. Measurements of omega and phi production in different decay channels give consistent results. New results for the omega are in agreement with previously published data and extend the measured p(T) coverage. The spectral shapes of all hadron transverse momentum distributions measured by PHENIX are well described by a Tsallis distribution functional form with only two parameters, n and T, determining the high-p(T) and characterizing the low-p(T) regions of the spectra, respectively. The values of these parameters are very similar for all analyzed meson spectra, but with a lower parameter T extracted for protons. The integrated invariant cross sections calculated from the fitted distributions are found to be consistent with existing measurements and with statistical model predictions.

Measurement of transverse single-spin asymmetries for J/psi production in polarized p plus p collisions at root s=200 GeV

We report the first measurement of transverse single-spin asymmetries in J/psi production from transversely polarized p + p collisions at root s = 200 GeV with data taken by the PHENIX experiment in 2006 and 2008. The measurement was performed over the rapidity ranges 1.2 < vertical bar y vertical bar < 2.2 and vertical bar y vertical bar < 0.35 for transverse momenta up to 6 GeV/c. J/psi production at the Relativistic Heavy Ion Collider is dominated by processes involving initial-state gluons, and transverse single-spin asymmetries of the J/psi can provide access to gluon dynamics within the nucleon. Such asymmetries may also shed light on the long-standing question in QCD of the J/psi production mechanism. Asymmetries were obtained as a function of J/psi transverse momentum and Feynman-x, with a value of -0.086 +/- 0.026(stat) +/- 0.003(syst) in the forward region. This result suggests possible nonzero trigluon correlation functions in transversely polarized protons and, if well defined in this reaction, a nonzero gluon Sivers distribution function.

High p(T) direct photon and pi(0) triggered azimuthal jet correlations and measurement of k(T) for isolated direct photons in p plus p collisions at root s=200 GeV

Correlations of charged hadrons of 1< p(T) < 10 Gev/c with high pT direct photons and pi(0) mesons in the range 5< p(T) < 15 Gev/c are used to study jet fragmentation in the gamma + jet and dijet channels, respectively. The magnitude of the partonic transverse momentum, k(T), is obtained by comparing to a model incorporating a Gaussian kT smearing. The sensitivity of the associated charged hadron spectra to the underlying fragmentation function is tested and the data are compared to calculations using recent global fit results. The shape of the direct photon-associated hadron spectrum as well as its charge asymmetry are found to be consistent with a sample dominated by quark-gluon Compton scattering. No significant evidence of fragmentation photon correlated production is observed within experimental uncertainties.

Detailed measurement of the e(+)e(-) pair continuum in p plus p and Au plus Au collisions at root s(NN)=200 GeV and implications for direct photon production

PHENIX has measured the e(+)e(-) pair continuum in root s(NN) = 200 GeV Au+Au and p+p collisions over a wide range of mass and transverse momenta. The e(+)e(-) yield is compared to the expectations from hadronic sources, based on PHENIX measurements. In the intermediate-mass region, between the masses of the phi and the J/psi meson, the yield is consistent with expectations from correlated c (c) over bar production, although other mechanisms are not ruled out. In the low-mass region, below the phi, the p+p inclusive mass spectrum is well described by known contributions from light meson decays. In contrast, the Au+Au minimum bias inclusive mass spectrum in this region shows an enhancement by a factor of 4.7 +/- 0.4(stat) +/- 1.5(syst) +/- 0.9(model). At low mass (m(ee) < 0.3 GeV/c(2)) and high p(T) (1 < p(T) < 5 GeV/c) an enhanced e(+)e(-) pair yield is observed that is consistent with production of virtual direct photons. This excess is used to infer the yield of real direct photons. In central Au+Au collisions, the excess of the direct photon yield over the p+p is exponential in p(T), with inverse slope T = 221 +/- 19(stat) +/- 19(syst) MeV. Hydrodynamical models with initial temperatures ranging from T(init) similar or equal to 300-600 MeV at times of 0.6-0.15 fm/c after the collision are in qualitative agreement with the direct photon data in Au+Au. For low p(T) < 1 GeV/c the low-mass region shows a further significant enhancement that increases with centrality and has an inverse slope of T similar or equal to 100 MeV. Theoretical models underpredict the low-mass, low-p(T) enhancement.

Measurement of Bottom Versus Charm as a Function of Transverse Momentum with Electron-Hadron Correlations in p plus p Collisions at s=200 GeV

The momentum distribution of electrons from semileptonic decays of charm and bottom quarks for midrapidity |y|< 0.35 in p+p collisions at s=200 GeV is measured by the PHENIX experiment at the Relativistic Heavy Ion Collider over the transverse momentum range 2 < p(T)< 7 GeV/c. The ratio of the yield of electrons from bottom to that from charm is presented. The ratio is determined using partial D/D -> e(+/-)K(-/+)X (K unidentified) reconstruction. It is found that the yield of electrons from bottom becomes significant above 4 GeV/c in p(T). A fixed-order-plus-next-to-leading-log perturbative quantum chromodynamics calculation agrees with the data within the theoretical and experimental uncertainties. The extracted total bottom production cross section at this energy is sigma(bb)=3.2(-1.1)(+1.2)(stat)(-1.3)(+1.4)(syst)mu b.

Event structure and double helicity asymmetry in jet production from polarized p plus p collisions at root s=200 GeV

We report on the event structure and double helicity asymmetry (A(LL)) of jet production in longitudinally polarized p + p collisions at root s = 200 GeV. Photons and charged particles were measured by the PHENIX experiment at midrapidity vertical bar eta vertical bar < 0.35 with the requirement of a high-momentum (> 2 GeV/c) photon in the event. Event structure, such as multiplicity, p(T) density and thrust in the PHENIX acceptance, were measured and compared with the results from the PYTHIA event generator and the GEANT detector simulation. The shape of jets and the underlying event were well reproduced at this collision energy. For the measurement of jet A(LL), photons and charged particles were clustered with a seed-cone algorithm to obtain the cluster pT sum (p(T)(reco)). The effect of detector response and the underlying events on p(T)(reco) was evaluated with the simulation. The production rate of reconstructed jets is satisfactorily reproduced with the next-to-leading-order and perturbative quantum chromodynamics jet production cross section. For 4< p(T)(reco) < 12 GeV/c with an average beam polarization of < P > = 49% we measured Lambda(LL) = -0.0014 +/- 0.0037(stat) at the lowest p(T)(reco) bin (4-5 GeV= c) and -0.0181 +/- 0.0282(stat) at the highest p(T)(reco) bin (10-12 GeV= c) with a beam polarization scale error of 9.4% and a pT scale error of 10%. Jets in the measured p(T)(reco) range arise primarily from hard-scattered gluons with momentum fraction 0: 02 < x < 0: 3 according to PYTHIA. The measured A(LL) is compared with predictions that assume various Delta G(x) distributions based on the Gluck-Reya-Stratmann-Vogelsang parameterization. The present result imposes the limit -a.1 < integral(0.3)(0.02) dx Delta G(x, mu(2) = GeV2) < 0.4 at 95% confidence level or integral(0.3)(0.002) dx Delta G(x, mu(2) = 1 GeV2) < 0.5 at 99% confidence level.

Direct Measurement of the Bose-Einstein Condensation Universality Class in NiCl(2)-4SC(NH(2))(2) at Ultralow Temperatures

In this work, we demonstrate field-induced Bose-Einstein condensation (BEC) in the organic compound NiCl(2)-4SC(NH(2))(2) using ac susceptibility measurements down to 1 mK. The Ni S=1 spins exhibit 3D XY antiferromagnetism between a lower critical field H(c1)similar to 2 T and a upper critical field H(c2)similar to 12 T. The results show a power-law temperature dependence of the phase transition line H(c1)(T)-H(c1)(0)=aT(alpha) with alpha=1.47 +/- 0.10 and H(c1)(0)=2.053 T, consistent with the 3D BEC universality class. Near H(c2), a kink was found in the phase boundary at approximately 150 mK.

Direct measurement of spin correlations using magnetostriction

We demonstrate that the short-range spin correlator < S(i)center dot S(j)>, a fundamental measure of the interaction between adjacent spins, can be directly measured in certain insulating magnets. We present magnetostriction data for the insulating organic compound NiCl(2)-4SC(NH(2))(2), and show that the magnetostriction as a function of field is proportional to the dominant short-range spin correlator. Furthermore, the constant of proportionality between the magnetostriction and the spin correlator gives information about the spin-lattice interaction. Combining these results with the measured Young's modulus, we are able to extract dJ/dz, the dependence of the superexchange constant J on the Ni interionic distance z.

Measurement of the Parity-Violating Longitudinal Single-Spin Asymmetry for W(+/-) Boson Production in Polarized Proton-Proton Collisions at root s=500 GeV

We report the first measurement of the parity-violating single-spin asymmetries for midrapidity decay positrons and electrons from W(+) and W(-) boson production in longitudinally polarized proton-proton collisions at root s = 500 GeV by the STAR experiment at RHIC. The measured asymmetries, A(L)(W+) = -0.27 +/- 0.10(stat.) +/- 0.02(syst.) +/- 0.03(norm.) and A(L)(W-) = 0.14 +/- 0.19(stat.) +/- 0.02(syst.) +/- 0.01(norm.), are consistent with theory predictions, which are large and of opposite sign. These predictions are based on polarized quark and antiquark distribution functions constrained by polarized deep-inelastic scattering measurements.

Measurement of the Bottom Quark Contribution to Nonphotonic Electron Production in p plus p Collisions at root s=200 GeV

The contribution of B meson decays to nonphotonic electrons, which are mainly produced by the semileptonic decays of heavy-flavor mesons, in p + p collisions at root s = 200 GeV has been measured using azimuthal correlations between nonphotonic electrons and hadrons. The extracted B decay contribution is approximately 50% at a transverse momentum of p(T) >= 5 GeV/c. These measurements constrain the nuclear modification factor for electrons from B and D meson decays. The result indicates that B meson production in heavy ion collisions is also suppressed at high p(T).

Measurement of D* mesons in jets from p plus p collisions at root s=200 GeV

We report the measurement of charged D* mesons in inclusive jets produced in proton-proton collisions at a center-of-mass energy root s = 200 GeV with the STAR experiment at the Relativistic Heavy Ion Collider. For D* mesons with fractional momenta 0.2< z< 0.5 in inclusive jets with 11.5 GeV mean transverse energy, the production rate is found to be N(D*(+) + D*(-))/N(jet) = 0.015 +/- 0.008(stat) +/- 0.007(sys). This rate is consistent with perturbative QCD evaluation of gluon splitting into a pair of charm quarks and subsequent hadronization.