R/get.better.R

In rxSeq: Combined Total and Allele Specific Reads Sequencing Study

`get.better` = 
function(lst, ind, trecase=TRUE){
  lls = Inf
  theta_i = iphi_i = inits = NULL
  for(item in lst){
    if(!is.null(item)){
      if(lls>item[ind]){
        lls = item[ind]
        inits = item[1:(ind-1)]
        iphi_i = item[ind + 1]
        if(trecase)theta_i = item[ind + 2]
      }
    }
  }
  if(lls<Inf){
    if(trecase){
      return(c(inits, lls, iphi_i, theta_i))  
    }else{
      return(c(inits, lls, iphi_i))  
    }
  }else{
    return(NULL)
  }
}