Creates a tidy data.table from the output of trace_asco().
tidy_trace(trace)a nested list output from trace_asco()
A tidy data.table of trace_asco() output.
if (FALSE) { # interactive()
Newmarracarra <-
read.csv(system.file("extdata",
"1998_Newmarracarra_weather_table.csv", package = "ascotraceR"))
station_data <-
system.file("extdata", "stat_dat.csv", package = "ascotraceR")
weather_dat <- format_weather(
x = Newmarracarra,
POSIXct_time = "Local.Time",
temp = "mean_daily_temp",
ws = "ws",
wd_sd = "wd_sd",
rain = "rain_mm",
wd = "wd",
station = "Location",
time_zone = "Australia/Perth",
lonlat_file = station_data)
traced <- trace_asco(
weather = weather_dat,
paddock_length = 20,
paddock_width = 20,
initial_infection = "1998-06-10",
sowing_date = as.POSIXct("1998-06-09"),
harvest_date = as.POSIXct("1998-06-09") + lubridate::ddays(100),
time_zone = "Australia/Perth",
primary_infection_foci = "centre")
tidied <- tidy_trace(traced)
# take a look at the infectious growing points on day 102
library(ggplot2)
ggplot(data = subset(tidied, i_day == 102),
aes(x = x, y = y, fill = infectious_gp)) +
geom_tile()
}