diff --git a/tp1.R b/tp1.R
index 3a665572add184254bbe1f60a8308254813e6ffd..5cc42e9223874b31cf894c9869c6f690760591da 100644
--- a/tp1.R
+++ b/tp1.R
@@ -1,7 +1,7 @@
 wdir="."
-dir.create(wdir, showWarnings = F, recursive = T)
+dir.create(<dir, showWarnings = F, recursive = T)
 setwd(wdir)
-  
+
 #library(devtools)
 
 if (!require("BiocManager", quietly = TRUE))
@@ -13,3 +13,11 @@ if (!require("snpStats", quietly = TRUE))
 if (!require("SNPRelate", quietly = TRUE))
   BiocManager::install("SNPRelate")
 
+# Les données analysées nécessitant beaucoup de RAM, nous allons sélectionner aléatoirement 250000 SNPs et réecrire des fichiers bed, bim, fam
+penncath_bed_path = "results/data/penncath.bed"
+penncath_bim_path = "results/data/penncath.bim"
+penncath_fam_path = "results/data/penncath.fam"
+geno <- snpStats::read.plink(penncath_bed_path, penncath_bim_path, penncath_fam_path, select.snps=sample(1:861473, 25000, replace = FALSE ), na.strings = ("-9"))
+
+
+