6.6 Metadata

We’ll also read in a metadata file, which tells us which population each individual is from.

metadata <- read.table("integrated_call_samples.txt",
                       header = TRUE)

head(metadata)
##    sample pop superpop    sex
## 1 HG00096 GBR      EUR   male
## 2 HG00097 GBR      EUR female
## 3 HG00099 GBR      EUR female
## 4 HG00100 GBR      EUR female
## 5 HG00101 GBR      EUR   male
## 6 HG00102 GBR      EUR female