Wbp.29=FVS%>%
filter(Species=="PIAL")%>% #restrict StdStk table to WbP species only
filter(Year==2029)%>% #Year 2029 is the SARA value year
filter(LiveTpa>5) #filter 2029 StdStk table by SARA value definition (>5 Tpa)
Wbp.35=FVS%>%
filter(Species=="PIAL")%>% #grab a 2035 snapshot for visualization
filter(Year==2035)
Wbp=FVS%>%
filter(Species=="PIAL")%>%
filter(Year>2019)%>%
filter(Year<2050)%>%
filter(StandID %in% Wbp.29$StandID)%>%
group_by(CaseID)%>%
mutate(dTpa=LiveTpa-LiveTpa[4])%>%
mutate(pDeltaTpa = dTpa/LiveTpa[4])
Wbp.S = Wbp %>%
group_by(MgmtID, Year, LiveTpa) %>%
summarise()
head(Wbp)%>%
kbl()%>%
kable_minimal()
MgmtID | StandID | Year | Species | DBHClass | LiveTpa | MrtTPA | HrvTPA | RsdTPA | LiveBA | MrtBA | HrvBA | RsdBA | LiveTCuFt | MrtTCuFt | HrvTCuFt | RsdTCuFt | LiveMCuFt | MrtMCuFt | HrvMCuFt | RsdMCuFt | LiveBdFt | MrtBdFt | HrvBdFt | RsdBdFt | RunTitle | CaseID | dTpa | pDeltaTpa |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
CMCC | 10513 | 2025 | PIAL | All | 23.30804 | 0.3538850 | 0 | 23.30804 | 20.85810 | 0.2799503 | 0 | 20.85810 | 455.5582 | 5.9491976 | 0 | 455.5582 | 450.2285 | 5.8767860 | 0 | 450.2285 | 2396.320 | 30.797997 | 0 | 2396.320 | Mechanical Clearcut (CMCC) 1 | 4511ec67-9af4-47d8-b35d-2bddb8397a98 | 0.4557459 | 0.0199431 |
CMCC | 10513 | 2029 | PIAL | All | 22.95957 | 0.3484694 | 0 | 22.95957 | 22.19158 | 0.3005502 | 0 | 22.19158 | 488.3684 | 6.4419498 | 0 | 488.3684 | 476.2163 | 6.2734580 | 0 | 476.2163 | 2606.046 | 33.824861 | 0 | 2606.046 | Mechanical Clearcut (CMCC) 1 | 4511ec67-9af4-47d8-b35d-2bddb8397a98 | 0.1072765 | 0.0046943 |
CMCC | 10513 | 2030 | PIAL | All | 22.87331 | 0.0862572 | 0 | 22.87331 | 22.49967 | 0.0757507 | 0 | 22.49967 | 495.8381 | 1.6253898 | 0 | 495.8381 | 483.4164 | 1.5798903 | 0 | 483.4164 | 2642.762 | 8.500325 | 0 | 2642.762 | Mechanical Clearcut (CMCC) 1 | 4511ec67-9af4-47d8-b35d-2bddb8397a98 | 0.0210194 | 0.0009198 |
CMCC | 10513 | 2031 | PIAL | All | 22.85229 | 0.0210194 | 0 | 22.85229 | 23.05260 | 0.0160744 | 0 | 23.05260 | 508.6741 | 0.3343615 | 0 | 508.6741 | 500.3529 | 0.3292578 | 0 | 500.3529 | 2724.836 | 1.726479 | 0 | 2724.836 | Mechanical Clearcut (CMCC) 1 | 4511ec67-9af4-47d8-b35d-2bddb8397a98 | 0.0000000 | 0.0000000 |
CMCC | 10513 | 2035 | PIAL | All | 22.76833 | 0.0839580 | 0 | 22.76833 | 25.26712 | 0.0714898 | 0 | 25.26712 | 560.4919 | 1.4955262 | 0 | 560.4919 | 546.1117 | 1.4490127 | 0 | 546.1117 | 3024.594 | 7.734826 | 0 | 3024.594 | Mechanical Clearcut (CMCC) 1 | 4511ec67-9af4-47d8-b35d-2bddb8397a98 | -0.0839580 | -0.0036739 |
CMCC | 10513 | 2036 | PIAL | All | 22.74728 | 0.0210532 | 0 | 22.74728 | 25.82911 | 0.0183441 | 0 | 25.82911 | 573.5971 | 0.3841275 | 0 | 573.5971 | 563.0073 | 0.3739107 | 0 | 563.0073 | 3138.175 | 1.999429 | 0 | 3138.175 | Mechanical Clearcut (CMCC) 1 | 4511ec67-9af4-47d8-b35d-2bddb8397a98 | -0.1050111 | -0.0045952 |
(2029 for modeling purposes)
asp.plot=plot_ly(data = Wbp.29, y = ~ LiveTpa, x = ~ MgmtID, type = 'box', color = ~ MgmtID)%>%
layout(title='Whitebark Pine Live Tpa Value Pre-Treatment (2029)')
asp.plot
w = Wbp%>%
filter(MgmtID=='CMCC')
asp.plot=plot_ly(data = w, y = ~ LiveTpa, x = ~ Year, type = 'box', color = ~ MgmtID)%>%
layout(title='Timeseries: Whitebark Pine Live Tpa Value for Clearcut Treatment')
asp.plot
asp.plot=plot_ly(data = w, y = ~ LiveBA, x = ~ Year, type = 'box', color = ~ MgmtID)%>%
layout(title='Timeseries: Whitebark Pine Live BA Value for Clearcut Treatment')
asp.plot
(Species=="PIAL")
Given the listed status of Whitebark pine, regulation prohibits the harvest of these species. Thus all FVS runs excluded whitebark pine from cut lists. The signal of all response functions are the result of mortality following treatment activities, and/or natural mortality (e.g. mortality following mastication even though whitebark pine trees themselves were not masticated).
Treatments show a positive affect for Live BA in WbP, but a negative effect in Live TPA. Given our static SARA value (2029) any loss in individual WbP trees (even due to natural mortality) in years 2030-2035 propagates as a loss in SARA value.
t = get_summary_stats(Wbp.S)
t %>%
kbl()%>%
kable_minimal()
MgmtID | Year | variable | n | min | max | median | q1 | q3 | iqr | mad | mean | sd | se | ci |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
CMCC | 2025 | LiveTpa | 985 | 0.000 | 1496.266 | 72.434 | 17.788 | 168.055 | 150.267 | 89.897 | 137.696 | 190.927 | 6.083 | 11.938 |
CMCC | 2029 | LiveTpa | 984 | 0.806 | 1377.098 | 65.353 | 16.805 | 160.218 | 143.413 | 82.067 | 129.032 | 180.088 | 5.741 | 11.266 |
CMCC | 2030 | LiveTpa | 984 | 0.679 | 1313.679 | 63.809 | 15.945 | 155.678 | 139.732 | 80.601 | 126.205 | 176.134 | 5.615 | 11.019 |
CMCC | 2031 | LiveTpa | 394 | 0.928 | 1030.397 | 65.460 | 17.273 | 159.088 | 141.814 | 86.935 | 127.627 | 167.737 | 8.450 | 16.614 |
CMCC | 2035 | LiveTpa | 394 | 0.928 | 982.915 | 65.253 | 17.220 | 158.603 | 141.383 | 86.691 | 125.668 | 161.409 | 8.132 | 15.987 |
CMCC | 2036 | LiveTpa | 394 | 0.928 | 965.724 | 65.201 | 17.206 | 158.482 | 141.276 | 86.629 | 124.951 | 159.186 | 8.020 | 15.767 |
CMCC | 2040 | LiveTpa | 394 | 0.920 | 925.823 | 64.993 | 17.152 | 157.997 | 140.845 | 86.384 | 122.762 | 152.756 | 7.696 | 15.130 |
CMCC | 2045 | LiveTpa | 394 | 0.902 | 782.697 | 64.673 | 17.084 | 157.393 | 140.309 | 85.987 | 119.715 | 144.102 | 7.260 | 14.273 |
CMUR | 2025 | LiveTpa | 340 | 5.359 | 1140.978 | 72.062 | 17.800 | 164.570 | 146.770 | 89.361 | 126.784 | 166.167 | 9.012 | 17.726 |
CMUR | 2029 | LiveTpa | 340 | 5.120 | 1069.748 | 65.313 | 17.557 | 152.381 | 134.824 | 80.506 | 119.864 | 154.832 | 8.397 | 16.517 |
CMUR | 2030 | LiveTpa | 340 | 5.018 | 1047.170 | 64.804 | 17.365 | 148.654 | 131.289 | 80.035 | 117.637 | 151.416 | 8.212 | 16.152 |
CMUR | 2031 | LiveTpa | 340 | 5.001 | 1032.135 | 64.356 | 17.330 | 148.440 | 131.110 | 79.422 | 116.998 | 150.081 | 8.139 | 16.010 |
CMUR | 2035 | LiveTpa | 340 | 4.932 | 986.688 | 63.939 | 17.191 | 147.839 | 130.648 | 78.581 | 114.868 | 144.810 | 7.853 | 15.448 |
CMUR | 2036 | LiveTpa | 340 | 4.914 | 949.740 | 63.831 | 17.155 | 147.726 | 130.570 | 78.475 | 113.953 | 142.265 | 7.715 | 15.176 |
CMUR | 2040 | LiveTpa | 340 | 4.844 | 823.330 | 63.396 | 16.970 | 147.272 | 130.302 | 77.709 | 111.219 | 135.470 | 7.347 | 14.451 |
CMUR | 2045 | LiveTpa | 340 | 4.674 | 778.793 | 62.846 | 16.591 | 146.009 | 129.418 | 76.986 | 107.513 | 127.928 | 6.938 | 13.647 |
HCTA | 2025 | LiveTpa | 744 | 0.000 | 1496.266 | 72.328 | 17.840 | 182.136 | 164.296 | 90.066 | 142.359 | 196.596 | 7.208 | 14.150 |
HCTA | 2029 | LiveTpa | 743 | 0.806 | 1377.098 | 69.208 | 17.353 | 167.547 | 150.194 | 86.669 | 134.812 | 185.761 | 6.815 | 13.379 |
HCTA | 2030 | LiveTpa | 743 | 0.679 | 1313.679 | 68.202 | 16.960 | 162.421 | 145.461 | 86.519 | 132.241 | 181.770 | 6.668 | 13.091 |
HCTA | 2031 | LiveTpa | 743 | 0.582 | 1272.637 | 65.497 | 16.788 | 161.042 | 144.254 | 84.634 | 130.250 | 178.966 | 6.566 | 12.889 |
HCTA | 2035 | LiveTpa | 742 | 0.000 | 1198.652 | 63.573 | 15.288 | 150.657 | 135.369 | 82.189 | 123.535 | 170.111 | 6.245 | 12.260 |
HCTA | 2036 | LiveTpa | 741 | 0.138 | 1155.259 | 62.294 | 14.910 | 148.979 | 134.069 | 79.949 | 121.382 | 166.852 | 6.129 | 12.033 |
HCTA | 2040 | LiveTpa | 741 | 0.000 | 1088.151 | 58.074 | 12.821 | 142.076 | 129.255 | 76.019 | 114.655 | 158.587 | 5.826 | 11.437 |
HCTA | 2045 | LiveTpa | 740 | 0.000 | 1005.496 | 52.268 | 11.436 | 129.095 | 117.659 | 68.619 | 106.260 | 147.933 | 5.438 | 10.676 |
HERB | 2025 | LiveTpa | 988 | 0.000 | 1496.266 | 72.393 | 17.789 | 168.278 | 150.489 | 89.821 | 138.447 | 193.036 | 6.141 | 12.052 |
HERB | 2029 | LiveTpa | 987 | 0.806 | 1377.098 | 65.176 | 16.883 | 160.379 | 143.496 | 81.814 | 129.734 | 182.016 | 5.794 | 11.369 |
HERB | 2030 | LiveTpa | 987 | 0.679 | 1313.679 | 63.745 | 16.052 | 155.864 | 139.813 | 80.297 | 126.891 | 178.007 | 5.666 | 11.119 |
HERB | 2031 | LiveTpa | 987 | 0.582 | 1272.637 | 62.364 | 15.316 | 152.200 | 136.883 | 78.797 | 124.679 | 175.131 | 5.574 | 10.939 |
HERB | 2035 | LiveTpa | 986 | 0.000 | 1198.652 | 56.307 | 12.068 | 146.130 | 134.062 | 73.839 | 117.343 | 166.238 | 5.294 | 10.389 |
HERB | 2036 | LiveTpa | 985 | 0.103 | 1155.259 | 54.126 | 11.866 | 141.967 | 130.101 | 71.309 | 114.947 | 162.790 | 5.187 | 10.179 |
HERB | 2040 | LiveTpa | 984 | 0.000 | 1088.151 | 50.934 | 11.200 | 134.820 | 123.619 | 66.942 | 108.129 | 154.501 | 4.925 | 9.665 |
HERB | 2045 | LiveTpa | 977 | 0.000 | 1005.496 | 46.200 | 9.893 | 124.277 | 114.384 | 60.121 | 100.198 | 144.084 | 4.610 | 9.046 |
MRCT | 2025 | LiveTpa | 618 | 0.000 | 1247.261 | 79.488 | 19.640 | 209.892 | 190.253 | 101.469 | 154.708 | 197.861 | 7.959 | 15.630 |
MRCT | 2029 | LiveTpa | 617 | 0.806 | 1186.833 | 72.475 | 17.176 | 193.559 | 176.384 | 93.749 | 142.480 | 186.688 | 7.516 | 14.760 |
MRCT | 2030 | LiveTpa | 617 | 0.679 | 1170.509 | 69.668 | 16.232 | 185.853 | 169.621 | 90.860 | 138.576 | 182.833 | 7.361 | 14.455 |
MRCT | 2031 | LiveTpa | 617 | 0.667 | 1160.492 | 69.287 | 15.631 | 185.438 | 169.807 | 90.679 | 137.680 | 181.806 | 7.319 | 14.374 |
MRCT | 2035 | LiveTpa | 617 | 0.549 | 1133.803 | 67.219 | 14.752 | 181.005 | 166.253 | 88.521 | 133.682 | 177.408 | 7.142 | 14.026 |
MRCT | 2036 | LiveTpa | 617 | 0.436 | 1123.942 | 65.799 | 14.130 | 177.540 | 163.410 | 86.741 | 131.761 | 175.235 | 7.055 | 13.854 |
MRCT | 2040 | LiveTpa | 616 | 0.000 | 1099.558 | 62.770 | 12.794 | 169.624 | 156.830 | 83.823 | 126.839 | 169.740 | 6.839 | 13.431 |
MRCT | 2045 | LiveTpa | 615 | 0.093 | 1069.090 | 58.375 | 11.331 | 159.208 | 147.877 | 78.105 | 120.543 | 162.427 | 6.550 | 12.862 |
MTTH | 2025 | LiveTpa | 173 | 5.359 | 1140.978 | 114.940 | 29.773 | 287.287 | 257.514 | 145.063 | 183.690 | 208.003 | 15.814 | 31.215 |
MTTH | 2029 | LiveTpa | 173 | 5.120 | 1069.748 | 106.013 | 29.602 | 263.817 | 234.215 | 132.474 | 171.543 | 193.466 | 14.709 | 29.033 |
MTTH | 2030 | LiveTpa | 173 | 5.018 | 1047.170 | 104.007 | 28.944 | 258.159 | 229.215 | 129.815 | 167.689 | 189.145 | 14.380 | 28.385 |
MTTH | 2031 | LiveTpa | 173 | 4.925 | 1030.902 | 103.833 | 28.036 | 256.224 | 228.188 | 129.702 | 166.370 | 187.330 | 14.242 | 28.112 |
MTTH | 2035 | LiveTpa | 173 | 4.672 | 977.602 | 102.189 | 26.047 | 252.464 | 226.417 | 127.726 | 160.019 | 178.565 | 13.576 | 26.797 |
MTTH | 2036 | LiveTpa | 173 | 4.581 | 964.674 | 101.313 | 25.755 | 252.104 | 226.348 | 126.632 | 157.822 | 175.975 | 13.379 | 26.408 |
MTTH | 2040 | LiveTpa | 173 | 4.374 | 848.656 | 97.096 | 24.064 | 241.430 | 217.366 | 121.663 | 150.686 | 165.915 | 12.614 | 24.899 |
MTTH | 2045 | LiveTpa | 173 | 3.032 | 771.767 | 89.121 | 22.913 | 228.216 | 205.303 | 113.562 | 141.042 | 154.927 | 11.779 | 23.250 |
RMGP | 2025 | LiveTpa | 385 | 5.359 | 1183.188 | 73.405 | 17.966 | 187.174 | 169.209 | 91.356 | 143.768 | 192.602 | 9.816 | 19.300 |
RMGP | 2029 | LiveTpa | 385 | 5.120 | 1069.748 | 71.523 | 17.933 | 168.283 | 150.350 | 88.760 | 136.287 | 179.758 | 9.161 | 18.013 |
RMGP | 2030 | LiveTpa | 385 | 5.018 | 1047.170 | 70.938 | 17.924 | 166.702 | 148.779 | 87.912 | 133.850 | 175.660 | 8.952 | 17.602 |
RMGP | 2031 | LiveTpa | 385 | 4.952 | 953.618 | 64.581 | 17.914 | 154.894 | 136.979 | 79.344 | 123.696 | 159.748 | 8.141 | 16.007 |
RMGP | 2035 | LiveTpa | 385 | 4.538 | 890.154 | 63.709 | 17.843 | 147.784 | 129.941 | 78.413 | 118.344 | 152.209 | 7.757 | 15.252 |
RMGP | 2036 | LiveTpa | 385 | 4.369 | 844.378 | 62.704 | 17.832 | 145.502 | 127.670 | 76.301 | 116.100 | 148.485 | 7.567 | 14.879 |
RMGP | 2040 | LiveTpa | 385 | 3.777 | 784.424 | 59.041 | 17.153 | 137.059 | 119.906 | 71.168 | 109.567 | 138.500 | 7.059 | 13.878 |
RMGP | 2045 | LiveTpa | 385 | 2.535 | 741.658 | 55.199 | 15.891 | 124.015 | 108.124 | 65.831 | 101.199 | 126.961 | 6.471 | 12.722 |
RMMA | 2025 | LiveTpa | 581 | 0.000 | 1247.261 | 67.027 | 16.831 | 150.612 | 133.781 | 82.398 | 122.405 | 175.745 | 7.291 | 14.320 |
RMMA | 2029 | LiveTpa | 580 | 0.806 | 1186.833 | 55.879 | 12.808 | 141.281 | 128.473 | 70.189 | 113.059 | 166.457 | 6.912 | 13.575 |
RMMA | 2030 | LiveTpa | 580 | 0.679 | 1170.509 | 52.480 | 12.309 | 137.825 | 125.516 | 67.082 | 110.088 | 163.072 | 6.771 | 13.299 |
RMMA | 2031 | LiveTpa | 580 | 0.667 | 1160.348 | 52.152 | 12.225 | 137.202 | 124.978 | 66.766 | 109.415 | 162.068 | 6.730 | 13.217 |
RMMA | 2035 | LiveTpa | 580 | 0.626 | 1135.341 | 51.159 | 11.942 | 134.391 | 122.448 | 65.771 | 107.408 | 158.524 | 6.582 | 12.928 |
RMMA | 2036 | LiveTpa | 580 | 0.598 | 1126.142 | 50.395 | 11.841 | 133.683 | 121.842 | 65.672 | 106.662 | 156.887 | 6.514 | 12.795 |
RMMA | 2040 | LiveTpa | 580 | 0.514 | 1103.287 | 49.096 | 11.587 | 131.101 | 119.513 | 64.002 | 104.464 | 153.007 | 6.353 | 12.478 |
RMMA | 2045 | LiveTpa | 580 | 0.268 | 1074.578 | 47.163 | 11.298 | 127.115 | 115.817 | 61.350 | 101.488 | 147.931 | 6.142 | 12.064 |
RMTF | 2025 | LiveTpa | 995 | 0.000 | 1496.266 | 72.352 | 17.775 | 168.501 | 150.726 | 89.794 | 138.607 | 193.065 | 6.121 | 12.011 |
RMTF | 2029 | LiveTpa | 996 | 0.806 | 1377.098 | 64.923 | 16.776 | 160.218 | 143.442 | 81.428 | 129.652 | 182.009 | 5.767 | 11.317 |
RMTF | 2030 | LiveTpa | 997 | 0.679 | 1313.679 | 63.439 | 15.811 | 156.238 | 140.427 | 79.844 | 127.674 | 180.025 | 5.701 | 11.188 |
RMTF | 2031 | LiveTpa | 997 | 0.582 | 1272.637 | 60.710 | 15.277 | 153.716 | 138.438 | 77.896 | 125.464 | 177.179 | 5.611 | 11.011 |
RMTF | 2035 | LiveTpa | 998 | 0.000 | 1198.652 | 55.294 | 11.986 | 146.130 | 134.143 | 71.820 | 117.929 | 168.272 | 5.327 | 10.453 |
RMTF | 2036 | LiveTpa | 998 | 0.100 | 1155.259 | 53.730 | 11.781 | 143.257 | 131.476 | 70.711 | 115.862 | 165.085 | 5.226 | 10.255 |
RMTF | 2040 | LiveTpa | 1001 | 0.000 | 1088.151 | 50.387 | 11.162 | 134.794 | 123.632 | 66.129 | 108.750 | 156.672 | 4.952 | 9.717 |
RMTF | 2045 | LiveTpa | 1003 | 0.000 | 1005.496 | 45.250 | 9.881 | 126.165 | 116.284 | 58.943 | 101.656 | 146.913 | 4.639 | 9.103 |
RXAI | 2025 | LiveTpa | 988 | 0.000 | 1496.266 | 72.393 | 17.789 | 168.278 | 150.489 | 89.821 | 138.447 | 193.036 | 6.141 | 12.052 |
RXAI | 2029 | LiveTpa | 987 | 0.806 | 1377.098 | 65.176 | 16.883 | 160.379 | 143.496 | 81.814 | 129.734 | 182.016 | 5.794 | 11.369 |
RXAI | 2030 | LiveTpa | 987 | 0.679 | 1313.679 | 63.745 | 16.052 | 155.864 | 139.813 | 80.297 | 126.891 | 178.007 | 5.666 | 11.119 |
RXAI | 2031 | LiveTpa | 987 | 0.320 | 612.719 | 31.011 | 8.449 | 72.661 | 64.212 | 37.729 | 60.253 | 82.537 | 2.627 | 5.156 |
RXAI | 2035 | LiveTpa | 987 | 0.319 | 605.557 | 29.907 | 8.082 | 71.290 | 63.208 | 36.505 | 59.379 | 81.804 | 2.604 | 5.110 |
RXAI | 2036 | LiveTpa | 987 | 0.318 | 603.777 | 29.785 | 7.915 | 70.898 | 62.983 | 36.408 | 59.047 | 81.551 | 2.596 | 5.094 |
RXAI | 2040 | LiveTpa | 987 | 0.187 | 596.705 | 29.060 | 7.594 | 69.448 | 61.854 | 35.895 | 57.825 | 80.336 | 2.557 | 5.018 |
RXAI | 2045 | LiveTpa | 987 | 0.089 | 592.513 | 27.864 | 7.178 | 66.285 | 59.108 | 34.597 | 56.096 | 78.614 | 2.502 | 4.910 |
RXGF | 2025 | LiveTpa | 969 | 0.000 | 1496.266 | 72.581 | 17.788 | 168.948 | 151.160 | 90.085 | 138.274 | 191.594 | 6.155 | 12.078 |
RXGF | 2029 | LiveTpa | 968 | 0.806 | 1377.098 | 65.353 | 16.805 | 161.035 | 144.230 | 82.172 | 129.579 | 180.703 | 5.808 | 11.398 |
RXGF | 2030 | LiveTpa | 968 | 0.679 | 1313.679 | 63.809 | 15.945 | 157.475 | 141.529 | 80.833 | 126.732 | 176.724 | 5.680 | 11.147 |
RXGF | 2031 | LiveTpa | 968 | 0.339 | 625.599 | 31.898 | 8.660 | 74.939 | 66.279 | 38.766 | 61.522 | 83.740 | 2.692 | 5.282 |
RXGF | 2035 | LiveTpa | 968 | 0.337 | 618.285 | 30.701 | 8.484 | 73.809 | 65.325 | 37.453 | 60.582 | 82.943 | 2.666 | 5.232 |
RXGF | 2036 | LiveTpa | 968 | 0.337 | 616.467 | 30.573 | 8.185 | 72.943 | 64.759 | 37.353 | 60.210 | 82.619 | 2.655 | 5.211 |
RXGF | 2040 | LiveTpa | 968 | 0.199 | 609.247 | 29.779 | 7.717 | 71.825 | 64.108 | 36.750 | 58.906 | 81.284 | 2.613 | 5.127 |
RXGF | 2045 | LiveTpa | 968 | 0.096 | 599.705 | 28.711 | 7.427 | 67.404 | 59.977 | 35.513 | 57.091 | 79.544 | 2.557 | 5.017 |
RF.Wbp = t %>%
group_by(MgmtID)%>%
mutate(Base.Value = mean[2])%>% #set SARA value as fixed at year 2029 (prior to treatment)
mutate(RF.Med = ((mean - Base.Value)/Base.Value))%>%
mutate(RF.rd=round(RF.Med, digits = 3)) #3 decimal places for modeled RFs
RF.Wbp.35 = RF.Wbp%>%
filter(Year == 2035) #only grab RF for 2035 (5 years post-treatment)
RF.Wbp.35%>%
select(MgmtID, Year, RF.Med, RF.rd)%>%
kbl()%>%
kable_minimal()
MgmtID | Year | RF.Med | RF.rd |
---|---|---|---|
CMCC | 2035 | -0.0260711 | -0.026 |
CMUR | 2035 | -0.0416806 | -0.042 |
HCTA | 2035 | -0.0836498 | -0.084 |
HERB | 2035 | -0.0955108 | -0.096 |
MRCT | 2035 | -0.0617490 | -0.062 |
MTTH | 2035 | -0.0671785 | -0.067 |
RMGP | 2035 | -0.1316560 | -0.132 |
RMMA | 2035 | -0.0499828 | -0.050 |
RMTF | 2035 | -0.0904190 | -0.090 |
RXAI | 2035 | -0.5423019 | -0.542 |
RXGF | 2035 | -0.5324705 | -0.532 |
MgmtID | REBA.Code | treatment.name | Year | variable | n | min | max | median | q1 | q3 | iqr | mad | mean | sd | se | ci | Base.Value | RF.Med | RF.rd | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
83 | RXGF | RXF-GFL | Rx Fire - Ground | 2030 | LiveTpa | 968 | 0.679 | 1313.679 | 63.809 | 15.945 | 157.475 | 141.529 | 80.833 | 126.732 | 176.724 | 5.680 | 11.147 | 129.579 | -0.0219712 | -0.022 |
84 | RXGF | RXF-GFL | Rx Fire - Ground | 2031 | LiveTpa | 968 | 0.339 | 625.599 | 31.898 | 8.660 | 74.939 | 66.279 | 38.766 | 61.522 | 83.740 | 2.692 | 5.282 | 129.579 | -0.5252163 | -0.525 |
85 | RXGF | RXF-GFL | Rx Fire - Ground | 2035 | LiveTpa | 968 | 0.337 | 618.285 | 30.701 | 8.484 | 73.809 | 65.325 | 37.453 | 60.582 | 82.943 | 2.666 | 5.232 | 129.579 | -0.5324705 | -0.532 |
86 | RXGF | RXF-GFL | Rx Fire - Ground | 2036 | LiveTpa | 968 | 0.337 | 616.467 | 30.573 | 8.185 | 72.943 | 64.759 | 37.353 | 60.210 | 82.619 | 2.655 | 5.211 | 129.579 | -0.5353414 | -0.535 |
87 | RXGF | RXF-GFL | Rx Fire - Ground | 2040 | LiveTpa | 968 | 0.199 | 609.247 | 29.779 | 7.717 | 71.825 | 64.108 | 36.750 | 58.906 | 81.284 | 2.613 | 5.127 | 129.579 | -0.5454047 | -0.545 |
88 | RXGF | RXF-GFL | Rx Fire - Ground | 2045 | LiveTpa | 968 | 0.096 | 599.705 | 28.711 | 7.427 | 67.404 | 59.977 | 35.513 | 57.091 | 79.544 | 2.557 | 5.017 | 129.579 | -0.5594116 | -0.559 |
These are not errors in FVS model runs, merely a function of how FVS calculate mortality (e.g. change in TPA in the year following a fire), or a function of our multi-step treatments (e.g. pile and burn or mastication in the year following fuels thinning)
FVS_Code
of NA
.
Pending revisions to previously released landscapes (e.g. SOFRC,
Trinity) these treatment response functions may be populated at a later
date.initial_method_code | initial_prescription_code | follow_up_method_code | follow_up_prescription_code | sara_id | spatial_predictor | predictor_value | value | name | objective | treatment_name |
---|---|---|---|---|---|---|---|---|---|---|
CMR | CC | NONE | NONE | 288 | GLOBAL | 0 | -0.026 | Whitebark Pine | bdv | Complex Mechanical Removal - Clearcut |
CMR | TH | NONE | NONE | 288 | GLOBAL | 0 | -0.067 | Whitebark Pine | bdv | Complex Mechanical Removal - Thinning |
CMR | URC | NONE | NONE | 288 | GLOBAL | 0 | -0.042 | Whitebark Pine | bdv | Complex Mechanical Removal - Urban Clearing |
HCD | BR | NONE | NONE | 288 | GLOBAL | 0 | 0.000 | Whitebark Pine | bdv | Herbicides - Broad |
HCD | TA | NONE | NONE | 288 | GLOBAL | 0 | -0.084 | Whitebark Pine | bdv | Herbicides - Targeted |
HRB | HRB | NONE | NONE | 288 | GLOBAL | 0 | -0.096 | Whitebark Pine | bdv | Herbivory - Herbivory |
MRM | CC | NONE | NONE | 288 | GLOBAL | 0 | -0.026 | Whitebark Pine | bdv | Mechanical Removal - Clearcut |
MRM | CT | NONE | NONE | 288 | GLOBAL | 0 | -0.062 | Whitebark Pine | bdv | Mechanical Removal - Commercial Thinning |
MRM | SD | NONE | NONE | 288 | GLOBAL | 0 | 0.000 | Whitebark Pine | bdv | Mechanical Removal - Small Diameter |
MRM | TH | NONE | NONE | 288 | GLOBAL | 0 | 0.000 | Whitebark Pine | bdv | Mechanical Removal - Thinning |
MRM | VDT | NONE | NONE | 288 | GLOBAL | 0 | 0.000 | Whitebark Pine | bdv | Mechanical Removal - Variable Density Thinning |
MTH | ISR | NONE | NONE | 288 | GLOBAL | 0 | -0.090 | Whitebark Pine | bdv | Manual - Invasives Removal |
MTH | TH | NONE | NONE | 288 | GLOBAL | 0 | -0.067 | Whitebark Pine | bdv | Manual - Thinning |
MTH | URC | NONE | NONE | 288 | GLOBAL | 0 | -0.042 | Whitebark Pine | bdv | Manual - Urban Clearing |
REV | A | NONE | NONE | 288 | GLOBAL | 0 | 0.000 | Whitebark Pine | bdv | Revegetation - Revegetation |
RGM | GMP | NONE | NONE | 288 | GLOBAL | 0 | -0.132 | Whitebark Pine | bdv | Mechanical Rearrangement - Grapple Pile |
RGM | MASL | NONE | NONE | 288 | GLOBAL | 0 | -0.050 | Whitebark Pine | bdv | Mechanical Rearrangement - Mastication |
RGM | TFF | NONE | NONE | 288 | GLOBAL | 0 | -0.090 | Whitebark Pine | bdv | Mechanical Rearrangement - Mowing |
RXF | AI | NONE | NONE | 288 | GLOBAL | 0 | -0.542 | Whitebark Pine | bdv | Rx Fire - Rx Fire - Aerial |
RXF | GFL | NONE | NONE | 288 | GLOBAL | 0 | -0.532 | Whitebark Pine | bdv | Rx Fire - Rx Fire - Ground |