Yoshi Nishikawa's Blog

〜医学・疫学・統計学編〜

ゲーム理論ことはじめ

まずは1-3章を読んだ。 ナッシュ均衡 相手が戦略を変えない限り、自分に戦略を変えるメリットがない戦略の組み合わせ 両性の戦い チキンゲーム 囚人のジレンマゲーム パレート最適 他者を犠牲にすることなく厚生を改善する余地が、誰にとってもなくなってい…

NAのある行を削除する

ある変数hogeの中に、NA(欠測値)があり、それを含む行だけを削除したい時。 data[complete.cases(data$hoge), ] でOK。

記述統計に立ち返る

記述統計にはいつもsummaryを使っていたが、 summary(data) を利用していたが、psych パッケージの describe(data) mean, sd, median, range, se, (歪度: skewness)https://en.wikipedia.org/wiki/Skewness, (尖度: Kurtosis)https://en.wikipedia.org/wiki/…

施設に関するCase studyを記す

Developing a methodological framework for organisational case studies: a rapid review and consensus development process | The EQUATOR Network 読んだ。 Case studyもヘルスサービスリサーチとして、研究のガイドラインが存在する。 本編はこちら チ…

Content analysisことはじめ

読んだ。意訳だが、各項目を列挙しておく。 歴史 その概念形成 その利用と推論 内容分析のロジック 利用 サンプリング 記録・コーディング データ言語 分析の構成 分析の代表的手法 コンピューターによる補助 信頼性 妥当性 実践ガイド 内容分析の歴史は古く…

ベイズ統計を学ぶ

ネット上にも、いろいろと教材が有るので、備忘録として載せておく。 Stan超初心者入門 階層モデルの分散パラメータの事前分布について 階層ベイズとWAIC Stanコードの書き方 中級編 Prior distributions for variance parameters in hierarchical models. A…

StanとRでベイズ統計モデリング

StanとRでベイズ統計モデリング (Wonderful R) 読んでいる。これは名著である。前半部分は、ベイズ統計関わらず、役立つ内容。とくに6章、モデリングの視点から確率分布に迫る箇所は、なかなかまとまったものがないので一見の価値あり。 1章 統計モデリング…

insetのある地図を作りたい

ArcGISでinsetのある地図を作る方法に関しては、ここで解説されている。 insetを入れるたびに探しているので、ここにリンクを貼り付けておく。 desktop.arcgis.com プレビュー画面と行き来しながら見ると良い。

csvファイルの読み込み

元データから読み込む際に、高頻度で出現する(していた)R上のエラー。対応策を備忘録として残しておく。 読み込みたいcsvを「UTF-8」形式にしておく。 私はmiというエディタを用いている。 “as.is=T"を入れる 文字列、数値が混在していてもOK。もし、数値…

5-FUによる高アンモニア血症では、異化産物のFBALとFMAが蓄積

5-FUによる高アンモニア血症の5-FU代謝物血中濃度分析論文についての解説記事 わたしたちのチームの研究で、5-FUによる高アンモニア血症を来した時にフルオロ・ベータ・アラニン(FBAL)とモノフルオロ酢酸(FMA)という異化産物の血中濃度が上昇していることが…

フォレストプロットを描く~JMP編

JMP JMPという統計ソフトを用いて学んだ。 www.jmp.com SAS社が開発したソフトウェア。 Forest plot 以下のページでforest plotというアドインを以下で取得。Rでいうpackageですね。 community.jmp.com 同ページ内においてあったRosaMeta.jmpで作成 同ページ…

Epidemiology and the Delivery of Health Care Services: Methods and Applications

ざっと読んだ。 PART1 1. An Epidemiologic Framework for the Delivery of Health Care Services 2. Measurement Issues in the Use of Epidemiologic Data 3. Descriptive Epidemiological Methods 4. Epidemiological Study Designs for Evaluating Healt…

ggplot2でグラフ作成 (3)

ggplot2 2016年に出版されたggplot2に関するHadley Wickhamによる著書。その名もggplot2、の続き。この本、最後まで読み解いて分かるが、本当におすすめ。データを扱うすべての人が読んだら良いと思う。 コンテンツ 大きく分けて3部で構成され、全12章にわか…

MCMC

はじめての統計データ分析 「はじめての統計データ分析」 に関して覚書。 1章の記事はこちら。 2章のハイライト 1) MCMC(Markov chain Monte Carlo method)により事後分布のデータを用い、 2) MCMCには、Metropolis-Hasting method, Gibbs sampling methods,…

survival aNalysis ~no no R more~

(またおっさんホイホイなタイトルをつけてしまった) ggplot2でsurvival analysisのプロットをしたいときは、 GGally パッケージに含まれる、ggsurv関数がある。 https://cran.r-project.org/web/packages/GGally/GGally.pdf 自分でggsurv関数を入れ込みた…

Pubmedで文献検索

文献検索には型がある ネット上で色々と教材が有るので、貼り付けておきます。 音羽病院廣江先生 20170424_PubMedの使い方@洛和会音羽病院 東京大学図書館 http://www.lib.m.u-tokyo.ac.jp/manual/pubmedmanual.pdf 慶應大学図書館 http://www.med.lib.keio.…

グラフ理論

読んだ。 1章 グラフの基礎 2章 最小全域木 3章 最短経路問題 4章 オイラー回路とハミルトン閉路 5章 グラフの彩色 6章 最大流問題 7章 マッチング グラフ理論の入門編としておすすめできる。 私が入門編の書物に求める要素 1. 読みやすい 具体例とともに簡…

ggraphでグラフ作成

ggraph/README.md at master · thomasp85/ggraph · GitHub ggplot2 extensions: ggraph ggplotのextensionはこんなところにも。

大腸がんスクリーニング (Colorectal Cancer Screening)

Colorectal Cancer Screening に関する書籍をざっと読んだ。大腸がんのスクリーニングが焦点。分子生物学的な観点も含めて記述してあって知識の確認になった。 Colorectal Cancer Pathways Petr Protiva Risk Factors and Screening for Colorectal Cancer J…

Rで21世紀の相関係数を算出する(MIC)

相関 2変数とも連続データで,正規分布に従っているならばPearsonの相関係数を用いることができる。 少なくとも1変数が非連続データの時にはノンパラメトリック検定のSpearmanやKendallの相関係数を用いる。これらは、実データでなく、順位付け(大小関係)…

データ解析チートシート cheat sheet

チートシート R Studioのウェブサイトにまとまっており、 www.rstudio.com 以下のように様々なチートシートがある。 Data Wrangling with dplyr and tidyr Cheat Sheet Data Visualization with ggplot2 Cheat Sheet これは便利だ。

Network Meta-Analysisについて

Network Meta-Analysis (NMA) 従来のメタアナリシスの課題 従来のメタアナリシスでは、2つの治療についての比較しか得られない。例えばうつ病に対する薬剤には多くの種類があり、有効性・安全性・薬価は様々である。 そこで、利用可能なすべての薬剤につい…

ポリファーマシー(Polypharmacy) [介護施設編]

ポリファーマシー(Polypharmacy)とは 単純に「4~6種類以上の薬を飲んでいる状態」をさすこともあるが、問題になるのは、不適切な多剤併用である。 これを受けて、Polypharmacyとは、往々にして「必要以上に薬を飲んでいて、薬による有害事象が起こっている…

可視セルのみをコピーする、隠れたA列を表示する

可視セルのみをコピーする ・Windows → 「Alt+;(セミコロン)」 ・Mac → 「shift + command + Z」 のショートカットキーで選択して貼り付ければOKです。 A列が表示されない場合 隠れたA列を表示するには、B列から左にドラッグして 2C と表示されてからマ…

qplotで日本語を使う

ggplotを手軽に ggplotのコードが覚えられない、忘れる。 何故かエラーになる、そんなことはありませんか? ggplotを手軽に行うための関数があります。 その名も、qplot(quickplot) qplot qplot. ggplot2 2.1.0 ここを参考に色々な事ができる。 plot関数に近…

GLMMを実装する

The glmmADMB package このパッケージを用いて、以下でinstallする。 install.packages("R2admb") install.packages("glmmADMB", repos=c("http://glmmadmb.r-forge.r-project.org/repos", getOption("repos")), type="source") これでzero inflationなどの…

google mapを用いてある地点のシェープファイルを生成する

シェープファイル(Shapefile)とは シェープファイル(Shapefile)とは、Esri 社の提唱したベクトル形式の業界標準フォーマットです。Esri 製品はもちろん、多くの GIS ソフトウェアで利用が可能です。シェープファイルは仕様が公開されています。シェープフ…

Mac bookのbootcampで起動したwindowsでスクリーンショットを撮る

ブートキャンプのウィンドウズでスクリーンショット Windowsをブートキャンプで起動することがあります。 私の場合、ArcGISを使う場合がほとんどです。 以下のサイトを参考にしました。 www.webantena.net 主に3つの方法があります。 osk (スクリーンキーボ…

UGT1A1*28*28と肺癌リスクについて

UGT1A1*28*28と肺癌リスクについての解説記事 最新の研究で、UGT1A1*28*28と肺癌の罹患との関係が明らかになりました。この研究は2016年11月10日のInternational Journal of Clinical Oncology誌(オンライン版)に掲載されました。本稿は、今年読んだ一番好…

circular statistics (角度統計)について

髄膜炎に関する論文 Seasonal dynamics of bacterial meningitis: a time-series analysis. The Lancet Global Health http://dx.doi.org/10.1016/S2214-109X(16)30064-X 以前読んだ論文。前回はWaveletについて述べた。今日は、角度統計(circular statistic…