require(randomForest)

require(randomForestExplainer)

randomF <- randomForest(max_intensity ~ nrgroups + GDPlog_lag + logPopulation + Polity_lag + Asia + Africa + MiddleEast + Europe + Americas, data=MAR_regressions, na.action=na.exclude)

plot(randomF, type="l", main= "Random Forest Protest and Rebellion")

varImpPlot(randomF, main="Variable Importance Random Forest Prot & Reb", col="blue")

plot_min_depth_distribution(randomF)