DyMEP_prediction_visualizer {DyMEP}

DyMEP_prediction_visualizer

Description

Visualizes the predictions of the DyMEP model.

Usage

DyMEP_prediction_visualizer(detailed_output)

Arguments

detailed_output

Output of the pheno_phase_prediction function with output_type = "detailed_information".

Value

A plot with one panel per phenology phase, showing the environmental covariate responses, the GLM prediction, and the phase prediction (points).

Fields

timestamp

Description of timestamp column.

Examples

phase_covariate_list <- best_DyMEP_model(env_covariates =
             c("tas","tasmin","VPD","SPI","global_radiation","tasmax","RH"),
           pheno_phases = c("sowing-emergence","emergence-jointing",
           "jointing-heading"),
                                        crop_abbrev = "WW",
                                        output_list_for_prediction = TRUE)

# Create dummy environmental data
environmental_data <- data.frame("DATE" =
                                       seq.Date(from = as.Date("2021-01-01"),
                                         to = as.Date("2023-12-31"), by = 1),
                      "tas" = runif(1095, min = -10, max = 40),
                       "RH" = runif(1095, min = 0, max = 100),
                       "tasmin" = runif(1095, min = -10, max = 40),
                       "tasmax" = runif(1095, min = 0, max = 40),
                       "VPD" = runif(1095, min = 0, max = 40),
                       "SPI" = runif(1095, min = -1, max = 4),
                       "global_radiation" = runif(1095, min = 0, max = 3500))

DyMEP_prediction_visualizer(detailed_output = pheno_phase_prediction(
                                 phase_covariate_list = phase_covariate_list,
                                 environmental_data = environmental_data,
                                 phase_starting_date = as.Date("2021-01-01"),
                                 crop_abbrev = "WW",
                                 output_type = "detailed_information"))