Nonparametric Hierarchical Clustering of Functional Data

In this paper, we deal with the problem of curves clustering. We propose a nonparametric method which partitions the curves into clusters and discretizes the dimensions of the curve points into intervals. The cross-product of these partitions forms a data-grid which is obtained using a Bayesian model selection approach while making no assumptions regarding the curves. Finally, a post-processing technique, aiming at reducing the number of clusters in order to improve the interpretability of the clustering, is proposed. It consists in optimally merging the clusters step by step, which corresponds to an agglomerative hierarchical classification whose dissimilarity measure is the variation of the criterion. Interestingly this measure is none other than the sum of the Kullback-Leibler divergences between clusters distributions before and after the merges. The practical interest of the approach for functional data exploratory analysis is presented and compared with an alternative approach on an artificial and a real world data set

Structurama: Bayesian Inference of Population Structure

Structurama is a program for inferring population structure. Specifically, the program calculates the posterior probability of assigning individuals to different populations. The program takes as input a file containing the allelic information at some number of loci sampled from a collection of individuals. After reading a data file into computer memory, Structurama uses a Gibbs algorithm to sample assignments of individuals to populations. The program implements four different models: The number of populations can be considered fixed or a random variable with a Dirichlet process prior; moreover, the genotypes of the individuals in the analysis can be considered to come from a single population (no admixture) or as coming from several different populations (admixture). The output is a file of partitions of individuals to populations that were sampled by the Markov chain Monte Carlo algorithm. The partitions are sampled in proportion to their posterior probabilities. The program implements a number of ways to summarize the sampled partitions, including calculation of the ‘mean’ partition—a partition of the individuals to populations that minimizes the squared distance to the sampled partitions

Professionalism, Golf Coaching and a Master of Science Degree: A commentary

As a point of reference I congratulate Simon Jenkins on tackling the issue of professionalism in coaching. As he points out coaching is not a profession, but this does not mean that coaching would not benefit from going through a professionalization process. As things stand I find that the stimulus article unpacks some critically important issues of professionalism, broadly within the context of golf coaching. However, I am not sure enough is made of understanding what professional (golf) coaching actually is nor how the development of a professional golf coach can be facilitated by a Master of Science Degree (M.Sc.). I will focus my commentary on these two issues

A measurement of the tau mass and the first CPT test with tau leptons

We measure the mass of the tau lepton to be 1775.1+-1.6(stat)+-1.0(syst.) MeV using tau pairs from Z0 decays. To test CPT invariance we compare the masses of the positively and negatively charged tau leptons. The relative mass difference is found to be smaller than 3.0 10^-3 at the 90% confidence level.Comment: 10 pages, 4 figures, Submitted to Phys. Letts.

First Measurement of Z/gamma* Production in Compton Scattering of Quasi-real Photons

We report the first observation of Z/gamma* production in Compton scattering of quasi-real photons. This is a subprocess of the reaction e+e- to e+e-Z/gamma*, where one of the final state electrons is undetected. Approximately 55 pb-1 of data collected in the year 1997 at an e+e- centre-of-mass energy of 183 GeV with the OPAL detector at LEP have been analysed. The Z/gamma* from Compton scattering has been detected in the hadronic decay channel. Within well defined kinematic bounds, we measure the product of cross-section and Z/gamma* branching ratio to hadrons to be (0.9+-0.3+-0.1) pb for events with a hadronic mass larger than 60 GeV, dominated by (e)eZ production. In the hadronic mass region between 5 GeV and 60 GeV, dominated by (e)egamma* production, this product is found to be (4.1+-1.6+-0.6) pb. Our results agree with the predictions of two Monte Carlo event generators, grc4f and PYTHIA.Comment: 18 pages, LaTeX, 5 eps figures included, submitted to Physics Letters

Search for Higgs Bosons in e+e- Collisions at 183 GeV

The data collected by the OPAL experiment at sqrts=183 GeV were used to search for Higgs bosons which are predicted by the Standard Model and various extensions, such as general models with two Higgs field doublets and the Minimal Supersymmetric Standard Model (MSSM). The data correspond to an integrated luminosity of approximately 54pb-1. None of the searches for neutral and charged Higgs bosons have revealed an excess of events beyond the expected background. This negative outcome, in combination with similar results from searches at lower energies, leads to new limits for the Higgs boson masses and other model parameters. In particular, the 95% confidence level lower limit for the mass of the Standard Model Higgs boson is 88.3 GeV. Charged Higgs bosons can be excluded for masses up to 59.5 GeV. In the MSSM, mh > 70.5 GeV and mA > 72.0 GeV are obtained for tan{beta}>1, no and maximal scalar top mixing and soft SUSY-breaking masses of 1 TeV. The range 0.8 < tanb < 1.9 is excluded for minimal scalar top mixing and m{top} < 175 GeV. More general scans of the MSSM parameter space are also considered.Comment: 49 pages. LaTeX, including 33 eps figures, submitted to European Physical Journal

A Measurement of the Product Branching Ratio f(b->Lambda_b).BR(Lambda_b->Lambda X) in Z0 Decays

The product branching ratio, f(b->Lambda_b).BR(Lambda_b->Lambda X), where Lambda_b denotes any weakly-decaying b-baryon, has been measured using the OPAL detector at LEP. Lambda_b are selected by the presence of energetic Lambda particles in bottom events tagged by the presence of displaced secondary vertices. A fit to the momenta of the Lambda particles separates signal from B meson and fragmentation backgrounds. The measured product branching ratio is f(b->Lambda_b).BR(Lambda_b->Lambda X) = (2.67+-0.38(stat)+0.67-0.60(sys))% Combined with a previous OPAL measurement, one obtains f(b->Lambda_b).BR(Lambda_b->Lambda X) = (3.50+-0.32(stat)+-0.35(sys))%.Comment: 16 pages, LaTeX, 3 eps figs included, submitted to the European Physical Journal

Time-integrated luminosity recorded by the BABAR detector at the PEP-II e+e- collider

This article is the Preprint version of the final published artcile which can be accessed at the link below.We describe a measurement of the time-integrated luminosity of the data collected by the BABAR experiment at the PEP-II asymmetric-energy e+e- collider at the ϒ(4S), ϒ(3S), and ϒ(2S) resonances and in a continuum region below each resonance. We measure the time-integrated luminosity by counting e+e-→e+e- and (for the ϒ(4S) only) e+e-→μ+μ- candidate events, allowing additional photons in the final state. We use data-corrected simulation to determine the cross-sections and reconstruction efficiencies for these processes, as well as the major backgrounds. Due to the large cross-sections of e+e-→e+e- and e+e-→μ+μ-, the statistical uncertainties of the measurement are substantially smaller than the systematic uncertainties. The dominant systematic uncertainties are due to observed differences between data and simulation, as well as uncertainties on the cross-sections. For data collected on the ϒ(3S) and ϒ(2S) resonances, an additional uncertainty arises due to ϒ→e+e-X background. For data collected off the ϒ resonances, we estimate an additional uncertainty due to time dependent efficiency variations, which can affect the short off-resonance runs. The relative uncertainties on the luminosities of the on-resonance (off-resonance) samples are 0.43% (0.43%) for the ϒ(4S), 0.58% (0.72%) for the ϒ(3S), and 0.68% (0.88%) for the ϒ(2S).This work is supported by the US Department of Energy and National Science Foundation, the Natural Sciences and Engineering Research Council (Canada), the Commissariat à l’Energie Atomique and Institut National de Physique Nucléaire et de Physiquedes Particules (France), the Bundesministerium für Bildung und Forschung and Deutsche Forschungsgemeinschaft (Germany), the Istituto Nazionale di Fisica Nucleare (Italy), the Foundation for Fundamental Research on Matter (The Netherlands), the Research Council of Norway, the Ministry of Education and Science of the Russian Federation, Ministerio de Ciencia e Innovación (Spain), and the Science and Technology Facilities Council (United Kingdom). Individuals have received support from the Marie-Curie IEF program (European Union) and the A.P. Sloan Foundation (USA)

Measurement of the Michel Parameters in Leptonic Tau Decays

The Michel parameters of the leptonic tau decays are measured using the OPAL detector at LEP. The Michel parameters are extracted from the energy spectra of the charged decay leptons and from their energy-energy correlations. A new method involving a global likelihood fit of Monte Carlo generated events with complete detector simulation and background treatment has been applied to the data recorded at center-of-mass energies close to sqrt(s) = M(Z) corresponding to an integrated luminosity of 155 pb-1 during the years 1990 to 1995. If e-mu universality is assumed and inferring the tau polarization from neutral current data, the measured Michel parameters are extracted. Limits on non-standard coupling constants and on the masses of new gauge bosons are obtained. The results are in agreement with the V-A prediction of the Standard Model.Comment: 32 pages, LaTeX, 9 eps figures included, submitted to the European Physical Journal

Measurement of the B0-anti-B0-Oscillation Frequency with Inclusive Dilepton Events

The $B^0$-$\bar B^0$ oscillation frequency has been measured with a sample of 23 million \B\bar B pairs collected with the BABAR detector at the PEP-II asymmetric B Factory at SLAC. In this sample, we select events in which both B mesons decay semileptonically and use the charge of the leptons to identify the flavor of each B meson. A simultaneous fit to the decay time difference distributions for opposite- and same-sign dilepton events gives $\Delta m_d = 0.493 \pm 0.012{(stat)}\pm 0.009{(syst)}$ ps$^{-1}$.Comment: 7 pages, 1 figure, submitted to Physical Review Letter