R/dlrs.R

In seandavi/Rpressa: A collection of useful code for bioinformatics

dlrs <-
 function(x) {
   nx <- length(x)
   if (nx<3) {
     stop("Vector length>2 needed for computation")
   }
   tmp <- embed(x,2)
   diffs <- tmp[,2]-tmp[,1]
   dlrs <- IQR(diffs,na.rm=TRUE)/(sqrt(2)*1.34)
   return(dlrs)
 }