distances <- french_cities_distances[1:15, 1:15] / 1000 ## convert to km
production <- log(french_cities$population[1:15])
attractiveness <- log(french_cities$area[1:15])
all_flows_log <- grid_blvim(
distances, production, c(1.1, 1.25, 1.5),
c(1, 2, 3, 4) / 500, attractiveness,
epsilon = 0.1,
bipartite = FALSE,
iter_max = 750
)
as.data.frame(all_flows_log, sim_column = "log flows")
Run the code above in your browser using DataLab