Carly J. Stevens
Lancaster Environment Centre
Lancaster University
Lancaster
LA1 4YQ, UK
E-mail: [email protected]

Eric M. Lind
Department of Ecology, Evolution, and Behavior
University of Minnesota
St. Paul
MN 55108 USA
E-mail: [email protected]

Yann Hautier
Department of Ecology, Evolution, and Behavior
University of Minnesota
St. Paul
MN 55108 USA
and
Department of Plant Sciences
University of Oxford
South Parks Road
Oxford
OX1 3RB, UK

W. Stanley Harpole
Department of Ecology, Evolution and Organismal Biology
Iowa State University
Ames
IA 50011 USA
and
German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig
Deutscher Platz 5e
D-04103 Leipzig
Germany
and
Department of Physiological Diversity
Helmholtz Center for Environmental Research – UFZ
Permoserstr. 15
04318 Leipzig
Germany
and
Institute of Biology, Martin Luther University Halle-Wittenberg
Am Kirchtor 1
06108 Halle (Saale)
Germany

Elizabeth T. Borer
Department of Ecology, Evolution, and Behavior
University of Minnesota
St. Paul
MN 55108 USA

Sarah Hobbie
Department of Ecology, Evolution, and Behavior
University of Minnesota
St. Paul
MN 55108 USA

Eric W. Seabloom
Department of Ecology, Evolution, and Behavior
University of Minnesota
St. Paul
MN 55108 USA

Laura Ladwig
Department of Biology
University of New Mexico
Albuquerque
NM 87131 USA

Jonathan D. Bakker
School of Environmental and Forest Sciences
University of Washington
Box 354115
Seattle
WA 98195-4115, USA

Chengjin Chu
State Key Laboratory of Grassland and Agro-Ecosystems
Research Station of Alpine Meadow and Wetland Ecosystems
School of Life Sciences
Lanzhou University
Lanzhou 730000
China

Scott Collins
Department of Biology
University of New Mexico
Albuquerque
NM 87131 USA

Kendi F. Davies
Department of Ecology and Evolutionary Biology
University of Colorado
Boulder
CO 80309, USA

Jennifer Firn
School of Earth, Environmental and Biological Sciences
Queensland University of Technology
Brisbane Qld, 4001
Australia

Helmut Hillebrand
Institute for Chemistry and Biology of the Marine Environment
Carl-von Ossietzky University
Oldenburg
Germany

Kimberly J. La Pierre
Department of Ecology and Evolutionary Biology
Yale University
New Haven
CT 06511, USA

Andrew MacDougall
Department of Integrative Biology
University of Guelph
Guelph
Ontario
N1G 2W1, Canada

Brett Melbourne
Department of Ecology and Evolutionary Biology
University of Colorado
Boulder
CO 80309, USA

Rebecca L. McCulley
Department of Plant & Soil Sciences
University of Kentucky
Lexington
KY 40546-0091, USA

John Morgan
Department of Botany
La Trobe University
Bundoora 3086
Victoria
Australia

John L. Orrock
Department of Zoology
University of Wisconsin
Madison
WI 53704, USA

Suzanne M. Prober
CSIRO Land and Water Flagship and the Great Western Woodlands Supersite
Private Bag 5
Wembley WA 6913
Australia

Anita C. Risch
Swiss Federal Institute for Forest
Snow and Landscape Research
Birmensdorf 8903
Switzerland

Martin Schuetz
Swiss Federal Institute for Forest
Snow and Landscape Research
Birmensdorf 8903
Switzerland

Peter D. Wragg
Department of Ecology, Evolution, and Behavior
University of Minnesota
St. Paul
MN 55108 USA

R_code.txt (MD5: b098d58a19eb2994276fc868c14660e7)

R script used to implement the hierarchical model used in analysis of the relationship between grassland primary productivity and potential driver variables in the paper ‘Anthropogenic nitrogen deposition predicts local grassland primary production worldwide’. A full description of the hierarchical model is given in Appendix A. The hierarchical model predicting plot-level aboveground biomass was constructed by combining plot-level and site-level predictors as follows. The natural log of the live biomass harvested at peak standing crop (an estimate of aboveground annual productivity) was the response variable. The following code was run in R v3.0.2 and JAGS 3.1.0 using R2Jags 0.03-11.