QTLsim: Simulate QTL + residual effects

In rramadeu/diaQTLSimulations: Simulate multi-population of polyploids for linkage analysis and QTL

Description

Simulates phenotypic data on top of the pedigreeSimR output and save then as diaQTL csv format

Usage

QTLsim(
  parents = 3,
  ploidy = 4,
  workingfolder = "PedigreeSimR_files",
  QTLmarker = NULL,
  QTLh2 = 0.3,
  run_diaQTL = FALSE
)

Arguments

parents	number of parents in the pedigree
ploidy	integer even number
workingfolder	string with folder name to load the input and write output files
QTLmarker	name of marker to be used as QTL, if NULL a marker is sampled randomly (default)
QTLh2	numeric between 0 and 1 with the QTL heritability
run_diaQTL	logical, if TRUE runs diaQTL scan1 and fitQTL functions in the data.

Value

if run_diaQTL=TRUE, return a list with diaQTL results. Otherwise, returns NULL.

Author(s)

Rodrigo R Amadeu, [e-mail]

Examples

## Not run: 
map = 1:100
haplotypes = fake_haplo(m=50,n=100,seed=1234)
pedigree = diallel_pedigree(parents=3,popsize=100)
pedigreesimR(map,haplotypes,pedigree,workingfolder="PedigreeSimR_files")
QTLsim(parents=3,
       ploidy=4,
       workingfolder="PedigreeSimR_files",
       QTLmarker=NULL,
       QTLh2=0,
       run_diaQTL=FALSE)

## End(Not run)

rramadeu/diaQTLSimulations documentation built on Feb. 11, 2022, 1:34 a.m.