Learn R Programming

mappoly (version 0.2.1)

calc_homoprob: Homolog probabilities

Description

Compute homolog probabilities for all individuals in the full-sib population given a map and conditional genotype probabilities.

Usage

calc_homoprob(input.genoprobs, verbose = TRUE)

Arguments

input.genoprobs

an object of class mappoly.genoprob

verbose

if TRUE (default), the current progress is shown; if FALSE, no output is produced

References

Mollinari M., Olukolu B. A., Pereira G. da S., Khan A., Gemenet D., Yencho G. C., Zeng Z-B. (2020), Unraveling the Hexaploid Sweetpotato Inheritance Using Ultra-Dense Multilocus Mapping, _G3: Genes, Genomes, Genetics_. https://doi.org/10.1534/g3.119.400620

Examples

Run this code
# NOT RUN {
   
# }
# NOT RUN {
     ## tetraploid example
     w1 <- calc_genoprob(solcap.dose.map[[1]])
     h.prob <- calc_homoprob(w1)
     print(h.prob)
     plot(h.prob, ind = 5, use.plotly = FALSE)
     ## using error modeling (removing noise)
     w2 <- calc_genoprob_error(solcap.err.map[[1]])
     h.prob2 <- calc_homoprob(w2)
     print(h.prob2)
     plot(h.prob2, ind = 5, use.plotly = FALSE)
  
# }
# NOT RUN {
# }

Run the code above in your browser using DataLab