inde <-  c("a","b","c","d","e","f","g")


AncestMatA <- c()
for (runIndex in inde)   {

load(file=paste("AncestMat",runIndex,".Rda",sep=""))
AncestMatA <- rbind(AncestMatA,AncestMat)

}

AncestMat <- AncestMatA
colnames(AncestMat) <- c("allCount1","allCount","ibdC","posAll","#SNPs","#Nea","#Ref","#Alt","#New","NeaR","RefR","AltR","NewR","xAltCon","xRefCon","mAltCon","mRefCon","NeaPosC","AltPosC","RefPosC")

save(AncestMat,file=paste("AncestMat.Rda",sep=""))



dim(AncestMat)

summary(AncestMat)


oo <- which(AncestMat[,8]==1.0)


perf <- AncestMat[oo,1:4]

perf

oo1 <- which(AncestMat[,4]>40)

perf1 <- AncestMat[oo1,1:4]

perf1
AncestMat[oo1,]


inde <-  c("a","b","c","d","e","f","g")

sum <- 0
for (runIndex in inde)   {

load(file=paste("Mat",runIndex,".Rda",sep=""))
print(dim(AncestMat))
sum <- sum + dim(AncestMat)[1]

}
sum

AncestMat[oo1,]

a <- c()
for (i in 1:7)   {

a <- c(a,length(which(PopulationMat[,4]<(i*100+1))))
}


a
> a
[1]  27723  52974  77365 101629 125134 149446 160588
diff(a)
> diff(a)
[1] 25251 24391 24264 23505 24312 11142

27723 25251 24391 24264 23505 24312 11142