交叉验证模式

library(spInterp)

data(TempBrazil) # Temperature for some poins of Brazil

X <- TempBrazil[, 1:2] %>% set_names(c("lon", "lat"))
Y <- TempBrazil[, 3] %>% as.matrix()  # Vector with observations in points
range <- c(-78, -34, -36, 5)
res = 1
r = kford_ml(X, Y, FUN = spInterp_adw, 
  range = range, res = res, cdd = 450)
#> 
  |                                                                            
  |                                                                      |   0%
  |                                                                            
  |==============                                                        |  20%
  |                                                                            
  |============================                                          |  40%
  |                                                                            
  |==========================================                            |  60%
  |                                                                            
  |========================================================              |  80%
  |                                                                            
  |======================================================================| 100%

print(r)
#> 
[38;5;246m# A tibble: 6 × 12
[39m
#>   kfold     R pvalue    R2   NSE   KGE  RMSE   MAE   Bias Bias_perc    AI n_sim
#>   
[3m
[38;5;246m<chr>
[39m
[23m 
[3m
[38;5;246m<dbl>
[39m
[23m  
[3m
[38;5;246m<dbl>
[39m
[23m 
[3m
[38;5;246m<dbl>
[39m
[23m 
[3m
[38;5;246m<dbl>
[39m
[23m 
[3m
[38;5;246m<dbl>
[39m
[23m 
[3m
[38;5;246m<dbl>
[39m
[23m 
[3m
[38;5;246m<dbl>
[39m
[23m  
[3m
[38;5;246m<dbl>
[39m
[23m     
[3m
[38;5;246m<dbl>
[39m
[23m 
[3m
[38;5;246m<dbl>
[39m
[23m 
[3m
[38;5;246m<int>
[39m
[23m
#> 
[38;5;250m1
[39m Fold1 0.805      0 0.649 0.608 0.79   1.25 0.934 -
[31m0
[39m
[31m.
[39m
[31m329
[39m    -
[31m0
[39m
[31m.
[39m
[31m0
[39m
[31m13
[39m 0.886    53
#> 
[38;5;250m2
[39m Fold2 0.714      0 0.51  0.434 0.704  1.29 0.998 -
[31m0
[39m
[31m.
[39m
[31m304
[39m    -
[31m0
[39m
[31m.
[39m
[31m0
[39m
[31m12
[39m 0.832    50
#> 
[38;5;250m3
[39m Fold3 0.787      0 0.619 0.602 0.599  1.45 1.1   -
[31m0
[39m
[31m.
[39m
[31m0
[39m
[31m69
[39m    -
[31m0
[39m
[31m.
[39m
[31m00
[39m
[31m3
[39m 0.841    51
#> 
[38;5;250m4
[39m Fold4 0.628      0 0.395 0.366 0.501  1.90 1.28   0.394     0.016 0.752    50
#> 
[38;5;250m5
[39m Fold5 0.834      0 0.695 0.695 0.778  1.17 0.854  0.023     0.001 0.906    52
#> 
[38;5;250m6
[39m all   0.742      0 0.55  0.546 0.671  1.43 1.03  -
[31m0
[39m
[31m.
[39m
[31m0
[39m
[31m6
[39m     -
[31m0
[39m
[31m.
[39m
[31m00
[39m
[31m2
[39m 0.846   256
#> 
#> Fold index:
#> List of 5
#>  $ Fold1: int [1:55] 1 10 16 17 24 27 30 37 43 54 ...
#>  $ Fold2: int [1:52] 3 7 13 19 21 22 23 26 35 38 ...
#>  $ Fold3: int [1:52] 2 6 8 14 34 40 41 45 49 56 ...
#>  $ Fold4: int [1:53] 4 12 15 20 25 31 33 46 51 57 ...
#>  $ Fold5: int [1:53] 5 9 11 18 28 29 32 36 39 42 ...
#> NULL