library(tidyverse)


fst = choose.files()



fst1 = read.table(fst, header = T)
head(fst1)
fst2 = fst1 %>% select(1,2,3,6) %>% top_frac(0.05,wt = MEAN_FST)
head(fst2)
write.table(fst2,"fst_vasuclar.txt",quote = F,row.names = F,sep = "\t")

indica_pi = choose.files()
pi1 = read.table(indica_pi, header = T)
head(pi1)
pi2 = pi1 %>% select(1,2,3,5) %>% slice_min(PI,prop = 0.05)
head(pi2)
write.table(pi2,"indica_pi2_vasuclar.txt",quote = F,row.names = F,sep = "\t")


japonica_pi = choose.files()
pi1 = read.table(japonica_pi, header = T)
head(pi1)
pi2 = pi1 %>% select(1,2,3,5) %>% slice_min(PI,prop = 0.05)
head(pi2)
write.table(pi2,"japonica_pi2_vasuclar.txt",quote = F,row.names = F,sep = "\t")