ITSAとは ITSデザインに適した状況は？ インパクトモデルを選ぶ 分析実行例 使用するデータ 散布図を描く Level Change Model Slope Change Model Level + Slope Change Model 自己相関を確認する おわりに 参考文献 ITSAとは 例えば新しい政策が実施される…

前回、「Joint Frailty-Copula Modelを使った生存時間解析」の理論的なことについてまとめました。 necostat.hatenablog.jp 読んだ本： Survival Analysis with Correlated Endpoints: Joint Frailty-Copula Models (JSS Research Series in Statistics)作者…

前回、パラメトリック生存時間解析を勉強してみて、ワイブル分布について学び直す必要性を痛感したのでまとめました。 necostat.hatenablog.jp 指数分布（exponential distribution） ワイブル分布（Weibull distribution） 形状パラメータ 尺度パラメータ …

前回、競合リスクデータの解析方針の1つとして、生存時間を対象としてモデル化することをあげたので、読んだ本の内容を整理してみました。 前回の記事： necostat.hatenablog.jp 読んだ本： Survival Analysis with Correlated Endpoints: Joint Frailty-Cop…

率と割合の1対1対応関係 競合リスクがあると困ること 率と割合の1対1対応関係が崩れる CIFと1-KMの違い 競合リスクデータに対する解析方針 1. ハザードを対象としたモデル 2. リスクを対象としたモデル 3. イベント発生までの時間を対象としたモデル おわり…

カッパ係数（Cohen's kappa statistic）とは 定義 解釈 重み付けカッパ係数（weighted kappa statistic） 定義 Rで計算してみる サンプルデータを作成する 方法1：irrパッケージを用いる 方法2：vcdパッケージを用いる 方法3：psychパッケージを用いる おわ…

はじめに GISデータをプロットしてみる 例1：plot( )を使う 例2：geom_sf( )を使う 例3：離散変数に従って塗り分ける 例4：連続変数に従って塗り分ける 例5：複数の図を重ねる おわりに 参考資料 はじめに Rのsfパッケージを使って、神戸市の校区境界と学校…

以前、Dropbox経由でZoteroを同期させる方法をまとめました。 necostat.hatenablog.jp この方法だとPDFしか同期できず、Webページやパワポは別の方法で管理しなくてはなりませんでした。 WebDAV（ウェブダブ）経由で同期させる方法だと解決したので、備忘録…

分散分析で満たすべき仮定として球面性の仮定（the assumption of sphericity）*1をよく目にしますが、きちんと勉強したことがなかったので（わかるところだけ）読んでみました。 分散分析に必要な仮定 球面性の仮定とは 2水準間の差の分散に基づいた定義方…

臨床現場では、対象者の年齢や症状、検査値などをもとにして計算される「何とかスコア」が使われることがよくあると思います。 臨床研究でも、既に収集された項目からcomposite scoreを計算し、変数に追加したいことはよくあります。 今回はtidyverseパッケ…

以前、データフレームの列（=変数）に対して、追加・名前変更・要約などの操作方法についてまとめました。 necostat.hatenablog.jp necostat.hatenablog.jp 最近は列の操作でacross関数を使うことが多くなってきたので、ここらへんで学び直しておこうと思い…

ロジスティック回帰モデルなどから予測されるアウトカム発生確率（=予測確率）をグラフで可視化したいとします。 1つの変数に対して予測確率がどのように変化していくかを示したければ、折れ線グラフにしたり、変数をカテゴリー化して棒グラフで表したりする…

分散分析（Analysis of Variance, ANOVA）を学ぶ目的でKutner先生のApplied Linear Statistical Models（5th edition）を拾い読みし始めました。 前回は二元配置分散分析で分散分析表を作るところまで確認しました。。 necostat.hatenablog.jp 今回は要因の…

分散分析（Analysis of Variance, ANOVA）を学ぶ目的でKutner先生のApplied Linear Statistical Models（5th edition）を拾い読みし始めました。今回は主に"Chapter 16： Single Factor Studies" から、分散分析の基本となる平方和・自由度の分割についてで…

以前、segmentedパッケージを使って変化点を探す方法を調べました。 necostat.hatenablog.jp segment( )はglmオブジェクト（=glm関数で当てはめたモデル）とlmeオブジェクト（=変量効果を含む線形回帰モデル）には対応しているんですが、lme( )は変量効果を…

アウトカムが「割合」、というか「試行数と成功数」として与えられているときに二項回帰モデルを当てはめる方法は以前まとめました。 今回は「オフセット項」を使った回帰モデルに当てはめる方法をまとめてみようと思います。 オフセット項（offset）とは 使…

構造方程式モデリング（Structural Equation Modeling, SEM）の初学者が、 タイトルの "Handbook" に誘われて買ってしまったHoyle先生の分厚い本を拾い読みしたメモです。 といっても、私にとって重すぎる内容は拾い上げられていません。 今回は "Chapter 22…

構造方程式モデリング（Structural Equation Modeling, SEM）の初学者が、 タイトルの "Handbook" に誘われて買ってしまったHoyle先生の分厚い本を拾い読みしたメモです。 といっても、私にとって重すぎる内容は拾い上げられていません。 今回は "Chapter 3…

構造方程式モデリング（Structural Equation Modeling, SEM）の初学者が、 タイトルの "Handbook" に誘われて買ってしまったHoyle先生の分厚い本を拾い読みしたメモです。 といっても、私にとって重すぎる内容は拾い上げられていません。 まずは "Chapter 1…

測定されたデータをもとにして母集団を比較するときには2種類の誤りが生じる可能性があります。 第1種過誤（αエラー）：本当は差がないのに「差がある」と判断してしまう誤り 第2種過誤（βエラー）：本当は差があるのに「差がない」と判断してしまう誤り エ…

この記事では回帰モデルが満たすべき前提条件のうち、「線形性」について評価する方法について学んでみるつもりです。回帰モデルの満たすべき仮定の概要については以下の記事を参考にしてください。 necostat.hatenablog.jp 線形性（linearity）の仮定とは C…

Rothman先生のModern Epidemiology（4th edition）の拾い読みシリーズとして"Chapter 21：Regression Analysis Part II: Model Fitting and Assessment" をまとめようと思ったんですが、もう少し統計学寄りの内容を勉強しようと思い、Applied linear statist…

Rothman先生のModern Epidemiology（4th edition）をパラパラめくって拾い読みしたメモです。 今回は "Chapter 6：Epidemiologic Study Design With Validity and Efficiency Considerations" から、研究デザインの型について。 Chapter 6では研究デザインの…

Rothman先生のModern Epidemiology（4th edition）をパラパラめくって拾い読みしたメモです。 今回は "Chapter 13：Measurement and Measurement Error" から、測定誤差について。 測定誤差（measurement error）とは 測定における誤差の分類 測定誤差の影響…

Rothman先生のModern Epidemiology, 4th edition（以下ME）をパラパラめくって拾い読みしたメモです。 今回は "Chapter 5：Measures of Effect and Measures of Association" から、効果・関連の指標について。 効果（effect） 因果効果（causal effect）と…

Rothman先生のModern Epidemiology, 4th edition（以下ME）をパラパラめくって拾い読みしたメモです。 今回は "Chapter 4：Measures of Occurrence" から、population（集団）について。 疫学とは 集団（population） 閉じた集団（closed population）と開い…

Fisher正確確率検定（Fisher's exact test） Fisher正確確率検定の問題点 固定されている周辺度数いろいろ 行も列も固定されているケース 行（あるいは列）のみ固定されているケース 合計数のみ固定されているケース Fisher正確確率検定の過剰な保守性に対処…

ランドマーク（LM）解析とは スーパーランドマーク解析とは 必要なパッケージとデータの読み込み 使用するスーパーランドマーク解析の概要 cutLM( )を使ってLM解析用データセットを準備する。 おわりに 参考資料 ランドマーク（LM）解析とは Immortal time b…

Tex( )でLaTeX表記を使う LaTeX表記の方法 おわりに 参考資料 文字や図版などの要素を紙面に配置するという作業を「組版」というが、TeXはこの作業を行うために開発された言語。 LaTeX（ラテフ, ラテク）はTeXを使いやすくした無料の文書処理システムで、数…

サブグループ解析とは、全体からある特性を持った一部の集団（例：男性）を取り出して解析をすること。ランダム化が崩れたり、解析対象集団のサイズが小さくなるから検出力が落ちたり、多重検定の問題が生じたり、と色々注意はある。 でもこの記事では、そこ…