diff --git a/HPC_analysis/output/Popgen/angsd/README.md b/HPC_analysis/output/Popgen/angsd/README.md index d79ee22..f550efa 100644 --- a/HPC_analysis/output/Popgen/angsd/README.md +++ b/HPC_analysis/output/Popgen/angsd/README.md @@ -267,15 +267,6 @@ vcftools --vcf out.5.recode.vcf --out out.5 --missing-indv library(readr) library(tidyverse) -out_5_imiss <- read_tsv("out.5.imiss") %>% - extract(INDV, "pop", "(\\w+)_", remove = FALSE) - -out_5_imiss %>% count(pop) -``` - -# Load the data in Rstudio -# open the out.5.imiss file in R -# note the read_tsv is in the package readr out_5_imiss <- read_tsv("out.5.imiss") %>% extract(INDV, "pop", "(\\w+)_", remove = FALSE)