Advertisement
Light1992

Indice fuori range

Sep 30th, 2019
120
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.64 KB | None | 0 0
  1. library(GOplot)
  2. EC$david <- read.delim("EC.david", header=TRUE)
  3. EC$genelist<- read.delim("EC.genelist", header=TRUE)
  4. EC$genelist[,1] <- toupper(EC$genelist[,1])
  5. circ <- circle_dat(EC$david, EC$genelist)
  6. EC$genes <- read.delim("ECgenes3.txt", header=TRUE)
  7. EC$genes[,1] <- toupper(EC$genes[,1])
  8. EC$processes <- c('protein heterodimerization activity','oxoacid metabolic process','response to lipid hydroperoxide','chloroplast stroma','nucleotide diphosphatase activity','oxidoreductase activity','cytoplasmic part','carboxylic acid biosynthetic process','fatty acid transmembrane transporter activity')
  9. chord <- chord_dat(circ, EC$genes, EC$processes)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement