A Fast Algorithm for Robust Regression with Penalised Trimmed Squares

The presence of groups containing high leverage outliers makes linear regression a difficult problem due to the masking effect. The available high breakdown estimators based on Least Trimmed Squares often do not succeed in detecting masked high leverage outliers in finite samples. An alternative to the LTS estimator, called Penalised Trimmed Squares (PTS) estimator, was introduced by the authors in \cite{ZiouAv:05,ZiAvPi:07} and it appears to be less sensitive to the masking problem. This estimator is defined by a Quadratic Mixed Integer Programming (QMIP) problem, where in the objective function a penalty cost for each observation is included which serves as an upper bound on the residual error for any feasible regression line. Since the PTS does not require presetting the number of outliers to delete from the data set, it has better efficiency with respect to other estimators. However, due to the high computational complexity of the resulting QMIP problem, exact solutions for moderately large regression problems is infeasible. In this paper we further establish the theoretical properties of the PTS estimator, such as high breakdown and efficiency, and propose an approximate algorithm called Fast-PTS to compute the PTS estimator for large data sets efficiently. Extensive computational experiments on sets of benchmark instances with varying degrees of outlier contamination, indicate that the proposed algorithm performs well in identifying groups of high leverage outliers in reasonable computational time.Comment: 27 page

Population Dynamics of Aedes aegypti and Dengue as Influenced by Weather and Human Behavior in San Juan, Puerto Rico

Previous studies on the influence of weather on Aedes aegypti dynamics in Puerto Rico suggested that rainfall was a significant driver of immature mosquito populations and dengue incidence, but mostly in the drier areas of the island. We conducted a longitudinal study of Ae. aegypti in two neighborhoods of the metropolitan area of San Juan city, Puerto Rico where rainfall is more uniformly distributed throughout the year. We assessed the impacts of rainfall, temperature, and human activities on the temporal dynamics of adult Ae. aegypti and oviposition. Changes in adult mosquitoes were monitored with BG-Sentinel traps and oviposition activity with CDC enhanced ovitraps. Pupal surveys were conducted during the drier and wetter parts of the year in both neighborhoods to determine the contribution of humans and rains to mosquito production. Mosquito dynamics in each neighborhood was compared with dengue incidence in their respective municipalities during the study. Our results showed that: 1. Most pupae were produced in containers managed by people, which explains the prevalence of adult mosquitoes at times when rainfall was scant; 2. Water meters were documented for the first time as productive habitats for Ae. aegypti; 3. Even though Puerto Rico has a reliable supply of tap water and an active tire recycling program, water storage containers and discarded tires were important mosquito producers; 4. Peaks in mosquito density preceded maximum dengue incidence; and 5. Ae. aegypti dynamics were driven by weather and human activity and oviposition was significantly correlated with dengue incidence

The Role of Extramembranous Cytoplasmic Termini in Assembly and Stability of the Tetrameric K+-Channel KcsA

Membrane-active alcohol 2,2,2-trifluoroethanol has been proven to be an attractive tool in the investigation of the intrinsic stability of integral membrane protein complexes by taking K+-channel KcsA as a suitable and representative ion channel. In the present study, the roles of both cytoplasmic N and C termini in channel assembly and stability of KcsA were determined. The N terminus (1–18 residues) slightly increased tetramer stability via electrostatic interactions in the presence of 30 mol.% acidic phosphatidylglycerol (PG) in phosphatidylcholine lipid bilayer. Furthermore, the N terminus was found to be potentially required for efficient channel (re)assembly. In contrast, truncation of the C terminus (125–160 residues) greatly facilitated channel reversibility from either a partially or a completely unfolded state, and this domain was substantially involved in stabilizing the tetramer in either the presence or absence of PG in lipid bilayer. These studies provide new insights into how extramembranous parts play their crucial roles in the assembly and stability of integral membrane protein complexes

Wegner on Hallucinations, Inconsistency, and the Illusion of Free Will. Some Critical Remarks

Wegner's argument on the illusory nature of conscious will, as developed in The Illusion of Conscious Will (2002) and other publications, has had major impact. Based on empirical data, he develops a theory of apparent mental causation in order to explain the occurrence of the illusion of conscious will. Part of the evidence for his argument is derived from a specific interpretation of the phenomenon of auditory verbal hallucinations as they may occur in schizophrenia. The aim of this paper is to evaluate the validity of the evidence on auditory verbal hallucinations as employed by Wegner. I conclude that auditory hallucinations do not provide solid evidence for Wegner's theory. Moreover, the phenomena in schizophrenia provide, in fact, an argument against part of Wegner's theory of apparent mental causation. © 2010 The Author(s)

Evaluation of Location-Specific Predictions by a Detailed Simulation Model of Aedes aegypti Populations

Skeeter Buster is a stochastic, spatially explicit simulation model of Aedes aegypti populations, designed to predict the outcome of vector population control methods. In this study, we apply the model to two specific locations, the cities of Iquitos, Peru, and Buenos Aires, Argentina. These two sites differ in the amount of field data that is available for location-specific customization. By comparing output from Skeeter Buster to field observations in these two cases we evaluate population dynamics predictions by Skeeter Buster with varying degrees of customization.Skeeter Buster was customized to the Iquitos location by simulating the layout of houses and the associated distribution of water-holding containers, based on extensive surveys of Ae. aegypti populations and larval habitats that have been conducted in Iquitos for over 10 years. The model is calibrated by adjusting the food input into various types of containers to match their observed pupal productivity in the field. We contrast the output of this customized model to the data collected from the natural population, comparing pupal numbers and spatial distribution of pupae in the population. Our results show that Skeeter Buster replicates specific population dynamics and spatial structure of Ae. aegypti in Iquitos. We then show how Skeeter Buster can be customized for Buenos Aires, where we only had Ae. aegypti abundance data that was averaged across all locations. In the Argentina case Skeeter Buster provides a satisfactory simulation of temporal population dynamics across seasons.This model can provide a faithful description of Ae. aegypti populations, through a process of location-specific customization that is contingent on the amount of data available from field collections. We discuss limitations presented by some specific components of the model such as the description of food dynamics and challenges that these limitations bring to model evaluation

Multiple treatment comparisons in epilepsy monotherapy trials

© 2007 Smith et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution Licens

Jet energy measurement with the ATLAS detector in proton-proton collisions at root s=7 TeV

The jet energy scale and its systematic uncertainty are determined for jets measured with the ATLAS detector at the LHC in proton-proton collision data at a centre-of-mass energy of √s = 7TeV corresponding to an integrated luminosity of 38 pb-1. Jets are reconstructed with the anti-kt algorithm with distance parameters R=0. 4 or R=0. 6. Jet energy and angle corrections are determined from Monte Carlo simulations to calibrate jets with transverse momenta pT≥20 GeV and pseudorapidities {pipe}η{pipe}<4. 5. The jet energy systematic uncertainty is estimated using the single isolated hadron response measured in situ and in test-beams, exploiting the transverse momentum balance between central and forward jets in events with dijet topologies and studying systematic variations in Monte Carlo simulations. The jet energy uncertainty is less than 2. 5 % in the central calorimeter region ({pipe}η{pipe}<0. 8) for jets with 60≤pT<800 GeV, and is maximally 14 % for pT<30 GeV in the most forward region 3. 2≤{pipe}η{pipe}<4. 5. The jet energy is validated for jet transverse momenta up to 1 TeV to the level of a few percent using several in situ techniques by comparing a well-known reference such as the recoiling photon pT, the sum of the transverse momenta of tracks associated to the jet, or a system of low-pT jets recoiling against a high-pT jet. More sophisticated jet calibration schemes are presented based on calorimeter cell energy density weighting or hadronic properties of jets, aiming for an improved jet energy resolution and a reduced flavour dependence of the jet response. The systematic uncertainty of the jet energy determined from a combination of in situ techniques is consistent with the one derived from single hadron response measurements over a wide kinematic range. The nominal corrections and uncertainties are derived for isolated jets in an inclusive sample of high-pT jets. Special cases such as event topologies with close-by jets, or selections of samples with an enhanced content of jets originating from light quarks, heavy quarks or gluons are also discussed and the corresponding uncertainties are determined. © 2013 CERN for the benefit of the ATLAS collaboration

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal

Observation of associated near-side and away-side long-range correlations in √sNN=5.02 TeV proton-lead collisions with the ATLAS detector

Two-particle correlations in relative azimuthal angle (Δϕ) and pseudorapidity (Δη) are measured in √sNN=5.02  TeV p+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using approximately 1  μb-1 of data as a function of transverse momentum (pT) and the transverse energy (ΣETPb) summed over 3.1<η<4.9 in the direction of the Pb beam. The correlation function, constructed from charged particles, exhibits a long-range (2<|Δη|<5) “near-side” (Δϕ∼0) correlation that grows rapidly with increasing ΣETPb. A long-range “away-side” (Δϕ∼π) correlation, obtained by subtracting the expected contributions from recoiling dijets and other sources estimated using events with small ΣETPb, is found to match the near-side correlation in magnitude, shape (in Δη and Δϕ) and ΣETPb dependence. The resultant Δϕ correlation is approximately symmetric about π/2, and is consistent with a dominant cos⁡2Δϕ modulation for all ΣETPb ranges and particle pT