統計シェアするTwitterGoogle+Pocketさて、今回からイデオロギー色の少ないブログ活動に戻ります。Analysis")」というテクニックを紹介しました。私の言葉足らずか、このテクニックに関して多くの方に誤解と混乱を招いたようなので今回改めて整理してみたいと思います。参考文献なぜ感度分析が有効なのか?その意義は?観察データを用いた分析から因果関係を言うためには、様々なバイアスに対処しなければいけません。数あるバイアスの中でも最もよくあるものが、交絡によるバイアスです。DAG上ではアウトカムと曝露因子の共通原因によって引き起こされるバイアスと整理されます。詳しくは以下の過去記事を参照ください。この交絡バイアスですが様々な対処方法があります。概要は上記の記事にまとめてあるので興味がある方は読んでみてください。一般的には、見たい因果関係を交絡する交絡因子に関するデータをすべて把握し、研究デザインや統計解析を通してその影響を制御する、というのが基本方針になります。あらゆる交絡因子の影響を制御した場合、観察データからであってもドンピシャ正確な因果関係が得られるのです。これは統計的因果推論の数学的背景を確認すればすぐわかります。問題は、「すべての交絡因子を把握、対処」という前提が成り立つことが実質不可能という点です。データがとれないことも多いでしょう。Confounders)の影響が、観察研究から得られた結果を解釈するうえで非常に重要なポイントになります。ランダム化をしていない以上、多かれ少なかれ未測定交絡因子は必ず存在します。Effect)とは異なる値になっています。因果関係を”正確に”推定できていないことになりますので、よくこの辺りを批判されます。統計学の世界で著名なフィッシャーが喫煙の健康被害について懐疑的であったという話は有名ですが、フィッシャーの批判も上記の点を問題視したものでした。彼は「遺伝要因」による交絡の可能性を指摘していたそうです。以下のようなイメージ。タバコを吸いやすくなる遺伝子って何だよ、とつっこっみたくなるのですが、まあ仮にそういうものが存在するした場合、たしかに交絡によって計算された因果効果がゆがめられていそうです。ここで考えるべきは「未測定交絡因子によってどのくらいの因果効果の推定がゆがめられているか(どの程度のバイアスがあるのか)」がありうるか、ということです。バイアスの程度に関する具体例はこちらを再読してみてください。未測定交絡因子によるバイアスの可能性を指摘するのであれば、どのような交絡因子が存在しうるのか、それによって観察されている関連は完全に説明できてしまうのか、逆に観察されている関連が完全に交絡によるものだというためにはどのような交絡因子が存在する必要があるのか、といった点をディスカッションしなければ建設的な批判とはいえません。交絡によるバイアスを考える際には、1バイアスの方向(Direction)と2その大きさ(Magnitude)の両方を検討する必要があります。この二つは1交絡因子とアウトカムの関連、2交絡因子と曝露因子の関連、3注目している集団における交絡因子の存在割合(Prevalence)という三つの要素によって決まります。これは超重要なポイントです。以下に詳細を説明します。交絡バイアスの向きバイアスにも向きがあります。喫煙が肺がんに与える”真の”因果効果がリスク比1.5だったとしましょう。未測定交絡因子によるバイアスがある場合、実際に観察される喫煙と肺がんの関連は1.5とは異なる値をとることになります。この時、観察されたリスク比が2.0だったとすると、真の因果効果よりも過大に因果効果を見積もっていることになります。bias)と呼ぶことにしましょう。bias)が生じていることになります。value;リスク比=1またはリスク差=0)から離れていくようなバイアスのほうが問題だと考えます。value)かもしれず注意が必要です。バイアスが下向きならば「真の因果効果は観察されているものよりももっと強い」といえるので「喫煙が肺がんのリスクをあげる」という結論自体は変わらなさそうです。やってバイアスの向きを予想するのか。交絡因子とアウトカム、交絡因子と曝露因子それぞれの関連の向きに注目すれば、交絡バイアスの向きが分かります。以下の図において、+は正の関連(片方が増えると、もう一方も増える)、-は負の関連(片方が増えると、もう一方は減る)を表すとしてください。要するに、交絡因子と曝露因子の関連、交絡因子とアウトカムの関連の向きが同じならば上向きバイアス、向きが逆なら下向きバイアスということになります。フィッシャーの例をとるならば、ある遺伝子を持っている人で肺がんのリスクが高い(+)かつタバコを吸いやすい(+)、もしくは肺がんのリスクが低い(-)かつタバコを吸わない(-)ならば上向きバイアスが生じることになります。仮に、交絡バイアスを引き起こす遺伝要因があったとしても、遺伝子と喫煙・肺がんそれぞれとの関連の向きが一致していない場合は、因果効果が”ない”可能性の指摘としては的外れ(むしろ因果関係の存在をサポートする)といえます。valueに近づくものであったとして、でははたして考えうる未測定交絡因子は観察されている関連を交絡バイアスだけで引き起こしてしまうほど強力なものなのか、どのような未測定交絡因子が存在すればそのような状況になるのかを考える必要があります。Analysis)があります。一般に感度分析という言葉は、統計モデリングのやり方を変えたり、異なる定義をつかって変数を定義したりなどすることで結果の頑健性を確認するテクニックのことを指しますが、今回は交絡バイアスの大きさの程度に目星をつける手法のことだと思ってください。Assumptions”とはこれより先は、前回の記事の中で引用していた以下の論文で紹介されている感度分析を中心に書きます。368.最初に結論だけ言っておくと、交絡バイアスの大きさは1交絡因子とアウトカムの関連の強さ2交絡因子と曝露因子の関連の強さ3集団における交絡因子の存在割合によって決まります。assumptions"とあります。”前提を置いていない”ということです。以下紹介する感度分析のテクニックは1と2の条件のみを使い、3の条件を考慮していないため、このように書かれています。*1実は注目している集団において交絡因子が珍しいほど(割合0%に近づくほど)、もしくは逆に非常にありふれたものであるほど(100%に近づくほど)、交絡バイアスは弱くなることが知られています。この交絡因子の存在割合に関する条件を指定していないため、以下の手法ではバイアスの大きさが過剰に見積もられています。言い換えると、12の条件下で起こりうる最大のバイアスを推定しています。交絡因子のprevalenceによっては、推定されたよりも実際のバイアスは小さい可能性もありますので、これからご紹介するテクニックは非常に保守的(Conservative)な方法であると言われています。アウトカムすべて二値の変数であると仮定してください。今回は省略しますが、連続変数の交絡因子に対する感度分析も基本的には同じアイデアです。以下、簡単な数学表現も使うので一度用語を定義しておきます。conditionsと呼ばれる古典的な条件は次のようなものです。”UA>RRかつUY>RRでないならば、Uだけで観察されている関連(RR)を説明することができない”言い換えると、RRがUAまたはUYの少なくともどちらか一方よりも大きい場合は、交絡バイアスだけで観察されている関連を説明できないということです。フィッシャー大先生の具体例を使いましょう。1958)に対して、遺伝子Uによる交絡の存在を問題視しました。10.73を満たす遺伝子Uでなければ喫煙と肺がんの因果関係は覆されないと言えます。Cornfield条件にはもう一つ条件があって、UA=10.73で観測されている関連を説明するためには、UとYの関連がPerfect(UY=∞)でなければいけないというものです。UY=10.73のときにはUA=∞です。このような遺伝子が存在するかどうかは、その他の研究結果を踏まえて考えるのですが、基本的には分野の常識から言ってありえない条件だろうと言えます。したがってフィッシャー大先生の指摘する遺伝要因による交絡バイアスだけでは、喫煙と肺がんの因果関係を否定することは難しそうです。factorをもとめるCornfield条件ではUA、UYの大きさから交絡バイアスで観察されている関連が説明されてしまうかが分かりました。実はUA、UYが分かっていればでバイアスの程度を直接求めることもできます。factorと呼んでいます。このBはバイアスの程度を表す数値であり、観測された関連をBで割ることで、交絡バイアスがない場合の真の因果効果を計算することができます。RR=10.73の場合を考えます。UY=10.73、UA=10.73だったとしましょう。value)なので、これほど強力な交絡因子であってもそれだけで観察された関連を説明することができないと分かります。-1)の式でUY→∞でB=UAUA→∞でB=UYとなりますのでCornfield条件と一致します。factorを解釈するうえで、注意すべきポイントが2点あります。maximum)さきほども書きましたが、交絡因子のPrevalenceに関する条件を指定していないので、B=5.63はUY=UA=10.73を満たすUが起こしうる最大のバイアスです。逆に言うと、UY=UA=10.73でも、実際のバイアスはもっと小さいことになります。これらの感度分析がConservativeと言われるゆえんです。confoundersこれも非常に重要なポイントです。元論文を見れば詳しく数学的背景も書いてあるのですが、実は上記の式に含まれているパラメータは”すでに調整されている交絡因子を考慮したうえでの”値を採択します。多くの研究では単純に曝露因子とアウトカムの関連を見たりせず、回帰分析などを使って性、年齢、学歴など主な交絡因子の影響を調整したうえで注目している二つの因子の関連を見ていると思います。Bの計算に使うUA、UYはこれら交絡因子の影響を除いたうえでの値になります。交絡因子が多く観測され、調整されているほどUA、UYが大きい値をとるのが難しい→Bが大きい値を取りにくくなる→未測定交絡因子によるバイアスだけで説明がつきにくくなるということです。どんな交絡因子ならば観測された関連を説明できるか実際にはUA、UYの値はわかりません(Uが未測定なので当然です)。ほかの研究から得られた知見からUAとUYのおおよその値を予想したり、さまざまな値を入れてみたときのBの一覧を表にするなどの方法がとられます。発想を逆転させて、「観察されている関連を完全に交絡バイアスだけで説明するためには、未測定の交絡因子はどれくらい強力である必要があるか」を求めることができます。それには次のような条件を使います。(RR-1)]*3よりも大きい場合は交絡バイアスだけで観測されている関連を説明できるが、これより”弱い”交絡では説明できない。20.95なのでUY=UA=20.95かそれ以上の交絡因子UならばRR=10.73をすべて交絡バイアスで説明できるが、これより”弱い”交絡では説明できないということです。これはUYとUAが必ず20.95以上でないといけないという意味ではありません。UYが20.95を下回る場合、UAが20.95よりも大きくなければUY=UA=20.95のUが起こす交絡バイアスよりも”弱い”交絡となってしまいます。UYが小さくなればUAがさらに大きい必要があるということです。factorと同じです。causality)と呼んでいます。factorとE-valueの関係をグラフにしました。使用したRコードは以下の通り。Rコードのうまい貼り付け方法ってあるんですかね?なんか味気ない、、、詳しい人いたら教えてください。cex=0.75)アウトプットが以下のグラフです。10.73(すなわち、観測されたリスク比をすべて説明できてしまう)の未測定交絡因子となるようなUY=UAの組み合わせをプロットしたものです。UY=UA=E-value=20.9ならば、交絡バイアスだけで説明がつくことがわかります。注目すべき点はUYが20.9より小さくなるほど(グラフの右のほうにくほど)、UAが大きくなっていることです。これが「UY=UA=20.9の交絡よりも”弱い”交絡では観測された関連を説明できない」の意味するところです。UY=∞でUA=10.73に収束していきますし、UA=∞でUY=10.73となりますのでCornfield条件が成立していることも確認できます。以上が未測定の交絡因子によるバイアスの程度を考える際に便利な感度分析の概要です。疫学の世界でも比較的新しい概念で、実際の論文で使われることはまだまだ稀です。これらのテクニックを使って個々の観察研究のバイアスの程度を考えることでより建設的な議論・批判的吟味が行われ、分野全体で研究の質が上がっていくのではと期待されています。analysisについて書かれています。Analysis)に関する本なので、その文脈での感度分析の話ですが、基本的な考え方はこの本でわかるはずです。関連記事2018-12-28回帰分析を使った因果推論の仮定:パラメトリックモデルを使うということお久しぶりです。Analysis)入門1...お久しぶりです。(3)スターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。