diff --git a/src/tp1.R b/src/tp1.R new file mode 100644 index 0000000000000000000000000000000000000000..1230310402f23cd574cf9be80fc7e2f330ac4689 --- /dev/null +++ b/src/tp1.R @@ -0,0 +1,17 @@ +options(repos = c(CRAN = "https://cloud.r-project.org")) + +if (!require("BiocManager", quietly = TRUE)) + install.packages("BiocManager") + +if (!require("snpStats", quietly = TRUE)) + BiocManager::install("snpStats") + +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")) +