Function for spatial fuzzy overlay.
Note
Independent variables in the data
provided to st_fuzzyoverlay()
must be discretized variables,
and dependent variable are continuous variable.
References
Yongze Song & Peng Wu (2021) An interactive detector for spatial associations, International Journal of Geographical Information Science, 35:8, 1676-1701, DOI:10.1080/13658816.2021.1882680
Author
Wenbo Lv lyu.geosocial@gmail.com
Examples
data('sim')
sim = sim %>%
dplyr::mutate(dplyr::across(4:6,\(.x) st_unidisc(.x,4,"quantile")))
fo1 = st_fuzzyoverlay(y~xa+xb+xc,data = sim, method = 'and')
fo2 = st_fuzzyoverlay(y~xa+xb+xc,data = sim, method = 'or')
fo1
#> 1 2 3 4 5 6 7 8 9 10 11
#> "xa_1" "xa_3" "xa_1" "xc_1" "xb_2" "xa_1" "xa_3" "xc_1" "xa_2" "xb_1" "xb_2"
#> 12 13 14 15 16 17 18 19 20 21 22
#> "xa_2" "xa_1" "xc_1" "xc_1" "xb_1" "xb_2" "xb_1" "xa_2" "xb_1" "xa_3" "xb_2"
#> 23 24 25 26 27 28 29 30 31 32 33
#> "xb_2" "xa_2" "xb_2" "xa_2" "xc_1" "xa_3" "xc_1" "xc_1" "xa_3" "xa_2" "xa_3"
#> 34 35 36 37 38 39 40 41 42 43 44
#> "xa_2" "xa_3" "xa_2" "xa_2" "xb_2" "xc_1" "xc_1" "xa_1" "xa_3" "xa_1" "xa_2"
#> 45 46 47 48 49 50 51 52 53 54 55
#> "xa_2" "xb_2" "xa_2" "xc_1" "xc_1" "xc_1" "xa_1" "xb_2" "xa_3" "xb_2" "xa_2"
#> 56 57 58 59 60 61 62 63 64 65 66
#> "xc_1" "xc_1" "xc_1" "xc_1" "xc_1" "xa_2" "xa_3" "xa_2" "xa_2" "xb_1" "xb_1"
#> 67 68 69 70 71 72 73 74 75 76 77
#> "xc_1" "xc_1" "xc_1" "xc_1" "xa_2" "xc_1" "xa_3" "xc_1" "xb_1" "xa_2" "xc_1"
#> 78 79 80
#> "xc_1" "xc_1" "xc_1"
fo2
#> 1 2 3 4 5 6 7 8 9 10 11
#> "xb_3" "xc_3" "xb_3" "xb_2" "xc_3" "xb_3" "xc_3" "xa_1" "xc_2" "xc_3" "xc_3"
#> 12 13 14 15 16 17 18 19 20 21 22
#> "xc_3" "xb_3" "xb_2" "xa_1" "xa_3" "xc_3" "xa_3" "xc_2" "xa_3" "xc_3" "xc_3"
#> 23 24 25 26 27 28 29 30 31 32 33
#> "xc_3" "xc_3" "xc_3" "xc_3" "xa_1" "xc_3" "xb_2" "xa_2" "xc_3" "xc_3" "xc_3"
#> 34 35 36 37 38 39 40 41 42 43 44
#> "xc_2" "xc_3" "xb_3" "xc_2" "xc_3" "xa_2" "xa_2" "xb_3" "xc_3" "xb_3" "xb_3"
#> 45 46 47 48 49 50 51 52 53 54 55
#> "xc_2" "xc_3" "xc_2" "xa_2" "xa_2" "xa_1" "xb_3" "xc_3" "xc_3" "xc_3" "xc_2"
#> 56 57 58 59 60 61 62 63 64 65 66
#> "xa_1" "xb_2" "xb_2" "xa_1" "xa_1" "xb_3" "xc_3" "xb_3" "xc_3" "xa_3" "xc_2"
#> 67 68 69 70 71 72 73 74 75 76 77
#> "xa_1" "xa_1" "xa_1" "xa_1" "xb_3" "xb_2" "xc_3" "xb_2" "xa_3" "xc_2" "xa_2"
#> 78 79 80
#> "xa_1" "xa_1" "xa_1"