- 第一棒:CNS里的泡泡图
- 第二棒:怎样用HOMER找motif?
- 第三棒:怎样用HOMER算出的P-value画出CNS级别的泡泡图?
这是小伙伴的命题作文,下面这个数据是曾老湿给出来的,其中header行竟然放在了最后一行,我手工把它挪到了第一行,里面还有一个数据缺失的,会导致读了数据后会有NA值,需要去除。
$ file summary_ATAC_seq_motif.xlsx
summary_ATAC_seq_motif.xlsx: Microsoft Excel 2007+
这是一个Excel文件,我们可以用readxl
来读:
require(readxl)
x <- read_xlsx("summary_ATAC_seq_motif.xlsx")