2026年8月– date –
-
P値依存症から脱却せよ！臨床試験のモロさを「人」の数で直感評価する Fragility Index (FI) 完全ガイド
「有名誌（NEJMやLancetなど）に掲載された最新のランダム化比較試験（RCT）で、$P = 0.042$ の有意差が出ている。新薬群は対照群に対して主要評価項目の発生率を統計学的に有意に抑制した！……よし、この新薬は間違いなく有効だから、明日からの日常診療に... -
「未測定交絡の定量的な感度分析を行え」と言われたら？E-value算出・95%信頼区間解釈と英文テンプレート
「各専門科の英文誌に論文を投稿したところ、査読者から『本研究は観察研究であり、喫煙歴や詳細な重症度などの未測定交絡（Unmeasured Confounding）の影響が否定できない。単にLimitationに可能性を書くだけでなく、定量的な感度分析（Sensitivity Analy... -
Lasso正則化はP値を出さない手法である！過学習を防ぐL1ペナルティの仕組みと lambda.1se による実務モデル構築
「査読者から『多変量モデルでのステップワイズ法による自動変数選択は過学習やP-hackingを招くため不可』と指摘されました。そこで機械学習で流行りのLasso（ラッソ）を試したのですが、Rの glmnet パッケージで解析したら……あれ？ P値も標準誤差（SE）も... -
ANOVAは「平均の比較」じゃない！回帰モデルの比較に anova() を使う本当の理由とあてはまり評価の実践
「先生！モデルに新しいリスク因子（$X_2$）を追加してみたのですが、この変数を追加したことでモデルのあてはまり（説明力）が本当に良くなったのか、客観的に評価する方法が分かりません」 「それなら、シンプルなモデル（model1）と変数を増やしたモデ... -
U字・J字型リスクを見逃さない！現代医学統計のデファクトスタンダード「Restricted Cubic Splines」によるCox回帰と連続HR可視化
「先生！BMIと死亡リスクの関係を見たくて多変量Cox比例ハザードモデルを組んだのですが、$P > 0.05$ で有意差が出ませんでした。でもデータをグラフで見ると、どう見てもU字型（低体重でも肥満でも死亡リスクが高い）なんです！直線モデルのまま論文に... -
順序ロジスティック回帰で比例オッズ性が崩れたらどうする？グラフ法による平行性確認と部分比例オッズモデルの実行
「先生！疾患の重症度（軽症/中等症/重症/最重症）を目的変数にして順序ロジスティック回帰モデルを組もうとし、Brant test（ブラント検定）をやったら $P < 0.05$ と判定されてしまいました！比例オッズ仮定が崩れているから順序ロジスティック回帰は... -
30項目以上の背景因子と非線形リスクをどう解析すべきか？Random Survival Forest（RSF）とrandomForestSRCパッケージによる解決法
「先生！多変量Cox比例ハザードモデルを組もうとしたら、比例ハザード性（PH仮定）が複数の変数で崩れてしまい、単一のハザード比（HR）が出せなくなりました！おまけに背景因子が30項目以上あって、ステップワイズ法を入れてもモデルが過学習で崩壊してし... -
同じ患者が何度も再入院したデータをどう解析すべきか？Cox比例ハザードモデルの拡張とAndersen-Gillモデルの実践手順
「先生！心不全の患者さんの再入院データを解析したくて、2回目や3回目の再入院データもたくさんあったのですが、解析方法がよく分からなかったので、とりあえず『初回再入院までの時間（Time-to-first event）』だけでCox比例ハザードモデルを実行しまし... -
なぜ他病死を打切りにするとイベント発生率が高く見積もられるのか？累積発生関数（CIF）とtidycmprskによる実践解決
「先生！がん患者さんの5年がん特異死亡率を調べたくて、追跡途中で他病死（老衰や心不全など）された患者さんを『中途打切り（Censoring）』として1-KM（1 - Kaplan-Meier）曲線を描きました！綺麗に死亡率が上昇していくグラフが描けました！」 ……ちょっ... -
P > 0.05だから非劣性」は査読で即却下！Cox比例ハザードモデルによる非劣性生存時間解析とICH-E9準拠の正攻法
「先生！低侵襲な新手術と標準手術の全生存期間（OS）を比較したところ、$P = 0.32$ で有意差なしでした！つまり、新手術の非劣性が証明されたということですよね？早速論文を書こうと思います！」 ……ちょっと待ってほしい。それは論文査読（Peer Review）... -
【R言語】カルテの「繰り返し測定データ」や「途中の治療変更」をそのままCox回帰に入れていませんか？時間依存共変量とtmerge関数の完全ガイド
研究者: 観察期間中に毎月測定した血圧の変化や、追跡途中で追加した抗がん剤の影響をそのままCox回帰モデルの共変量（1/0）として投入しても大丈夫でしょうか？ 統計ER: ちょっと待ってください！何の気なしに扱おうとしているそのデータは、統計学で『時... -
【R言語】Cox回帰で比例ハザード性が怪しい時の対処法！RMST（制限付き平均生存時間）の使い方と層別化・期間分割の使い分け
研究者: カプランマイヤー曲線を描いたら、途中で2本の線が交差（Cross）してしまいました……。このままCox回帰を行ってハザード比（HR）を出しても大丈夫でしょうか？ 統計ER: 曲線が交差している（＝ハザード比が時間で変動している）場合、単一のハザー... -
【R言語】層別化Cox回帰の完全ガイド！strata()でハザード比が消える理由から論文の記述作法まで
研究者: 査読者から「施設間差を層別化Coxモデルで調整しなさい」と言われ、strata() を使って解析したのですが……指定した変数のハザード比（HR）が出力から消えてしまいました！エラーでしょうか？ 統計ER: ご安心ください、それはエラーではなく層別化Co... -
3群以上のログランク検定で差が出たらどうする？pairwise_survdiff()とggsurvplot()で解く多群解析の作法
研究者: 治療法が3群（A薬・B薬・C薬）ある生存時間解析を行っています。全体のログランク検定で有意差が出たのですが、どの群同士に差があるのかの調べ方や、Cox回帰の基準群の設定方法、グラフの綺麗な描き方が分かりません……。 統計ER: 3群以上の比較で... -
傾向スコア重み付けでSMDが下がらない！RのIPTW解析でバランスを劇的に改善する実践テクニック
研究者: 傾向スコア逆確率重み付け（IPTW）を行ったのですが、一部の主要な共変量で標準化平均差（SMD）が 0.1 未満にならずバランスが取れません……。傾向スコアのモデルには変数を投入しているのですが、どう解決すればいいでしょうか？ 統計ER: 単純な二... -
査読者（Reviewer）のツッコミを未然に防ぐ！IPTW解析における重み分布の確認と箱ひげ図（Boxplot）の描き方
研究者: 傾向スコア逆確率重み付け（IPTW）解析を行っているのですが、査読者や指導医から「重みの分布を箱ひげ図（Boxplot）で確認したか？」と聞かれました。なぜ箱ひげ図が必要で、どう描いて外れ値を処置すればいいのでしょうか……？ 統計ER: IPTWはた... -
【R言語】3群以上の生存時間データをIPTW解析する方法！重み付きKM曲線描画からCox回帰（HR算出）まで完全攻略
研究者: 3群以上（例：「標準治療 vs 治療A vs 治療B」）ある治療群の生存時間データで交絡調整を行いたいのですが、傾向スコア重み付け（IPTW）のコードや重み付きカプランマイヤー曲線の描き方が分かりません……。 統計ER: 3群以上の多群（Multi-treatmen... -
【R言語】回帰分析の結果をExcelへきれいに貼り付ける方法！手入力卒業で転記ミスゼロへ
研究者: Rでロジスティック回帰を実行したのですが、コンソール画面の summary() の数値を1つずつExcelに手入力するのが大変です。直接コピペすると1つのセルに固まって崩れてしまうのですが、簡単に綺麗な表にする方法はありますか？ 統計ER: summary() ... -
Rで臨床データ前処理！文字列操作（結合・抽出・置換・桁揃え）鉄板テクニック
【研究者からの質問】 > 電子カルテやアンケートから抽出したcsvデータをRに読み込ませたら、エラーが出て解析が止まってしまいました……。「120 mmHg」などの単位付きデータはどう処理すればいいですか？ 【統計ERからの回答】 > 臨床データ解析に費や... -
【R言語】因子水準（Factor Levels）を逆順にする鉄板テクニック：回帰分析の基準群変更とggplot2の軸順反転
研究者: Rでロジスティック回帰分析を行ったら、対照群（Control）ではなく介入群が基準になってオッズ比が 0.2 と表示されてしまいました……。基準を入れ替えて見やすくする方法はありますか？ 統計ER: それはRの「因子水準（Factor Levels）の順序ルール...
1