I have an Rscript that constructs plots for a number of genes. I have hardcoded the gene name and color assignments in the script itself with the following code:
grid.col = c(gene0059="aquamarine",gene0086="bisque",gene0107="blue",gene1658="brown",
gene0171="brown1",gene0246="burlywood4",gene0404="cadetblue",gene0439="cadetblue1",
gene0518="chartreuse",gene0716="chartreuse4",gene1099="chocolate",gene1525="chocolate4",
gene1527="coral",gene1621="coral4",gene1661="cyan",gene1701="cyan3",gene1907="darkgoldenrod",
gene2223="purple",gene2244="green",gene2552="green4",gene2584="honeydew",
gene2692="hotpink",gene2874="indianred1",gene2884="khaki1",gene2888="red",
gene0016="lightblue1",gene0096="lightgoldenrod3",gene1959="lightpink",gene2271="lightsalmon",
gene2804="lightskyblue",gene2885="magenta",gene3000="blue",gene1689="mediumorchid3",
gene1776="mediumpurple1",gene2628="olivedrab2",gene0963="orange",gene1137="orangered",
gene0621="palegreen1",gene2733="royalblue1",gene2806="seagreen1",gene1702="purple1")
When plotting, I can call the grid.col for the colors. Specifically I am using this with the chordDiagram command in the circlize package.
When I import this as a list (example pasted next) the colors do not get assigned properly to the gene IDs, so the plots do not use the right colors.
gene0059="aquamarine"
...
gene1702="purple1"
How do I both format the file and import it to correctly assign the colors to the gene labels?
2 answers
Ahh, apologies for the confusion. My original code worked fine. But I realized that if when wanted to use another set of genes and colors, I'd have to open the file and change the code. That's fine for me, but I want to share this code with people who haven't coded in R ever. So I wanted to find a way to have the colors and gene names provided as a command line argument so that it was easy to change for different experiments.
So both versions work, the only difference is my original was hard coded and the updated version requires a command line input.
The problem was the original script contained the gene and color designations within the code. My question was how do replicate this action but through importing files through the command line and parse them so the graph would print properly. I did not know how to do this and that lead to my initial question. I had figured how (below) how to correctly import and parse the files to replicate the hardcoded part of the script used before. So now I do not have to edit the code every time I want to change the gene IDs or the associated color, just import a new file.
I got this to work and want to put it here case anyone finds this in the future.
I imported two data.frames, one with the gene names (only_names.txt) and one with the colors (only_colors.txt). Then I converted both into vectors and combined them into a single named vector:
names <- read.table(file = "only_names.txt", header = F)
nameVec <- unlist(names, use.names = F)
colors <- read.table(file = "only_colors.txt", header = F)
colorsVec <- unlist(colors, use.names = F)
my_vector <- setNames(colorsVec, nameVec)
When I used my_vector in the chordDiagram command it worked to get the colors correct:
chordDiagram(matrix1, annotationTrack = "grid", transparency = 0.1, grid.col = my_vector, preAllocateTracks = list(track.height = max(strwidth(unlist(dimnames(matrix1))))))
You could add a single .csv file as a second command line arg, read it in as a data.frame and make a named vector using the information from the first two columns. That is a lot more user friendly than inputting two separate .txt files that are not obviously paired.
Yes, thanks! Here is the code I use to import a single two column text file, then split it up. Much cleaner!
# import name_color_file and set vector for colors to use
name_color_file <- read.csv(file, header = F)
# split file by columns into two vectors
nameVec <- unlist(name_color_file[[1]], use.names = F)
colorVec <- unlist(name_color_file[[2]], use.names = F)
# make a named vector for gene labes and colors
grid.col <- setNames(colorVec, nameVec)
Log in to answer this question.
Please show us your full plotting code.
Sure thing. Wasn't sure if it was too much to post:
The input file is just a tab delimited, three column file that has the first gene ID (geneID_1), the connecting gene ID (geneID_2), and value that is used for drawing the chord diagrams (Count).
The manual uses a matrix (wide form of the long form dataset you're using). See: https://jokergoo.github.io/circlize_book/book/the-chorddiagram-function.html#basic-usage-of-making-chord-diagram
Try getting your data in that format, that might solve the problem.
Thanks for the reply. I use
acastto convert the DF into the format used by chordDiagram. So I don't think the input data format is the issue, but I will certainly take a look.