实验记录 | 6/5

最新推荐文章于 2024-03-18 16:52:53 发布

今天也是个妖精头子呀

最新推荐文章于 2024-03-18 16:52:53 发布

阅读量1.2k

点赞数

分类专栏：谱系追踪

本文链接：https://blog.csdn.net/weixin_40640700/article/details/117585671

版权

谱系追踪专栏收录该内容

82 篇文章 11 订阅

订阅专栏

今天早上（8：57）过来，发现R已经安装完成了，非常棒。
./Rscript

Usage:
/path/to/Rscript [–options] [-e expr [-e expr2 …] | file] [args]
–options accepted are
–help Print usage and exit
–version Print version and exit
–verbose Print information on progress
–default-packages=list
Where ‘list’ is a comma-separated set
of package names, or ‘NULL’
or options to R, in addition to --slave --no-restore, such as
–save Do save workspace at the end of the session
–no-environ Don’t read the site and user environment files
–no-site-file Don’t read the site-wide Rprofile
–no-init-file Don’t read the user R profile
–restore Do restore previously saved objects at startup
–vanilla Combine --no-save, --no-restore, --no-site-file
–no-init-file and --no-environ

接着要做的事情，就是检查一下新安装到的包会存放在哪里？权限是否有？
在目录下找到了lib64的文件夹（就是不知道这个文件夹是否是存放R包的那个文件夹），修改其权限为可读可写可操作。

来重新看一下，Rscript的路径：

/home/xxzhang/miniconda3/bin/Rscript

现在，在这边的准备工作基本完成。接下来，就是去window那边，修改somatic.pl的文件，然后重新运行代码文件。
修改记录：

178行：$annovar_path="/home/xxzhang/workplace/QBRC/annovar/";
Rscript添加路径
samtools的路径更换为bcftools
495行添加strelka的路径

以上，保存后，重新上传。
接下来，重新运行代码，看是否有报错。
切换至Linux系统。

程序正在运行中。不知道会不会报错。（9：28）
现在就是等待，在等待的过程中，自己可以看一下，cnv.pl的代码，看一下，这个文件的作用是什么？
（9：50）运行结束，大概用了20分钟，新产生的文件比之前产生的文件多了更多的内容。
简直太给我面子了，有那么多的错误，这个程序居然还可以如此坚强的运行下去。前进的每一步都是艰难的，但是前进的每一步都是有意义的。

出错1：
/home/xxzhang/miniconda3/bin/bcftools mpileup -d 7000 -I -f ./geneome/hg19/hg19.fa ./output/tumor/tumor.bam > ./output/tumor/tumor.mpileup

[mpileup] 1 samples in 1 input files
Error: mpileup has failed!

其实，上一步骤就出现了一些问题，产生的index.out几个文件都是0，tumor.bam确实是有内容的。为什么使用gatk做这一步呢？这个过程，对于我们的体系而言是有意义的吗？必须要有结果才“正确”吗？

 # indel realignment
    if ($pdx=~/mouse/) {$known=" -known ".$resource_dbsnp." ";} else {$known=" -known ".$resource_mills." -known ".$resource_1000g." ";}
    system_call("java -Djava.io.tmpdir=".$type_output."/tmp -jar ".$gatk." -T RealignerTargetCreator -R ".$index." --num_threads ".$thread.
      $known." -o ".$type_output."/".$type."_intervals.list -I ".$type_output."/dupmark.bam > ".$type_output."/index.out");

    system_call("java -Djava.io.tmpdir=".$type_output."/tmp -jar ".$gatk." -T IndelRealigner ".
      " --filter_bases_not_stored --disable_auto_index_creation_and_locking_when_reading_rods -R ".$index.$known.
      " -targetIntervals ".$type_output."/".$type."_intervals.list -I ".
      $type_output."/dupmark.bam -o ".$type_output."/realigned.bam >".$type_output."/".$type."_realign.out");

    unlink_file($type_output."/dupmark.bam");
    unlink_file($type_output."/dupmark.bam.bai");

    # base recalibration
    if ($skip_recal==0)
    {
      if ($pdx=~/mouse/) {$known=" ";} else {$known=" -knownSites ".$resource_mills." ";}
      system_call("java -Djava.io.tmpdir=".$type_output."/tmp -jar ".$gatk." -T BaseRecalibrator -R ".$index." ".
        " -knownSites ".$resource_dbsnp.$known." -I ".$type_output."/realigned.bam -o ".$type_output."/".
        $type."_bqsr > ".$type_output."/table.out");
      system_call("java -Djava.io.tmpdir=".$type_output."/tmp -jar ".$gatk." -T PrintReads -rf NotPrimaryAlignment -R ".
        $index." -I ".$type_output."/realigned.bam -BQSR ".$type_output."/".$type."_bqsr -o ".
        $type_output."/".$type.".bam > ".$type_output."/".$type."_recal.out");

      unlink_file($type_output."/realigned.bam");
      unlink_file($type_output."/realigned.bai");
      unlink_file($type_output."/".$type."_bqsr");
      unlink_file($type_output."/".$type."_intervals.list");

这个出错，我其实有些没有办法理解，是因为只有一个文件吗？

这个错误是我最不知道怎么办的了。

继续回到查看运行指令，发现这个流程，光会跑通是远远不够的，我希望自己能够理解他每一步为什么会这样做？
就有点像孔子学琴的故事：

孔子学古琴师襄子，十日不进。师襄子曰：“可以益矣。”孔子曰：“丘已习其曲矣，未得其数也。”有间，曰：“已习其数，可以益矣。”孔子曰：“丘未得其志也。”有间，曰：“已习其志，可以益矣。”孔子曰：“丘未得其为人也。”有间，有所穆然深思焉，有所怡然高望而远志焉。曰：“丘得其为人，黯然而黑，几然而长，眼如望羊，如王四国，非文王其谁能为此也！”师襄子辟席再拜，曰：“师盖云《文王操》也。”

学习，有一个逐渐深化的过程，我现在的理解还远远不够。
看到昨天的使用指南，我想，或许我可以去掉原始文件中的-I指令，因为在源文件中，根本就没有说有这个option。
不对，又确认了一下，是有的。

-I, --skip-indels do not perform indel calling

意思就是说，跳过indel，不进行indel calling。
通过解释，这行代码的意思是：
/home/xxzhang/miniconda3/bin/bcftools mpileup -d 7000 -I -f ./geneome/hg19/hg19.fa ./output/tumor/tumor.bam > ./output/tumor/tumor.mpileup

-d的参数是7000，即为了避免对内存的过度使用，规定每一个文件的最大的深度（depth）为7000。对于我们这个没有多大的意义吧。

单独的运行一行代码，查看运行结果是否顺利。
/home/xxzhang/miniconda3/bin/bcftools mpileup -f ./geneome/hg19/hg19.fa /home/xxzhang/workplace/output/dupmark.bam > /home/xxzhang/workplace/output/tumor.mpileup

运行成功，生成了结果文档，那就说明是我们的指令本身出现了问题，所以，程序才运行不起来。我先让他转起来，至于具体的生物学意义，我们再进行调整参数。

另外，加不加-I这个参数，对于结果文件还是有一定的差异的。前后，进行两次的实验操作。生成的结果文件如下所示：

tumor.mpileup 131,804,995
tumor_2.mpileup 131,702,722

所以，还是加上去这个参数更好（起码会有影响）。

最终修改代码为：
/home/xxzhang/miniconda3/bin/bcftools mpileup -I -f ./geneome/hg19/hg19.fa ./output/tumor/tumor.bam > ./output/tumor/tumor.mpileup

出错2：
/home/xxzhang/workplace/software/strelka/bin/configureStrelkaGermlineWorkflow.py --bam /home/xxzhang/workplace/QBRC/output/tumor/tumor.bam --referenceFasta ./geneome/hg19/hg19.fa --runDir ./output/strelka/ --exome

Can’t exec “/home/xxzhang/workplace/software/strelka/bin/configureStrelkaGermlineWorkflow.py”: Permission denied at somatic.pl line 523.

出错原因是：拒绝访问。接下来，就是修改这个文件夹的使用权限。这个并不难。

解决方案：

chmod -R 777 bin/
ll

total 2
drwxrwxrwx. 2 xxzhang Zhu_lab 6 6月 3 12:13 bin
drwxr-xr-x. 3 xxzhang Zhu_lab 1 6月 3 12:13 lib
drwxr-xr-x. 2 xxzhang Zhu_lab 29 6月 3 12:14 libexec
drwxr-xr-x. 5 xxzhang Zhu_lab 6 6月 3 12:14 share

cd bin
ll

total 26
-rwxrwxrwx. 1 xxzhang Zhu_lab 8537 6月 2 21:18 configureStrelkaGermlineWorkflow.py
-rwxrwxrwx. 1 xxzhang Zhu_lab 606 6月 2 21:18 configureStrelkaGermlineWorkflow.py.ini
-rwxrwxrwx. 1 xxzhang Zhu_lab 6085 6月 2 21:18 configureStrelkaSomaticWorkflow.py
-rwxrwxrwx. 1 xxzhang Zhu_lab 2580 6月 2 21:18 configureStrelkaSomaticWorkflow.py.ini
-rwxrwxrwx. 1 xxzhang Zhu_lab 3465 6月 2 21:18 runStrelkaGermlineWorkflowDemo.bash
-rwxrwxrwx. 1 xxzhang Zhu_lab 3384 6月 2 21:18 runStrelkaSomaticWorkflowDemo.bash

修改完成。

出错3：

installing source package ‘statmod’ …
** package ‘statmod’ successfully unpacked and MD5 sums checked
** libs
x86_64-conda_cos6-linux-gnu-gfortran -fpic -fopenmp -march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -pipe -I/home/xxzhang/miniconda3/include -fdebug-prefix-map==/usr/local/src/conda/- -fdebug-prefix-map==/usr/local/src/conda-prefix -c gaussq2.f -o gaussq2.o
make: x86_64-conda_cos6-linux-gnu-gfortran: Command not found
make: *** [gaussq2.o] Error 127
ERROR: compilation failed for package ‘statmod’
removing ‘/home/xxzhang/miniconda3/lib/R/library/statmod’

安装包出错，主要原因是在运行x86_64-conda_cos6-linux-gnu-gfortran这行指令的过程中，显示没有这个指令，导致这个包的编译过程出了问题。

解决方案：

查了一下，系统中确实有这个指令的存在，但是，链接不到这里，估计并没有存放绝对路径。
whereis x86_64-conda_cos6-linux-gnu-gfortran

x86_64-conda_cos6-linux-gnu-gfortran: /opt/anaconda3/bin/x86_64-conda_cos6-linux-gnu-gfortran /opt/anaconda3/bin/x86_64-conda_cos6-linux-gnu-gfortran.bin

那么，我们只需要按照老办法，将其放在我们原始代码文件中的绝对路径即可（这种方法的平台依赖很大，甚至只适用于我们的系统，但是能用就行，以后能够跑起来的话，就可以尝试扩大它的平台的灵活性）。

/opt/anaconda3/bin/x86_64-conda_cos6-linux-gnu-gfortran

这个好像暂时还不行。
我找不到这个程序运行，调用的是哪里的？系统自身的吗？
明明系统中是有这个插件的，为什么检测不到呢？我们可以先把包安装完成，然后到时候已经存在，是不是就可以不用它自己安装了?

出错4：
/home/xxzhang/workplace/QBRC/annovar/table_annovar.pl ./output/germline_mutations.txt /home/xxzhang/workplace/QBRC/annovar//humandb/ -buildver hg19 -out ./output/germline_mutations -remove -protocol refGene,ljb26_all,cosmic70,esp6500siv2_all,exac03,1000g2015aug_all -operation g,f,f,f,f,f -nastring

NOTICE: the --polish argument is set ON automatically (use --nopolish to change this behavior)
Error: the required database file /home/xxzhang/workplace/QBRC/annovar/humandb/hg19_ljb26_all.txt does not exist.

是什么问题呢？这里，我检查了一下我的数据库，发现就是hg19_ljb26_all.txt这个文件没有下载完全，显示的是filepart,还是我可以重新将其命名为txt.
因为查看实验记录：

(base) zxx@zxx-Lenovo-Yoga710-14ISK:/media/zxx/TOSHIBA/QBRC/annovar$ perl annotate_variation.pl -buildver hg19 -downdb -webfrom annovar ljb26_all humandb/
NOTICE: Web-based checking to see whether ANNOVAR new version is available … Done
NOTICE: Downloading annotation database http://www.openbioinformatics.org/annovar/download/hg19_ljb26_all.txt.gz … OK
NOTICE: Downloading annotation database http://www.openbioinformatics.org/annovar/download/hg19_ljb26_all.txt.idx.gz … OK
NOTICE: Uncompressing downloaded files
NOTICE: Finished downloading annotation files for hg19 build version, with files saved at the ‘humandb’ directory

我明白了，并不是我没有下载完全，是我将我的这个文件，上传到服务器的时候，并没有上传完全。实际上，这个文件有11GB，而我这边只有8GB，我实在是太蠢笨了，计算机要比我严谨，敏锐的多。我知道怎么处理了。

最后，将本次操作运行产生的结果文件，转移到备份文件夹中，好有一个结果的跟踪记录。

转移到window界面下，继续编辑指令，上传缺失的文件。

~~上传缺失文件。~~
上传完成。
~~修改bcftools指令。~~
~~添加指令的位置。~~
选择手动安装R包。
install.packages("statmod")

回到Linux平台下。重新运行指令，又是一堆的错误。

/home/xxzhang/workplace/software/strelka/bin/configureStrelkaGermlineWorkflow.py --bam /home/xxzhang/workplace/QBRC/output/tumor/tumor.bam --referenceFasta ./geneome/hg19/hg19.fa --runDir ./output/strelka/ --exome

bin/sh: /home/xxzhang/workplace/software/strelka/libexec/htsfile: Permission denied

修改权限：chmod -R 777 strelka/
权限修改完成。

还有的那个，就是老问题。

x86_64-conda_cos6-linux-gnu-gfortran -fpic -fopenmp -march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -pipe -I/home/xxzhang/miniconda3/include -fdebug-prefix-map==/usr/local/src/conda/- -fdebug-prefix-map==/usr/local/src/conda-prefix -c gaussq2.f -o gaussq2.o
make: x86_64-conda_cos6-linux-gnu-gfortran: Command not found
make: *** [gaussq2.o] Error 127
ERROR: compilation failed for package ‘statmod’

removing ‘/home/xxzhang/miniconda3/lib/R/library/statmod’
The downloaded source packages are in
‘/tmp/RtmpQAtxv3/downloaded_packages’
Updating HTML index of packages in ‘.Library’
Making ‘packages.html’ … done
Warning message:
In install.packages(name_pkg[bool_nopkg]) :
installation of package ‘statmod’ had non-zero exit status
Error in FUN(X[[i]], …) : there is no package called ‘statmod’
Calls: lapply -> FUN
Execution halted

我来修改指令，重新运行。

遇到了一个问题，挺严重的，一时不知道该如何解决。
我来简单的描述一下这个问题：
（1）安装某个非常重要的R包的时候，显示编译不成功。主要原因是缺少“依赖”——

make: x86_64-conda_cos6-linux-gnu-gfortran: Command not found

在编译的过程中，出现错误。显示缺乏这个指令。我不知道R环境中，R包一般编译的位置是在哪里？我去网上提问：R的编译器放在什么位置？不是自己的电脑实在是太麻烦了。
所以，我选择在外面安装并编译好R包，然后移动到文件夹中，必要的话，可以修改原始代码。
刚开始的时候，选择的是同济的安装镜像，安装不了。很困惑。后来，选择清华的镜像，就安装成功了。
options("repos" = c(CRAN="http://mirrors.tuna.tsinghua.edu.cn/CRAN/"))
install.packages("statmod")

trying URL ‘http://mirrors.tuna.tsinghua.edu.cn/CRAN/src/contrib/statmod_1.4.36.tar.gz’
Content type ‘application/x-gzip’ length 93897 bytes (91 KB)
==================================================
downloaded 91 KB

installing source package ‘statmod’ …
** package ‘statmod’ successfully unpacked and MD5 sums checked
** using staged installation
** libs
x86_64-conda_cos6-linux-gnu-gfortran -fno-optimize-sibling-calls -fpic -fopenmp -march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -ffunction-sections -pipe -isystem /home/zxx/miniconda3/include -fdebug-prefix-map=/tmp/build/80754af9/r-base_1589917437985/work=/usr/local/src/conda/r-base-3.6.1 -fdebug-prefix-map=/home/zxx/miniconda3=/usr/local/src/conda-prefix -c gaussq2.f -o gaussq2.o
x86_64-conda_cos6-linux-gnu-cc -I"/home/zxx/miniconda3/lib/R/include" -DNDEBUG -DNDEBUG -D_FORTIFY_SOURCE=2 -O2 -isystem /home/zxx/miniconda3/include -I/home/zxx/miniconda3/include -Wl,-rpath-link,/home/zxx/miniconda3/lib -fpic -march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -ffunction-sections -pipe -isystem /home/zxx/miniconda3/include -fdebug-prefix-map=/tmp/build/80754af9/r-base_1589917437985/work=/usr/local/src/conda/r-base-3.6.1 -fdebug-prefix-map=/home/zxx/miniconda3=/usr/local/src/conda-prefix -c init.c -o init.o
x86_64-conda_cos6-linux-gnu-cc -shared -L/home/zxx/miniconda3/lib/R/lib -Wl,-O2 -Wl,–sort-common -Wl,–as-needed -Wl,-z,relro -Wl,-z,now -Wl,–disable-new-dtags -Wl,–gc-sections -Wl,-rpath,/home/zxx/miniconda3/lib -Wl,-rpath-link,/home/zxx/miniconda3/lib -L/home/zxx/miniconda3/lib -L/home/zxx/miniconda3/lib -Wl,-rpath-link,/home/zxx/miniconda3/lib -o statmod.so gaussq2.o init.o -lgfortran -lm -lgomp -lquadmath -lpthread -L/home/zxx/miniconda3/lib/R/lib -lR
installing to /home/zxx/miniconda3/lib/R/library/00LOCK-statmod/00new/statmod/libs
** R
** data
** inst
** byte-compile and prepare package for lazy loading
** help
*** installing help indices
** building package indices
** testing if installed package can be loaded from temporary location
** checking absolute paths in shared objects and dynamic libraries
** testing if installed package can be loaded from final location
** testing if installed package keeps a record of temporary installation path
DONE (statmod)
The downloaded source packages are in
‘/tmp/Rtmp2nx94T/downloaded_packages’
Updating HTML index of packages in ‘.Library’
Making ‘packages.html’ … done

安装好了之后，我将其移到服务器中指定的文件夹中。

虽然能够加载，但是出现了新的问题：
library(statmod)

Error: package or namespace load failed for ‘statmod’ in rbind(info, getNamespaceInfo(env, “S3methods”)):
number of columns of matrices must match (see arg 2)
In addition: Warning message:
package ‘statmod’ was built under R version 3.6.1

所以，这件事情，最核心的问题应该是什么？就是筛选的这个指令无法运行对吧？影响大吗？

perl somatic.pl NA NA ./example/example_dataset/sequencing/SRR7246238_1.fastq.gz ./example/example_dataset/sequencing/SRR7246238_2.fastq.gz 32 hg19 ./geneome/hg19/hg19.fa /home/xxzhang/workplace/software/java/jdk1.7.0_80/bin/java ./output human 1 ./disambiguate_pipeline

环境的问题，很难搭建好。

重新运行之后，到strelka这一步都是没有问题的。接下来，就是加载R包的问题。
/home/xxzhang/miniconda3/bin/Rscript /home/xxzhang/workplace/QBRC//somatic_script/filter_vcf.R ./output NA

Error: package or namespace load failed for ‘statmod’ in rbind(info, getNamespaceInfo(env, “S3methods”)):
number of columns of matrices must match (see arg 2)
In addition: Warning message:
package ‘statmod’ was built under R version 3.6.1
Execution halted

下面的报错，都是源于，缺少这一步产生的结果文件。所以，问题的根源还是在于：
（1）如何更新R版本？
（2）错误加载的原因是什么？

参考链接：https://www.it1352.com/1721382.html
出现这个错误，是不是说明，这个包是损坏的？（或者我这样“武力”的移植，本身就是不契合其实际的操作环境的？）
最本质的方法，还是跟R有关。

还是要想办法说，怎样更新我的R版本到（3.6.1），以及怎样解决编译的过程中，依赖库缺失的问题（明明系统上是存在的，但是自动编译的时候，却显示缺失，我一直想不明白，大概是没有把这个放在绝对的根目录中）。

如果，我能够在九点之前解决这个问题，今天就是成功的。
（19：06）

 ./conda create -n R3.6
 ./conda activate R3.6

但是，当我想要更新R环境的时候，系统说这个shell版本太老了。需要更新。

CommandNotFoundError: Your shell has not been properly configured to use ‘conda activate’.
To initialize your shell, run
$ conda init <SHELL_NAME>
Currently supported shells are:

bash
fish
tcsh
xonsh
zsh
powershell
See ‘conda init --help’ for more information and options.
IMPORTANT: You may need to close and restart your shell after running ‘conda init’.

那我就更新吧。
./conda init bash

-bash-4.2$ ./conda init bash
no change /home/xxzhang/miniconda3/condabin/conda
no change /home/xxzhang/miniconda3/bin/conda
no change /home/xxzhang/miniconda3/bin/conda-env
no change /home/xxzhang/miniconda3/bin/activate
no change /home/xxzhang/miniconda3/bin/deactivate
no change /home/xxzhang/miniconda3/etc/profile.d/conda.sh
no change /home/xxzhang/miniconda3/etc/fish/conf.d/conda.fish
no change /home/xxzhang/miniconda3/shell/condabin/Conda.psm1
no change /home/xxzhang/miniconda3/shell/condabin/conda-hook.ps1
no change /home/xxzhang/miniconda3/lib/python3.6/site-packages/xontrib/conda.xsh
no change /home/xxzhang/miniconda3/etc/profile.d/conda.csh
modified /home/xxzhang/.bashrc
> For changes to take effect, close and re-open your current shell. <

接下来，就是重启shell,但是我不知道如何重启shell。

1 hour later ……
尝试了好多方法都不行。

20:32
决定下载R。重新安装全新的R版本，不用那么麻烦了，要有壮士断腕的决心。
./conda uninstall R
因为之前的那个R版本是3.5,不再适用于我们的这个statmod这个包。

出现了新的问题：

Found conflicts! Looking for incompatible packages.

我好像最近尝试这个那个变得一团糟，难道还有从新的安装conda吗？心累。

20：50

我需要休息了。

今天也是个妖精头子呀

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
实验记录 | 6/5

今天早上（8：57）过来，发现R已经安装完成了，非常棒。./RscriptUsage:/path/to/Rscript [–options] [-e expr [-e expr2 …] | file] [args]–options accepted are–help Print usage and exit–version Print version and exit–verbose Print information on
复制链接

扫一扫