diff --git a/R/write_rqtl.R b/R/write_rqtl.R index 0926c5c3a44d009a46244ee5d52835afb1309978..ff77d827447e2bc92cc0f7a60c317db946a7e179 100755 --- a/R/write_rqtl.R +++ b/R/write_rqtl.R @@ -59,6 +59,9 @@ write_rqtl <- function(geno,pheno,tab,ref,par1,par2,prefix,pos,path=NA){ #merge genotypes with parents geno <- left_join(geno,ref,by=c("marker"="marker")) + #remove snps with no position + geno <- geno %>% filter(is.na(chr)==FALSE) + #recode "-" in "N" in geno file geno <- geno %>% mutate(allele_1 = recode(allele_1, "-" = "N"))