Learn R Programming

DataFusionGDM (version 1.3.2)

simulate_genetic_distances: Simulate genetic distances using realistic population structure

Description

Generates a synthetic genetic distance matrix and metadata using hierarchical population structure, admixture and bottleneck options.

Usage

simulate_genetic_distances(
  n_pops = 50,
  n_major_groups = 5,
  n_subgroups = 12,
  geo_dims = 2,
  genetic_dims = 2,
  group_separation = 15,
  subgroup_separation = 5,
  pop_dispersion = 0.5,
  isolation_factor = 0.8,
  admixture_prob = 0.1,
  bottleneck_prob = 0.05,
  noise_level = 0.1,
  nonlinear_factor = 0.7,
  use_subgroups = TRUE,
  use_genetic_dims = TRUE,
  use_admixture = TRUE,
  use_bottlenecks = TRUE,
  use_isolation_by_distance = TRUE,
  use_nonlinear = TRUE,
  use_noise = TRUE,
  seed = NULL,
  verbose = TRUE
)

Value

A list with distance_matrix, population_info, position_matrix, and parameters.

Arguments

n_pops

Number of populations

n_major_groups

Number of major groups

n_subgroups

Number of subgroups

geo_dims

Geographic dimensions

genetic_dims

Additional genetic drift dimensions

group_separation

Separation between major groups

subgroup_separation

Separation between subgroups

pop_dispersion

Within-subgroup dispersion

isolation_factor

Weight for geography in isolation-by-distance model (0-1)

admixture_prob

Proportion of admixed populations

bottleneck_prob

Proportion of bottlenecked populations

noise_level

Noise level in transformation

nonlinear_factor

Nonlinearity factor in transformation

use_subgroups

Whether to create subgroups

use_genetic_dims

Whether to include genetic dimensions

use_admixture

Whether to include admixture

use_bottlenecks

Whether to include bottlenecks

use_isolation_by_distance

Whether to weight geographic distance

use_nonlinear

Whether to apply nonlinear transformation

use_noise

Whether to add noise

seed

Optional seed for reproducibility (NULL leaves the RNG state unchanged)

verbose

Print diagnostics