Skip to contents

Retrieves QTL haplotype data

Usage

pullQtlHaplo(
  pop,
  trait = 1,
  haplo = "all",
  chr = NULL,
  asRaw = FALSE,
  simParam = NULL
)

Arguments

pop

an object of Pop-class

trait

an integer. Indicates which trait's QTL haplotypes to retrieve.

haplo

either "all" for all haplotypes or an integer for a single set of haplotypes. Use a value of 1 for female haplotypes and a value of 2 for male haplotypes in diploids.

chr

a vector of chromosomes to retrieve. If NULL, all chromosome are retrieved.

asRaw

return in raw (byte) format

simParam

an object of SimParam

Value

Returns a matrix of QTL haplotypes.

Examples

#Create founder haplotypes
founderPop = quickHaplo(nInd=10, nChr=1, segSites=15)

#Set simulation parameters
SP = SimParam$new(founderPop)
SP$addTraitA(10)
SP$addSnpChip(5)

#Create population
pop = newPop(founderPop, simParam=SP)
pullQtlHaplo(pop, simParam=SP)
#>      1_1 1_3 1_7 1_8 1_10 1_11 1_12 1_13 1_14 1_15
#> 1_1    1   1   0   0    0    0    1    0    1    0
#> 1_2    1   1   1   1    1    1    0    0    1    1
#> 2_1    0   1   1   1    0    0    0    0    0    0
#> 2_2    0   0   1   0    0    0    1    1    1    1
#> 3_1    0   0   0   0    0    0    1    1    0    1
#> 3_2    0   1   0   0    0    1    0    0    0    1
#> 4_1    1   1   0   0    1    1    0    0    0    1
#> 4_2    0   0   0   0    0    0    1    1    1    0
#> 5_1    0   1   0   0    1    0    1    1    1    1
#> 5_2    0   0   1   0    0    1    1    1    0    0
#> 6_1    0   1   0   0    0    1    0    0    0    1
#> 6_2    1   1   0   0    0    0    1    1    0    0
#> 7_1    1   1   0   0    1    1    0    1    1    1
#> 7_2    1   0   1   0    1    1    0    0    1    0
#> 8_1    0   1   0   1    1    1    0    1    1    0
#> 8_2    1   0   1   1    0    0    0    0    0    1
#> 9_1    1   0   0   0    0    1    1    0    1    1
#> 9_2    1   0   1   1    0    1    1    1    0    1
#> 10_1   0   1   0   0    1    0    1    0    1    1
#> 10_2   1   0   1   0    1    0    0    0    1    1