devAskNewPage(ask = TRUE)

for (i in 1:length(mergedIBD)) {

sampEU <- mergedIBD[[i]]$noSamples        
ibb <- mergedIBD[[i]]$noSnps 
ibb <- as.integer(sort.int(as.integer(unique(ibb))))


samp <-  EUeurop[sampEU]


labelsEU <- indi[samp,2]
labelsNA <- indi[samp,1]
labels_ALL <- paste(labelsNA,labelsEU,sep="_")
  
Lout <- readSamplesSpfabia(X=paste("ALL.chr1.merged_beagle_mach.20101123.snps_indels_svs.genotypes",pRange,"_mat",sep=""),samples=samp,lowerB=0,upperB=1000.0)


plotRange(Lout=Lout,ibb=ibb,annot=annot,colRamp=12,val=c(0.0,2.0,1.0),chrom=1,labelsNA=labels_ALL)

}

devAskNewPage(ask = FALSE)