From fa4bc9a68b37bceaaacf6d00e69413ba829b3e77 Mon Sep 17 00:00:00 2001
From: Alexander Kabui
Date: Wed, 13 Oct 2021 13:25:45 +0300
Subject: perform the ctl scan

---
 scripts/ctl_analysis.R | 10 ++++++++++
 1 file changed, 10 insertions(+)

(limited to 'scripts/ctl_analysis.R')

diff --git a/scripts/ctl_analysis.R b/scripts/ctl_analysis.R
index bfb6da4..04b55f9 100644
--- a/scripts/ctl_analysis.R
+++ b/scripts/ctl_analysis.R
@@ -4,3 +4,13 @@ library(ctl)
 genotypes <- read.csv("genotypes.csv",row.names=1, header=FALSE, sep="\t")
 # The phenotypes.csv file containing individuals (rows) x traits (columns) measurements:
 traits <- read.csv("phenotypes.csv",row.names=1, header=FALSE, sep="\t")
+
+
+ctls <- CTLscan(geno,traits,strategy=input$strategy,
+	nperm=input$nperms,parametric =input$parametric,
+	nthreads=6,verbose=TRUE)
+
+
+# plots and ctl networks
+
+
-- 
cgit v1.2.3