```
library(GenomeAdmixR)
library(ggplot2)
packageVersion("GenomeAdmixR")
```

`## [1] '2.1.1'`

To demonstrate the applicability and versatility of joyplots, we create a population, and select on this population on two separate locations, for two different ancestors.

```
<- matrix(ncol = 5, nrow = 2)
select_matrix <- 0.1
s 1, ] <- c(0.23, 1.0, 1 + 0.5 * s, 1 + s, 0)
select_matrix[2, ] <- c(0.27, 1.0, 1 + 0.5 * s, 1 + s, 1)
select_matrix[
<- seq(from = 0.2, to = 0.3, length.out = 100) markers
```

We track frequencies of all ancestors in the range [0.2, 0.3], using 100 markers.

```
<- simulate_admixture(
selected_pop module = ancestry_module(number_of_founders = 2,
morgan = 1,
markers = markers),
pop_size = 100,
total_runtime = 1001,
select_matrix = select_matrix)
```

`## Found a selection matrix, performing simulation including selection`

We can now plot all the ancestors over time, we choose here to plot them only for a select subset of times.

```
<- seq(from = 0, to = 1000, by = 100)
time_points plot_joyplot_frequencies(selected_pop$frequencies,
time_points,picked_ancestor = "ALL")
```

Alternatively, we can plot them separately.

```
<- plot_joyplot_frequencies(selected_pop$frequencies,
p1
time_points,picked_ancestor = 0)
+ ggplot2::geom_vline(xintercept = 0.23, lty = 2) p1
```

And for the other ancestor:

```
<- plot_joyplot_frequencies(selected_pop$frequencies,
p2
time_points,picked_ancestor = 1)
+ ggplot2::geom_vline(xintercept = 0.27, lty = 2) p2
```