R

Heatmap demystified part II: Using heatmap to represent continuous values

In my last blog post, I showed an example to use heatmap to represent discrete values. I am going to continue the theme to introduce heatmap to represent…

Heatmap demystified part I: use heatmap to represent discrete values

In many Genomic papers, you will see heatmaps. Heatmaps are of no mystery. It is a way to visualize the data a.k.a. using colors to represent values.…

bring the power of R to command line

If you ever wants to use R on the command line, now you have several opitions.

Breakpoints clustering for structural variants

When structural variants are represented in bedpe format, each structural variant is represented as two linked breakpoints. Each breakpoint is represented…

read in all the files in a folder to R

see a post\ I downloaded some broadpeak data from UCSC and put them in a folder /Users/mtang1/projects/ENCODE_ChIPseq/data/\ I cut only first four columns…

Upset plot for overlapping ChIP-seq peaks: an alternative to Venn diagram

I am going to demonstrate how to make Upset plot using UpSetR package for ChIP-seq peaks. After running a peak-calling method for 3 samples, you want to…

How to make a box-plot with jittered points for multiple groups

Let me demonstrate it with an example. Use the ToothGrowth data set in the ggplot2 library

I promised a post with heatmaps, and I delivered more with PCA, MDS and clustering

I am back from the 1st summer Institute in statistics for big data held in the university of Washington. I have learned so much!

get all the promoter sequences of human hg19 genome

One of my former friends (a biologist who does not know much about computer) wants me to help her to get all the promoter sequences from the human genome.…

simulation of distribution of means draw from exponential distribution

“central limit theorem (CLT) states that, given certain conditions, the arithmetic mean of a sufficiently large number of iterates of independent random…