Helper function that combines objects of class cv_ammi
,
cv_ammif
or cv_blup
. It is useful when looking for a boxplot
containing the RMSPD values of those cross-validation procedures.
Arguments
- ...
Input objects of class
cv_ammi
,cv_ammif
orcv_blup
.- bind
What data should be used? To plot the RMSPD, use 'boot' (default). Use
bind = 'means'
to return the RMSPD mean for each model.- sort
Used to sort the RMSPD mean in ascending order.
Value
An object of class cv_ammif
. The results will depend on the
argument bind
. If bind = 'boot'
then the RMSPD of all models
in ...
will be bind to a unique data frame. If bind = 'means'
then the RMSPD mean of all models in ...
will be bind to an unique
data frame.
Author
Tiago Olivoto tiagoolivoto@gmail.com
Examples
# \donttest{
library(metan)
# Two examples with only 5 resampling procedures
AMMI <- cv_ammi(data_ge,
resp = GY,
gen = GEN,
env = ENV,
rep = REP,
nboot = 5)
#> Validating 1 of 5 sets |======== | 20% 00:00:01
Validating 2 of 5 sets |================ | 40% 00:00:03
Validating 3 of 5 sets |========================= | 60% 00:00:04
Validating 4 of 5 sets |================================= | 80% 00:00:06
Validating 5 of 5 sets |=========================================| 100% 00:00:07
BLUP <- cv_blup(data_ge,
resp = GY,
gen = GEN,
env = ENV,
rep = REP,
nboot = 5)
#> Validating 1 of 5 sets |======== | 20% 00:00:01
Validating 2 of 5 sets |================ | 40% 00:00:02
Validating 3 of 5 sets |========================= | 60% 00:00:04
Validating 4 of 5 sets |================================= | 80% 00:00:05
Validating 5 of 5 sets |=========================================| 100% 00:00:06
bind_data <- bind_cv(AMMI, BLUP)
plot(bind_data)
print(bind_cv(AMMI, BLUP, bind = 'means'))
#> $RMSPD
#> MODEL mean sd se Q2.5 Q97.5
#> 1 BLUP_g_RCBD 0.4117647 0.018842733 0.008426726 0.3835109 0.4289504
#> 2 AMMI2 0.4300084 0.009145797 0.004090125 0.4181219 0.4415320
#>
#> attr(,"class")
#> [1] "cvalidation"
# }