startRun <- 0 endRun <- 639 runIndex="a" lengthAll <- 3201157 shift <- 5000 intervallAll <- 10000 over <- intervallAll%/%shift N1 <- lengthAll%/%shift endRunAll <- (N1-over+1) # 639 source("/seppdata/sepp/linkage/release/funct.R") load(file=paste("resDups.Rda",sep="")) res <- list() allCount <- 0 allCount1 <- 0 for (posAll in (startRun+1):(endRun+1)) { start <- (posAll-1)*shift end <- start + intervallAll if (end > lengthAll) { end <- lengthAll } pRange <- paste("_",format(start,scientific=FALSE),"_",format(end,scientific=FALSE),sep="") load(file=paste("resIBD_chr1",pRange,".Rda",sep="")) noIBD <- length(mergedIBD) if (noIBD > 0) { for (ibdC in 1:noIBD) { allCount <- allCount + 1 if (!dups[allCount]) { allCount1 <- allCount1 + 1 res[[allCount1]] <- c(posAll,ibdC,allCount,allCount1) } } } } rr <- unlist(res) l <- 4 IndexMat <- matrix(rr,nrow=allCount1,ncol=l,byrow=TRUE) colnames(IndexMat) <- c("posAll","ibdC","allCount","allCount1") write.table(IndexMat,file=paste("IndexMat.txt",sep="")) save(IndexMat,file=paste("IndexMat.Rda",sep=""))