aboutsummaryrefslogtreecommitdiff
path: root/scripts/ctl_analysis.R
blob: 04b55f99ab2e253230d7c3e79b38d9105729f519 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
library(ctl)

# The genotypes.csv file containing the genotype matrix is stored individuals (rows) x genetic marker (columns):
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