scanone.perm: Perform a genome scan.

Description Usage Arguments Value Author(s) See Also Examples

Description

This function is the main QTL mapping function for point mapping at each marker. The user must supply phenotypes, genotype probabilities and marker locations. Optional kinship, additive and interactive covariates may be passed in. scanone.perm regresses the phenotype on the genotype probabilities and produces a LOD score and founder allele effects at each marker.

Usage

1
2
  scanone.perm(pheno, pheno.col = 1, probs, addcovar, intcovar, snps, model = 
  c("additive", "full"), path = ".", nperm = 1000)

Arguments

pheno

data.frame containing phenotype data. Required. rownames must contain sample IDs and there must be a column labelled 'sex' to perform correct mapping on the X chromosome.

pheno.col

numeric or character vector: Either a vector of number that indicate columns to use or a set of column names in pheno.

probs

3D numeric array, containing the founder haplotype contributions or genotype probabilities. The sample IDs, founder letter codes and markers IDs must be in dimnames.

addcovar

data.frame or numeric matrix, containing any additive covariates. Sample IDs must be in rownames.

intcovar

data.frame or numeric matrix, containing any covariates that interact with the QTL. Sample IDs must be in rownames.

snps

data.frame containing 4 columns with marker location information. SNP ID, chr, Mb, cM in columns 1 through 4, respectively.

model

Character string, containing one of "additive" or "full", indicating the type of model to fit.

path

Character string containing the location where permutation results should be written.

nperm

Numeric, containing the number of permutations to run.

Value

A list containing two elements:

lod

Data.frame with seven columns containing the marker locations and LOD scores for each phenotype.

coef

Numeric matrix containing the founder allele effects at each locus.

Author(s)

Daniel Gatti

See Also

scanone.perm

Examples

1
  ## Not run: scanone.perm(pheno, pheno.col = 1, probs)

DOQTL documentation built on May 6, 2019, 3:09 a.m.