SIMQTL: QTL detection by SIM

SIMQTLR Documentation

QTL detection by SIM

Description

Perform QTL detection by simple interval mapping

Usage

SIMQTL(
  cross,
  numeric.chr.format = TRUE,
  response.in.cross = TRUE,
  pheno.col = "y",
  response = NULL,
  method = "em",
  geno.joinmap = NULL,
  phase,
  threshold = NA,
  nperm = 100,
  alpha = c(0.05, NA),
  plot = FALSE,
  QTL_position = NULL,
  verbose = 0
)

Arguments

cross

object

numeric.chr.format

logical to indicate if chromosome names are numeric

response.in.cross

logical to indicate if response studied is in cross$pheno, default is TRUE

pheno.col

character indicating column to study in cross$pheno

response

named numeric or vector for response if not in cross$pheno, default is NULL

method

method to detect QTL in qtl::scanone, default is "em"

geno.joinmap

genotypes at markers in the JoinMap format, if NULL (default), no estimation of allelic effects is given

phase

marker phases

threshold

genomewide significance LOD threshold, if NA (default), is found by permutations (with nperm parameter).

nperm

number of permutations to be done in qtl::scanone, default is 100

alpha

vector of length 1 or 2 (optional) with thresholds for (1) the significance of QTL presence (based on permutations) and (2) the significance of linear regression on QTL effect, if NA, no threshold is applied and all allelic effects are kept.

plot

logical, default is FALSE.

QTL_position

matrix with genetic.distance and linkage.group as columns indicating QTL positions for plotting

verbose

verbosity level (0/1/2)

Value

list of 3 elements: qtl.df is a data frame with QTL informations (linkage.group, position, LOD, interval.inf and interval.sup) / selected markers is a character vector for markers inside confidence interval /allelic effetcs is a data frame with a column predictor and a column effect with estimated allelic effects.

Author(s)

Charlotte Brault [aut], Timothee Flutre [ctb]

See Also

MIMQTL


timflutre/rutilstimflutre documentation built on Feb. 7, 2024, 8:17 a.m.