Extract list with tree sequence metadata saved by SLiM
Examples
check_dependencies(python = TRUE) # make sure dependencies are present
init_env()
#> The interface to all required Python modules has been activated.
# load an example model with an already simulated tree sequence
slendr_ts <- system.file("extdata/models/introgression.trees", package = "slendr")
model <- read_model(path = system.file("extdata/models/introgression", package = "slendr"))
# load the tree-sequence object from disk
ts <- ts_load(slendr_ts, model)
# extract the list of metadata information from the tree sequence
ts_metadata(ts)
#> $version
#> [1] "slendr_0.4.0.9000"
#>
#> $description
#> [1] ""
#>
#> $sampling
#> # A tibble: 13 × 3
#> name time pop
#> <chr> <int> <chr>
#> 1 NEA_1 70000 NEA
#> 2 NEA_2 40000 NEA
#> 3 AFR_1 0 AFR
#> 4 AFR_2 0 AFR
#> 5 AFR_3 0 AFR
#> 6 AFR_4 0 AFR
#> 7 AFR_5 0 AFR
#> 8 CH_1 0 CH
#> 9 EUR_1 0 EUR
#> 10 EUR_2 0 EUR
#> 11 EUR_3 0 EUR
#> 12 EUR_4 0 EUR
#> 13 EUR_5 0 EUR
#>
#> $map
#> NULL
#>
#> $arguments
#> $arguments$BURNIN_LENGTH
#> [1] 0
#>
#> $arguments$MAX_ATTEMPTS
#> [1] 1
#>
#> $arguments$RECOMB_RATE
#> [1] 1e-08
#>
#> $arguments$SEED
#> [1] 314159
#>
#> $arguments$SEQUENCE_LENGTH
#> [1] 500000
#>
#> $arguments$SIMULATION_LENGTH
#> [1] 216667
#>
#>