diff options
| author | Alexander Kabui | 2021-10-13 13:25:45 +0300 |
|---|---|---|
| committer | BonfaceKilz | 2022-01-22 09:23:14 +0300 |
| commit | fa4bc9a68b37bceaaacf6d00e69413ba829b3e77 (patch) | |
| tree | 5b233c5530aaa00f73634c0d642d43d4d4952b28 /scripts/ctl_analysis.R | |
| parent | dbb88497846253a23e028fb627f27cec82047be3 (diff) | |
| download | genenetwork3-fa4bc9a68b37bceaaacf6d00e69413ba829b3e77.tar.gz | |
perform the ctl scan
Diffstat (limited to 'scripts/ctl_analysis.R')
| -rw-r--r-- | scripts/ctl_analysis.R | 10 |
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 + + |