formatファイルをダウンロードして、Rで読み込む方法を解説します。アノテーション情報を含むVCFファイルをダウンロードする解析結果の画面(前回の内容参照)で、設定部分の右側"Download"(下図)から、"VCF"をクリックします。ダウンロードされるファイルを別名で保存します。ここでは"MHC_test4.var.annotated.vcf"とします。このファイルをCotEditorなどで開くとわかりますが、このファイルは改行コードがCR/LFになっているので、Terminal等で開くことも考えて、LFに変換しておきます。file"のような絵(上図)をクリックすると、.csv形式でファイルが保存され、Excelで開くことができます(Excelを使う場合、これも悪くない)。アノテーション付き".vcf"ファイルをRで読み込むダウンロードしたVCFファイル"MHC_test4.var.annotated.vcf"は、Rで読み込むことができます。Rに読み込む前に一手間かける必要があります。"MHC_test4.var.annotated.vcf"をCotEditorで開きます。「検索」メニューから「検索...」を選択します。(3)「検索と置換」ウィンドウの「正規表現」をチェックし、検索文字列の窓(上)に"\t"と入力、置換文字列の窓(下)に"|"と入力します(下図)。「すべて置換」をクリックして、"MHC_test4.var.annotated.vcf"のすべてのtabを"|"に変換します。ファイルを上書き保存します。改行コードがLFになっていることを確認。「作業ディレクトリの変更...」で、"MHC_test4.var.annotated.vcf"のあるディレクトリを選択します。以下のコマンドで、"MHC_test4.var.annotated.vcf"の内容をRに読み込みます。sep="|")エラーが出なければ、以下のコマンドで読み込まれていることを確認します。10)-データを読み込めれば、Rを使ってTableの必要な部分のみを抜き出したり、特定の遺伝子や特定のタイプの変異を選択的に得ることができます。V9=="missense_variant")とします。条件に一致するSNP数をカウントするにはlength(HLA_DPB1_missense[,1])とします。(4)月別アーカイブはてなブログをはじめよう!hashiyukiさんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。