diff options
| author | Alexander_Kabui | 2024-11-11 12:07:45 +0300 |
|---|---|---|
| committer | Alexander_Kabui | 2024-11-11 12:07:45 +0300 |
| commit | 006f3f73c20431754b8b6e5c2fc50202b2b7d926 (patch) | |
| tree | a6d6e71c1c8290e4bfa86ad71bd39c14eabf8032 | |
| parent | 0a5c1f7d4ccc88d3555e44b50beec8ee91bc43fc (diff) | |
| download | genenetwork3-006f3f73c20431754b8b6e5c2fc50202b2b7d926.tar.gz | |
refactor: adding pseudomarkers o genetic probabilites map.
| -rw-r--r-- | scripts/rqtl2_wrapper.R | 11 |
1 files changed, 4 insertions, 7 deletions
diff --git a/scripts/rqtl2_wrapper.R b/scripts/rqtl2_wrapper.R index 6581ed4..d025935 100644 --- a/scripts/rqtl2_wrapper.R +++ b/scripts/rqtl2_wrapper.R @@ -178,16 +178,13 @@ perform_genetic_pr <- function(cross, cat("Finding the genetic Probabilities\n") if (use_pseudomarkers){ + cat("Using pseudo markers for genetic probabilites\n") map <- insert_pseudomarkers(cross$gmap, step=step) + } return(calc_genoprob(cross, map=map, error_prob=error_prob, map_function=map_function, - quiet=FALSE, cores=cores)) - } - else { - return (calc_genoprob(cross, map=map, error_prob=error_prob, - quiet = FALSE, map_function =map_function, - cores = cores)) - }} + quiet=FALSE, cores=cores)) +} Pr = perform_genetic_pr(dataset) cat("Summaries on the genetic probabilites \n") |
