The lack of scientists equally trained and prepared to understand both mathematics and biology/medicine hampers the development and application of computer simulation methods in biology and ...neurogastrobiology. Currently, there are no texts for navigating the extensive and intricate field of mathematical and computational modeling in neurogastrobiology. This book bridges the gap between mathematicians, computer scientists and biologists, and thus assists in the study and analysis of complex biological phenomena that cannot be done through traditional in vivo and in vitro experimental approaches.
This paper concerns robust inference on average treatment effects following model selection. Under selection on observables, we construct confidence intervals using a doubly-robust estimator that are ...robust to model selection errors and prove their uniform validity over a large class of models that allows for multivalued treatments with heterogeneous effects and selection amongst (possibly) more covariates than observations. The semiparametric efficiency bound is attained under appropriate conditions. Precise conditions are given for any model selector to yield these results, and we specifically propose the group lasso, which is apt for treatment effects, and derive new results for high-dimensional, sparse multinomial logistic regression. Both a simulation study and revisiting the National Supported Work demonstration show our estimator performs well in finite samples.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
This paper develops a new standard-error estimator for linear panel data models. The proposed estimator is robust to heteroskedasticity, serial correlation, and cross-sectional correlation of unknown ...forms. The serial correlation is controlled by the Newey–West method. To control for cross-sectional correlations, we propose to use the thresholding method, without assuming the clusters to be known. We establish the consistency of the proposed estimator. Monte Carlo simulations show the method works well. An empirical application is considered.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
The marginal distribution of count data processes rarely follows a simple Poisson model in practice. Instead, one commonly observes deviations such as overdispersion or zero inflation. To express the ...extend of such deviations from a Poisson model, one can compute an appropriately defined dispersion index or zero index. In this article, we develop several tests based on such indexes, including joint tests being based on an index combination. The asymptotic distribution of the resulting test statistics under the null hypothesis of a Poisson INAR(1) model is derived, and the finite-sample performance of the resulting tests is analyzed. Real data examples illustrate the application of these tests in practice.
Paleotopographic reconstructions of the Tibetan Plateau based on stable isotope paleoaltimetry methods conclude that most of the Plateau's current elevation was already reached by the Eocene, ~40 ...million years ago. However, changes in atmospheric and hydrological dynamics affect oxygen stable isotopes in precipitation and may thus bias such reconstructions. We used an isotope-equipped general circulation model to assess the influence of changing Eocene paleogeography and climate on paleoelevation estimates. Our simulations indicate that stable isotope paleoaltimetry methods are not applicable in Eocene Asia because of a combination of increased convective precipitation, mixture of air masses, and widespread aridity. Rather, a model-data comparison suggests that the Tibetan Plateau only reached low to moderate (less than 3000 meters) elevations during the Eocene, reconciling oxygen isotope data with other proxies.
Explosive volcanic eruptions are one of the most important driver of climate variability. Yet, we still lack a fundamental understanding of how climate change may affect future eruptions. Here, we ...use an ensemble of simulations by 1‐D and 3‐D volcanic plume models spanning a large range of eruption source and atmospheric conditions to assess changes in the dynamics of future eruptive columns. Our results shed new light on differences between the predictions of 1‐D and 3‐D plume models. Furthermore, both models suggest that as a result of ongoing climate change, for tropical eruptions, (i) higher eruption intensities will be required for plumes to reach the upper troposphere/lower stratosphere and (ii) the height of plumes currently reaching the upper troposphere/lower stratosphere or above will increase. We discuss the implications of these results for the climatic impacts of future eruptions. Our simulations can directly inform climate model experiments on climate‐volcano feedback.
Key Points
We compare the impacts of climate change on the dynamics of eruptive columns, as predicted by 1‐D and 3‐D plume models
Both models agree that higher eruption intensities will be required to inject sulfur into the tropical stratosphere
Eruptive column‐climate interactions are key to understand the climatic impacts of future eruptions
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
We describe Global Atmosphere 6.0 and Global Land 6.0 (GA6.0/GL6.0): the latest science configurations of the Met Office Unified Model and JULES (Joint UK Land Environment Simulator) land surface ...model developed for use across all timescales. Global Atmosphere 6.0 includes the ENDGame (Even Newer Dynamics for General atmospheric modelling of the environment) dynamical core, which significantly increases mid-latitude variability improving a known model bias. Alongside developments of the model's physical parametrisations, ENDGame also increases variability in the tropics, which leads to an improved representation of tropical cyclones and other tropical phenomena. Further developments of the atmospheric and land surface parametrisations improve other aspects of model performance, including the forecasting of surface weather phenomena. We also describe GA6.1/GL6.1, which includes a small number of long-standing differences from our main trunk configurations that we continue to require for operational global weather prediction. Since July 2014, GA6.1/GL6.1 has been used by the Met Office for operational global numerical weather prediction, whilst GA6.0/GL6.0 was implemented in its remaining global prediction systems over the following year.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
Extensions of linear models are very commonly used in the analysis of biological data. Whereas goodness of fit measures such as the coefficient of determination (R2) or the adjusted R2 are well ...established for linear models, it is not obvious how such measures should be defined for generalized linear and mixed models. There are by now several proposals but no consensus has yet emerged as to the best unified approach in these settings. In particular, it is an open question how to best account for heteroscedasticity and for covariance among observations present in residual error or induced by random effects. This paper proposes a new approach that addresses this issue and is universally applicable for arbitrary variance‐covariance structures including spatial models and repeated measures. It is exemplified using three biological examples.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
Multivariate count time series models are an important tool for analyzing and predicting the spread of infectious disease. We consider the endemic-epidemic framework, a class of autoregressive models ...for infectious disease surveillance counts, and replace the default autoregression on counts from the previous time period with more flexible weighting schemes inspired by discrete-time serial interval distributions. We employ three different parametric formulations, each with an additional unknown weighting parameter estimated via a profile likelihood approach, and compare them to an unrestricted nonparametric approach. The new methods are illustrated in a univariate analysis of dengue fever incidence in San Juan, Puerto Rico, and a spatiotemporal study of viral gastroenteritis in the 12 districts of Berlin. We assess the predictive performance of the suggested models and several reference models at various forecast horizons. In both applications, the performance of the endemic-epidemic models is considerably improved by the proposed weighting schemes.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP