atac下游分析整理得到所有样本的.Narrowpeak文件后,依次进行peak merge/insertion count,即可得到所有样本peak的矩阵。行为各个peak区域的起始位点在染色体上的坐标,列是各个样本。
注意:下面两个输入路径都是nf-core中atac pipeline的输出结果
peak mergescript:
12345678910Rscript ./createIterativeOverlapPeakSet.R \ --metadata ./file.txt \ --macs2dir /data2/kunpeng/ATAC/integrate_for_heatmap_update/results/bwa/merged_library/macs2/narrow_peak/ \ --outdir ./peak_merging \ --suffix _summits.bed \ --blacklist ./hg38.blacklist.bed \ --genome hg38 \ --spm 0 \ --rule "(n+1)/2 ...
七月早起的晨总是不思考地生活是不行的,之前自己为何不愿整日沉湎于思考中?
理由嘛,一是琐事缠身,尽是不思考也能做的简单事,不思考也无妨;二是久违地想体验单纯地快乐,不思考、自由玩乐的快乐,过去总为各种各样的事情烦恼,我也想过一段不思考也没关系的生活;三是,孤独,阅读、电影、游戏、动画,这些能富足我精神的又恰恰全是自我封闭的活动,我希望不再由这些故事支撑我的精神,我希望在现实中找到精神的支柱。所以,我刻意去减少自我封闭的活动,去现实中寻找朋友,哪怕是酒肉朋友也无妨,哪怕仅仅是聊天喝酒打牌的关系也无妨,只要是能帮助我面对孤独的现实中的一切都可以。但是,即便如此,孤独仍时常伴随着我,我开始怀疑,这个年纪,“伴侣”是否孤独唯一的解药。我开始渴望伴侣。
嘛,总之回归到以往,把思考当成一种习惯吧,不能让kai失望。
昨天推完了药师线,最近很喜欢八方旅人,甚至可以说是向往,在这份向往的驱使下,我开始了音乐的学习,首先是英式八孔竖笛,之后可能还会接触爱尔兰哨笛和吉他。
睡不着的夜晚睡不着的夜晚就写写博客好了,写博客和夜晚意外得很搭配呢。为什么写博客?我希望记录一些自己的生活以及学到的知识,顺便积攒点睡意。而且能锻炼锻炼思考和交谈的能力,还能分享自己的喜欢的游戏好笼统的回答
最近的生活总觉得自己这个暑假会度过得十分充实呢。无拘无束,好不自在晚上几点睡都行,吃点烧烤、喝着啤酒、在这悠哉游哉地写博客,明明玩游戏也好不过现在没有兴致。每天就算不做科研也很充实,学竖笛、练长跑、打球,想做的事情多到做不完。最近习惯了不思考,无聊的时候打打lol,刷刷手机,刷手机的时间长了,总觉得这样不好,浪费时间不说,隐隐约约觉得自己思考和交谈的能力都在减弱。尽管不思考的日子很快活,但这样并不适合我,还是回归曾经的样子好了。最近的日子怎么过也有头绪,算法要学,rust要学,正好项目不忙,这一次不能半途而废了。竖笛、跑步、打球、之后还有吉他,这些带着练习。忙碌了一天之后就写写博客,挺好,这样自由的日子真是久违了。




