美文网首页
无生物学重复的差异分析—Gfold

无生物学重复的差异分析—Gfold

作者: 白米饭睡不醒 | 来源:发表于2021-03-19 11:24 被阅读0次

一、安装gfold

Gfold软件是做没有生物学重复样本首选的软件https://zhanglab.tongji.edu.cn/softwares/GFOLD/index.html

#安装gsl最新版
##以下所有路径均为解压后gsl内的路径
wget http://mirrors.ocf.berkeley.edu/gnu/gsl/gsl-latest.tar.gz
tar -zxv -f gsl-latest.tar.gz
cd gsl-2.6/
./configure --prefix=/home/fanzhiyu_31946012_cyl_150/chayifenxi/gsl-2.6
make
make check(选做)
make install

#安装gfold
cd ..
wget https://zhanglab.tongji.edu.cn/softwares/GFOLD/gfold.V1.1.4.tar.gz
tar -zxf gfold.V1.1.4.tar.gz 
cd gfold.V1.1.4/
make
#发现报错,输入以下命令
export CXXFLAGS="-g -O3 -I/home/fanzhiyu_31946012_cyl_150/chayifenxi/gsl-2.6/include -L/home/fanzhiyu_31946012_cyl_150/chayifenxi/gsl-2.6/lib"
export LD_LIBRARY_PATH="/home/fanzhiyu_31946012_cyl_150/chayifenxi/gsl-2.6/lib:"$LD_LIBRARY_PATH
g++ -O3 -Wall -g main.cc -o gfold -lgsl -lgslcblas -I/home/fanzhiyu_31946012_cyl_150/chayifenxi/gsl-2.6/include -L/home/fanzhiyu_31946012_cyl_150/chayifenxi/gsl-2.6/lib
#有warnning不用理会
./gfold -h

出现此 安装成功

二、差异分析

1、把实验组和对照组的表达值分别建立两个csv文件,各包含如图五列,保证第一列顺序相同,如果不需要计算后两列的值,后两列可以随便填值
1.png
2、将自己生成的csv文件,再在 R 中统一一下格式(下面的代码是在简书作者没有猫但是有猫饼文章中粘的)
#新生成的文件名我一般不改,就是 raw_data_Vec , raw_data_LSH 反正后面还会改掉
file1 = "ck.csv" 
raw_data_ck <- read.csv(file1, stringsAsFactors=FALSE) 
file2 = "VI.csv" 
raw_data_vi <- read.csv(file2, stringsAsFactors=FALSE)
#统一列名别出错 
colnames(raw_data_ck) <- c("GeneSymbol", "GeneName", "Read Count", "Gene exon length", "RPKM") 
colnames(raw_data_vi) <- c("GeneSymbol", "GeneName", "Read Count", "Gene exon length", "RPKM") 
#导出文件,这个文件名我也不改,不然在 Linux 里改好麻烦
write.table(raw_data_ck, file="C:/Users/Administrator/Desktop/GFOLD/Sample1ck.read_cnt", row.names=F, col.names=F, quote=F, sep="\t") 
write.table(raw_data_vi, file="C:/Users/Administrator/Desktop/GFOLD/Sample2vi.read_cnt", row.names=F, col.names=F, quote=F, sep="\

3、将文件上传到服务器,运行如下命令
./gfold diff -s1 /home/fanzhiyu_31946012_cyl_150/chayifenxi/gfold.V1.1.4/Sample1ck -s2 /home/fanzhiyu_31946012_cyl_150/chayifenxi/gfold.V1.1.4/Sample2vi -suf .read_cnt -o sample1VSsample2.diff 

输入后出现此页面
4.筛选差异基因
#大于0为上调,小于0下调(第3列是gfold值所在列,第4列是E值)
awk '{if($3>0 && $4=1) print $0}' sample1VSsample2.diff OFS='\t' > up_diff.gene
awk '{if($3<0 && $4=1) print $0}' sample1VSsample2.diff  OFS='\t' > down_diff.gene

#筛选两倍差异基因,即大于1和小于-1的
awk '{if($3>1 && $4=1) print $0}' sample1VSsample2.diff OFS='\t' > up_diff.gene-2
awk '{if($3<-1 && $4=1) print $0}' sample1VSsample2.diff  OFS='\t' > down_diff.gene-2

相关文章

  • 无生物学重复的差异分析—Gfold

    一、安装gfold Gfold软件是做没有生物学重复样本首选的软件https://zhanglab.tongji....

  • 无重复RNAseq样本Gfold值的计算

    使用Gfold进行无重复RNAseq数据基因差异分析 1. 准备资料 输入的资料为read_count,要对两个样...

  • 在线作图|在线做有/无生物学重复的差异分析

    有/无生物学重复的差异分析 差异分析在转录组分析中运用广泛,它是计算两组间差异表达的基因的分析,常包含差异表达水平...

  • 无生物重复样本的差异分析

    目前由于测序的价格,还有样本自身的珍贵稀缺性,部分实验设计仍然是没有生物学重复的。无重复的数据做差异分析是一件很麻...

  • 差异分析+火山图+COX模型构建

    简书:没有生物学重复的转录组数据怎么进行差异分析? edgeR需要的数据是reads数,可以设置BCV值,做单样本...

  • gfold 分析

    之前分析遇到2个问题:1.样本只有一个,没有replicates,是不是就没有办法进行分析?2.一个样本有3个重复...

  • 转录组数据分析—clusterProfiler富集分析

    差异分析之后,为了了解我们获取的差异基因的功能及参与的生物学进程,我们需要对差异基因进行功能富集分析。基因富集分析...

  • 几种常用的差异分析方法简介

    几种常用的差异分析方法简介 如今在生物学研究中,差异分析越来越普遍,也有许多做差异分析的方法可供选择。但是在实际应...

  • DEseq2差异表达分析

    转录组数据绕不过差异分析。为什么选择这两个包呢?DEseq2针对有生物学重复的样本。(一般情况下应该是都需要生物学...

  • 生物学重复与技术重复

    生物学重复与技术重复 在利用基因芯片或者RNA-seq做基因表达分析的时候,经常听说生物学重复和技术重复,这篇文章...

网友评论

      本文标题:无生物学重复的差异分析—Gfold

      本文链接:https://www.haomeiwen.com/subject/rfescltx.html