Rank-Based Analysis of Linear Models Using R

It is well-known that Wilcoxon procedures out perform least squares procedures when the data deviate from normality and/or contain outliers. These procedures can be generalized by introducing weights; yielding so-called weighted Wilcoxon (WW) techniques. In this paper we demonstrate how WW-estimates can be calculated using an L1 regression routine. More importantly, we present a collection of functions that can be used to implement a robust analysis of a linear model based on WW-estimates. For instance, estimation, tests of linear hypotheses, residual analyses, and diagnostics to detect differences in fits for various weighting schemes are discussed. We analyze a regression model, designed experiment, and autoregressive time series model for the sake of illustration. We have chosen to implement the suite of functions using the R statistical software package. Because R is freely available and runs on multiple platforms, WW-estimation and associated inference is now universally accessible.

Resolving on 100 pc scales the UV-continuum in Lyman-α\alpha emitters between redshift 2 to 3 with gravitational lensing

We present a study of seventeen LAEs at redshift 2$<z<$3 gravitationally lensed by massive early-type galaxies (ETGs) at a mean redshift of approximately 0.5. Using a fully Bayesian grid-based technique, we model the gravitational lens mass distributions with elliptical power-law profiles and reconstruct the UV-continuum surface brightness distributions of the background sources using pixellated source models. We find that the deflectors are close to, but not consistent with isothermal models in almost all cases, at the $2\sigma$-level. We take advantage of the lensing magnification (typically $\mu\simeq$ 20) to characterise the physical and morphological properties of these LAE galaxies. From reconstructing the ultra-violet continuum emission, we find that the star-formation rates range from 0.3 to 8.5 M$_{\odot}$ yr$^{-1}$ and that the galaxies are typically composed of several compact and diffuse components, separated by 0.4 to 4 kpc. Moreover, they have peak star-formation rate intensities that range from 2.1 to 54.1 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$. These galaxies tend to be extended with major axis ranging from 0.2 to 1.8 kpc (median 561 pc), and with a median ellipticity of 0.49. This morphology is consistent with disk-like structures of star-formation for more than half of the sample. However, for at least two sources, we also find off-axis components that may be associated with mergers. Resolved kinematical information will be needed to confirm the disk-like nature and possible merger scenario for the LAEs in the sample.Comment: 19 pages, 7 figures, accepted for publication on MNRA

Rank-Based Analysis of Linear Models Using R

It is well-known that Wilcoxon procedures out perform least squares procedures when the data deviate from normality and/or contain outliers. These procedures can be generalized by introducing weights; yielding so-called weighted Wilcoxon (WW) techniques. In this paper we demonstrate how WW-estimates can be calculated using an L1 regression routine. More importantly, we present a collection of functions that can be used to implement a robust analysis of a linear model based on WW-estimates. For instance, estimation, tests of linear hypotheses, residual analyses, and diagnostics to detect differences in fits for various weighting schemes are discussed. We analyze a regression model, designed experiment, and autoregressive time series model for the sake of illustration. We have chosen to implement the suite of functions using the R statistical software package. Because R is freely available and runs on multiple platforms, WW-estimation and associated inference is now universally accessible

Robust General Linear Models and Graphics via a User Interface (Web RGLM)

Rank-based procedures provide superior estimation and testing techniques when the data deviate from normality or contain gross outliers. However, these robust techniques are rarely incorporated in a nonparametric statistics or methods courses due to the lack of computational tools. One reason for this is the existence of certain unavoidable complexities in the numerical methods due to the absence of a closedform solution for the rank estimation problem. This article introduces a user interface, Web RGLM, which may be used to perform rank-based analyses of linear models across the World Wide Web. These models include simple location problems to complicated ANOVA and ANCOVA designs with multiple comparison procedures. The robust and least squares analyses are presented side-by-side for immediate comparisons. Web RGLM meets many of the computational demands of the classroom as well as the computational demands of quantitative researchers. Several illustrative examples are provided

The Cosmic Lens All-Sky Survey parent population - I. Sample selection and number counts

We present the selection of the Jodrell Bank Flat-spectrum (JBF) radio source sample, which is designed to reduce the uncertainties in the Cosmic Lens All-Sky Survey (CLASS) gravitational lensing statistics arising from the lack of knowledge about the parent population luminosity function. From observations at 4.86 GHz with the Very Large Array, we have selected a sample of 117 flat-spectrum radio sources with flux densities greater than 5 mJy. These sources were selected in a similar manner to the CLASS complete sample and are therefore representative of the parent population at low flux densities. The vast majority (~90 per cent) of the JBF sample are found to be compact on the arcsecond scales probed here and show little evidence of any extended radio jet emission. Using the JBF and CLASS complete samples we find the differential number counts slope of the parent population above and below the CLASS 30 mJy flux density limit to be -2.07+/-0.02 and -1.96+/-0.12, respectively.Comment: 10 pages, 4 figures, accepted for publication in MNRA

Time-Series Intervention Analysis Using ITSACORR: Fatal Flaws

The ITSACORR method (Crosbie, 1993, 1995) is evaluated for the analysis of two-phase interrupted time-series designs. It is shown that each component of the ITSACORR framework (including the structural model, the design matrix, the autocorrelation estimator, the ultimate parameter estimation scheme, and the inferential method) contains fatal flaws

Partial survival and inelastic collapse for a randomly accelerated particle

We present an exact derivation of the survival probability of a randomly accelerated particle subject to partial absorption at the origin. We determine the persistence exponent and the amplitude associated to the decay of the survival probability at large times. For the problem of inelastic reflection at the origin, with coefficient of restitution $r$, we give a new derivation of the condition for inelastic collapse, $r<r_c=e^{-\pi/\sqrt{3}}$, and determine the persistence exponent exactly.Comment: 6 page

Drug Residues in Food Animals

During the past several years, drug residues in food animal products have become a serious problem fr livestock producers and veterinarians. The reasons for concern are threefold: first, increased sensitivity of testing methods; second, percentage of product containing residues, and third, restrictions on potential carcinogens dictated by the Delaney Amendment. The Federal government monitors foods for residues in order to provide the American people with food that is safe and unadulterated by exogenous chemicals. The use of drugs in the livestock industry today has become widespread, both as fee additives and therapeutic agents

Low Odor, High Yield Kraft Pulping

In laboratory cooks pure oxygen was profiled into the circulation line of a batch digester during two periods of the cooking cycle: The first injection occurred during the heating steps for the purpose of in-situ generation of polysulfide. This chip treatment was studied to explore stabilization against alkaline induced carbohydrate peeling and to increase pulp yield. Under optimum conditions small amounts of polysulfide were produced with yield increase of about 0.5% These increases fell below earlier reports suggesting that unknown differences in liquor composition may influence the relative amounts of polysulfide and thiosulfate generated during the oxidation. Consequently, further studies are required to understand the factors that influence the ratios of those two sulfur species