@@ -3,7 +3,7 @@ title: "Synergistic-Unique-Redundant Decomposition of Causality (SURD)"
33author : " Wenbo Lyu"
44date : |
55 | Last update: 2026-05-05
6- | Last run: 2026-05-05
6+ | Last run: 2026-05-26
77output : rmarkdown::html_vignette
88vignette : >
99 %\VignetteIndexEntry{surd}
@@ -152,26 +152,26 @@ tibble::as_tibble(res_cvds)
152152# # # A tibble: 20 × 3
153153# # vars types values
154154# # <chr> <chr> <dbl>
155- # # 1 V1 U 0.00219
156- # # 2 V2 U 0.0155
157- # # 3 V3 U 0.00252
158- # # 4 V1_V2 R 0.00345
159- # # 5 V1_V4 R 0.000711
160- # # 6 V1_V2_V3 R 0.00506
161- # # 7 V1_V2_V4 R 0.00991
162- # # 8 V1_V2_V3_V4 R 0.0269
163- # # 9 V1_V2 S 0.0205
164- # # 10 V1_V3 S 0.0224
165- # # 11 V1_V4 S 0.0451
166- # # 12 V2_V3 S 0.0256
167- # # 13 V2_V4 S 0.0352
168- # # 14 V3_V4 S 0.0366
169- # # 15 V1_V2_V3 S 0.0350
170- # # 16 V1_V2_V4 S 0.104
155+ # # 1 V1 U 0.00220
156+ # # 2 V2 U 0.0156
157+ # # 3 V3 U 0.00253
158+ # # 4 V1_V2 R 0.00347
159+ # # 5 V1_V4 R 0.000714
160+ # # 6 V1_V2_V3 R 0.00508
161+ # # 7 V1_V2_V4 R 0.00995
162+ # # 8 V1_V2_V3_V4 R 0.0232
163+ # # 9 V1_V2 S 0.0206
164+ # # 10 V1_V3 S 0.0225
165+ # # 11 V1_V4 S 0.0453
166+ # # 12 V2_V3 S 0.0257
167+ # # 13 V2_V4 S 0.0353
168+ # # 14 V3_V4 S 0.0367
169+ # # 15 V1_V2_V3 S 0.0352
170+ # # 16 V1_V2_V4 S 0.105
171171# # 17 V1_V3_V4 S 0.151
172- # # 18 V2_V3_V4 S 0.0467
173- # # 19 V1_V2_V3_V4 S 0.404
174- # # 20 InfoLeak InfoLeak 0.535
172+ # # 18 V2_V3_V4 S 0.0468
173+ # # 19 V1_V2_V3_V4 S 0.406
174+ # # 20 InfoLeak InfoLeak 0.530
175175```
176176
177177The SURD results are shown in the figure below:
@@ -188,8 +188,8 @@ utils_plot_surd(res_cvds)
188188
189189``` r
190190popd_nb = spdep :: read.gal(system.file(" case/popd_nb.gal" ,package = " spEDM" ))
191- # # Warning in spdep::read.gal(system.file("case/popd_nb.gal", package = "spEDM")): neighbour object has 4
192- # # sub-graphs
191+ # # Warning in spdep::read.gal(system.file("case/popd_nb.gal", package =
192+ # # "spEDM")): neighbour object has 4 sub-graphs
193193popd = readr :: read_csv(system.file(" case/popd.csv" ,package = " spEDM" ))
194194popd_sf = sf :: st_as_sf(popd , coords = c(" lon" ," lat" ), crs = 4326 )
195195popd_sf
@@ -256,15 +256,15 @@ utils_plot_surd(res_popd)
256256``` r
257257npp = terra :: rast(system.file(" case/npp.tif" , package = " spEDM" ))
258258npp
259- # # class : SpatRaster
259+ # # class : SpatRaster
260260# # size : 404, 483, 5 (nrow, ncol, nlyr)
261261# # resolution : 10000, 10000 (x, y)
262262# # extent : -2625763, 2204237, 1877078, 5917078 (xmin, xmax, ymin, ymax)
263- # # coord. ref. : CGCS2000_Albers
264- # # source : npp.tif
265- # # names : npp, pre, tem, elev, hfp
266- # # min values : 164.00 , 384.3409 , -47.8194 , -122.2004 , 0.03390418
267- # # max values : 16606.33 , 23878.3555 , 263.6938 , 5350.4902 , 44.90312195
263+ # # coord. ref. : CGCS2000_Albers
264+ # # source : npp.tif
265+ # # names : npp, pre, tem, elev, hfp
266+ # # min values : 164, 384.340942 , -47.819405 , -122.200386 , 0.033904
267+ # # max values : 16606.333984 , 23878.355469 , 263.693787 , 5350.490234 , 44.903122
268268```
269269
270270
0 commit comments