tp1

src/tp1.R

@@ -2,7 +2,7 @@ output_dir = "results/tp1"
 dir.create(output_dir, showWarnings = F, recursive = T)
 
 #if (!require("BiocManager", quietly = TRUE))
-  install.packages("BiocManager")
+ # install.packages("BiocManager")
 
 #if (!require("snpStats", quietly = TRUE))
  # BiocManager::install("snpStats")
@@ -28,17 +28,16 @@ write.plink(plink_base, snps=geno$genotypes, pedigree=geno$fam[,1], id=geno$fam[
 
 genoBim<-geno$map
 colnames(genoBim)<-c("chr", "SNP", "gen.dist", "position", "A1", "A2")
-head(genoBim)
+#head(genoBim)
 
 genotype<-geno$genotype
-print(genotype)
+#print(genotype)
 
 genoFam<-geno$fam
-head(genoFam)
+#head(genoFam)
 
 # On commence par libérer de l'espace
 rm(geno)
 
-
 rdata_path = file.path(output_dir, "TP1_asbvg.RData")
 save.image(rdata_path)