Adaptive trimmed likelihood estimation in regression

In this paper we derive an asymptotic normality result for an adaptive trimmed likelihood estimator of regression starting from initial high breakdownpoint robust regression estimates. The approach leads to quickly and easily computed robust and efficient estimates for regression. A highlight of the method is that it tends automatically in one algorithm to expose the outliers and give least squares estimates with the outliers removed. The idea is to begin with a rapidly computed consistent robust estimator such as the least median of squares (LMS) or least trimmed squares (LTS) or for example the more recent MM estimators of Yohai. Such estimators are now standard in statistics computing packages, for example as in SPLUS or R. In addition to the asymptotics we provide data analyses supporting the new adaptive approach. This approach appears to work well on a number of data sets and is quicker than the related brute force adaptive regression approach described in Clarke (2000). This current approach builds on the work of Bednarski and Clarke (2002) which considered the asymptotics for the location estimator only

Crime Script Sequencing: An optimal forensic combination for cold case analysis

Criminal cases go cold when investigative leads or forensic testing does not lead to a successful arrest. In these cases, investigators are often keen to use novel methods to derive fresh ideas or insights. Recently, academics from a range of fields, including Psychology, Criminology, and Forensic Sciences have developed a range of new methods and tests to assist with police investigations. The current paper outlines a novel approach to assisting with police cold case investigations: Crime Script Sequencing. The new method combines two leading temporal methods, Crime Script Analysis and Behaviour Sequence Analysis. A real-world cold case, the bombing of Canadian Pacific Airlines Fight 21, is presented and analysed using Crime Script Sequencing to offer readers a guide of how to use the method for other investigations. Impacts, insights, and potential future developments of the method are outlined

Relationship of milk expected progeny differences (EPD\u27s)to milk production and calf weaning weight

Spring-calving Angus (n=37) and Simmental (n=28) cows were used to evaluate the relationship between milk expected progeny difference (EPD), milk production, total milk energy, and weaning weight. There were positive correlations between milk EPD\u27s and milk produced during lactation, total milk energy per lactation, and calf weaning weight. On average, a 1 lb increase in milk EPD predicted a 1.8 lb increase in calf weaning weight and an increase of 56.6 (Angus) and 70.2 (Simmental) lb in total milk per lactation. Our results suggest that producers can use milk EPD\u27s as selection tools to influence the milk production of their cows and subsequent weaning weights of their calves

Determination of step--edge barriers to interlayer transport from surface morphology during the initial stages of homoepitaxial growth

We use analytic formulae obtained from a simple model of crystal growth by molecular--beam epitaxy to determine step--edge barriers to interlayer transport. The method is based on information about the surface morphology at the onset of nucleation on top of first--layer islands in the submonolayer coverage regime of homoepitaxial growth. The formulae are tested using kinetic Monte Carlo simulations of a solid--on--solid model and applied to estimate step--edge barriers from scanning--tunneling microscopy data on initial stages of Fe(001), Pt(111), and Ag(111) homoepitaxy.Comment: 4 pages, a Postscript file, uuencoded and compressed. Physical Review B, Rapid Communications, in press

Debris disk size distributions: steady state collisional evolution with P-R drag and other loss processes

We present a new scheme for determining the shape of the size distribution, and its evolution, for collisional cascades of planetesimals undergoing destructive collisions and loss processes like Poynting-Robertson drag. The scheme treats the steady state portion of the cascade by equating mass loss and gain in each size bin; the smallest particles are expected to reach steady state on their collision timescale, while larger particles retain their primordial distribution. For collision-dominated disks, steady state means that mass loss rates in logarithmic size bins are independent of size. This prescription reproduces the expected two phase size distribution, with ripples above the blow-out size, and above the transition to gravity-dominated planetesimal strength. The scheme also reproduces the expected evolution of disk mass, and of dust mass, but is computationally much faster than evolving distributions forward in time. For low-mass disks, P-R drag causes a turnover at small sizes to a size distribution that is set by the redistribution function (the mass distribution of fragments produced in collisions). Thus information about the redistribution function may be recovered by measuring the size distribution of particles undergoing loss by P-R drag, such as that traced by particles accreted onto Earth. Although cross-sectional area drops with 1/age^2 in the PR-dominated regime, dust mass falls as 1/age^2.8, underlining the importance of understanding which particle sizes contribute to an observation when considering how disk detectability evolves. Other loss processes are readily incorporated; we also discuss generalised power law loss rates, dynamical depletion, realistic radiation forces and stellar wind drag.Comment: Accepted for publication by Celestial Mechanics and Dynamical Astronomy (special issue on EXOPLANETS

Longitudinal double-spin asymmetry and cross section for inclusive neutral pion production at midrapidity in polarized proton collisions at sqrt(s) = 200 GeV

We report a measurement of the longitudinal double-spin asymmetry A_LL and the differential cross section for inclusive Pi0 production at midrapidity in polarized proton collisions at sqrt(s) = 200 GeV. The cross section was measured over a transverse momentum range of 1 < p_T < 17 GeV/c and found to be in good agreement with a next-to-leading order perturbative QCD calculation. The longitudinal double-spin asymmetry was measured in the range of 3.7 < p_T < 11 GeV/c and excludes a maximal positive gluon polarization in the proton. The mean transverse momentum fraction of Pi0's in their parent jets was found to be around 0.7 for electromagnetically triggered events.Comment: 6 pages, 3 figures, submitted to Phys. Rev. D (RC

Longitudinal scaling property of the charge balance function in Au + Au collisions at 200 GeV

We present measurements of the charge balance function, from the charged particles, for diverse pseudorapidity and transverse momentum ranges in Au + Au collisions at 200 GeV using the STAR detector at RHIC. We observe that the balance function is boost-invariant within the pseudorapidity coverage [-1.3, 1.3]. The balance function properly scaled by the width of the observed pseudorapidity window does not depend on the position or size of the pseudorapidity window. This scaling property also holds for particles in different transverse momentum ranges. In addition, we find that the width of the balance function decreases monotonically with increasing transverse momentum for all centrality classes.Comment: 6 pages, 3 figure

Energy and system size dependence of \phi meson production in Cu+Cu and Au+Au collisions

We study the beam-energy and system-size dependence of \phi meson production (using the hadronic decay mode \phi -- K+K-) by comparing the new results from Cu+Cu collisions and previously reported Au+Au collisions at \sqrt{s_NN} = 62.4 and 200 GeV measured in the STAR experiment at RHIC. Data presented are from mid-rapidity (|y|<0.5) for 0.4 < pT < 5 GeV/c. At a given beam energy, the transverse momentum distributions for \phi mesons are observed to be similar in yield and shape for Cu+Cu and Au+Au colliding systems with similar average numbers of participating nucleons. The \phi meson yields in nucleus-nucleus collisions, normalised by the average number of participating nucleons, are found to be enhanced relative to those from p+p collisions with a different trend compared to strange baryons. The enhancement for \phi mesons is observed to be higher at \sqrt{s_NN} = 200 GeV compared to 62.4 GeV. These observations for the produced \phi(s\bar{s}) mesons clearly suggest that, at these collision energies, the source of enhancement of strange hadrons is related to the formation of a dense partonic medium in high energy nucleus-nucleus collisions and cannot be alone due to canonical suppression of their production in smaller systems.Comment: 20 pages and 5 figure

Measurement of the Bottom contribution to non-photonic electron production in p+pp+p collisions at s\sqrt{s} =200 GeV

The contribution of $B$ meson decays to non-photonic electrons, which are mainly produced by the semi-leptonic decays of heavy flavor mesons, in $p+p$ collisions at $\sqrt{s} =$ 200 GeV has been measured using azimuthal correlations between non-photonic electrons and hadrons. The extracted $B$ decay contribution is approximately 50% at a transverse momentum of $p_{T} \geq 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}$.Comment: 6 pages, 4 figures, accepted by PR