bcRobMVNmtQTL {rmtQTL} | R Documentation |
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.
bcRobMVNmtQTL(inData, phenoCols=1, Chrs, Beta=0)
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. |
It will return a dataframe with the values of chromosome number, marker position (in cM) and robust LOD statistic for multi-trait QTL mapping.
Md. Jahangir Alam <jahangir_statru63@yahoo.com>
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)