Libraries

library(regionaldrivers)
library(magrittr)

Data loading

SSCT_data and SSCT_labelled_data are loaded with the package. Here is what 10 random rows of each looks like:

SSCT_labelled_data
CatAreaSqKm WsAreaSqKm CHYD CCHEM CSED CCONN CTEMP CHABT ICI WHYD WCHEM WSED WCONN WTEMP WHABT IWI PctCarbResidCat PctNonCarbResidCat PctAlkIntruVolCat PctSilicicCat PctExtruVolCat PctColluvSedCat PctGlacTilClayCat PctGlacTilLoamCat PctGlacTilCrsCat PctGlacLakeCrsCat PctGlacLakeFineCat PctHydricCat PctEolCrsCat PctEolFineCat PctSalLakeCat PctAlluvCoastCat PctCoastCrsCat PctWaterCat PctCarbResidWs PctNonCarbResidWs PctAlkIntruVolWs PctSilicicWs PctExtruVolWs PctColluvSedWs PctGlacTilClayWs PctGlacTilLoamWs PctGlacTilCrsWs PctGlacLakeCrsWs PctGlacLakeFineWs PctHydricWs PctEolCrsWs PctEolFineWs PctSalLakeWs PctAlluvCoastWs PctCoastCrsWs PctWaterWs MineDensCat MineDensWs MineDensCatRp100 MineDensWsRp100 PctOw2011Cat PctIce2011Cat PctUrbOp2011Cat PctUrbLo2011Cat PctUrbMd2011Cat PctUrbHi2011Cat PctBl2011Cat PctDecid2011Cat PctConif2011Cat PctMxFst2011Cat PctShrb2011Cat PctGrs2011Cat PctHay2011Cat PctCrop2011Cat PctWdWet2011Cat PctHbWet2011Cat PctOw2011Ws PctIce2011Ws PctUrbOp2011Ws PctUrbLo2011Ws PctUrbMd2011Ws PctUrbHi2011Ws PctBl2011Ws PctDecid2011Ws PctConif2011Ws PctMxFst2011Ws PctShrb2011Ws PctGrs2011Ws PctHay2011Ws PctCrop2011Ws PctWdWet2011Ws PctHbWet2011Ws PctOw2011CatRp100 PctIce2011CatRp100 PctUrbOp2011CatRp100 PctUrbLo2011CatRp100 PctUrbMd2011CatRp100 PctUrbHi2011CatRp100 PctBl2011CatRp100 PctDecid2011CatRp100 PctConif2011CatRp100 PctMxFst2011CatRp100 PctShrb2011CatRp100 PctGrs2011CatRp100 PctHay2011CatRp100 PctCrop2011CatRp100 PctWdWet2011CatRp100 PctHbWet2011CatRp100 PctOw2011WsRp100 PctIce2011WsRp100 PctUrbOp2011WsRp100 PctUrbLo2011WsRp100 PctUrbMd2011WsRp100 PctUrbHi2011WsRp100 PctBl2011WsRp100 PctDecid2011WsRp100 PctConif2011WsRp100 PctMxFst2011WsRp100 PctShrb2011WsRp100 PctGrs2011WsRp100 PctHay2011WsRp100 PctCrop2011WsRp100 PctWdWet2011WsRp100 PctHbWet2011WsRp100 Precip8110Cat Tmax8110Cat Tmean8110Cat Tmin8110Cat Precip8110Ws Tmax8110Ws Tmean8110Ws Tmin8110Ws RunoffCat RunoffWs ClayCat SandCat ClayWs SandWs OmCat PermCat RckDepCat WtDepCat OmWs PermWs RckDepWs WtDepWs SLOPE CONFINEMEN RUSLE SO LDD aspect_max.rstr aspect_mean.rstr aspect_median.rstr aspect_min.rstr aspect_sd.rstr aspect_skew.rstr curvplan_max.rstr curvplan_mean.rstr curvplan_median.rstr curvplan_min.rstr curvplan_sd.rstr curvplan_skew.rstr curvprof_max.rstr curvprof_mean.rstr curvprof_median.rstr curvprof_min.rstr curvprof_sd.rstr curvprof_skew.rstr flowdir_max.rstr flowdir_mean.rstr flowdir_median.rstr flowdir_min.rstr flowdir_sd.rstr flowdir_skew.rstr layer_max.rstr layer_mean.rstr layer_median.rstr layer_min.rstr layer_sd.rstr layer_skew.rstr roughness_max.rstr roughness_mean.rstr roughness_median.rstr roughness_min.rstr roughness_sd.rstr roughness_skew.rstr slope_max.rstr slope_mean.rstr slope_median.rstr slope_min.rstr slope_sd.rstr slope_skew.rstr tpi_max.rstr tpi_mean.rstr tpi_median.rstr tpi_min.rstr tpi_sd.rstr tpi_skew.rstr tri_max.rstr tri_mean.rstr tri_median.rstr tri_min.rstr tri_sd.rstr tri_skew.rstr aspect_max.nrch aspect_mean.nrch aspect_median.nrch aspect_min.nrch aspect_sd.nrch aspect_skew.nrch curvplan_max.nrch curvplan_mean.nrch curvplan_median.nrch curvplan_min.nrch curvplan_sd.nrch curvplan_skew.nrch curvprof_max.nrch curvprof_mean.nrch curvprof_median.nrch curvprof_min.nrch curvprof_sd.nrch curvprof_skew.nrch flowdir_max.nrch flowdir_mean.nrch flowdir_median.nrch flowdir_min.nrch flowdir_sd.nrch flowdir_skew.nrch layer_max.nrch layer_mean.nrch layer_median.nrch layer_min.nrch layer_sd.nrch layer_skew.nrch roughness_max.nrch roughness_mean.nrch roughness_median.nrch roughness_min.nrch roughness_sd.nrch roughness_skew.nrch slope_max.nrch slope_mean.nrch slope_median.nrch slope_min.nrch slope_sd.nrch slope_skew.nrch tpi_max.nrch tpi_mean.nrch tpi_median.nrch tpi_min.nrch tpi_sd.nrch tpi_skew.nrch tri_max.nrch tri_mean.nrch tri_median.nrch tri_min.nrch tri_sd.nrch tri_skew.nrch H.640 H.960 H.1280 H.1600 H.1920 H.2240 H.2560 H.2880 H.3200 H.3840 H.4480 H.5120 H.5760 H.6400 H.7680 H.8960 H.10240 H.11520 H.12800 H.15360 H.17920 H.20480 H.23040 H.25600 H.30720 H.35840 H.40960 H.46080 H.51200 H.61440 H.71680 H.81920
201241 4.892 361.055 1.000 0.995 0.993 1.000 1.000 0.990 0.977 0.995 0.994 0.991 0.992 0.994 0.987 0.954 0 0 0 0.000 0.000 0 0 0 0 0 0 0 0 0 0 100.000 0 0 26.998 34.981 0 32.728 0.000 0 0 0 0.000 0 0 0 0 0 0 5.293 0 0.000 0 0 0 0 0.000 0 0.629 0 0.0 0 0.000 0.000 0.143 9.871 66.094 19.800 0.200 0.000 3.262 0.000 0.047 0 2.390 0.020 0.006 0.000 0.197 5.649 22.570 6.016 41.264 21.493 0.189 0.060 0.099 0.001 0 0 1.838 7.353 0 0 18.015 1.103 0.000 0.000 0.000 29.044 1.103 8.824 32.721 0.000 0.136 0.012 2.369 0.378 0.194 0.041 1.199 1.449 45.596 1.434 22.462 13.855 4.535 4.568 0.312 1.459 789.413 23.998 16.835 9.666 1049.550 20.324 13.615 6.902 521.00 280.846 33.016 26.206 25.453 32.828 0.580 1.173 139.832 181.074 1.206 5.379 70.159 182.221 0.005 34 1.748 4 0.729 6.281 2.568 2.215 0.001 1.759 0.354 0.023 0.000 0.000 -0.025 0.003 -0.420 0.022 -0.002 0.000 -0.046 0.007 -1.559 128 8.426 1 1 11.871 4.874 1588.414 1520.509 1517.787 1484.161 24.691 0.410 17.716 8.032 7.956 0.034 3.273 -0.040 0.855 0.348 0.354 0.000 0.148 -0.046 1.411 -0.041 0.000 -1.834 0.309 -1.145 5.441 2.341 2.320 0.007 0.955 0.004 5.117 4.508 4.892 2.695 0.635 -0.692 0.003 0.000 0.000 -0.007 0.003 -0.299 -0.004 -0.013 -0.012 -0.026 0.007 -0.626 58.667 17.333 16.000 11.000 9.387 NA 1493.268 1490.205 1490.015 1488.065 1.497 0.352 7.944 5.379 5.392 2.426 1.450 -0.256 0.391 0.238 0.242 0.066 0.087 -0.288 -0.034 -0.411 -0.366 -0.896 0.253 -0.435 2.548 1.640 1.605 0.792 0.465 0.033 0.879 0.914 0.932 0.882 0.896 0.862 0.823 0.865 0.831 0.837 0.771 0.804 0.760 0.728 0.656 0.722 0.636 0.689 0.635 0.648 0.607 0.583 0.574 0.549 0.559 0.574 0.576 0.588 0.567 0.486 0.497 0.533
106305 3.163 454.925 0.992 0.968 0.968 0.988 0.984 0.957 0.865 0.977 0.967 0.970 0.971 0.966 0.963 0.827 0 0 0 100.000 0.000 0 0 0 0 0 0 0 0 0 0 0.000 0 0 0.000 0.000 0 100.000 0.000 0 0 0 0.000 0 0 0 0 0 0 0.000 0 0.000 0 0 0 0 0.000 0 0.000 0 0.0 0 0.000 0.000 92.583 0.000 7.278 0.139 0.000 0.000 0.000 0.000 0.000 0 0.000 0.000 0.000 0.000 0.000 0.000 92.583 0.000 7.278 0.139 0.000 0.000 0.000 0.000 0 0 0.000 0.000 0 0 0.000 0.000 1.429 0.000 78.571 20.000 0.000 0.000 0.000 0.000 0.000 0.000 0.104 0.007 0.000 0.000 1.116 0.000 36.124 0.000 52.597 5.869 0.000 0.000 0.000 4.183 1170.733 14.109 7.702 1.290 1170.733 14.109 7.702 1.290 292.41 292.410 14.656 55.949 14.656 55.949 0.628 9.536 123.600 182.880 0.628 9.536 123.600 182.880 0.016 57 5.519 5 0.874 6.281 3.052 3.025 0.000 2.073 -0.019 0.065 0.000 0.001 -0.064 0.010 -0.834 0.035 -0.002 -0.002 -0.043 0.010 -0.891 128 10.033 1 1 15.447 4.946 432.050 307.860 297.054 250.146 47.468 0.463 28.007 12.160 12.481 0.596 6.164 0.026 1.248 0.504 0.518 0.003 0.259 0.030 4.300 -0.093 -0.056 -4.079 0.729 -0.569 9.912 3.794 3.759 0.180 1.966 0.157 3.156 1.804 1.394 0.746 1.025 0.169 0.006 0.000 0.000 -0.008 0.004 -0.417 -0.004 -0.013 -0.013 -0.025 0.005 -0.165 85.667 8.781 1.000 1.000 21.476 NA 258.434 253.570 253.797 250.720 2.096 0.446 12.658 6.901 6.907 1.550 3.486 -0.004 0.544 0.293 0.278 0.060 0.155 0.057 -0.117 -0.666 -0.693 -1.332 0.352 -0.350 4.212 2.251 2.138 0.473 1.178 0.108 0.971 0.919 0.912 0.898 0.849 0.792 0.854 0.820 0.776 0.739 0.664 0.633 0.698 0.608 0.621 0.561 0.533 0.591 0.556 0.516 0.524 0.406 0.487 0.395 0.443 0.362 0.405 0.398 0.449 0.404 0.368 0.362
467423 0.031 11.259 1.000 0.990 0.994 1.000 1.000 0.992 0.976 1.000 0.988 0.997 1.000 0.999 0.997 0.982 0 0 0 100.000 0.000 0 0 0 0 0 0 0 0 0 0 0.000 0 0 0.000 0.000 0 100.000 0.000 0 0 0 0.000 0 0 0 0 0 0 0.000 0 0.000 0 0 0 0 0.000 0 0.000 0 0.0 0 0.000 0.000 0.128 0.000 30.393 69.479 0.000 0.000 0.000 0.000 0.000 0 0.000 0.000 0.000 0.000 0.000 0.000 0.128 0.000 30.393 69.479 0.000 0.000 0.000 0.000 0 0 0.000 0.000 0 0 0.000 0.000 20.000 0.000 80.000 0.000 0.000 0.000 0.000 0.000 0.019 0.000 2.072 0.337 0.125 0.046 0.112 0.046 68.305 0.055 20.554 4.839 1.230 0.132 0.056 2.072 1029.869 21.757 15.874 9.987 1029.869 21.757 15.874 9.987 848.00 848.000 20.360 35.030 20.360 35.030 1.390 6.730 48.750 182.880 1.390 6.730 48.750 182.880 0.135 0 1.392 2 0.962 6.281 3.643 3.748 0.001 1.991 -0.313 0.040 0.002 0.002 -0.113 0.009 -2.688 0.061 -0.003 -0.001 -0.103 0.014 -1.915 64 10.776 16 1 7.624 0.058 400.136 288.385 279.215 207.182 46.196 0.422 34.737 15.579 15.341 1.956 5.152 0.218 1.413 0.653 0.648 0.040 0.215 0.024 3.611 -0.055 0.047 -4.784 0.926 -0.944 11.291 4.911 4.812 0.590 1.600 0.284 5.419 3.579 3.974 1.554 1.691 -0.165 0.007 -0.006 -0.004 -0.027 0.012 -0.474 0.003 -0.018 -0.012 -0.056 0.020 -0.524 22.000 10.526 11.083 2.167 7.459 0.099 241.109 233.435 235.316 223.608 6.064 -0.284 17.593 11.661 11.881 4.338 3.895 -0.169 0.766 0.478 0.480 0.118 0.197 -0.197 0.396 -1.023 -0.719 -3.366 1.204 -0.514 5.605 3.707 3.718 1.341 1.254 -0.157 0.859 0.858 0.883 0.773 0.787 0.762 0.847 0.744 0.695 0.677 0.631 0.809 0.726 0.684 0.788 0.607 0.724 0.591 0.766 0.000 0.744 0.468 0.598 0.680 0.000 0.551 0.626 0.491 0.750 0.563 0.648 0.654
106898 3.391 75.776 1.000 0.978 0.987 0.975 0.967 0.981 0.892 1.000 0.980 0.989 0.980 0.975 0.983 0.910 0 0 0 40.972 0.000 0 0 0 0 0 0 0 0 0 0 59.028 0 0 0.000 0.000 0 95.079 0.000 0 0 0 0.000 0 0 0 0 0 0 4.921 0 0.000 0 0 0 0 0.000 0 0.000 0 0.0 0 0.000 0.000 15.865 0.000 78.466 1.286 0.000 0.000 0.000 4.383 0.007 0 0.000 0.000 0.000 0.000 1.985 0.000 46.502 0.000 45.544 5.097 0.000 0.000 0.000 0.866 0 0 0.000 0.000 0 0 0.000 0.000 0.000 0.000 100.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.349 0.000 98.651 0.000 0.000 0.000 0.000 0.000 561.389 14.134 6.153 -1.833 484.746 11.495 3.908 -3.684 106.00 108.153 16.850 56.820 6.440 71.718 1.010 13.840 80.250 182.880 1.681 25.854 66.136 181.420 0.054 6 4.944 3 1.001 6.277 3.188 3.039 0.007 1.764 -0.086 0.020 0.000 0.001 -0.050 0.005 -2.708 0.059 -0.002 0.000 -0.091 0.015 -1.834 128 9.197 8 1 10.755 4.412 404.347 345.802 341.377 310.978 22.211 0.624 31.097 7.671 6.719 0.911 4.736 1.398 1.247 0.352 0.311 0.005 0.220 1.160 3.361 -0.042 0.018 -3.683 0.653 -1.349 8.289 2.341 2.085 0.287 1.429 1.172 3.961 1.930 1.949 0.902 0.903 0.342 0.009 0.001 0.001 -0.014 0.005 -0.805 0.026 -0.010 -0.007 -0.059 0.023 -0.515 48.500 5.326 1.000 1.000 11.470 NA 325.718 318.516 318.467 313.572 3.779 0.505 16.530 10.037 9.810 4.496 3.274 0.276 0.811 0.464 0.436 0.132 0.180 0.238 1.543 -0.254 -0.085 -2.478 1.053 -0.714 5.526 3.192 3.017 1.311 1.149 0.370 0.913 0.841 0.786 0.724 0.745 0.685 0.620 0.628 0.577 0.543 0.483 0.484 0.431 0.466 0.440 0.455 0.367 0.380 0.409 0.409 0.451 0.406 0.404 0.395 0.466 0.362 0.405 0.513 0.449 0.404 0.368 0.362
71212 3.090 15.234 1.000 0.979 0.985 0.989 0.985 0.983 0.923 1.000 0.980 0.989 0.994 0.995 0.990 0.949 0 0 0 100.000 0.000 0 0 0 0 0 0 0 0 0 0 0.000 0 0 0.000 0.000 0 100.000 0.000 0 0 0 0.000 0 0 0 0 0 0 0.000 0 0.000 0 0 0 0 0.000 0 0.000 0 0.0 0 1.885 0.000 26.873 0.000 47.849 22.910 0.000 0.000 0.000 0.483 0.000 0 0.000 0.000 0.000 0.000 1.885 0.000 26.873 0.000 47.849 22.910 0.000 0.000 0.000 0.483 0 0 0.000 0.000 0 0 0.000 0.000 0.000 0.000 1.724 87.931 0.000 0.000 0.000 10.345 0.140 0.000 1.183 0.140 0.009 0.000 1.622 0.021 25.500 0.020 35.210 23.034 5.517 4.328 0.179 3.096 416.166 16.435 8.518 0.597 416.166 16.435 8.518 0.597 67.00 67.000 36.047 22.123 36.047 22.123 0.554 1.317 84.623 182.880 0.554 1.317 84.623 182.880 0.052 7 17.985 2 0.672 6.229 3.610 3.119 0.992 1.489 0.107 0.026 0.001 0.002 -0.068 0.008 -2.257 0.027 -0.003 -0.001 -0.069 0.010 -2.462 64 8.393 4 1 7.524 0.217 813.296 697.510 688.268 633.075 36.338 0.796 26.005 10.702 10.854 2.412 3.471 0.112 1.071 0.472 0.479 0.025 0.166 -0.051 1.972 -0.060 0.011 -3.549 0.563 -1.562 7.112 3.220 3.216 0.712 1.074 0.265 -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA 653.921 649.534 649.436 645.154 2.731 0.039 -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA 0.880 0.851 0.845 0.800 0.846 0.744 0.746 0.732 0.746 0.737 0.714 0.709 0.720 0.651 0.685 0.715 0.617 0.638 0.615 0.604 0.587 0.534 0.551 0.504 0.478 0.451 0.430 0.412 0.395 0.377 0.332 0.328
104640 4.808 12.181 0.999 0.981 0.988 0.981 0.976 0.982 0.911 0.999 0.982 0.989 0.991 0.990 0.987 0.940 0 0 0 100.000 0.000 0 0 0 0 0 0 0 0 0 0 0.000 0 0 0.000 0.000 0 100.000 0.000 0 0 0 0.000 0 0 0 0 0 0 0.000 0 0.000 0 0 0 0 0.000 0 0.000 0 0.0 0 0.000 0.000 80.567 15.789 3.036 0.607 0.000 0.000 0.000 0.000 0.000 0 0.000 0.000 0.000 0.000 0.000 0.000 80.567 15.789 3.036 0.607 0.000 0.000 0.000 0.000 0 0 2.489 0.000 0 0 0.000 0.226 66.742 30.317 0.226 0.000 0.000 0.000 0.000 0.000 0.000 0.000 2.489 0.000 0.000 0.000 0.000 0.226 66.742 30.317 0.226 0.000 0.000 0.000 0.000 0.000 3239.750 17.588 11.979 6.365 3239.750 17.588 11.979 6.365 1335.00 1335.000 23.240 31.510 23.240 31.510 1.690 2.720 110.360 182.880 1.690 2.720 110.360 182.880 0.026 960 38.944 2 0.592 6.282 2.827 1.529 0.002 2.233 0.552 0.014 0.000 0.000 -0.026 0.002 -3.456 0.026 -0.001 0.000 -0.043 0.004 -3.309 128 8.807 1 1 16.890 5.009 76.720 42.920 40.844 27.716 10.314 1.018 15.020 2.462 2.077 0.121 2.004 2.308 0.729 0.114 0.094 0.004 0.097 2.224 0.952 -0.025 -0.005 -1.894 0.182 -3.507 4.709 0.760 0.628 0.026 0.607 2.309 3.287 1.884 1.491 0.729 0.862 -0.105 0.005 -0.001 0.000 -0.006 0.003 -0.606 0.009 -0.003 -0.003 -0.012 0.005 0.205 43.333 9.500 6.000 1.000 13.944 NA 37.561 35.502 35.862 32.784 1.551 -0.256 2.518 1.610 1.515 0.992 0.448 -0.012 0.120 0.070 0.066 0.028 0.023 -0.231 0.376 -0.111 -0.117 -0.453 0.212 0.077 0.837 0.506 0.462 0.338 0.137 0.349 0.695 0.887 0.884 0.891 0.926 0.899 0.822 0.828 0.902 0.813 0.741 0.746 0.716 0.616 0.630 0.631 0.637 0.602 0.575 0.509 0.492 0.507 0.448 0.468 0.419 0.419 0.397 0.385 0.366 0.464 0.420 0.474
200288 6.838 43.591 0.961 0.975 0.973 0.934 0.961 0.956 0.783 0.988 0.989 0.981 0.988 0.987 0.972 0.907 0 0 0 42.227 57.773 0 0 0 0 0 0 0 0 0 0 0.000 0 0 0.000 0.000 0 54.932 29.111 0 0 0 0.000 0 0 0 0 0 0 15.478 0 0.479 0 0 0 0 3.628 0 0.000 0 0.0 0 0.909 0.000 38.922 0.000 29.975 25.069 0.728 0.002 0.000 0.767 1.156 0 0.833 0.074 0.017 0.001 1.150 0.010 35.619 0.004 34.848 18.542 3.845 2.753 0.066 1.083 0 0 0.000 0.000 0 0 0.000 0.000 100.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.726 78.125 1.166 18.983 0.000 0.000 0.000 0.000 0.000 434.153 17.102 9.833 2.559 474.029 15.277 7.987 0.693 250.00 91.505 29.938 25.784 27.102 32.325 1.334 2.836 89.788 182.446 1.118 4.077 97.365 171.960 0.001 3320 0.035 2 0.888 6.283 1.999 1.571 0.000 1.637 0.836 0.001 0.000 0.000 -0.002 0.000 -0.008 0.001 0.000 0.000 -0.002 0.000 0.112 128 21.000 1 1 38.572 2.124 1284.345 1280.846 1280.542 1279.920 0.801 1.879 0.547 0.166 0.123 0.000 0.137 0.624 0.023 0.007 0.005 0.000 0.006 0.572 0.083 -0.001 0.000 -0.099 0.022 0.045 0.176 0.050 0.043 0.000 0.041 0.564 -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA 1280.469 1280.406 1280.397 1280.361 NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA 0.000 0.000 0.500 0.500 0.522 0.500 0.470 0.500 0.500 0.832 0.483 0.688 0.560 0.432 0.636 0.491 0.721 0.687 0.668 0.719 0.700 0.688 0.735 0.577 0.742 0.621 0.683 0.677 0.569 0.691 0.497 0.533
474515 2.885 5.252 1.000 0.971 0.997 0.999 1.000 0.987 0.955 1.000 0.971 0.998 0.996 1.000 0.992 0.958 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0.000 54.559 0 39.504 0.000 0 0 0 0.000 0 0 0 0 0 0 5.937 0 0.000 0 0 0 0 0.000 0 0.000 0 62.5 0 37.500 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.621 0 13.319 7.989 4.213 0.159 0.181 0.015 1.551 1.710 55.367 12.602 1.916 0.000 0.334 0.023 0 0 22.472 1.966 0 0 0.000 0.000 21.910 1.685 51.966 0.000 0.000 0.000 0.000 0.000 0.000 0.000 22.472 1.966 0.000 0.000 0.000 0.000 21.910 1.685 51.966 0.000 0.000 0.000 0.000 0.000 370.256 18.479 16.061 13.635 510.190 23.697 17.630 11.558 NA 145.247 NA NA 24.475 27.587 NA NA NA NA 0.925 4.370 60.185 175.831 0.037 295 0.320 2 1.027 6.283 3.580 4.495 0.002 2.281 -0.478 0.014 0.000 0.000 -0.015 0.002 -0.502 0.015 0.000 0.000 -0.016 0.002 -0.795 128 21.010 16 1 30.107 2.791 2769.465 2741.920 2741.289 2732.860 7.337 0.545 9.352 1.761 1.580 0.000 1.298 1.198 0.391 0.073 0.067 0.000 0.055 1.226 0.867 -0.019 -0.001 -0.870 0.135 -0.526 3.089 0.540 0.504 0.000 0.400 1.400 -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA 2739.646 2737.482 2737.479 2735.674 NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA -Inf NA NA Inf NA NA 0.782 0.924 0.970 0.981 0.948 0.934 0.897 0.843 0.859 0.909 0.849 0.802 0.813 0.818 0.732 0.799 0.687 0.687 0.675 0.725 0.516 0.677 0.589 0.609 0.586 0.482 0.495 0.528 0.498 0.473 0.504 0.493
157401 3.230 118.786 0.935 0.932 0.943 0.934 0.907 0.926 0.644 0.961 0.961 0.964 0.960 0.952 0.953 0.776 0 0 0 100.000 0.000 0 0 0 0 0 0 0 0 0 0 0.000 0 0 0.181 0.000 0 84.002 0.000 0 0 0 15.817 0 0 0 0 0 0 0.000 0 0.000 0 0 0 0 0.000 0 1.031 0 0.0 0 0.000 0.000 83.505 0.000 15.464 0.000 0.000 0.000 0.000 0.000 0.000 0 0.261 0.019 0.000 0.000 0.000 0.000 75.222 0.000 22.986 0.631 0.000 0.000 0.000 0.881 0 0 5.708 0.000 0 0 0.000 0.000 68.950 0.000 25.342 0.000 0.000 0.000 0.000 0.000 0.000 0.000 5.708 0.000 0.000 0.000 0.000 0.000 68.950 0.000 25.342 0.000 0.000 0.000 0.000 0.000 1454.515 17.191 9.930 2.665 1392.046 14.852 8.434 2.010 803.00 800.087 15.957 43.967 10.114 55.211 1.930 5.011 123.535 180.324 2.017 9.962 151.288 148.726 0.004 269 0.524 3 0.734 6.244 3.734 3.384 0.003 1.454 -0.051 0.055 0.001 0.000 -0.046 0.007 0.859 0.050 -0.002 0.000 -0.042 0.007 -0.491 128 9.518 16 1 10.842 4.983 117.220 60.774 55.472 47.111 14.565 1.674 22.274 4.312 2.182 0.068 4.668 1.739 0.956 0.182 0.089 0.001 0.198 1.714 2.946 -0.017 -0.005 -2.045 0.434 0.729 6.655 1.271 0.635 0.014 1.377 1.748 2.783 2.260 2.207 1.991 0.203 0.639 0.001 0.000 0.000 -0.001 0.001 0.348 0.003 -0.003 -0.004 -0.006 0.002 0.795 8.500 2.266 1.000 1.000 2.606 NA 50.361 49.203 49.106 48.498 0.548 0.543 2.087 1.374 1.481 0.527 0.452 -0.346 0.094 0.057 0.062 0.019 0.022 -0.278 0.103 -0.086 -0.120 -0.207 0.088 0.620 0.605 0.385 0.415 0.133 0.142 -0.280 0.728 0.717 0.784 0.755 0.704 0.719 0.629 0.709 0.760 0.670 0.651 0.596 0.550 0.665 0.548 0.489 0.577 0.507 0.606 0.606 0.500 0.601 0.518 0.538 0.541 0.522 0.526 0.505 0.348 0.478 0.408 0.418
90261 4.736 958.816 0.999 0.986 0.998 0.999 0.999 0.998 0.980 0.987 0.975 0.981 0.984 0.983 0.978 0.892 0 0 0 95.531 0.000 0 0 0 0 0 0 0 0 0 0 4.469 0 0 0.000 0.000 0 92.051 0.000 0 0 0 0.000 0 0 0 0 0 0 7.949 0 0.000 0 0 0 0 0.000 0 7.402 0 0.0 0 0.000 1.816 6.844 0.000 4.888 78.212 0.000 0.000 0.000 0.838 0.010 0 0.330 0.000 0.000 0.000 0.518 7.577 45.108 0.202 22.064 24.139 0.000 0.000 0.021 0.031 0 0 0.000 0.000 0 0 0.000 0.000 100.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 2.583 0.000 0.000 0.000 0.296 0.000 94.623 0.000 2.498 0.000 0.000 0.000 0.000 0.000 512.333 24.893 17.035 9.172 782.823 19.753 13.468 7.179 285.00 286.404 13.700 60.255 16.704 54.322 0.580 8.570 61.350 181.380 0.652 7.696 107.035 182.639 0.001 108 1.750 4 0.808 6.281 4.019 4.936 0.006 1.731 -0.174 0.020 0.000 0.000 -0.015 0.003 0.864 0.029 -0.001 -0.001 -0.047 0.005 -1.834 128 9.762 16 1 11.946 5.587 128.507 65.343 55.893 45.160 21.523 0.960 23.192 5.338 5.324 0.042 4.248 0.657 0.964 0.226 0.218 0.000 0.183 0.640 1.547 -0.036 -0.020 -2.898 0.276 -1.808 7.581 1.576 1.608 0.009 1.264 0.802 5.636 3.250 3.544 0.645 1.403 0.468 0.001 -0.001 0.000 -0.004 0.001 -1.463 0.001 -0.002 -0.002 -0.005 0.002 0.394 90.667 13.947 6.000 1.000 27.319 2.127 45.737 45.407 45.330 45.185 0.190 0.661 0.883 0.483 0.488 0.161 0.178 0.138 0.040 0.016 0.015 0.001 0.010 0.373 0.025 -0.099 -0.099 -0.199 0.060 0.418 0.276 0.158 0.155 0.047 0.056 0.118 0.873 0.948 0.927 0.842 0.901 0.809 0.784 0.743 0.673 0.704 0.735 0.632 0.613 0.670 0.644 0.592 0.678 0.633 0.705 0.614 0.641 0.560 0.466 0.563 0.525 0.508 0.524 0.437 0.421 0.397 0.420 0.424
SSCT_data
SiteID POINT_Y POINT_X ward
NCC_PH_1_770 40.153 -123.997 8
SAC_WS_3 39.915 -122.532 8
SC_CL_4_59606 32.753 -116.452 10
SAC_LSR_41 41.208 -122.271 1
SCC_CH_5_34910 36.390 -121.595 10
SFE_2018_1640 39.983 -123.822 4
SC_UL_2_16853 34.442 -119.712 5
SAC_WS_505PS0174 40.999 -122.060 1
SC_PL_1_28522 34.099 -118.715 4
SC_PL_4_28521 34.097 -118.717 9

Data transformation

We transform the data to retain what we need for the analysis:

drivers_data <- cbind(SSCT_labelled_data, SSCT_data) %>%
    dplyr::select(c("layer_sd.rstr", "slope_mean.rstr", "SLOPE", "RUSLE", "ward"))  %>%
    # dplyr::select(c("CONFINEMEN", "layer_sd.rstr", "H.960", "slope_mean.rstr", "SLOPE", "RUSLE", "ward")) %>%
    dplyr::rename(ward.grp = ward) %>%
    dplyr::mutate(RUSLE = log10(RUSLE)) %>%
    reshape2::melt(id.vars = "ward.grp")
head(drivers_data)
#>   ward.grp      variable    value
#> 1        8 layer_sd.rstr 38.79167
#> 2        5 layer_sd.rstr 34.09278
#> 3        8 layer_sd.rstr 29.79244
#> 4        5 layer_sd.rstr 16.15452
#> 5        9 layer_sd.rstr 35.78888
#> 6       10 layer_sd.rstr 50.37402

Statistical testing and visualization

We can now visualize the results from a pair-wise statistical comparison with significative_diff_boxplot().

Box plots

Dunn_results <- significative_diff_boxplot(drivers_data, type = "Dunn", padjmeth = "bonferroni", pvalue = 0.10, no_label = TRUE)
Dunn_results$p

Matrices

As the boxplots might be difficult to read, the p-value can be recasted into a matrix form with plot_pairwise_pvalue().

plotly::ggplotly(
    plot_pairwise_pvalue(Dunn_results$stats),
    tooltip = c("group1", "group2", "p.value")
)   

Lollipop charts

We can repeat the above procedure for all data and visualize the results at the number of significant pairwise comparisons. This time we select most data, removing the data with near-zero variance and the Indices of Catchment Integrity. significative_count_lollipop() produces a lollipop chart counting the number of significant pairwise comparison. If normalize = TRUE, the counts are normalized by the number of pairwise comparisons.

all_drivers_data <- cbind(SSCT_labelled_data, SSCT_data %>% dplyr::select("ward")) %>%
    dplyr::select(-caret::nearZeroVar(., names = TRUE)) %>%
    dplyr::select(-c("CHYD", "CCHEM", "CSED", "CCONN", "CTEMP", "CHABT", "ICI", "WHYD", "WCHEM", "WSED", "WCONN", "WTEMP", "WHABT")) %>%
    dplyr::select(-dplyr::contains("H.")) %>%
    dplyr::rename(ward.grp = ward) %>%
    dplyr::mutate(RUSLE = log10(RUSLE)) %>%
    reshape2::melt(id.vars = "ward.grp") %>%
    na.omit()
Dunn_results <- significative_diff_boxplot(all_drivers_data, type = "Dunn", padjmeth = "bonferroni", pvalue = 0.10)
significative_count_lollipop(Dunn_results, pvalue = 0.10, first = 15)$p %>% plotly::ggplotly() 
significative_count_lollipop(Dunn_results, pvalue = 0.10, normalize = TRUE, first = 15)$p %>% plotly::ggplotly()