newNuclearFamily: Simulates genotypes for 1 family with 1 offspring

Description Usage Arguments Value Examples

Description

Simulates genotypes for 1 family with 1 offspring

Usage

1
newNuclearFamily(family_id)

Arguments

family_id

What will be the family_id (for example: 100)

Value

family structure object

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
library("sim1000G")

examples_dir = system.file("examples", package = "sim1000G")
vcf_file = file.path(examples_dir, "region.vcf.gz")
vcf = readVCF( vcf_file, maxNumberOfVariants = 100 ,
   min_maf = 0.12 ,max_maf = NA)

genetic_map_of_region = system.file("examples",
   "chr4-geneticmap.txt",
   package = "sim1000G")
readGeneticMapFromFile(genetic_map_of_region)

startSimulation(vcf, totalNumberOfIndividuals = 1200)
fam1 = newNuclearFamily(1)
fam2 = newNuclearFamily(2)

# See also the documentation on our github page

sim1000G documentation built on June 10, 2019, 1:01 a.m.