統計解析ソフトR
Rに読み込んだVCFに、アノテーションを付加しましょう。 staffblog.amelieff.jp VariantAnnotationの公式説明書はこちらです。 アノテーションの付加には、Bioconducterのアノテーションデータパッケージを利用します。 パッケージ一覧はこちらから見られま…
VCF(Variant Call Format)をRで読みましょう！ Rパッケージ "VariantAnnotation" をご紹介します。 VariantAnnotationはVCFを読み込み、更にアノテーションパッケージを利用してアノテーションを付与することが可能です。 列の分類やPOS表記に少しクセがあり…
SNP解析、GWASの結果の図示に広く使われるマンハッタンプロット。 これをカラフルにする方法、色の例をご紹介します。 マンハッタンプロットの作図には、統計解析ソフトRの、CRANパッケージ qqmanを利用すると簡便です。 install.packages("qqman") #パッケ…
kimoton です。 皆様休日は何をされてますかね。 わたくしkimoton、最近料理をしているのですが、書店でこのようなものを見つけてしまいました。 週末2時間で作れる! 1週間分のおかず (ei cooking)作者: 結城寿美江,ei cooking編集部出版社/メーカー: エイ出…
hr-kです。今回はRでvcfファイルを扱う時の話題で書きたいと思います。 Rでvcfファイルを扱うソースを書いており、気付いたことが一点。 Rでvcfを開くとき、vcfR(説明書のリンク)というライブラリにあるread.vcfRを使うのがすごく便利です。 具体的な使用方…
Rで、ベクトルなどに含まれる重複した値をユニークにしたいとき、unique()を使用します。 > a [1] "AAA" "BBB" "CCC" "AAA" > unique(a) [1] "AAA" "BBB" "CCC" linuxコマンドの、隣り合っている重複行をユニークにするuniqコマンドとよく混同して、uniqと書…
Bioconductorのseq2pathwayパッケージを使ってパスウェイのアノテーションを行います。 ※seq2pathwayのインストールはこちらの記事をご覧ください。 seq2pathwayパッケージに含まれているテストデータを使ってテストしてみましょう。 library(seq2pathway) d…
ChIP-seqやRIP-seqでピークが得られた後は、遺伝子をアノテーションしたり、それらの遺伝子がどんなパスウェイに関連するか調べたりしたいですね。 Bioconductorのパッケージseq2pathwayを使って、ゲノム領域に遺伝子をアノテーションしたり、遺伝子にパスウ…
バイオインフォには関係ないのですが、少しつまづいたことです。 ReactomePAを使っています。 hatさんがかかれた紹介記事では、enrichPathway()でエンリッチされているパスウェイを取得した結果を棒グラフに描画していますが、summary()することでデータフレ…
短い間ですが、弊社のトレーニングのTAを担当させていただいていた時期があります。 RのTAをやるなかで、とても便利なのに、初めてRに触れる方が戸惑われることが多いと感じたポイントのひとつに paste 関数があります。 今回は paste 関数をあれこれいじく…
heatmap関数のRowSideColors/ColSideColorsオプションの紹介です。 gplotsパッケージのheatmap.2やheatmap3パッケージでも同様の機能があります。 ヒートマップのデンドログラムの横に、カラフルなサイドバーを付ける機能です。 テスト用のマトリクスを作り…
Rで解析をしていると、データフレーム（やマトリクスやベクトル）から、条件を満たすデータだけを抜き出したいことがよくあります。 そんなときはsubset関数を使います。 例として、みんな大好きirisを使います。 > head(iris) Sepal.Length Sepal.Width Pet…
Rにはデフォルトでもヒートマップを作成するheatmap関数がありますが、私はheatmap関数より多機能なgplotsパッケージのheatmap.2関数をよく使います。 基本的な使い方はheatmap関数と同じです。 gplotsパッケージをインストールしていない場合はインストール…
ある環境で使用していたとあるBioconductorパッケージを、別の環境でも使うためインストールしようとしたらエラーが出ました。 具体的にはggbioとbiovizBase、他にもhatさんご紹介のReactomePAもダメでした。 エラー文を見てみると、「RSQLite」パッケージの…
以前、Rのgrep()についての記事で、grep()は文字列のベクトルにのみ使えるとご紹介しましたが、今回はマトリクスやデータフレームではどうしたらよいか？という記事です。 例として以下のような文字列のマトリクスを用意しました。 1列目が季節、2列目が英語…
Rのオブジェクトから、特定のキーワードを含む要素だけ抜き出したいときはどうしたらいいでしょうか。 ベクトルの場合はgrep()関数が使用できます（grepの仲間にgrepl、regexpr、gregexpr、regexec、検索後に置換を行うsub、gsubがあります）。 このgrepは文…
前回はbiomaRtパッケージを使って、「ADH1B」「ALDH2」遺伝子にGeneOntologyのアノテーションをつけました。 今回は、RamiGOパッケージを使って、それらのGeneOntologyIDがどのような関係にあるか、図示してみたいと思います RamiGOパッケージをインストール…
＿★ ＿▲ ▲▲▲ R Advent Calendar 2013に参加させていただいております。 ＿■ Bioconductorは、バイオ解析関連のRパッケージを集めたサイトです。 750個以上（2013年12月時点）の有用なパッケージが登録されています。 複数の遺伝子の特徴の傾向を見たい時に、…
ベクトルやリストを並び替えるときはsort()関数を使いますが 行列（マトリクス）やデータフレームを、特定の列をキーにして 並び替える場合はorder()という関数を使います。 同じジュースが、250ml・100円、350ml・120円、180ml・80円で 売られています。 1…
家系解析を行っている方は、家系図を描く機会が多いのでは ないでしょうか。 kinship2というパッケージを使うと、Rで家系図を描くことができます。 （１）家系図情報ファイルの作成 以下の内容のテキストファイルを作成します。 id father mother sex name 1…
Rの非常に有名なテストデータの一つ「iris」は、 アヤメの花弁の長さ・幅、ガクの長さ・幅のデータです。 ※アヤメのガクは、黄色い筋が入っていて外に垂れ下がっている部分です。 setosa（セトサ）種、versicolor（バージカラー）種、 virginica（バージニカ…
毎日寒いですね。 家ではこたつから出られなくなって、こたつの中で生活している方もいらっしゃるのではないでしょうか。 今日は、RでNCBI SRAのデータを検索するRパッケージ「SRAdb」をご紹介します。 最近出たばかりのSRAdbの論文はこちらです。 Zhu Y, St…
hatです。 先週土曜日、2012年度統計数理研究所共同研究集会「データ解析環境Rの整備と利用」の「国内ユーザーによる報告」に参加してきました。 普段はバイオの話ばかりなので、疫学やTwitter解析や感染症シミュレーションなど、他分野のお話がたくさん伺え…
明日、理化学研究所で開催される “Rでつなぐ次世代オミックス情報統合解析研究会” に参加してきます！ （詳細はコチラ） NGSの解析に特化したRを学べる事は、大変嬉しいです。 パワーアップをめざします！
箱ひげ図の概念に対し、ご指摘いただきました。 誠に申し訳ありません。訂正し、御詫び申し上げます。 またご指摘いただきまして、誠にありがとうございます。 下記に訂正、追加説明を行います。 箱ひげ図の概念図（訂正版）をFigure1に示します。 Figure1 …
昨日は、中央値について記述しました。 本日は、Rを用いた中央値の算出方法と箱ひげ図の作成について記述します。 よろしくお願いいたします。 例としまして 21歳, 23歳, 23歳, 25歳, 22歳, 60歳の方々の中央値を算出します。 中央値を算出する際は、median(…
最近鼻血がよく出るタノです。 先日Rスクリプトを作成するためのツール“Tinn-R”を紹介しました（参照）。 本日は、“Tinn-R”とRを用いてちょっと簡単なグラフ作成をしてみたいと思います。 よろしくお願いします。 突然ですが、夏といえば・・・・ダイエット…
こんにちは。 本日より、R漬けの日々を過ごします、タノです。 学生時代にやっていたとはいえ、大量のデータに四苦八苦しております。というわけで、本日は、私が愛用しているRスクリプトを書くためのR code editor、“Tinn-R”を紹介します。 Tinn-Rは、R言語…
R 「関数textとlegend」では、文字列と凡例を書き込みました。 書きこんだ文字列が見にくいようなので、位置を調節します。 太字の部分を書き足しました。 # デバイス領域 png("110131_iris.png", height=700, width=500) # 作図領域 par(mar=c(4,4,2,6),xpd…
前回のおさらいとして、関数textで「○」をプロット領域に書きます。 次に、プロット領域の外に凡例を書き込みます。 前回のプログラムに、太字の部分だけ追加しました。 # デバイス領域 png("110114_iris.png", height=700, width=500) # 作図領域 par(mar=c…
