二更~corrplot我觉得有一点不好就是改变不了图的纵横比,其他一流 1. 构建相关系数矩阵 2. 构建P值矩阵 3. corrplot绘制相关系数热图 最终结果
二更~corrplot我觉得有一点不好就是改变不了图的纵横比,其他一流 1. 构建相关系数矩阵 2. 构建P值矩阵 3. corrplot绘制相关系数热图 最终结果
参考:https://www.jianshu.com/p/3b9eee8c31cc[https://www.jianshu.com/p/3b9eee8c31cc]优雅R:「W...
@5f1b9fff33d0 rMATs,LINUX系统下运行
RBP AS生信流程(N)RBP与AS相关性
以ZFAND4为例:在https://www.ncbi.nlm.nih.gov/[https://www.ncbi.nlm.nih.gov/]里选择Gene,搜索ZFAND4...
cytoscape就可以~
可变剪接文章套路(一)PSI值:(再翻墙找一下答案) PSI ( Percent-spliced-in,拼接百分比)可以针对isoform, exon, ASE 进行计算,对于ASE来说, PSI...
输入文件:1.RBP与AS相关性文件allcorResult.txt: 2.节点属性文件nodeType.txt(构建方法在脚本提取节点名.R里): 方法:将allcorRe...
我用的是GSE46224这个数据集,readLength是101
RBP AS生信流程( 四)rMATsrMATs rMATs运行: 在rMATstest这个文件夹里有官网下载好的2x2测试文件: 用samtools查看bam文件: 可以看到第三列为1,代表1号染色体(这里是1...
感谢大佬~
说说R语言的转置在群里看到了一个问题,他用R语言的转置函数t()对下面形式的数据进行了转置 转置之后的数据结构如下 他发现转置之后多了一行,也就是有了geneid这一行,和他预期的不符合。按...
在GSE46224文件夹里,定量以后,获得(15+8)个fc.txt文件,这个是RNA-seq测序数据的表达矩阵,RNA-seq的count数据符合负二项分布,因此用SVA包...
用jimmy老师的方法定量:
bioMart包是一个连接bioMart数据库的R语言接口,能通过这个软件包自由连接到bioMart数据库这个·包可以做以下几个工作1.查找某个基因在染色体上的位置。反之,给...
从数据库下载的文件为PSI_download_STAD.xlsx: 第一列是基因名字,第二列是数据库给的id,第三列为可变剪接类型,后续把三列统一为一列。4-10列为基因可变...
featurecounts的输出文件 查看featurecounts输出文件: 第一行是运行代码,第一列为基因名,第七列为每个基因上的count数,现在要舍去第一行,并只保留...
samtools 尝试了各种方法后,最后用conda install samtools安装成功 Piranha 安装参考http://smithlabresearch.org...
RSEM软件 会生成基于基因的定量结果_rsem.gene.results和基于转录本的定量结果_rsem.isoform.results,isoforms.results中...
STAR比对 比对原理参考生信媛的文章:https://www.jianshu.com/p/16c938fd3bd6[https://www.jianshu.com/p/16...
1.查看cpu信息 model name : Intel(R) Xeon(R) CPU E5-2683 v3 @ 2.00GHzmodel name : Intel(R)...