knitr::opts_chunk$set( collapse = TRUE, comment = "#>", fig.path = "man/figures/README-", out.width = "100%" )
A fork of {cropsim} (Hijmans et al. 2009) designed to make using the EPIRICE model (Savary et al. 2012) for rice diseases easier to use.
This version provides easy to use functions to fetch weather data from NASA POWER, via the {nasapower} package (Sparks 2018, Sparks 2020) and predict disease intensity of five rice diseases using a generic SEIR model (Zadoks 1971) function, SEIR()
.
The original EPIRICE manuscript, Savary et al. (2012), which details the model and results of its use to model global epidemics of rice diseases, was published in Crop Protection detailing global unmanaged disease risk of bacterial blight, brown spot, leaf blast, sheath blight and rice tungro, which are included in this package.
You can easily simulate any of the five diseases for rice grown anywhere in the world for years from 1983 to near current using get_wth()
to fetch data from the NASA POWER web API.
Alternatively, you can supply your own weather data for any time period as long as it fits the model’s requirements.
{epicrop} is not yet on CRAN. You can install it this way.
if (!require("remotes")) install.packages("remotes") remotes::install_github("adamhsparks/epicrop", build_vignettes = TRUE )
First you need to provide weather data for the model; {epicrop} provides the get_wth()
function to do this.
Using it you can fetch weather data for any place in the world from 1983 to near present by providing the and latitude and dates or length of rice growing season as shown below.
library("epicrop") # Fetch weather for year 2000 wet season for a 120 day rice variety at the IRRI # Zeigler Experiment Station wth <- get_wth( lonlat = c(121.25562, 14.6774), dates = "2000-07-01", duration = 120 ) wth
Once you have the weather data, run the model for any of the five rice diseases by providing the emergence or crop establishment date for transplanted rice.
bb <- predict_bacterial_blight(wth, emergence = "2000-07-01") bb
Lastly, you can visualise the result of the model run.
library("ggplot2") ggplot(data = bb, aes(x = dates, y = intensity)) + labs(y = "Intensity", x = "Date") + geom_line() + geom_point() + theme_classic()
Please report any issues or bugs.
License: GPL-3
To cite {epicrop}, please use the output from citation(package = "epicrop")
.
Please note that the epicrop project is released with a Contributor Code of Conduct. By contributing to this project, you agree to abide by its terms.
Robert J. Hijmans, Serge Savary, Rene Pangga and Jorrel Aunario. (2009). Simulation modeling of crops and their diseases. R package version 0.2-6.
Serge Savary, Andrew Nelson, Laetitia Willocquet, Ireneo Pangga and Jorrel Aunario.(2012). Modeling and mapping potential epidemics of rice diseases globally. Crop Protection, Volume 34, Pages 6-17, ISSN 0261-2194 DOI: 10.1016/j.cropro.2011.11.009.
Serge Savary, Stacia Stetkiewicz, François Brun, and Laetitia Willocquet. Modelling and Mapping Potential Epidemics of Wheat Diseases-Examples on Leaf Rust and Septoria Tritici Blotch Using EPIWHEAT. (2015). European Journal of Plant Pathology 142, no. 4:771–90. DOI: 10.1007/s10658-015-0650-7.
Adam Sparks (2018). nasapower: A NASA POWER Global Meteorology, Surface Solar Energy and Climatology Data Client for R. Journal of Open Source Software, 3(30), 1035, DOI: 10.21105/joss.01035.
Adam Sparks (2021). nasapower: NASA-POWER Data from R. DOI: 10.5281/zenodo.1040727, R package version 4.0.0, URL: https://CRAN.R-project.org/package=nasapower.
Jan C. Zadoks. (1971). Systems Analysis and the Dynamics of Epidemics. Phytopathology 61:600. DOI: 10.1094/Phyto-61-600.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.