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)