摘要：
2025年03月05日 安装 # add before submit job # source /home/zz950/softwares/miniconda3/bin/activate /home/zhl595/softwares/miniconda3/envs/RNAseq source /ho 阅读全文
摘要：
在我申请之前（6-8月），合作就已经开始了，估计是2022年1月，Cliff（Shirley组）在开发MIRA，已经到投稿阶段，而Nialy也已经在测序multi-omics了，就差人分析了。 multi-omics看似就是多了个ATAC-seq，但复杂度最少比单纯的RNA-seq多了5倍，Alle 阅读全文
摘要：
2023年01月12日 注意：要用conda进入r4p3进行开发。 发现了pkgdown可能更适配R包，而readthedocs则适合shell和Python。 用pkgdown来开发R包的网站，非常方便。 https://pkgdown.r-lib.org/ 用pkgdown为R程序包建立网站 准 阅读全文
摘要：
2023年01月19日 如何下载mm10的MT的SNP数据？ https://www.informatics.jax.org/snp - 这个没有mm10（GRCm38）的，只有GRCm39（mm39） Mouse Genomes Project: https://www.mousegenomes. 阅读全文
摘要：
We quantified this time dependent increase in basal-luminal mixing using entropy of cell-type probabilities, which robustly agrees with our per-sample 阅读全文
摘要：
2023年01月28日 Moreover, genomics often cannot distinguish “passenger” mutations from “drivers” that promote tumorigenesis and might be therapeutic targe 阅读全文
摘要：
除了刚接触10x的那会儿，还真没怎么亲自倒腾过fastq的制作。 正常从测序商那里拿到的应该是bcl的原始数据，需要自己做一步bcl2fastq。 后面大家都觉得这一步太麻烦了，没必要，所以大部分测序商也就自己做了。 这也是对的，生信分析的起点是fastq，不应该是bcl，数据保存的最原始格式也应该 阅读全文
摘要：
有的包在升级后函数和API就完全变了，甚至删除了一切原来的功能函数，这导致后面分析数据非常的畸形。 但没办法，接下了这样的项目，含着泪也得做出点东西出来，于是我就想同时把两个包装上。 python是肯定不行的，网上有一些教程，但都不能正常工作。 于是曲线救国，那我就把包下载到本地，然后用sublim 阅读全文
摘要：
transcription factor perturbation gene-sets 目的比较直接：我有个DEG基因集，我想要知道那个TF最有可能调控这个基因集。 Transcription factor perturbation gene ontology analysis (EnrichR) 阅读全文
摘要：
2023年08月28日 从我第一次接触motif分析开始，就意识到motif数据库的不完善，那么多的已知的TF，有较为明确motif的却不到一半，已知的部分质量也不高。 而motif其实至关重要，它是Protein与DNA互作的pattern和证据，蛋白与DNA的互作也决定了基因的转录，染色质的表观 阅读全文
摘要：
目前对peak的数据处理上，发现测序深度对peak的数量有很大影响，即使做了normalization也没办法，所以这里希望从原始的bam文件开始做downsampling。 参考一：Downsample BAM file to specific amount of reads input_dir= 阅读全文
摘要：
2025年01月15日 一个更直接的Heatmap出图方式，直接从bigwig读取数据。 教程：https://crazyhottommy.github.io/reproduce_genomics_paper_figures/07_figure1_i_j_k.html https://crazyho 阅读全文
摘要：
常在河边走，哪能不湿鞋。 作为生信工作者，尤其是计算生物学，遇到BUG真的是不能太常见，能否顺利DEBUG是此行业工作者的基本功。 Mac MacMini M1芯片则有些包装不了，在osx-64上的包，无法在osx-arm64上安装。 注意：千万不要在M1的Mac上装M1的conda！！！【直接装x 阅读全文
摘要：
2023年05月04日 在搞懂DiffBind后就几乎没有用过deeptools，因为R可以用jupyter notebook，代码和结果方便撰写和保存。 但是，最近被怼了，质疑了我的peak calling，一个核心的问题就是macs call出来的peak太多了，有6万个左右，我觉得这个比较正常 阅读全文
摘要：
一般在传输过程中中断，fastq就会truncate，但仍然可以用zless或者zcat打开查看。 真到万不得已的时候还是要截取可以利用的部分，然后用bbmap修复一下。 bbmap repair.sh -Xmx14g in1=1.fq.gz in2=2.fq.gz out1=1_1.fq.gz o 阅读全文
摘要：
10x是最成功的单细胞公司，目前最成功的商业平台，已经有很多成熟的单细胞产品线，不管是新手还是进阶的高手都在玩他们平台产生的数据。 这里试图对10x的核心技术和数据来做一个剖析，试图理解其核心原理，加速数据分析进程，辅助科研发现。 我目前玩过的10x数据： scRNA-seq (antibody h 阅读全文
摘要：
这个要严谨一点，众所周知，小鼠是小写，人是大写，以前为了方便都是直接一个toupper函数完成转换，但这样做实在是太粗糙了，大概有三分之一的基因会丢失。 我简单统计了一下： 人鼠同源的有23242个基因注释，其中有7695个基因是无法简单对应的，所以你之前的分析会漏掉很多基因。 建议使用ensemb 阅读全文
摘要：
2023年04月24日 参考：decoupleR：丰富你的通路富集分析结果 2022年12月28日 没想到碰到了一个神器：PROGENy。可以给常见的14种signaling pathway活性打分。 参考：看着14个肿瘤相关通路活性就够了 buik RNA-seq: http://localhos 阅读全文
摘要：
搜索历史文件真的是最基本的需求，听说everything是个好东西，但Mac平台自带的splotlight search已经足够好用，直达文件。 Linux我则不太会搜索，find命令有点复杂，最近无意搜索Linux平台是否有类似everything的工具，还真发现了一个非常活跃的fzf工具，非常好 阅读全文
摘要：
好久没有继续我的学习/写作/科普大业了，这个positive feedback loop是非常powerful的。 这一期的主题，Apc - Wnt - CRC 这也是我今后的研究主题。 先抛出问题： Wnt pathway的组成成分有哪些，都控制了哪些重要的细胞功能？ 正常情况下，Apc是如何与W 阅读全文
