読者です 読者をやめる 読者になる 読者になる

疫学と医療統計学と遺伝学と時々、大学院生活

疫学を専門とする大学院生の研究に関する備忘録的ページ。

R package"VennDiagram"でベン図を書く

今日はベン図を描きたいと思います。 (今回は例として、東海三県での多い名字30位についてベン図を書いてみようと思う) パッケージの読み込み library(VennDiagram) まずはベン図で書きたいそれぞれのカテゴリに含まれる要素を書き出す aichi <- c("鈴木","…

R for beginners vol.4 「データの要約とビジュアライゼーション」

今回はデータの整頓から少し分析に近いことを始めます。その中でデータを要約し、図示することが解析の一歩かと思い今回の内容にしています。R for beginners vol.1 「Rの紹介と基本的なコマンド」 jojoshin.hatenablog.comR for beginners vol.2 「データの…

R 役に立つ小技5選 vol.3

今回は第3回目となる『R役に立つ小技5選』についての記事です。 前回と前々回の記事は下にリンクがあるので、ご参照ください。 No.1 Macを使用して、Rでグラフを日本語で書くと、文字化けするのを解消する方法 par(family="Osaka") hist(データ名, xlab="x軸…

R package 'estatapi'を使用して、公開されている政府統計を取得、解析する

今回は先日公開されたばかりのRパッケージ'estatapi'を使用して、政府統計を取得し、可視(グラフ)化する。 ちなみに参照する統計表は様々な行動の時間などを収録している社会生活基本調査ということのみ決めておく。 (注意)横に大きなデータは見づらくて…

R package 'ggmap'を用いたマッピング

今回は'ggmap'と呼ばれるパッケージを使用して、愛知県にあるナゴヤドームへの経路(最寄駅(大曽根駅)から)を図示することにする。 1. 早速、パッケージのインストールと読み出す install.packages("ggmap") library(ggmap) 2. 地名で緯度経度をグーグル…

QCの実行

今回は先日書いたGWASのQCに関して、『Primer to Analysis of Genomic Data using R』の第3章のサンプルデータを使用して実施する。 はじめに ここで使用するデータはヒトではなく、83頭の羊54,977SNPsを使用してトレーニングする。今回のSNPsはイルミナ社の…

R package 'maptools'と'spsurvey'を利用して、カラーマップを作成する。

今回は様々な指標をカラーマップを使用して、地図上に見やすく可視化するコードについてまとめる。 (特別なパッケージとしてはmaptoolsとspsurveyが挙げられる。) 今回は事例として、愛知県の男女の平均寿命を地図上にプロットすることを最終目的とする。 …

R package 'qqman'の基本的な使用法

今回はRのqqmanというパッケージの基本的な使い方をまとめる。 【ページ下部に追記あり(20160922)】 概要 現在、ゲノム全体10万から50万ほどのSNPsと疾患との関連を探索するGWASが可能になっている。そこでその結果を可視化する際に用いるのがマンハッタン…

R package 'fitbitScraper'を使用したVisualization

今回はウェアラブルデバイス業界で先端を走るfitbitのfitbit flexのデータをRのパッケージ'fitbitScraper'を使用して解析する。まずはじめに必要になるパッケージを読み込む。 (インストールしてあることを前提とする。) library(fitbitScraper) library(l…