gravatar for marvel479

2 hours ago by

Forgive my confusion, but I can't figure out the correct usage of circlize to draw a histogram. I want to draw a histogram for one of its tracks and would like to use my read counts for Y-axis, and the x should be the genomic position, however wide or narrow. The syntax however is:

circos.trackHist(factors, x, ....)

that doesn't work though. How should I imply that my x should have chromosome no., start and stop position information, the bin height being no. of reads, for each row?

I am assuming I can write:

circos.trackHist(factors = mydataframe$reads, x = x)

This is what my file looks like: head(dataframe)

 Chrom     Start       End Reads
1  chr1  42938268  42938392     2
2  chr1  63224349  63225050     2
3  chr1  71646172  71646228     2
4  chr1  72904979  72905055     2
5  chr1 138065763 138065817     2
6  chr1 164635940 164635994     2

Source link