about summary refs log tree commit diff
path: root/scripts/rqtl2_wrapper.R
diff options
context:
space:
mode:
authorAlexander_Kabui2025-02-21 16:47:25 +0300
committerAlexander_Kabui2025-02-21 16:48:01 +0300
commitcd8e9c46341e878907e99dbc8baea2e1cf089d51 (patch)
tree2dac39005aa0c78b2a845b5670a1d2208769d8fd /scripts/rqtl2_wrapper.R
parente21226d606eb09a5f528524bde9339b947cdbf16 (diff)
downloadgenenetwork3-cd8e9c46341e878907e99dbc8baea2e1cf089d51.tar.gz
feat: Temporary disable performing genome scan with covariates for RQTL2.
Diffstat (limited to 'scripts/rqtl2_wrapper.R')
-rw-r--r--scripts/rqtl2_wrapper.R13
1 files changed, 7 insertions, 6 deletions
diff --git a/scripts/rqtl2_wrapper.R b/scripts/rqtl2_wrapper.R
index 5e4622c..59467c6 100644
--- a/scripts/rqtl2_wrapper.R
+++ b/scripts/rqtl2_wrapper.R
@@ -268,13 +268,14 @@ perform_genome_scan <- function(cross,
   if (method == "LMM") {
     # provide parameters for this
     cat("Performing scan1 using Linear mixed model\n")
+  
     out  <- scan1(
       genome_prob,
       cross$pheno,
       kinship = kinship,
-      addcovar = covar,
-      Xcovar = Xcovar,
-      intcovar = intcovar,
+      # addcovar = covar,
+      # Xcovar = Xcovar,
+      # intcovar = intcovar,
       model = model,
       cores = NO_OF_CORES
     )
@@ -284,10 +285,10 @@ perform_genome_scan <- function(cross,
       genome_prob,
       cross$pheno,
       kinship = kinship,
-      addcovar = covar,
-      intcovar = intcovar,
+      # addcovar = covar,
+      # intcovar = intcovar,
       model = model,
-      Xcovar = Xcovar,
+      # Xcovar = Xcovar,
       cores = NO_OF_CORES
     )
   }