From eab69afb5a4685cd448a5b361c3294065cd6c3bc Mon Sep 17 00:00:00 2001 From: Kham Chaimae <khamchaimae@MacBook-Pro-de-Kham.local> Date: Wed, 23 Oct 2024 10:54:40 +0200 Subject: [PATCH] commit tp1.R --- tp1.R | 12 ++++++++++-- 1 file changed, 10 insertions(+), 2 deletions(-) diff --git a/tp1.R b/tp1.R index 3a66557..5cc42e9 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")) + + + -- GitLab