bcMVNmtQTL {rmtQTL} | R Documentation |
This function calculates the LOD (log of odds) for multi-trait QTL analysis using multivariate normal distribution.
bcMVNmtQTL(inData, phenoCols=1, Chrs)
inData |
Input data set along with genotype probabilities for multi-trait QTL analysis. |
phenoCols |
Column numbers in the phenotype matrix which should be used in the analysis. |
Chrs |
Vector of chromosome number(s) for which LOD scores should be calculated. |
It will return a dataframe with the values of chromosome number, marker position (in cM) and 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<-bcMVNmtQTL(inData=multGenoProb, phenoCols=c(1:3), Chrs=c(1:13))