aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-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
+
+