Skip to content
Snippets Groups Projects
Commit 36148f0f authored by OTT Oceane's avatar OTT Oceane
Browse files

script

parent 2eaa0572
Branches
Tags
No related merge requests found
wdir="."
dir.create(wdir, showWarnings = F, recursive = T)
setwd(wdir)
#library(devtools)
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"))
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment