![[Toho University]](image/CLogo_Kihon_process.jpg "[Toho University]")
|
ノート/DRY解析教本_2017-02-12http://pepper.is.sci.toho-u.ac.jp/pepper/index.php?%A5%CE%A1%BC%A5%C8%2FDRY%B2%F2%C0%CF%B6%B5%CB%DC_2017-02-12 |
![[YYLAB]](image/yamanouc2.png "[YYLAB]")
|
|
|
ノート/DRY解析教本_2017-02-12http://pepper.is.sci.toho-u.ac.jp/pepper/index.php?%A5%CE%A1%BC%A5%C8%2FDRY%B2%F2%C0%CF%B6%B5%CB%DC_2017-02-12 |
|
ノート
マイニング
├マイニング
├バスケット解析
├バスケット解析をＲで
├図書貸出をＲで
├モデルクラスタリング
└
テキストマイニング
├テキストマイニングTM
├TMとシソーラス
├PubMedをTM
├TMとMeSH
├TMとNLTK
├テキストマイニングとtagger
├医薬品添付文書ＤＢ
├論文の処理1
├剽窃１
├Bing検索
└
分子進化学
├分子進化学
└
Pythonと論文アクセス
├Pythonを使ってみる
├PythonでPubMed
├Pythonで...続き
└
CUDA
├GPU
├ProgrammingGuide
├Selandメモ
├数値積分
└
ACS、PHPからLDAP
├ACS導入
├新規ホスト移行
├再度やり直し
├OpenLDAP導入
└
Linuxでビデオ
├Linuxでビデオデータを作る
├AVIフォーマット
├x264ソースを見る
├Macでフレーム抜き
└
MPIと並列化
├並列処理の実用化
├Fedora10でMPICH2
├EPMRを見る
├EPMRを見る２
├EPMRを見る〜Ｃ＋＋２
├EPMRを見る〜Ｃ＋＋２
├EPMRを見る〜Ｃ＋＋３
├EPMRを見る〜Ｃ＋＋４
├CASP7と日本人
└
Macメモ
├Script
自信ないので
├Winのバッチ
│├スクリプト入門
│├バッチファイル
│├SETコマンド
│├自分のメモ
│└
└
ノート/ノート
423 2017-02-12 (日) 17:03:19
DRY解析教本の練習のコメント
ginger環境で実験
yum install R
$ wget https://download2.rstudio.org/rstudio-server-rhel-1.0.136-x86_64.rpm $ sudo yum install --nogpgcheck rstudio-server-rhel-1.0.136-x86_64.rpm
sudo rstudio-server {start, stop, restart}cd ~/src wget http://ccb.jhu.edu/software/tophat/downloads/tophat-2.1.1.Linux_x86_64.tar.gz tar xvf tophat-2.1.1.Linux_x86_64.tar.gz sudo ln -s ~/src/tophat-2.1.1.Linux_x86_64/tophat2 /usr/local/bin/tophat2
cd ~/src wget http://cole-trapnell-lab.github.io/cufflinks/assets/downloads/cufflinks-2.2.1.Linux_x86_64.tar.gz tar xvf cufflinks-2.2.1.Linux_x86_64.tar.gz sudo ln -s ~/src/cufflinks-2.1.1.Linux_x86_64/cufflinks /usr/local/bin/cufflinks
source("http://bioconductor.org/biocLite.R")
biocLite("cummeRbund")sudo ln -s ~/src/FastQC/fastqc /usr/local/bin/fastqc
レファレンスデータのダウンロード��
wget -r http://shujunsha.com/NGS_DAT/Lv2_2/iGenome/Homo_sapiens/
はうまくいかなかった。問題はrobot.txtがあるためらしく、/etc/wgetrc上で robots=off を指定 するとOKとなった。
105ページからの記述は、別のデータを必要としている。
wget ftp://ftp.sra.ebi.ac.uk/vol1/fastq/ERR266/ERR266335/ERR266335.fastq.gz wget ftp://ftp.sra.ebi.ac.uk/vol1/fastq/ERR266/ERR266349/ERR266349.fastq.gz wget ftp://ftp.sra.ebi.ac.uk/vol1/fastq/ERR266/ERR266337/ERR266337.fastq.gz wget ftp://ftp.sra.ebi.ac.uk/vol1/fastq/ERR266/ERR266351/ERR266351.fastq.gz wget ftp://ftp.sra.ebi.ac.uk/vol1/fastq/ERR266/ERR266338/ERR266338.fastq.gz wget ftp://ftp.sra.ebi.ac.uk/vol1/fastq/ERR266/ERR266347/ERR266347.fastq.gz
106ページにある、SRAtoolkitのダウンロード。これはSRAファイルをダウンロードした場合のfastqファイルを取り出すツール。
本ではHomebrewでインストールしている
NCBI SRA Toolkit
NCBI SRA Toolkitの使い方
[[Linux用パッケージをダウンロード:
wget https://ftp-trace.ncbi.nlm.nih.gov/sra/sdk/2.8.1-3/sratoolkit.2.8.1-3-centos_linux64.tar.gzw
展開した後、ディレクトリノードを/usr/local/sratoolkitにmvしておいて、~/.bash_profile中に
PATH=/usr/local/sratoolkit:$PATH export PATH
まず、出力を置くためのディレクトリを作っておく。./FastQCとする。解析は、
fastqc --nogroup -o ./FastQC ERR266335.fastq
まず、fastx_toolkitをインストール。
Linux用パッケージをダウンロード。
展開するとbinディレクトリができるので、、そのbinノードを /usr/local/fastx_toolkit/bin にmvしておいて、~/.bash_profile中に
PATH=/usr/local/fastx_toolkit:$PATH export PATH
トリミングは
fastq_quality_trimmer -Q 33 -t 20 -l 30 -i ERR266335.fastq -o ERR266335_trim.fastq
tophat -p 32 -G Homo_sapiens/NCBI/build37.2/Annotation/Archives/archive-2014-06-02-13-47-29/Genes/genes.gtf -o tophat_results/ERR266335_P0 ./ERR266335_trim.fastq
