table
rbinom
genetic_diff
d()
Published with bookdown & the OTTR template)
Style adapted from: rstudio4edu-book (CC-BY 2.0)
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