day32 失败:EPIC-Seq流程

由于对cfDNA的兴趣,查到22年这篇在Nature biotechnology上的文章Inferring gene expression from cell-free DNA fragmentation profiles.想尝试跑一下EPIC-seq的流程。感谢文章作者提供的这么全面的资料。但,仍然很多不懂。

无论用conda安装R,还是下载gz文件本地安装,都是失败。
太难过了。
一、 Rscript 这个是在Linux里面使用R要调用的命令。

Rscript runEPIC.R –-bamdir  D:/work/research/ssDNA/result/ssDNA_1stsequence/bamdir –-tssinfo D:/work/research/ssDNA/result/ssDNA_1stsequence/code/priordata/all.tss.genes.canonical.ensembl75.txt  –-panelbed D:/work/research/ssDNA/result/ssDNA_1stsequence/code/priordata/all.tss.genes.canonical.ensembl75.txt.noheader –-targeted no  –-outdir D:/work/research/ssDNA/result/ssDNA_1stsequence/outdir  –-mapq 30 –-groupref <character> –-groupsize <double, number of genes per group> –-skipHist <binary, 0 or 1> –-skipOCFFrags <binary, 0 or 1>

二、先建立虚拟环境
conda create -n R_cfDNApipe python=3.8.3
输入:conda activate R_cfDNApipe,报错。
于是先输入:source activate,再输入 conda activate R_cfDNApipe

image.png

查看所有虚拟环境 conda env list
image.png

退出虚拟环境: conda deactivate
参考教程: 虚拟环境
虚拟环境:https://blog.csdn.net/robot8me/article/details/109471568
conda的配置:https://www.bilibili.com/read/cv8956636
conda瘦身:
conda clean --packages --tarballs 删除缓存包,重新下载
conda clean -i清除缓存,重新加载condo地址
conda clean -p //删除没有用的包
conda clean -t //删除tar包
conda clean -y --all //删除所有的安装包及cache
三、在Linux里直接安装R(直接安装失败,就是到PCRE那一步怎么都不成功)
参考教程:
https://www.136.la/jingpin/show-30899.html
https://blog.csdn.net/tandelin/article/details/1183937531
https://www.cnblogs.com/jessepeng/p/10984983.html
下载文件 R-4.1.1.tar.gz,放到software目录下
2,解压:tar -zxvf R-4.1.1.tar.gz
3,进入解压后的文件夹
进行configure,make, make install
./configure --prefix=/data/zds209/software/R --enable-R-shlib
make
make install
4,configure遇到各种报错:

  1. zlib version >= 1.2.5需要安装
    由于公共software里有,只需直接.bashrc添加:
    export PATH=/data/software/zlib1.22.11/bin:$PATH
    export CFLAGS="-fPIC -I/data/software/zlib1.22.11/include"
    export LDFLAGS="-fPIC -L/data/software/zlib1.22.11/lib"

  2. bzip2需要安装
    下载https://sourceforge.net/projects/bzip2/files/latest/download

 tar xzvf bzip2-1.0.6.tar.gz
 cd bzip2-1.0.6
 make -f Makefile-libbz2_so
 make clean
 make
 make install PREFIX=/data/zds209/software/bzip2-1.0.6
  1. lzma version >= 5.0.3
    下载http://tukaani.org/xz/xz-5.2.3.tar.gz
tar xzvf xz-5.2.3.tar.gz./configure --prefix=5. 
cd tar xzvf xz-5.2.3.tar.gz
./configure --prefix=/data/zds209/software/xz-5.2.3
make -j3
make install
  1. PCRE >= 8.32
    下载 https://nchc.dl.sourceforge.net/project/pcre/pcre/8.42/pcre-8.42.tar.gz
tar -zxvf pcre-8.42.tar.gz
cd pcre-8.42

./configure --enable-utf8 --prefix=/data/zds209/software/pcre-8.42
make -j3
make install

到.bashrc中

export LC_ALL="zh_CN.UTF-8"
export LANG="zh_CN.UTF-8"
export PATH=/data/zds209/software/pcre-8.42/bin:$PATH
export LD_LIBRARY_PATH=/data/software/pcre-8.42/bin:$LD_LIBRARY_PATH

查看一下pcre的版本:
zds209 07:14:44 ~pcre-config --version 8.44  查看UTF-8: zds209 07:15:09 ~locale
LANG=zh_CN.UTF-8
LC_CTYPE="zh_CN.UTF-8"
LC_NUMERIC="zh_CN.UTF-8"
LC_TIME="zh_CN.UTF-8"
LC_COLLATE="zh_CN.UTF-8"
LC_MONETARY="zh_CN.UTF-8"
LC_MESSAGES="zh_CN.UTF-8"
LC_PAPER="zh_CN.UTF-8"
LC_NAME="zh_CN.UTF-8"
LC_ADDRESS="zh_CN.UTF-8"
LC_TELEPHONE="zh_CN.UTF-8"
LC_MEASUREMENT="zh_CN.UTF-8"
LC_IDENTIFICATION="zh_CN.UTF-8"
LC_ALL=zh_CN.UTF-8

用which查位置,发现不对劲,怎么conda里面也有。赶紧到.bashrc里面把conda的PATH注释掉。也不行。

image.png

于是进入到conda虚拟环境中,conda安装PCRE。

到R-4.1.1目录下。

./configure --prefix=/data/zds209/software/R  --enable-R-shlib --with-readline=yes --with-libpng=yes --with-blas --with-pcre1 LDFLAGS="-L /data/zds209/software/bzip2-1.0.6/lib -L /data/zds209/software/xz-5.2.3/lib -L /data/zds209/software/bzip2-1.0.6/include -L /data/zds209/software/xz-5.2.3/include -L /data/zds209/software/pcre/lib -L /data/zds209/software/pcre/include" 

四、虚拟环境中安装R
https://www.jianshu.com/p/b9eb874fc8f4

conda clean -i #清缓存
conda create -n R_EPIC-seq python=3.8.3 #新建一个虚拟环境
source activate
conda activate R_EPIC-seq
conda search R
conda install r-base=4.0#总是在solving environment

仍然失败。。。。。



教程:
https://yunwei365.blog.csdn.net/article/details/108126550?spm=1001.2101.3001.6650.12&utm_medium=distribute.pc_relevant.none-task-blog-2~default~OPENSEARCH~default-12.no_search_link&depth_1-utm_source=distribute.pc_relevant.none-task-blog-2~default~OPENSEARCH~default-12.no_search_link&utm_relevant_index=18








©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 219,366评论 6 508
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,521评论 3 395
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 165,689评论 0 356
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,925评论 1 295
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,942评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,727评论 1 305
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,447评论 3 420
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,349评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,820评论 1 317
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,990评论 3 337
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,127评论 1 351
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,812评论 5 346
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,471评论 3 331
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,017评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,142评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,388评论 3 373
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,066评论 2 355

推荐阅读更多精彩内容