bcRobMVNmtQTL {rmtQTL}R Documentation

Multivariate normal (MVN) distribution based multi-trait QTL analysis.

Description

This function calculates the robust LOD (log of odds) scores for multi-trait QTL analysis using robust estimation of the parameters of multivariate normal distribution.

Usage

bcRobMVNmtQTL(inData, phenoCols=1, Chrs, Beta=0)

Arguments

inData

Input data set for multi-trait QTL analysis.

phenoCols

A numeric value of column number or a vector of column numbers in the phenotype matrix which should be used in the analysis. Default value is phenoCols=1

Chrs

Vector of chromosome number(s) for which LOD scores should be calculated.

Beta

Tuning parameter. Default value is 0 for which the estimates reduce to the classical estimates.

Value

It will return a dataframe with the values of chromosome number, marker position (in cM) and robust LOD statistic for multi-trait QTL mapping.

Author(s)

Md. Jahangir Alam <jahangir_statru63@yahoo.com>

Examples

library(qtl)
library(rmtQTL)
data(bcSimData)
multGenoProb <- calc.genoprob(bcSimData, step=1, stepwidth="max")
LODstat<-bcRobMVNmtQTL(inData=multGenoProb, phenoCols=c(1:3), Chrs=c(1:13), Beta=0.001)

[Package rmtQTL version 0.01 Index]