label on map: add color for a given factor

[priviere]

Example with cluster from bayesian gxe model

data_cluster = out_model_bh_GxE_vec_variables$out_parameter_groups_t$clust$clust
data_cluster$parameter = rownames(data_cluster)

vec_loc = vec_lat = vec_long = NULL
for(i in 1:nrow(data_cluster)){
  loc = unlist(strsplit(sub("theta\\[", "", as.character(data_cluster[i, "parameter"])), ":"))[1]
  vec_loc = c(vec_loc , loc)
  lat = data[which(data$location == loc), "lat"][1]
  vec_lat = c(vec_lat, lat)
  long = data[which(data$location == loc), "long"][1]
  vec_long = c(vec_long, long)
}
data_cluster$location = vec_loc
data_cluster$lat = vec_lat
data_cluster$long = vec_long

p = plot(data, plot_type = "map")
p = p$map + geom_nodelabel_repel(data = data_cluster, aes(x=as.numeric(long), y=as.numeric(lat), fill = clust, label=location))
p