Question: Suggestion for circlize package in R for plotting the datasets
0
gravatar for sinha.puja
7 months ago by
sinha.puja10
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 • 217 views
ADD COMMENTlink written 7 months ago by sinha.puja10

Cross-posted on Bioconductor: https://support.bioconductor.org/p/128388/

ADD REPLYlink written 7 months ago by Kevin Blighe65k
1
gravatar for Arup Ghosh
7 months ago by
Arup Ghosh2.7k
India
Arup Ghosh2.7k wrote:

Check this tutorial by the developer of the circlize package. https://jokergoo.github.io/blog/html/segmentation_WGBS_dmr.html

ADD COMMENTlink written 7 months ago by Arup Ghosh2.7k

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

ADD REPLYlink written 7 months 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 7 months ago by Arup Ghosh2.7k

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 7 months ago by sinha.puja10
1

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 7 months ago by Kevin Blighe65k

The code is below:

circos.initializeWithIdeogram(species = "rn6")
circos.par("track.height"=0.20)
circos.genomicTrackPlotRegion(data = yc,ylim = NULL, numeric.column = 4, 
+                               panel.fun = 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).

Thanks.

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

Help
Access

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