Skip to contents

Create a complete ggplot for a spatial interaction models data frame

Source: R/sim_df_autoplot.R
autoplot.sim_df.Rd

This function uses a tile plot from ggplot2 to display a single value for each of the parameter pairs used to produce the collection of spatial interaction models.

Usage

# S3 method for class 'sim_df'
autoplot(object, value, inverse = TRUE, ...)

Arguments

object

a data frame of spatial interaction models, an object of class sim_df

value

the value to display, defaults to diversity if unspecified

inverse

whether to use the cost scale parameter \(1/\beta\) (default, TRUE) or the inverse cost scale parameter \(\beta\) (FALSE) on the horizontal axis

...

additional parameters (not used currently)

Value

a ggplot object

Details

The value to display is specified via an expression evaluated in the context of the data frame. It defaults to the diversity as computed by diversity(). See below for examples of use.

The horizontal axis is used by default for the cost scale parameter, that is \(1/\beta\). This is generally easier to read than the inverse cost scale. The inverse parameter can be used to turn off this feature. The vertical axis is used by default for the return to scale parameter.

See also

sim_df(), diversity(), grid_autoplot(), grid_var_autoplot()

Examples

positions <- as.matrix(french_cities[1:10, c("th_longitude", "th_latitude")])
distances <- french_cities_distances[1:10, 1:10] / 1000 ## convert to km
production <- rep(1, 10)
attractiveness <- log(french_cities$area[1:10])
all_flows <- grid_blvim(distances, production, seq(1.05, 1.45, by = 0.1),
  seq(1, 3, by = 0.5) / 400,
  attractiveness,
  bipartite = FALSE,
  epsilon = 0.1, iter_max = 1000
)
all_flows_df <- sim_df(all_flows)

## default display: Shannon diversity
ggplot2::autoplot(all_flows_df)

## we leverage non-standard evaluation to compute a different diversity
ggplot2::autoplot(all_flows_df, diversity(sim, "RW"))

## or to refer to columns of the data frame, either default ones
ggplot2::autoplot(all_flows_df, converged)

ggplot2::autoplot(all_flows_df, iterations)

ggplot2::autoplot(all_flows_df, potential)

## or added ones
all_flows_df["Nystuen & Dacey"] <- diversity(sim_column(all_flows_df), "ND")
ggplot2::autoplot(all_flows_df, `Nystuen & Dacey`)