Hi,
I am quite new to R programming and could really appreciate your help to resolve this issue. I'm stuck with it, even after having carefully read the Gviz user's guide.
I have a datatrack for one sample : on the X axis the genomic positions on a particular chromosome, on the Y axis the coverage of SNPs. Besides the vectors containing the position and the coverage, I have a vector for the color of each point (color_point : eg all are "grey" and 1 in particular is "red"), and a vector for the size of each point (size_point : eg all are "1" and 1 in particular is "5").
When I create the datatrack to show the value at each genomic position, how come the point related to my SNP of interest is not highlighted? Every point is grey and of a size of 1. And there's no error reporting. I also changed all the "grey" in "blue", and this works, so the issue is not changing the color of my dots, but highlighting one.
Here is my code :
dtrack <- DataTrack(data = coverage,
start = coords[-length(coords)],
end = coords[-1],
chromosome = chrom,
genome = "hg19",
col.symbol=color_point,
cex.symbol=size_point
)
plotTracks(dtrack,
type = c("p"),
from = lim[1], to = lim[2])
And here is the plot : https://www.dropbox.com/s/08co2cobh3ghm2i/Rplot01.jpeg?dl=0
Thank you very much if you can show me my mistake.
Best,
Cedric