epidemiology
データの可視化に関する書籍＋論文を簡単なレビューとともに紹介したいと思います（2021年9月1日現在）。
gtsummaryパッケージのtbl_merge関数についてまとめています。
gtsummaryパッケージのtbl_stack関数についてまとめています。
前回は生存解析においてベースラインだけでなく、縦断的に測定している（繰り返し測定、反復測定値ともいう）値の影響も組み込んだ「時間依存性共変量」についてまとめました。今回は、さらに進化した（個人的見解です）Joint modelについて、簡単な紹介とR…
Rのsurvivalパッケージを使用した生存解析について、時間依存性共変量を組み込んだモデルと通常モデルの比較を行なっています。
今回は、前回に引き続いて「Epigenetic clock」についてUCLAのDr. Hovarthによるレビューの後半をまとめていきたいと思います。主な内容は、「年齢的な変化との関連」や、「組織・細胞レベルとの関連性」、「介入による若返りは可能？」、「この方法における…
今回は、「Epigenetic clock」についてUCLAのDr. Hovarthによるレビューをまとめておきたいと思います。この研究がどのような未来をもたらし得るかみてみましょう。その前編として、「Epigenetic clock」はどのような学術的背景から注目されたのか、そしてど…
かなり久方ぶりの記事になりますが、今回はJournal of Clinical Epidemiologyに先日アクセプトになった”A systematic review shows no performance benefit of machine learning over logistic regression for clinical prediction models”について、簡単に…
今回は、遺伝疫学でも近年注目されているmendelian randomization（以下：MR）について、情報を集約し簡単なまとめをしておきます。MRは、経済学や社会調査の解析において利用されてきた操作変数（instrumental variable）を導入した解析方法の一種でありま…
今回も先ほど公開した記事と同様に、『入門機械学習による異常検知-Rによる実践ガイド-』（コロナ社、井出剛著、2015）の中にある「7.1 近傍法による異常検知」ということで、この内容をなぞり、最後に時系列データ解析の際に重要になる窓幅（どれだけの地点…
今回は、『入門機械学習による異常検知-Rによる実践ガイド-』（コロナ社、井出剛著、2015）の「6.3 リッジ回帰と異常検知（P165-172）」の内容をもとに政府・官公庁データをもとにして作成した都道府県別の10万人あたりの自殺者数を予測するモデルを作成し、…
今回は、文献の紹介をするとともに、個人の遺伝子検査サービス前後で消費者の意識がどのように変わっているかを紹介したい。論文のタイトルはDirect-to-Consumer Genetic Testing: User Motivations, Decision Making, and Perceived Utility of Resultsであ…
あけましておめでとうございます（疫学会や予防早期医療創成センターのワークショップ等への参加もあり、一ヶ月ぶりの更新です）。 今回はRのパッケージ"BLR"を用いて、ベイジアンラッソを実行する。ちなみに参考にするのは、Cedric Gondroらによる『Genome …
今回はSpringerから出ている『Genome-Wide Association Studies and Genomic Prediction』という本からデータの読み込みについて記述する。GWASなどのゲノムデータは非常に膨大であり、読み込みだけでもかなりの時間を要する。そこで下記のTipsで幾分改善さ…
今回は、9月22日の記事に加えて、共変量で調整する（その他の細かい設定もする）GWASについてコマンドを記しておく。 jojoshin.hatenablog.com 最も簡単なコマンド plink --bfile データ名 --assoc --out アウトプットするデータ名 共変量を加えたロジスティ…
今回はレアバリアント解析についてまとめた資料を公開しておきます。 このスライドの内容 内容は、Lee S, et al (2014) Rare-variant association analysis: Study designs and statistical tests. Am J Hum Genet. をもとにしています。他にも日本語の書籍*…
今日は朝型人間（Morningness）を規定する遺伝子の紹介です。（結果が膨大なので、交互作用やpathway解析の部分は省略していますので、もし気になる方は本文をご参照下さい。また、かなり学術的な解説なので、あまりお勧めはしませんが、簡単な解説をご希望…
今日はベン図を描きたいと思います。 (今回は例として、東海三県での多い名字30位についてベン図を書いてみようと思う) パッケージの読み込み library(VennDiagram) まずはベン図で書きたいそれぞれのカテゴリに含まれる要素を書き出す aichi <- c("鈴木","…
今日はリストを比較する時に役にたつコマンドを紹介します。 これらは膨大なリストを扱う際に、「共通のものを抽出する」「片方に特有のものを抽出する」などに非常に役に立ちます。 早速、実践 下記のようにAとBのリストが与えられているとする。 #AとBそれ…
こんにちは。今回は異なる環境や場所で行った同様の実験や研究によって導かれたp値を統合しようという試みがあった時の対処法について、Fisher's methodを取り上げて、今回は説明します。 FIsher's methodとは i個あるp値のlogをとった値を足し合わせて、そ…
今回はPLINKで典型的なSinge-SNP GWAS（アレルベースのカイ2乗検定）を実施する。元々、plinkが使えるディレクトリ内にbim,bed,famファイルが生成されていることが前提条件である。その他の方法についてはまた別の機会に書くことにする。 コマンド plink --b…
こんばんわ。今日は以前の投稿でも取り扱ったGWASのcase-control studyで問題となり得る「集団の構造化（以下PF）」について書きます。今回は『Japanese Population Structure, Based on SNP Genotypes from 7003 Individuals Compared to Other Ethnic Grou…
こんばんは。前回のSKATの記事の内容に少し付け加え（Kernelの種類を増やして検討しただけ）をします。基本的な解析は前回の記事をご参照下さい。 jojoshin.hatenablog.com はじめに 前回の記事にも記載した通り、SKAT()は二値のアウトカムには適しておらず…
だいぶ更新の時期が空きましたが、研究の方が順調に進んでいるという証拠だと個人的にはポジティブに捉えているところです。今回は文献紹介です。特に学術的な「面白さ」というよりは社会的な「面白さ」で題材を選びました。GWASというのはこれまでもご説明…
MultivariateもMultivariableってどっちも「多変量」? こんな疑問を持ったことのある人もいるのではないかと思います。確かに、"Multivariate"と"Multivariable"のどちらも辞書では「多変量」となります。疫学や公衆衛生の分野だけでなく、共変量で調整した…
久しぶりの更新になりました。6月末の北海道でのISEE-ISES AC2016への参加とその後は論文執筆に集中していました。ISEE-ISESのアジア支部総会では、pre-conference workshopにも参加し、distributed lag non linear model（DLNM）という時系列分析の一種を学…
今回はデータの整頓から少し分析に近いことを始めます。その中でデータを要約し、図示することが解析の一歩かと思い今回の内容にしています。R for beginners vol.1 「Rの紹介と基本的なコマンド」 jojoshin.hatenablog.comR for beginners vol.2 「データの…
今日は以前のエントリで解説したインピュテーションされたGWASデータ（調整項目を指定したロジスティック回帰）を解析するmach2datによって出力されるデータの解釈について、記述する。（確か、前回はREAD ME読んでくださいとか言って、解説をしていませんで…
今回は疫学などヒトを対象とした研究で対象者の特性として、主な変数の要約することが多くありますが、その表1（table1）を描き、csvとして出力できる便利な関数の紹介です。 その便利な関数はCreateTableOne()関数です。（パッケージはtableoneと言うもので…
Rの例としてよく使用されるirisデータを用いて、残差を簡単に求める。 特に今回は例として、Sepal.Length（がくの長さ）とSepal.Width（がくの幅）について、残差を求める head(iris) ## Sepal.Length Sepal.Width Petal.Length Petal.Width Species ## 1 5.…
