aboutsummaryrefslogtreecommitdiff
path: root/scripts/ctl_analysis.R
diff options
context:
space:
mode:
authorAlexander Kabui2021-10-13 13:25:45 +0300
committerBonfaceKilz2022-01-22 09:23:14 +0300
commitfa4bc9a68b37bceaaacf6d00e69413ba829b3e77 (patch)
tree5b233c5530aaa00f73634c0d642d43d4d4952b28 /scripts/ctl_analysis.R
parentdbb88497846253a23e028fb627f27cec82047be3 (diff)
downloadgenenetwork3-fa4bc9a68b37bceaaacf6d00e69413ba829b3e77.tar.gz
perform the ctl scan
Diffstat (limited to 'scripts/ctl_analysis.R')
-rw-r--r--scripts/ctl_analysis.R10
1 files changed, 10 insertions, 0 deletions
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
+
+