R/label2jaq.R

In IBDLabels: Convert Between Different IBD-State Labelling Schemes

Defines functions label2jaq

Documented in label2jaq

label2jaq <-
function( label, phased = TRUE){
   
    out <- allJaq()

    jaq <- numeric( length( label )) 
    sel <- ifelse( phased, 2,3 )
    for( i in 1:length(label)){
        jaq[i] <- out[ out[,1]== label[i], sel]
        }
        
    return( jaq ) 
}

Try the IBDLabels package in your browser