问题最初产生于conda安装了java-jdk和fastqc后,fastqc仍然无法使用,用java --version确认之后发现java的VM存在问题。于是手动下载了JA...
问题最初产生于conda安装了java-jdk和fastqc后,fastqc仍然无法使用,用java --version确认之后发现java的VM存在问题。于是手动下载了JA...
1、创建账号将数据上传到GEO[https://www.plob.org/tag/geo]数据库,首先要创建并登陆NCBI帐号, 然后进入提交的网址:https://www....
发现awk比较的时候可以用-代替某一个文件,这样让整个代码更灵活: 取代第一个文件的位置: 这样在提取bam文件的信息的时候比较方面,比如:
python默认输出的pdf格式是type3,这种不能在 Adobe Illustrator中编辑文字等信息,需要改成type42. 或者直接更改matplotlib的默认设...
加载所需R包 设置工作路径 构建测试数据集 将热图结果按聚类后的顺序输出
不好意思 忘记改了!
ChIP-seq分析的一般流程方法现在ChIP-seq的数据基本是最常见的测序数据类型之一,主要有Transcription factor ChIP-seq和Histone ChIP-seq。前者是看转录因子...
z-score计算方法为:Z =(x-μ)/ σμ为均值,σ为标准差。 以下是R中将z-score转为p.value的方法: pnorm(q, mean = 0, sd = ...
有时ggplot作图时希望某些字变成斜体,可以导入Unicode包实现: 可以实现p为斜体: 注意保存的时候,要使用cairo_pdf()而不是pdf(),否则p有可能显示不...
来自https://github.com/r-dbi/RPostgres/issues/110[https://github.com/r-dbi/RPostgres/issu...
在本地pyplot画图可以运行,但是在服务器显示以下错误: 其实这是因为matplotlib是默认画图backend是TkAgg,这个需要有GUI的图形界面。只需要指定不需要...
bigwig转bed文件首先需要将bigwig转为wig文件 利用BEDOPS的 wig2bed]将wig文件转 bed文件 参考:https://bioinformatic...
作者 | Arno审稿 | 童蒙编辑 | amethyst 上一期我们介绍了ATAC-seq相关的背景知识。ATAC-Seq能帮助我们从全基因组范围内推测可能的开放染色质位点...
构建一个显示的数值的函数,将plt.pie中的autopct=该函数即可。 代码: 生成图:
有一个叫group_g的data.frame,想按group那一列转为list,代码如下:
以自己的seurat对象的分群来鉴定细胞类群 参考:https://github.com/ggjlab/scHCL[https://github.com/ggjlab/scHCL]
此功能使用样本中的信息通过指定的道具对每个分子的读数进行下采样。然后,它基于具有非零读取计数的分子构造一个UMI计数矩阵。目的是消除技术噪声中的差异,这些差异可以按批次进行聚...
参考:https://stackoverflow.com/questions/59101791/seurat-dimplot-highlight-specific-group...
linux下安装nbextensions,jupyter_nbextensions_configurator。 重启jupyter,在根目录下,点击Nbextensions标...
cellranger mkfastq Illumina测序下机后的数据为 原始数据(raw base call )BCL文件,拿到BCL文件之后,第一步是使用cellrang...