diff --git a/src/tp1.R b/src/tp1.R index 2c57e92c5cd9f79842898e7a8e0c93215874da3f..904a363a5343d3f96bc2a0184070e65963a5eae1 100644 --- a/src/tp1.R +++ b/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)