knitr::opts_chunk$set( collapse = TRUE, comment = "#>", fig.path = "man/figures/README-", out.width = "100%" )
The raceland package implements a computational framework for a pattern-based, zoneless analysis, and visualization of (ethno)racial topography (Dmowska et al., 2020). It is a reimagined approach for analyzing residential segregation and racial diversity based on the concept of ‘landscape’ used in the domain of landscape ecology. A racial landscape, represented by a high-resolution raster grid with each cell containing only inhabitants of a single race, is quantified by two metrics (entropy and mutual information) derived from Information Theory concept (IT). Entropy is the measure of racial diversity and mutual information measures racial segregation.
Racial landscape method is based on the raster gridded data, and unlike the previous methods, does not depend on the division of specific zones (census tract, census block, etc.). Calculation of racial diversity (entropy) and racial segregation (mutual information) can be performed for the whole area of interests (i.e., metropolitan area) without introducing any arbitrary divisions. Racial landscape method also allows for performing calculations at different spatial scales.
You can install the released version of raceland from CRAN with:
install.packages("raceland")
You can install the development version from GitHub with:
# install.packages("remotes") remotes::install_github("Nowosad/raceland")
library(raceland) library(terra)
# Plot the input data race_raster = rast(system.file("extdata/race_raster.tif", package = "raceland")) plot(race_raster)
# Construct racial landscape real_raster = create_realizations(x = race_raster, n = 100) race_colors = c("#F16667", "#6EBE44", "#7E69AF", "#C77213","#F8DF1D") plot(real_raster, col = race_colors, maxnl = 9)
# Plot racial ladnscape plot_realization(x = real_raster[[1]], y = race_raster, hex = race_colors)
# Calculate local subpopulation densities dens_raster = create_densities(real_raster, race_raster, window_size = 10) plot(dens_raster, maxnl = 9)
# Calculate IT-metrics metr_df = calculate_metrics(x = real_raster, w = dens_raster, neighbourhood = 4, fun = "mean", size = NULL, threshold = 1) head(metr_df)
# Summarize IT metrics summary(metr_df[, c("ent", "mutinf")])
Contributions to this package are welcome. The preferred method of contribution is through a GitHub pull request. Feel free to contact us by creating an issue.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.