Pn3_Pn5_manuscript_graphs_code/helical_wheel_PN3_PN5.R at main · Menukabh/Pn3_Pn5_manuscript_graphs_code · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
##make the helical wheel of the antimicrobial peptides

setwd("C:/Users/bhandari.72/Desktop/Hayes_poster/PN3 and PN5 results")
library("helixvis")
library(gridExtra)

##you can save the plot in the higher resolution using the dev.off function in R

data<-read.table('PN3_PN5_seq_helical_wheel.txt', header = TRUE, stringsAsFactors = FALSE)


tiff('pn5.tiff', units="in", width=5, height=6, res=300)

## Make helical wheel plot using the default color function
plot<-draw_wheel(data$Seq[1], labels = TRUE, label.col = "black",
                 fixed= TRUE,
                 legend=TRUE)


##If you would like to change the color of the circle
plot<-draw_wheel(data$Seq[1], col = c("purple",  "orange", "red", "lightgreen"),
           labels = TRUE, label.col = "black",
           fixed= TRUE,
           legend=TRUE)
plot


plot_pn5<-draw_wheel(data$Seq[2], col = c("purple",  "orange", "red", "lightgreen"),
                 labels = TRUE, label.col = "black",
                 fixed= TRUE,
                 legend=TRUE)
plot_pn5

##to save graph
dev.off

## tovisualize graph in R
dev.new()