Variable Selection and Model Averaging in Semiparametric Overdispersed Generalized Linear Models

We express the mean and variance terms in a double exponential regression model as additive functions of the predictors and use Bayesian variable selection to determine which predictors enter the model, and whether they enter linearly or flexibly. When the variance term is null we obtain a generalized additive model, which becomes a generalized linear model if the predictors enter the mean linearly. The model is estimated using Markov chain Monte Carlo simulation and the methodology is illustrated using real and simulated data sets.Comment: 8 graphs 35 page

Stellar Hydrodynamics in Radiative Regions

We present an analysis of the response of a radiative region to waves generated by a convective region of the star; this wave treatment of the classical problem of ``overshooting'' gives extra mixing relative to the treatment traditionally used in stellar evolutionary codes. The interface between convectively stable and unstable regions is dynamic and nonspherical, so that the nonturbulent material is driven into motion, even in the absence of ``penetrative overshoot.'' These motions may be described by the theory of nonspherical stellar pulsations, and are related to motion measured by helioseismology. Multi-dimensional numerical simulations of convective flow show puzzling features which we explain by this simplified physical model. Gravity waves generated at the interface are dissipated, resulting in slow circulation and mixing seen outside the formal convection zone. The approach may be extended to deal with rotation and composition gradients. Tests of this description in the stellar evolution code TYCHO produce carbon stars on the asymptotic giant branch (AGB), an isochrone age for the Hyades and three young clusters with lithium depletion ages from brown dwarfs, and lithium and beryllium depletion consistent with observations of the Hyades and Pleiades, all without tuning parameters. The insight into the different contributions of rotational and hydrodynamic mixing processes could have important implications for realistic simulation of supernovae and other questions in stellar evolution.Comment: 27 pages, 5 figures, accepted to the Astrophysical Journa

Observations from Preliminary Experiments on Spatial and Temporal Pressure Measurements from Near-Field Free Air Explosions

It is self-evident that a crucial step in analysing the performance of protective structures is to be able to accurately quantify the blast load arising from a high explosive detonation. For structures located near to the source of a high explosive detonation, the resulting pressure is extremely high in magnitude and highly non-uniform over the face of the target. There exists very little direct measurement of blast parameters in the nearfield, mainly attributed to the lack of instrumentation sufficiently robust to survive extreme loading events yet sensitive enough to capture salient features of the blast. Instead literature guidance is informed largely by early numerical analyses and parametric studies. Furthermore, the lack of an accurate, reliable data set has prevented subsequent numerical analyses from being validated against experimental trials. This paper presents an experimental methodology that has been developed in part to enable such experimental data to be gathered. The experimental apparatus comprises an array of Hopkinson pressure bars, fitted through holes in a target, with the loaded faces of the bars flush with the target face. Thus, the bars are exposed to the normally or obliquely reflected shocks from the impingement of the blast wave with the target. Pressure-time recordings are presented along with associated Arbitary-Langrangian-Eulerian modelling using the LS-DYNA explicit numerical code. Experimental results are corrected for the effects of dispersion of the propagating waves in the pressure bars, enabling accurate characterisation of the peak pressures and impulses from these loadings. The combined results are used to make comments on the mechanism of the pressure load for very near-field blast events

LBT and Spitzer Spectroscopy of Star-Forming Galaxies at 1 < z < 3: Extinction and Star Formation Rate Indicators

We present spectroscopic observations in the rest-frame optical and near- to mid-infrared wavelengths of four gravitationally lensed infrared (IR) luminous star-forming galaxies at redshift 1 < z < 3 from the LUCIFER instrument on the Large Binocular Telescope and the Infrared Spectrograph on Spitzer. The sample was selected to represent pure, actively star-forming systems, absent of active galactic nuclei. The large lensing magnifications result in high signal-to-noise spectra that can probe faint IR recombination lines, including Pa-alpha and Br-alpha at high redshifts. The sample was augmented by three lensed galaxies with similar suites of unpublished data and observations from the literature, resulting in the final sample of seven galaxies. We use the IR recombination lines in conjunction with H-alpha observations to probe the extinction, Av, of these systems, as well as testing star formation rate (SFR) indicators against the SFR measured by fitting spectral energy distributions to far-IR photometry. Our galaxies occupy a range of Av from ~0 to 5.9 mag, larger than previously known for a similar range of IR luminosities at these redshifts. Thus, estimates of SFR even at z ~ 2 must take careful count of extinction in the most IR luminous galaxies. We also measure extinction by comparing SFR estimates from optical emission lines with those from far-IR measurements. The comparison of results from these two independent methods indicates a large variety of dust distribution scenarios at 1 < z < 3. Without correcting for dust extinction, the H-alpha SFR indicator underestimates the SFR; the size of the necessary correction depends on the IR luminosity and dust distribution scenario. Individual SFR estimates based on the 6.2 micron PAH emission line luminosity do not show a systematic discrepancy with extinction, although a considerable, ~0.2 dex scatter is observed.Comment: Accepted for publication in The Astrophysical Journal; 14 pages, 8 figure

Large Binocular Telescope and Sptizer Spectroscopy of Star-forming Galaxies at 1 < z < 3: Extinction and Star Formation Rate Indicators

We present spectroscopic observations in the rest-frame optical and near- to mid-infrared wavelengths of four gravitationally lensed infrared (IR) luminous star-forming galaxies at redshift 1 < z < 3 from the LUCIFER instrument on the Large Binocular Telescope and the Infrared Spectrograph on Spitzer. The sample was selected to represent pure, actively star-forming systems, absent of active galactic nuclei. The large lensing magnifications result in high signal-to-noise spectra that can probe faint IR recombination lines, including Pa and Br at high redshifts. The sample was augmented by three lensed galaxies with similar suites of unpublished data and observations from the literature, resulting in the final sample of seven galaxies. We use the IR recombination lines in conjunction with H observations to probe the extinction, Av, of these systems, as well as testing star formation rate (SFR) indicators against the SFR measured by fitting spectral energy distributions to far-IR photometry. Our galaxies occupy a range of Av from 0 to 5.9 mag, larger than previously known for a similar range of IR luminosities at these redshifts. Thus, estimates of SFR even at z 2 must take careful count of extinction in the most IR luminous galaxies.We also measure extinction by comparing SFR estimates from optical emission lines with those from far- IR measurements. The comparison of results from these two independent methods indicates a large variety of dust distribution scenarios at 1 < z < 3. Without correcting for dust extinction, the H SFR indicator underestimates the SFR; the size of the necessary correction depends on the IR luminosity and dust distribution scenario. Individual SFR estimates based on the 6.2m polycyclic aromatic hydrocarbon emission line luminosity do not show a systematic discrepancy with extinction, although a considerable, 0.2 dex, scatter is observed

Identification of novel subgroup a variants with enhanced receptor binding and replicative capacity in primary isolates of anaemogenic strains of feline leukaemia virus

&lt;b&gt;BACKGROUND:&lt;/b&gt; The development of anaemia in feline leukaemia virus (FeLV)-infected cats is associated with the emergence of a novel viral subgroup, FeLV-C. FeLV-C arises from the subgroup that is transmitted, FeLV-A, through alterations in the amino acid sequence of the receptor binding domain (RBD) of the envelope glycoprotein that result in a shift in the receptor usage and the cell tropism of the virus. The factors that influence the transition from subgroup A to subgroup C remain unclear, one possibility is that a selective pressure in the host drives the acquisition of mutations in the RBD, creating A/C intermediates with enhanced abilities to interact with the FeLV-C receptor, FLVCR. In order to understand further the emergence of FeLV-C in the infected cat, we examined primary isolates of FeLV-C for evidence of FeLV-A variants that bore mutations consistent with a gradual evolution from FeLV-A to FeLV-C.&lt;p&gt;&lt;/p&gt; &lt;b&gt;RESULTS:&lt;/b&gt; Within each isolate of FeLV-C, we identified variants that were ostensibly subgroup A by nucleic acid sequence comparisons, but which bore mutations in the RBD. One such mutation, N91D, was present in multiple isolates and when engineered into a molecular clone of the prototypic FeLV-A (Glasgow-1), enhanced replication was noted in feline cells. Expression of the N91D Env on murine leukaemia virus (MLV) pseudotypes enhanced viral entry mediated by the FeLV-A receptor THTR1 while soluble FeLV-A Env bearing the N91D mutation bound more efficiently to mouse or guinea pig cells bearing the FeLV-A and -C receptors. Long-term in vitro culture of variants bearing the N91D substitution in the presence of anti-FeLV gp70 antibodies did not result in the emergence of FeLV-C variants, suggesting that additional selective pressures in the infected cat may drive the subsequent evolution from subgroup A to subgroup C.&lt;p&gt;&lt;/p&gt; &lt;b&gt;CONCLUSIONS:&lt;/b&gt; Our data support a model in which variants of FeLV-A, bearing subtle differences in the RBD of Env, may be predisposed towards enhanced replication in vivo and subsequent conversion to FeLV-C. The selection pressures in vivo that drive the emergence of FeLV-C in a proportion of infected cats remain to be established

Far Infrared Source Counts at 70 and 160 microns in Spitzer Deep Surveys

We derive galaxy source counts at 70 and 160 microns using the Multiband Imaging Photometer for Spitzer (MIPS) to map the Chandra Deep Field South (CDFS) and other fields. At 70 microns, our observations extend upwards about 2 orders of magnitude in flux density from a threshold of 15 mJy, and at 160 microns they extend about an order of magnitude upward from 50 mJy. The counts are consistent with previous observations on the bright end. Significant evolution is detected at the faint end of the counts in both bands, by factors of 2-3 over no-evolution models. This evolution agrees well with models that indicate most ofthe faint galaxies lie at redshifts between 0.7 and 0.9. The new Spitzer data already resolve about 23% of the Cosmic Far Infrared Background at 70 microns and about 7% at 160 microns.Comment: Small modifications to match printed version. Models in Differential Counts plots were changed. MIPS Source Counts are available at: http://lully.as.arizona.edu/GTODeep/Counts/ . Accepted for Publication in ApJS Special Issue on Spitze