R/mqmprepare.R

In kbroman/qtl: Tools for Analyzing QTL Experiments

#####################################################################
#
# mqmprepare.R
#
# Copyright (c) 2009-2010, Danny Arends
#
# Modified by Pjotr Prins and Karl Broman
#
#
# first written Februari 2009
# last modified April 2010
#
#     This program is free software; you can redistribute it and/or
#     modify it under the terms of the GNU General Public License,
#     version 3, as published by the Free Software Foundation.
#
#     This program is distributed in the hope that it will be useful,
#     but without any warranty; without even the implied warranty of
#     merchantability or fitness for a particular purpose.  See the GNU
#     General Public License, version 3, for more details.
#
#     A copy of the GNU General Public License, version 3, is available
#     at http://www.r-project.org/Licenses/GPL-3
#
# Part of the R/qtl package
# Contains: mqmfind.marker
#
#
#####################################################################

mqmfind.marker <- function(cross,mqmscan=NULL,perm=NULL,alpha=0.05,verbose=FALSE){
    thesum <- summary(mqmscan[,1:3],alpha=alpha,perms=perm,pvalues=FALSE)
    chr <- thesum$'chr'
    pos <- thesum$'pos'
    if(verbose) cat("INFO: Found",length(chr),"markers with alpha <",alpha,".\n")
    ret <- NULL
    for(i in 1:length(chr)){
        ret <- rbind(ret,cbind(find.marker(cross,chr=chr[i],pos=pos[i]),as.integer(chr[i]),as.double(pos[i])))
    }
    colnames(ret) <- c("marker","chr","pos (cM)")
    ret
}



# end of mqmprepare.R