Question: Suggestion for circlize package in R for plotting the datasets
gravatar for sinha.puja
6 weeks ago by
sinha.puja10 wrote:

Hi All,

I have large data set with methylation values ranging from 0-1. The data contains around 3 million rows of values corresponding to the chromosomal locations. I want to split/sort my data for plot using circos plot. I thought of doing like sorting the value range 0.2-0.4 as hypomethylation category and range 0.7-1.0 as hypermethylation category. However, with this approach the plot again looks cluttered. Another way is that I can split the data into two parts, one part with values lower than median, and another part with values larger than median.

Can anyone suggest me which approach is better or if there are other approaches I can implement in my datasets to make data possible for circos plot using circlize package.

Thanks in advance.

circlize • 112 views
ADD COMMENTlink written 6 weeks ago by sinha.puja10

Cross-posted on Bioconductor:

ADD REPLYlink written 6 weeks ago by Kevin Blighe56k
gravatar for Arup Ghosh
6 weeks ago by
Arup Ghosh2.3k
Arup Ghosh2.3k wrote:

Check this tutorial by the developer of the circlize package.

ADD COMMENTlink written 6 weeks ago by Arup Ghosh2.3k

Thank you for the tutorial but still I did not get the appropriate answer.

ADD REPLYlink written 6 weeks ago by sinha.puja10

It's hard to provide a better solution without knowing the number(large data is a very ambiguous term) of data sets you are trying accommodate in a plot and how the Circos plot outcome looks with your current configuration.

ADD REPLYlink written 6 weeks ago by Arup Ghosh2.3k

It's very dense and I don't see the individual dots except two separate bands in the track based on my two range of data sets.

ADD REPLYlink written 6 weeks ago by sinha.puja10

sinha.puja, to be fair, if you want help with something like this, then you need to post the exact code that you have been using, and also paste samples of the data that you currently have. You have not even provided a screenshot of what you are seeing on your screen.

ADD REPLYlink written 6 weeks ago by Kevin Blighe56k

The code is below:

circos.initializeWithIdeogram(species = "rn6")
circos.genomicTrackPlotRegion(data = yc,ylim = NULL, numeric.column = 4, 
+                      = function(region,value,...) {
+                                 cond <- value[,1] < 0.7
+                        circos.genomicPoints(region[cond,], value[cond,], pch = ".", cex = 0.1,
+                                                      col = "red")
+                      circos.genomicPoints(region[!cond,], value[!cond,], pch = ".", cex = 0.1,
+                                                      col = "blue")
+                               })

The snapshot of image link is below:

Circos Plot

About my data I already described above consisting of 4 columns( Chromosome, start, end and methylation values).


ADD REPLYlink modified 5 weeks ago • written 5 weeks ago by sinha.puja10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1246 users visited in the last hour