A Project Based Approach to Statistics and Data Science

In an increasingly data-driven world, facility with statistics is more important than ever for our students. At institutions without a statistician, it often falls to the mathematics faculty to teach statistics courses. This paper presents a model that a mathematician asked to teach statistics can follow. This model entails connecting with faculty from numerous departments on campus to develop a list of topics, building a repository of real-world datasets from these faculty, and creating projects where students interface with these datasets to write lab reports aimed at consumers of statistics in other disciplines. The end result is students who are well prepared for interdisciplinary research, who are accustomed to coping with the idiosyncrasies of real data, and who have sharpened their technical writing and speaking skills

Urban transportation: Perspectives on mobility and choice

A study of urban transportation systems are presented characterized by intensive scrutiny of many ideas, philosophies, and academic perspectives. This report is intended to communicate some dimensions of the urban transportation problem to the general public

GRB 050408: An Atypical Gamma-Ray Burst as a Probe of an Atypical Galactic Environment

The bright GRB 050408 was localized by HETE-II near local midnight, enabling an impressive ground-based followup effort as well as space-based followup from Swift. The Swift data from the X-Ray Telescope (XRT) and our own optical photometry and spectrum of the afterglow provide the cornerstone for our analysis. Under the traditional assumption that the visible waveband was above the peak synchrotron frequency and below the cooling frequency, the optical photometry from 0.03 to 5.03 days show an afterglow decay corresponding to an electron energy index of p_lc = 2.05 +/- 0.04, without a jet break as suggested by others. A break is seen in the X-ray data at early times (at ~12600 sec after the GRB). The spectral slope of the optical spectrum is consistent with p_lc assuming a host-galaxy extinction of A_V = 1.18 mag. The optical-NIR broadband spectrum is also consistent with p = 2.05, but prefers A_V = 0.57 mag. The X-ray afterglow shows a break at 1.26 x 10^4 sec, which may be the result of a refreshed shock. This burst stands out in that the optical and X-ray data suggest a large H I column density of N_HI ~ 10^22 cm^-2; it is very likely a damped Lyman alpha system and so the faintness of the host galaxy (M_V > -18 mag) is noteworthy. Moreover, we detect extraordinarily strong Ti II absorption lines with a column density through the GRB host that exceeds the largest values observed for the Milky Way by an order of magnitude. Furthermore, the Ti II equivalent width is in the top 1% of Mg II absorption-selected QSOs. This suggests that the large-scale environment of GRB 050408 has significantly lower Ti depletion than the Milky Way and a large velocity width (delta v > 200 km/s).Comment: ApJ submitte

Boolean delay equations on networks: An application to economic damage propagation

We introduce economic models based on Boolean Delay Equations: this formalism makes easier to take into account the complexity of the interactions between firms and is particularly appropriate for studying the propagation of an initial damage due to a catastrophe. Here we concentrate on simple cases, which allow to understand the effects of multiple concurrent production paths as well as the presence of stochasticity in the path time lengths or in the network structure. In absence of flexibility, the shortening of production of a single firm in an isolated network with multiple connections usually ends up by attaining a finite fraction of the firms or the whole economy, whereas the interactions with the outside allow a partial recovering of the activity, giving rise to periodic solutions with waves of damage which propagate across the structure. The damage propagation speed is strongly dependent upon the topology. The existence of multiple concurrent production paths does not necessarily imply a slowing down of the propagation, which can be as fast as the shortest path.Comment: Latex, 52 pages with 22 eps figure

George Wharton Rice Correspondence

Entry is typed resume with biographical information

Cosmological Constraints from calibrated Yonetoku and Amati relation implies Fundamental plane of Gamma-ray bursts

We consider two empirical relations using data only from the prompt emission of Gamma-Ray Bursts (GRBs), peak energy ($E_p$) - peak luminosity ($L_p$) relation (so called Yonetoku relation) and $E_p$-isotropic energy ($E_{\rm iso}$) relation (so called Amati relation). We first suggest the independence of the two relations although they have been considered similar and dependent. From this viewpoint, we compare constraints on cosmological parameters, $\Omega_m$ and $\Omega_{\Lambda}$, from the Yonetoku and Amati relations calibrated by low-redshift GRBs with $z < 1.8$. We found that they are different in 1-$\sigma$ level, although they are still consistent in 2-$\sigma$ level. This and the fact that both Amati and Yonetoku relations have systematic errors larger than statistical errors suggest the existence of a hidden parameter of GRBs. We introduce the luminosity time $T_L$ defined by $T_L\equiv E_{\rm iso}/L_p$ as a hidden parameter to obtain a generalized Yonetoku relation as $(L_p/{10^{52} \rm{erg s^{-1}}}) = 10^{-3.88\pm0.09}(E_p/{\rm{keV}})^{1.84\pm0.04} (T_L/{\rm{s}})^{-0.34\pm0.04}$. The new relation has much smaller systematic error, 30%, and can be regarded as "Fundamental plane" of GRBs. We show a possible radiation model for this new relation. Finally we apply the new relation for high-redshift GRBs with $1.8 < z < 5.6$ to obtain $(\Omega_m,\Omega_{\Lambda}) = (0.16^{+0.04}_{-0.06},1.20^{+0.03}_{-0.09})$, which is consistent with the concordance cosmological model within 2-$\sigma$ level.Comment: 5 pages, 6 figures, published in JCA

White-tailed deer (Odocoileus virginianus) positively affect the growth of mature northern red oak (Quercus rubra) trees

Understanding and predicting the effects of deer (Cervidae) on forest ecosystems present significant challenges in ecosystem ecology. Deer herbivory can cause large changes in the biomass and species composition of forest understory plant communities, including early life-cycle trees (i.e., seedlings and saplings). Such changes can impact juvenile to adult transitions and the future age structure and species composition of mature forests. Changes to understory vegetation also impact flow of energy and nutrients in forest ecosystems. Studies examining the influence of deer on mature trees, however, are rare and rely on extrapolating effects from early life cycle stages of trees. We tested the hypothesis that the absence of deer would result in an increase in the growth rate of mature trees by examining the impact of white-tailed deer (Odocoileus virginianus) on mature canopy trees. We examined incremental growth in individuals of Quercus rubra, an important component of temperate deciduous forests in North America, inside and outside 16-year deer exclosures in eastern U.S. deciduous forests. We found that adult trees inside exclosures grew less than those directly exposed to deer. Our findings highlight the indirect effects of white-tailed deer on the growth of adult individuals of Q. rubra in a way opposite of what would be expected from previous studies based on immature or understory tree populations. We suggest the increased growth of adult trees in the presence of deer may be explained by increased nutrient inputs through deer fecal and urine deposits and the alteration of the competitive environment belowground through the reduction of understory vegetation by browsing. Underscoring the ecological and demographic importance of adult trees in forest ecosystems, results from this study suggest the direct and indirect effects of deer on mature trees should not be overlooked

Electrical signature in polar night cloud base variations

Layer clouds are globally extensive. Their lower edges are charged negatively by the fair weather atmospheric electricity current flowing vertically through them. Using polar winter surface meteorological data from Sodankyla ̈ (Finland) and Halley (Antarctica), we find that when meteorological diurnal variations are weak, an appreciable diurnal cycle, on average, persists in the cloud base heights, detected using a laser ceilometer. The diurnal cloud base heights from both sites correlate more closely with the Carnegie curve of global atmospheric electricity than with local meteorological measurements. The cloud base sensitivities are indistinguishable between the northern and southern hemispheres, averaging a (4.0 ± 0.5) m rise for a 1% change in the fair weather electric current density. This suggests that the global fair weather current, which is affected by space weather, cosmic rays and the El Nin ̃o Southern Oscillation, is linked with layer cloud properties

Directly Measured Denitrification Reveals Oyster Aquaculture and Restored Oyster Reefs Remove Nitrogen at Comparable High Rates

Coastal systems are increasingly impacted by over-enrichment of nutrients, which has cascading effects for ecosystem functioning. Oyster restoration and aquaculture are both hypothesized to mitigate excessive nitrogen (N) loads via benthic denitrification. The degree to which these management activities perform similar functions for removing N, however, has not been extensively examined in New England, a place where nutrient runoff is high and increasing oyster (Crassostrea virginica) restoration and aquaculture activity is taking place. Here, we use a novel in situ methodology to directly measure net N2 and O2 fluxes across the sediment-water interface in a shallow (~1 m) coastal pond in southern Rhode Island. We collected data seasonally during 2013 and 2014 at restored oyster reefs, oyster aquaculture, oyster cultch (shell), and bare sediment. Restored oyster reefs and aquaculture had the highest mean (±SE) denitrification rates, 581.9 (±164.2) and 346 (±168.6) μmol N2−N m−2 h−1, respectively, and are among the highest recorded for oyster-dominated environments. Denitrification rates at sites with oyster cultch were 60.9 (±44.3) μmol N2−N m−2 h−1, which is substantially less than the sites with active oysters but still more than 50% higher than denitrification rates measured in bare sediment (24.4 ± 10.1 μmol N2–N m−2 h−1). The increase in denitrification rates at treatments, however, varied by season and the greatest rates for restored reefs were in the fall. Overall, the greatest aggregate denitrification rates occurred in the fall. Sediment oxygen demand (SOD) followed similar patterns but with greater overall rates in the summer, and displayed a strong linear relationship with denitrification (R2 = 0.93). Our results demonstrate that habitats associated with live oysters have higher net denitrification rates and that oyster reef restoration and oyster aquaculture may provide similar benefits to the ecosystem in terms of N removal. However, gas fluxes may also be affected where three-dimensional structure is introduced via oyster shell cultch and this appears to be seasonally-dependent. These data will be important for managers as they incorporate oysters into nutrient reduction strategies and consider system-level trade-offs in services provided by oyster reef restoration and aquaculture activities

Directly Measured Denitrification Reveals Oyster Aquaculture and Restored Oyster Reefs Remove Nitrogen at Comparable High Rates

Coastal systems are increasingly impacted by over-enrichment of nutrients, which has cascading effects for ecosystem functioning. Oyster restoration and aquaculture are both hypothesized to mitigate excessive nitrogen (N) loads via benthic denitrification. The degree to which these management activities perform similar functions for removing N, however, has not been extensively examined in New England, a place where nutrient runoff is high and increasing oyster (Crassostrea virginica) restoration and aquaculture activity is taking place. Here, we use a novel in situ methodology to directly measure net N2 and O2 fluxes across the sediment-water interface in a shallow (~1 m) coastal pond in southern Rhode Island. We collected data seasonally during 2013 and 2014 at restored oyster reefs, oyster aquaculture, oyster cultch (shell), and bare sediment. Restored oyster reefs and aquaculture had the highest mean (±SE) denitrification rates, 581.9 (±164.2) and 346 (±168.6) μmol N2−N m−2 h−1, respectively, and are among the highest recorded for oyster-dominated environments. Denitrification rates at sites with oyster cultch were 60.9 (±44.3) μmol N2−N m−2 h−1, which is substantially less than the sites with active oysters but still more than 50% higher than denitrification rates measured in bare sediment (24.4 ± 10.1 μmol N2–N m−2 h−1). The increase in denitrification rates at treatments, however, varied by season and the greatest rates for restored reefs were in the fall. Overall, the greatest aggregate denitrification rates occurred in the fall. Sediment oxygen demand (SOD) followed similar patterns but with greater overall rates in the summer, and displayed a strong linear relationship with denitrification (R2 = 0.93). Our results demonstrate that habitats associated with live oysters have higher net denitrification rates and that oyster reef restoration and oyster aquaculture may provide similar benefits to the ecosystem in terms of N removal. However, gas fluxes may also be affected where three-dimensional structure is introduced via oyster shell cultch and this appears to be seasonally-dependent. These data will be important for managers as they incorporate oysters into nutrient reduction strategies and consider system-level trade-offs in services provided by oyster reef restoration and aquaculture activities