2026年10月– date –
-
メタアナリシス
metaとmetaforでメタ回帰の結果が違う?原因は設定にあった【R】
研究者: 同じBCGワクチンのデータで、metafor と meta の両方でメタ回帰をしてみたんです。そうしたら、異質性 $\tau^2$ が 0.0504 と 0.3378 で全然違う結果になって……。meta のメタ回帰は $\tau^2$ を推定し直していない、ということでしょうか? 統計ER... -
ロバスト分散
ロバスト分散を用いたGLMM:サンドイッチ推定の考え方とRでの実装(lme4・merDeriv)
研究者: 反復測定の2値データでGLMM(glmer)を当てはめました。査読対策に「ロバスト分散」も出しておこうと思って sandwich::vcovCL() を使ったら、エラーで止まってしまいました……。そもそも、ロバスト分散を使えばGLMMの仮定が崩れていても大丈夫にな... -
多重代入法
多重代入×IPTWのTable 1は何を平均する?Nは?mice・MatchThem・gtsummaryで作る背景比較表【R】
研究者: 欠損値のある臨床データで多重代入(mice)とIPTW(傾向スコア逆確率重み付け)を両方行う必要があるのですが、代入データをどう統合して患者背景比較表(Table 1)を作成し、Excelに貼ればいいのかわかりません……。 統計ER: mice + MatchThem(We... -
検定の基礎知識
「関連は認められなかった」と書く前に—p値関数で防ぐ”有意差なし”の誤読
若手医師:先生、例のコホート研究なんですが、重症度スコアと年齢で調整したら、ハザード比が1.38、95%信頼区間が0.94〜2.02になりました。p = 0.096です。Discussionには「曝露とイベントとの間に関連は認められなかった」と書けばいいですよね? 統計家... -
検定の基礎知識
非劣性は「差がゼロ」ではない!「少し劣っていても合格」の本当の意味
本記事は前回の記事(有意差なしと同等性は異なる、をテーマにした記事)の続編である。 1. はじめに:「有意差なし=同じ」と言えないなら、何を示せばよいのか ある日の研究室での会話である。 後輩「先生、新しい薬を標準薬と比べたら、有意差なしでし... -
検定の基礎知識
「有意差なし」は「差なし」ではない――約30年続く誤用の正体
新薬Aと標準治療Bを比べた小さな臨床研究の解析が終わった。結果は次のとおりである。 先輩「P値は0.23。有意差なしだから、AとBは同等と書いておこう」後輩「え、でも……『差がない』と『差があるとは言えなかった』は、同じ意味なんですか?」 後輩の違和... -
回帰分析・線形回帰・重回帰
Shapiro-Wilk検定で有意になった…重回帰は諦めるべき?正規性の誤解を解く
「先生、目的変数をShapiro-Wilk検定にかけたら、P値が0.001未満でした。説明変数もほとんど正規分布ではありません。重回帰分析は諦めるしかないでしょうか…」 統計の相談を受けていると、こうした声を本当によく聞く。ここで、結論から述べておく。 絶望... -
検定の基礎知識
検定結果を見てからストーリーを決めてはいけない理由
後輩「先生、解析が終わりました。主要評価項目は有意になりませんでした……。でも、年齢×治療の交互作用だけp < 0.05でした!」先輩「それはいい。じゃあ『高齢者でだけ効く治療』というストーリーで論文を書こう」 一見、よくある研究室の風景だ。しか...
1
