load("dups.Rda") load("analyzeResult.Rda") library(hapFabia) library(fabia) #source("prepareAncientGenome.R") #source("plotAncientGenomes.R") shiftSize <- 5000 segmentSize <- 10000 fileName="KPGP1000_chr1" # without type posAll <- 442 ibdC <- 51 posAll <- 442 ibdC <- 52 posAll <- 443 ibdC <- 73 start <- (posAll-1)*shiftSize end <- start + segmentSize pRange <- paste("_",format(start,scientific=FALSE),"_",format(end,scientific=FALSE),sep="") load(file=paste(fileName,pRange,"_resAnno",".Rda",sep="")) IBD <- resHapFabia$mergedIBD plotAncientGenomes(IBD,filename=paste(fileName,pRange,"_mat",sep=""),ibdC) name="korDenis1b" pdf(paste(name,".pdf",sep=""),width=12) plotAncientGenomes(IBD,filename=paste(fileName,pRange,"_mat",sep=""),ibdC) dev.off() xfig(paste(name,".fig",sep=""),width=12) plotAncientGenomes(IBD,filename=paste(fileName,pRange,"_mat",sep=""),ibdC) dev.off()