tp3.R

#load data
load("results/tp2/TP2_asbvg.RData")

#fichier de sortie
output_dir = "results/tp3"
dir.create(output_dir, showWarnings = F, recursive = T)

#library
if (!require("BiocManager", quietly = TRUE))
  install.packages("BiocManager")

if (!require("snpStats", quietly = TRUE))
  BiocManager::install("snpStats")

#if (!require("doParallel", quietly = TRUE))
  #BiocManager::install("doParallel")

library(snpStats)
library(plyr)
library(GenABEL)
#library(doParallel)

rdata_path = file.path(output_dir,"TP3_asbvg.RData")
save.image(rdata_path)