getIncMatrix {epinetr}

Incidence matrix retrieval.

Description

Get an incidence matrix from a Population.

Usage

getIncMatrix(pop)

Arguments

pop

An object of class 'Population' which has an EpiNet object attached.

Details

getIncMatrix() retrieves the incidence matrix used in epistatic interactions within the given Population object. This is most useful for copying the network structure to a new Population object.

Value

An incidence matrix representing the epistatic network within the given Population object.

Author(s)

Dion Detterer, Paul Kwan, Cedric Gondro

See Also

attachEpiNet

Examples

# Create population
pop <- Population(
  popSize = 200, map = map100snp, QTL = 20,
  alleleFrequencies = runif(100),
  broadH2 = 0.9, narrowh2 = 0, traitVar = 40
)

# Attach random epistatic network and retrieve incidence matrix
pop <- attachEpiNet(pop)
inc <- getIncMatrix(pop)

# Create second population
pop2 <- Population(
  popSize = 200, map = map100snp, QTL = 20,
  alleleFrequencies = runif(100),
  broadH2 = 0.8, narrowh2 = 0.6, traitVar = 40
)

# Attach epistatic network to second population
# using incidence matrix from first
pop2 <- attachEpiNet(pop2, incmat = inc)