sim_qtl_pheno: 'Simulate QTL phenotype with normal distribution'

In jyanglab/g3tools: For Genomics, quantGen, and popGen studies, especially for crop species.

Description

Simulate phenotype using genotype data given number of QTLs and heritability.

Usage

sim_qtl_pheno(geno, h2, alpha, nqtl, distribution = "norm")

Arguments

geno	genotype data, col1=genoid, col2 and after=snpid, coding: 0,1,2 (no missing data allowed) [data.frame].
h2	Broad sense heritability of the trait [numeric(0-1)].
alpha	Alpha.
nqtl	number of QTL [interger].
distribution	[character=norm]

Value

return A list of many values.

Examples

geno <- read.table("data/geno.txt", header=TRUE)
pheno <- sim_qtl_pheno(geno, h2=0.7, alpha=0.5, nqtl=10, distribution="norm")
y <- pheno[['y']]

jyanglab/g3tools documentation built on May 20, 2019, 6:27 a.m.