allCount <- 0 for (posAll in startRun:endRun) { start <- (posAll-1)*shiftSize end <- start + segmentSize if (end > snps) { end <- snps } pRange <- paste("_",format(start,scientific=FALSE),"_",format(end,scientific=FALSE),sep="") load(file=paste(fileName,pRange,"_resAnno.Rda",sep="")) mergedIBD <- resHapFabia$mergedIBD noIBD <- length(mergedIBD) allCount <- allCount+noIBD } allCount