Spatially Variant Carbon Response Functions

carbon.RF %>%
  knitr::kable("html")%>%
  kable_material(c("striped", "hover")) %>%
  kable_styling(position = "center", full_width = F, fixed_thead = T)

Year	MgmtID	mean.C	RF.mean	RF.rd	spatial_predictor	predictor_value
2045	FIC1	55.083273	0.4233719	0.423	FRG	1
2045	FIC2	48.738645	0.2594244	0.259	FRG	1
2045	FIC3	44.398379	0.1472703	0.147	FRG	1
2045	FIC4	46.485945	0.0922138	0.092	FRG	1
2045	FIC5	24.661991	-0.4717252	-0.472	FRG	1
2045	FIC6	6.338624	-0.8362076	-0.836	FRG	1
2045	FIC1	30.729190	0.6689646	0.669	FRG	2
2045	FIC2	26.723818	0.4514248	0.451	FRG	2
2045	FIC3	24.076983	0.3076698	0.308	FRG	2
2045	FIC4	25.219283	0.2607104	0.261	FRG	2
2045	FIC5	13.527258	-0.3743072	-0.374	FRG	2
2045	FIC6	3.377306	-0.8165717	-0.817	FRG	2
2045	FIC1	56.476176	0.3998820	0.400	FRG	3
2045	FIC2	49.468151	0.2261732	0.226	FRG	3
2045	FIC3	44.798312	0.1104214	0.110	FRG	3
2045	FIC4	46.886683	0.0531861	0.053	FRG	3
2045	FIC5	24.851877	-0.4929934	-0.493	FRG	3
2045	FIC6	6.247995	-0.8451302	-0.845	FRG	3
2045	FIC1	37.527060	0.2352837	0.235	FRG	4
2045	FIC2	30.843008	0.0117015	0.012	FRG	4
2045	FIC3	27.836651	-0.0888614	-0.089	FRG	4
2045	FIC4	29.150006	-0.1539296	-0.154	FRG	4
2045	FIC5	15.820670	-0.5997971	-0.600	FRG	4
2045	FIC6	4.940831	-0.8547287	-0.855	FRG	4
2045	FIC1	30.725233	0.1634964	0.163	FRG	5
2045	FIC2	26.719476	0.0118073	0.012	FRG	5
2045	FIC3	24.072794	-0.0884167	-0.088	FRG	5
2045	FIC4	25.214902	-0.1541676	-0.154	FRG	5
2045	FIC5	13.525141	-0.5968329	-0.597	FRG	5
2045	FIC6	3.377033	-0.8721192	-0.872	FRG	5

plot = plot_ly(data=carbon.RF, x = ~MgmtID, y =~RF.rd, type = 'bar', color = ~as.factor(predictor_value))
plot <- plot %>% 
  layout(title = "Carbon Response Functions by FRG",
         legend = list(
           title = list(
             text = 'Fire Regime Group'
           )
         )
  )

plot