about summary refs log tree commit diff
path: root/scripts
diff options
context:
space:
mode:
authorAlexander_Kabui2025-02-07 14:17:17 +0300
committerAlexander_Kabui2025-02-07 14:17:17 +0300
commitfebd6ebebfba2b585d98c3e5e35c6713c9781036 (patch)
treeea14a276c555f034967ab81328783d9508d4311c /scripts
parent8a6c3966e63dd2e381892cd0f73deee88ccd2060 (diff)
downloadgenenetwork3-febd6ebebfba2b585d98c3e5e35c6713c9781036.tar.gz
fix: pass the correct genome probabilites.
Diffstat (limited to 'scripts')
-rw-r--r--scripts/rqtl2_wrapper.R4
1 files changed, 2 insertions, 2 deletions
diff --git a/scripts/rqtl2_wrapper.R b/scripts/rqtl2_wrapper.R
index 5ad812e..7cad7d5 100644
--- a/scripts/rqtl2_wrapper.R
+++ b/scripts/rqtl2_wrapper.R
@@ -226,9 +226,9 @@ print(Xcovar)
 #  Function to calculate the kinship 
 get_kinship <- function(probability, method="LMM"){
 if (opt$method == "LMM"){
-    kinship = calc_kinship(genome_prob)
+    kinship = calc_kinship(probability)
 } else if (opt$method == "LOCO"){
-    kinship = calc_kinship(genome_prob, "loco")
+    kinship = calc_kinship(probability, "loco")
 }else {
  kinship = NULL
 }