skip_if_not_installed("ggridges")
library(ggridges)
test_that(
desc = "ggridges basic ridgelines",
code = {
# simple ridge plot
data <- data.frame(x = 1:5, y = rep(1, 5), height = c(0, 1, 3, 4, 2))
p <- ggplot(data, aes(x, y, height = height)) + geom_ridgeline()
p2 <- ggplotly(p)
expect_doppelganger(p2, 'basic_ridgeline')
# Negative height
data <- data.frame(x = 1:5, y = rep(1, 5), height = c(0, 1, -1, 3, 2))
plot_base <- ggplot(data, aes(x, y, height = height))
## Negative height cut off
p <- plot_base + geom_ridgeline()
p2 <- ggplotly(p)
expect_doppelganger(p2, 'negative_height_cut')
## Negative height allowed
p <- plot_base + geom_ridgeline(min_height = -2)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'negative_height_retained')
# Multiple ridgelines at same time
d <- data.frame(
x = rep(1:5, 3),
y = c(rep(0, 5), rep(1, 5), rep(2, 5)),
height = c(0, 1, 3, 4, 0, 1, 2, 3, 5, 4, 0, 5, 4, 4, 1)
)
p <- ggplot(d, aes(x, y, height = height, group = y)) +
geom_ridgeline(fill = "lightblue")
p2 <- ggplotly(p)
expect_doppelganger(p2, 'multiple_ridgelines')
# stat = identity (works)
p <- ggplot(d, aes(x, y, height = height, group = y)) +
geom_density_ridges(stat = "identity", scale = 1)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'stat_identity')
}
)
test_that(
desc = "ggridges density_ridgeline",
code = {
# Density ridgeline plots
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
geom_density_ridges()
p2 <- ggplotly(p)
expect_doppelganger(p2, 'density_ridgeline')
# geom_density_ridges2 (closed polygon)
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) + geom_density_ridges2()
p2 <- ggplotly(p)
expect_doppelganger(p2, 'density_ridgeline2')
# Grouping aesthetic
# modified dataset that represents species as a number
iris_num <- transform(iris, Species_num = as.numeric(Species))
p <- ggplot(iris_num,
aes(x = Sepal.Length,
y = Species_num,
group = Species_num)) +
geom_density_ridges()
p2 <- ggplotly(p)
expect_doppelganger(p2, 'numeric_grouping')
# Cutting trailing tails (works)
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
geom_density_ridges(rel_min_height = 0.01)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'cutting_tails')
# Non-overlapping ridges (Works)
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
geom_density_ridges(scale = 0.9)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'overlapping_none')
# Exactly touching (Works)
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
geom_density_ridges(scale = 1)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'overlapping_touching')
# scale = 5, substantial overlap (Works)
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
geom_density_ridges(scale = 5)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'overlapping_lot')
# Panel scaling (Works)
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
geom_density_ridges(scale = 1) + facet_wrap(~Species)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'overlapping_facet_touching')
}
)
test_that(
desc = "ggridges fill colours",
code = {
# Varying fill colors along the x axis
# Example 1 (Works, but extra legend that is not shown in ggridges)
d <- data.frame(
x = rep(1:5, 3) + c(rep(0, 5), rep(0.3, 5), rep(0.6, 5)),
y = c(rep(0, 5), rep(1, 5), rep(3, 5)),
height = c(0, 1, 3, 4, 0, 1, 2, 3, 5, 4, 0, 5, 4, 4, 1))
p <- ggplot(d, aes(x, y, height = height, group = y, fill = factor(x+y))) +
geom_ridgeline_gradient() +
scale_fill_viridis_d(direction = -1, guide = "none")
p2 <- ggplotly(p)
expect_doppelganger(p2, 'varying_fill_colours')
# geom_density_ridges_gradient (Doesn't work)
# p <- ggplot(lincoln_weather, aes(x = `Mean Temperature [F]`, y = Month, fill = stat(x))) +
# geom_density_ridges_gradient(scale = 3, rel_min_height = 0.01) +
# scale_fill_viridis_c(name = "Temp. [F]", option = "C") +
# labs(title = 'Temperatures in Lincoln NE in 2016')
# ggplotly(p) # gets stuck
# Stats
## Quantile lines and coloring by quantiles or probabilities (Works)
# quantile multiple lines
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
stat_density_ridges(quantile_lines = TRUE)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'quantile_lines_multi')
# quantile single line
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
stat_density_ridges(quantile_lines = TRUE, quantiles = 2)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'quantile_lines_1')
# quantile by cut points
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
stat_density_ridges(quantile_lines = TRUE,
quantiles = c(0.025, 0.975),
alpha = 0.7)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'quantile_cut_points')
## Colour by quantile
# warning since ggridges uses stat(quantile)
suppressWarnings(
p <- ggplot(iris, aes(x=Sepal.Length, y=Species, fill = factor(stat(quantile)))) +
stat_density_ridges(
geom = "density_ridges_gradient", calc_ecdf = TRUE,
quantiles = 4, quantile_lines = TRUE
) +
scale_fill_viridis_d(name = "Quartiles")
)
suppressWarnings(
p2 <- ggplotly(p)
)
expect_doppelganger(p2, 'quantile_colouring')
# highglight tails of distributions (works)
p <- ggplot(iris, aes(x = Sepal.Length, y = Species, fill = factor(stat(quantile)))) +
stat_density_ridges(
geom = "density_ridges_gradient",
calc_ecdf = TRUE,
quantiles = c(0.025, 0.975)
) +
scale_fill_manual(
name = "Probability", values = c("#FF0000A0", "#A0A0A0A0", "#0000FFA0"),
labels = c("(0, 0.025]", "(0.025, 0.975]", "(0.975, 1]")
)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'quantile_colouring_tails_only')
# mapping prob onto colour (doesn't work)
# p <- ggplot(iris, aes(x = Sepal.Length, y = Species, fill = 0.5 - abs(0.5 - stat(ecdf)))) +
# stat_density_ridges(geom = "density_ridges_gradient", calc_ecdf = TRUE) +
# scale_fill_viridis_c(name = "Tail probability", direction = -1)
# ggplotly(p)
}
)
test_that(
desc = "ggridges points",
code = {
set.seed(123) # make jittering reproducible
# jittering points (works)
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
geom_density_ridges(jittered_points = TRUE)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'jittering points')
# raincloud effect (works)
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
geom_density_ridges(
jittered_points = TRUE, position = "raincloud",
alpha = 0.7, scale = 0.9
)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'raincloud_effect')
# rug effect (doesn't work, point shape not taken into account)
# p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
# geom_density_ridges(
# jittered_points = TRUE,
# position = position_points_jitter(width = 0.05, height = 0),
# point_shape = '|', point_size = 3, point_alpha = 1, alpha = 0.7,
# )
# styling points
p <- ggplot(iris, aes(x = Sepal.Length, y = Species, fill = Species)) +
geom_density_ridges(
aes(point_color = Species, point_fill = Species, point_shape = Species),
alpha = .2, point_alpha = 1, jittered_points = TRUE
) +
scale_point_color_hue(l = 40) +
scale_discrete_manual(aesthetics = "point_shape", values = c(21, 22, 23))
p2 <- ggplotly(p)
expect_doppelganger(p2, 'styling_points')
# styling points 2
p <- ggplot(iris, aes(x = Sepal.Length, y = Species, fill = Species)) +
geom_density_ridges(
aes(point_shape = Species, point_fill = Species, point_size = Petal.Length),
alpha = .2, point_alpha = 1, jittered_points = TRUE
) +
scale_point_color_hue(l = 40) + scale_point_size_continuous(range = c(0.5, 4)) +
scale_discrete_manual(aesthetics = "point_shape", values = c(21, 22, 23))
p2 <- ggplotly(p)
expect_doppelganger(p2, 'styling_points2')
# aesthetics for vertical line (works) (might need to check line on top of points)
p <- ggplot(iris, aes(x = Sepal.Length, y = Species)) +
geom_density_ridges(
jittered_points = TRUE, quantile_lines = TRUE, scale = 0.9, alpha = 0.7,
vline_size = 1, vline_color = "red",
point_size = 0.4, point_alpha = 1,
position = position_raincloud(adjust_vlines = TRUE)
)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'raincloud_vertical_line_points')
}
)
test_that(
desc = "ggridges alternate stats",
code = {
## stat_density_ridges (works)
suppressWarnings({
p <- ggplot(iris, aes(x = Sepal.Length, y = Species, height = stat(density))) +
geom_density_ridges(stat = "density")
p2 <- ggplotly(p)
})
expect_doppelganger(p2, 'stat_density')
skip_if_not_installed("dplyr")
iris_densities <- iris %>%
dplyr::group_by(Species) %>%
dplyr::group_modify(~ ggplot2:::compute_density(.x$Sepal.Length, NULL)) %>%
dplyr::rename(Sepal.Length = x)
p <- ggplot(iris_densities, aes(x = Sepal.Length, y = Species, height = density)) +
geom_density_ridges(stat = "identity")
p2 <- ggplotly(p)
expect_doppelganger(p2, 'manual_densities_stat_identity')
## histograms (works)
p <- ggplot(iris, aes(x = Sepal.Length, y = Species, height = stat(density))) +
geom_density_ridges(stat = "binline", bins = 20, scale = 0.95, draw_baseline = FALSE)
p2 <- ggplotly(p)
expect_doppelganger(p2, 'histogram_ridges')
}
)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.