
Changelog
epicrop 1.0.9000
Major Changes
Version now 1.0.0.9000, almost ready for first major release
SEIR()
is nowseir()
All helper functions drop the
predict_
prefix for simplicity, e.g.,predict_bacterial_blight()
becomes justbacterial_blight()
Addition of functions to predict rice blast and sheath blight based on Kim et al. 2015.
Adds the RcRH parameter to
seir()
to calculate the relative humidity correction factor for the disease progress curve as published in Savary et al. 2012.Adds the RcRH parameter to all helper functions.
New Functionality
-
Adds functions from Kim et al. 2015 for:
- Leaf blast in South Korean varieties
- Sheath blight in South Korean varieties
Adds functionality missing in original {cropsim} to handle complex leaf wetness calculations using daily RH to calculate hourly RH and then modify the value of
RcW
using.fn_Rc()
as already existing withRcT
andRcA
by setting the argument,simple_wetness == FALSE
inseir()
Adds functionality to use a matrix in
RRG
and modify it with.fn_Rc
over time as withRcA
andRcT
as found in Kim et al. 2015 (though not implemented in Kim et al. 2015)Allow pass through of more arguments from helper functions to
seir()
via...
argumentStrict checking of user-inputs with informative messages now provided if inputs are invalid
epicrop 0.0.3
Minor changes
- Extend the RcA curve for bacterial blight, the original version of EPIRICE had
...0.43, 0.41, 0.42, 0.41
, previous versions of {epicrop} had...0.43, 0.41, 0.41, 0.41
, now we have...0.43, 0.42, 0.42, 0.41...
.
Bug fixes
- Corrected documentation for
get_wth()
, {nasapower} is only client used.
epicrop 0.0.2
epicrop 0.0.1
First stable release on GitHub.
Fixes function title formatting.
Replaces outdated {tidyverse} functions found in vignettes with modern versions.
epicrop 0.0.0.9017
Update README R output.
Polish map output in the mapping vignette.
Polish names for figure in multiple locations figure.
Remove references and functionality called from {epifitter}; AUDPC is calculated internally now.
epicrop 0.0.0.9015
epicrop 0.0.0.9013
Bug fixes
Fixes a bug introduced into SEIR() with the last round of bug fixes. The calculations for latent and infectious sites were still calculating for the next, not current, day in the
for()
loop.Corrects values in the example for
SEIR()
to run the model for brown spot.
Minor changes
Updates input values for
Rc
curves in each of the EPIRICE models such that the values used to compute the curve are all supplied, none are computed before the interpolation is performed. This leads to longer, but clearer lines of code.All of the help files now reference NASA POWER and CHIRPS both as sources of possible weather data for use in the model.
epicrop 0.0.0.9012
Bug fixes
- Fixes bug in
SEIR()
whereRc
value was calculated incorrectly that was introduced with commit https://codeberg.org/r-forge/cropsim/commit/ba093b252deb222f1ffe54b3a7f0d14b6192f18c#diff-17fd120d02f173b87db7e09b563e2057b95391bef823eee29cbff30cb3d80287. See Lines 124:126 that were changed to Lines 118:120 where
Rc[day+1] <- baseRc * AFGen(ageRc, day) * AFGen(tmpRc, wthsub$tavg[day+1]) * RHCoef[day+1]
should have been
Rc[day+1] <- baseRc * AFGen(ageRc, day+1) * AFGen(tmpRc, wthsub$tavg[day+1]) * RHCoef[day+1]
to calculate the basic infection rate, Rc
, for day + 1
. Commit: 248df5512b5897e950dd58a0a8b3fd57922add47.
- Fixes typos in
RcA
values forpredict_bacterial_blight()
andpredict_sheath_blight()
.- In
predict_bacterial_blight()
the eighth entry was improperly encoded as 0.42, it has been corrected to 0.41, commit 8b2c36b61d4f949385dcac4c03ff435b8400e335. - In
predict_sheath_blight()
the seventh entry was improperly encoded as 0.83, it has been corrected to 0.84, Commit: 69b89ce24dd41c16d4627225c028355f76a64b80.
- In
Major changes
Implements a new function to calculate the correction values for crop age,
RcA
, and temperature,RcT
, usingapprox()
in a custom function that is faster than theselect_mod_value()
function.Implements new functionality for
get_wth()
to optionally get data from the CHIRPS/CHIRTS APIs in addition to the default NASA POWER API.Uses one-indexed
for()
loop, in place of the zero-indexedfor()
loop found in {cropsim}.
epicrop 0.0.0.9011
Minor changes
Fix title in vignette entry.
Remove
lazyData
field from DESCRIPTION.Remove redundant
maintainer
field from DESCRIPTION.Remove dates (date when knit) from vignettes.
Fix/clean up note on parallel processing in “Multiples” vignette.
Better, more consistent documentation.
Update
get_wth()
to work with the new (unreleased) version of {nasapower} that works with the new version of the POWER API.
epicrop 0.0.0.9010
Major changes
BREAKING CHANGE
intensity
column, inSEIR()
is now a proportion, i.e. values range between 0 - 1. If you wish to have the same values as previously reported,df$intensity*100
will give you these values.BREAKING CHANGE Removes column,
severity
, inSEIR()
output.Adds a new vignette illustrating how to run multiple simulations and use parallel processing to reduce run time for these simulations.
epicrop 0.0.0.9009
epicrop 0.0.0.9007
Minor changes
Only include Lat/Lon values if the
wth
object provides them.Better documentation about the Lat/Lon values in the output of
SEIR()
Parameters for
predict_bacterial_blight()
are used as basis for all of the otherpredict_*()
family of functions.Add DOI to CITATION and DESCRIPTION files.
Add Lifecycle and DOI badges to README files.
epicrop 0.0.0.9005
Bug fixes
- Convert
wth
(weather input object) to adata.table
type object internally if it is not already one. Prior, if the object was not adata.table
,SEIR()
would fail with a message that the dates did not align. This should fix that issue and anydata.frame
type object including atibble
can be provided now. Thanks to Jean Fabrice Adanve for helping me find this bug.
epicrop 0.0.0.9004
Bug fixes
Fixes bug where the relative humidity checks in
SEIR()
only checked if the daily RH value was equal to (==
) not equal to or greater than (>=
) the set parameter forrhlim
(default is 90%).Example for
SEIR()
in {roxygen} section now works properly when executed by the user.
Major changes
Any default parameter values are moved from
SEIR()
to thepredict_()
functions themselves, so any calls directly toSEIR()
must specify all parameters.Deleted
inst/alt_versions/tungrov2.R
file.Deleted
inst/workflows/spatsim.R
file.
Minor changes
Edit documentation for better clarity
Better commenting in code for self and others’ reference
-
SEIR()
is simplified and further optimised- In some cases in the
for()
loop, the value ofday + 1
was repeatedly calculated and assigned to a new object. This has been corrected with a single object,d1
being created at the beginning of each loop instance. - In some cases where RH, TEMP or RAIN are repeatedly checked against, vectors of these three values are created outside the loop to save time when checking or extracting values.
- Redundant code and other unused objects are cleaned up and removed.
- In some cases in the
Standardise library loading calls to use standard evaluation in documentation.
Standardise italics to use “_” rather than “*”.
Use
ggplot2::theme_classic()
for example figures in README and {epicrop} vignette.