inde <- c("a","b","c","d","e","f","g") DenisovaMatA <- c() for (runIndex in inde) { load(file=paste("DenisovaMat",runIndex,".Rda",sep="")) DenisovaMatA <- rbind(DenisovaMatA,DenisovaMat) } DenisovaMat <- DenisovaMatA colnames(DenisovaMat) <- c("allCount1","allCount","ibdC","posAll","#SNPs","#Nea","#Ref","#Alt","#New","NeaR","RefR","AltR","NewR","xAltCon","xRefCon","mAltCon","mRefCon","NeaPosC","AltPosC","RefPosC") save(DenisovaMat,file=paste("DenisovaMat.Rda",sep="")) dim(DenisovaMat) summary(DenisovaMat) oo <- which(DenisovaMat[,8]==1.0) perf <- DenisovaMat[oo,1:4] perf oo1 <- which(DenisovaMat[,4]>40) perf1 <- DenisovaMat[oo1,1:4] perf1 DenisovaMat[oo1,]