Universality Class of One-Dimensional Directed Sandpile Models

A general n-state directed `sandpile' model is introduced. The stationary properties of the n-state model are derived for n < infty, and analytical arguments based on a central limit theorem show that the model belongs to the universality class of the totally asymmetric Oslo model, with a crossover to uncorrelated branching process behavior for small system sizes. Hence, the central limit theorem allows us to identify the existence of a large universality class of one-dimensional directed sandpile models.Comment: 4 pages, 2 figure

Critical behavior of the three-dimensional bond-diluted Ising spin glass: Finite-size scaling functions and Universality

We study the three-dimensional (3D) bond-diluted Edwards-Anderson (EA) model with binary interactions at a bond occupation of 45% by Monte Carlo (MC) simulations. Using an efficient cluster MC algorithm we are able to determine the universal finite-size scaling (FSS) functions and the critical exponents with high statistical accuracy. We observe small corrections to scaling for the measured observables. The critical quantities and the FSS functions indicate clearly that the bond-diluted model for dilutions above the critical dilution p*, at which a spin glass (SG) phase appears, lies in the same universality class as the 3D undiluted EA model with binary interactions. A comparison with the FSS functions of the 3D site-diluted EA model with Gaussian interactions at a site occupation of 62.5% gives very strong evidence for the universality of the SG transition in the 3D EA model.Comment: Revised version. 10 pages, 9 figures, 2 table

An extended class of minimax generalized Bayes estimators of regression coefficients

We derive minimax generalized Bayes estimators of regression coefficients in the general linear model with spherically symmetric errors under invariant quadratic loss for the case of unknown scale. The class of estimators generalizes the class considered in Maruyama and Strawderman (2005) to include non-monotone shrinkage functions

Slip-velocity of large neutrally-buoyant particles in turbulent flows

We discuss possible definitions for a stochastic slip velocity that describes the relative motion between large particles and a turbulent flow. This definition is necessary because the slip velocity used in the standard drag model fails when particle size falls within the inertial subrange of ambient turbulence. We propose two definitions, selected in part due to their simplicity: they do not require filtration of the fluid phase velocity field, nor do they require the construction of conditional averages on particle locations. A key benefit of this simplicity is that the stochastic slip velocity proposed here can be calculated equally well for laboratory, field, and numerical experiments. The stochastic slip velocity allows the definition of a Reynolds number that should indicate whether large particles in turbulent flow behave (a) as passive tracers; (b) as a linear filter of the velocity field; or (c) as a nonlinear filter to the velocity field. We calculate the value of stochastic slip for ellipsoidal and spherical particles (the size of the Taylor microscale) measured in laboratory homogeneous isotropic turbulence. The resulting Reynolds number is significantly higher than 1 for both particle shapes, and velocity statistics show that particle motion is a complex non-linear function of the fluid velocity. We further investigate the nonlinear relationship by comparing the probability distribution of fluctuating velocities for particle and fluid phases

The Generic, Incommensurate Transition in the two-dimensional Boson Hubbard Model

The generic transition in the boson Hubbard model, occurring at an incommensurate chemical potential, is studied in the link-current representation using the recently developed directed geometrical worm algorithm. We find clear evidence for a multi-peak structure in the energy distribution for finite lattices, usually indicative of a first order phase transition. However, this multi-peak structure is shown to disappear in the thermodynamic limit revealing that the true phase transition is second order. These findings cast doubts over the conclusion drawn in a number of previous works considering the relevance of disorder at this transition.Comment: 13 pages, 10 figure

On the future of astrostatistics: statistical foundations and statistical practice

This paper summarizes a presentation for a panel discussion on "The Future of Astrostatistics" held at the Statistical Challenges in Modern Astronomy V conference at Pennsylvania State University in June 2011. I argue that the emerging needs of astrostatistics may both motivate and benefit from fundamental developments in statistics. I highlight some recent work within statistics on fundamental topics relevant to astrostatistical practice, including the Bayesian/frequentist debate (and ideas for a synthesis), multilevel models, and multiple testing. As an important direction for future work in statistics, I emphasize that astronomers need a statistical framework that explicitly supports unfolding chains of discovery, with acquisition, cataloging, and modeling of data not seen as isolated tasks, but rather as parts of an ongoing, integrated sequence of analyses, with information and uncertainty propagating forward and backward through the chain. A prototypical example is surveying of astronomical populations, where source detection, demographic modeling, and the design of survey instruments and strategies all interact.Comment: 8 pp, 2 figures. To appear in "Statistical Challenges in Modern Astronomy V," (Lecture Notes in Statistics, Vol. 209), ed. Eric D. Feigelson and G. Jogesh Babu; publication planned for Sep 2012; see http://www.springer.com/statistics/book/978-1-4614-3519-

Majority Dynamics and Aggregation of Information in Social Networks

Consider n individuals who, by popular vote, choose among q >= 2 alternatives, one of which is "better" than the others. Assume that each individual votes independently at random, and that the probability of voting for the better alternative is larger than the probability of voting for any other. It follows from the law of large numbers that a plurality vote among the n individuals would result in the correct outcome, with probability approaching one exponentially quickly as n tends to infinity. Our interest in this paper is in a variant of the process above where, after forming their initial opinions, the voters update their decisions based on some interaction with their neighbors in a social network. Our main example is "majority dynamics", in which each voter adopts the most popular opinion among its friends. The interaction repeats for some number of rounds and is then followed by a population-wide plurality vote. The question we tackle is that of "efficient aggregation of information": in which cases is the better alternative chosen with probability approaching one as n tends to infinity? Conversely, for which sequences of growing graphs does aggregation fail, so that the wrong alternative gets chosen with probability bounded away from zero? We construct a family of examples in which interaction prevents efficient aggregation of information, and give a condition on the social network which ensures that aggregation occurs. For the case of majority dynamics we also investigate the question of unanimity in the limit. In particular, if the voters' social network is an expander graph, we show that if the initial population is sufficiently biased towards a particular alternative then that alternative will eventually become the unanimous preference of the entire population.Comment: 22 page

Design strategies for optimizing holographic optical tweezers setups

We provide a detailed account of the construction of a system of holographic optical tweezers. While much information is available on the design, alignment and calibration of other optical trapping configurations, those based on holography are relatively poorly described. Inclusion of a spatial light modulator in the setup gives rise to particular design trade-offs and constraints, and the system benefits from specific optimization strategies, which we discuss.Comment: 16 pages, 15 figure

One-Dimensional Directed Sandpile Models and the Area under a Brownian Curve

We derive the steady state properties of a general directed ``sandpile'' model in one dimension. Using a central limit theorem for dependent random variables we find the precise conditions for the model to belong to the universality class of the Totally Asymmetric Oslo model, thereby identifying a large universality class of directed sandpiles. We map the avalanche size to the area under a Brownian curve with an absorbing boundary at the origin, motivating us to solve this Brownian curve problem. Thus, we are able to determine the moment generating function for the avalanche-size probability in this universality class, explicitly calculating amplitudes of the leading order terms.Comment: 24 pages, 5 figure

Incidence of keratitis of varying severity among contact lens wearers

AIM: To determine the incidence of non-severe keratitis (NSK) and severe keratitis (SK) among wearers of current generation contact lenses. METHODS: A 12 month, prospective, hospital based epidemiological study was conducted by examining all contact lens wearers presenting with a corneal infiltrate/ulcer to a hospital centre in Manchester. A clinical severity matrix was used to differentiate between NSK and SK, based on the severity of signs and symptoms. The size of the hospital catchment population and the wearing modalities (daily wear (DW) or extended wear (EW)) and lens types being used were estimated from relevant demographic and market data. RESULTS: During the survey period, 80 and 38 patients presented with NSK and SK, respectively. The annual incidences (cases per 10,000 wearers) for each wearing modality and lens type were: DW rigid--NSK 5.7, SK 2.9; DW hydrogel daily disposable--NSK 9.1, SK 4.9; DW hydrogel (excluding daily disposable)--NSK 14.1, SK 6.4; DW silicone hydrogel--NSK 55.9, SK 0.0; EW rigid--NSK 0.0, SK 0.0; EW hydrogel--NSK 48.2, SK 96.4; EW silicone hydrogel--NSK 98.8, SK 19.8. The difference in SK between EW hydrogel and EW silicone hydrogel was significant (p = 0.04). CONCLUSIONS: A clinical severity matrix has considerable utility in assessing contact lens related keratitis. There is a significantly higher incidence of SK in wearers who sleep in contact lenses compared with those who only use lenses during the waking hours. Those who choose to sleep in lenses should be advised to wear silicone hydrogel lenses, which carry a five times decreased risk of SK for extended wear compared with hydrogel lenses