Dear Ron,
Thanks for the help and the code. I found the cut option and I used a quantile method to class the data and compare the results from other method. Now, I have 3 Durations from 3 methods combined in shp
duration1= 0-6642
duration2=0-2089
duration3=0-4367
This is what I did
a<-cut(shp$Duration1, breaks=data.frame(classIntervals(shp$Duration1, n=5,method="quantile")[2])[,1],include.lowest=T)
f1= ggplot(shp) + geom_sf(aes(colour=a))
b<-cut(shp$Duration2, breaks=data.frame(classIntervals(shp$Duration2, n=5,method="quantile")[2])[,1],include.lowest=T)
f2= ggplot(shp) + geom_sf(aes(colour=b))
c<-cut(shp$Duration3, breaks=data.frame(classIntervals(shp$Duration3, n=5,method="quantile")[2])[,1],include.lowest=T)
f2= ggplot(shp) + geom_sf(aes(colour=c))
Now, I need to create a common legend and plot in one page. i tried grid_arrenge_shared_legend
grid_arrange_shared_legend(f1, f2, f3, ncol = 3, nrow = 1, position = "bottom")
But this is not creating a common legend, only display for the first plot (f1). any idea? thanks again for your help and time.
Best regards,
Solomon