過剰相対リスク (Excess Relative Risk) について

放影研のLSSコホート論文を読もうと思い、いろいろ下調べ。 最初は、「過剰相対リスク (ERR: Excess Relative Risk)」 について。ググってみると、長崎大の「放射線生命科学国際コンソーシアム」の用語解説のページが見つかる。
一般的な疫学研究において2つ以上の集団でリスクを比べる時、非曝露群と比べて曝露群が「何倍」のリスクがあるのかをみるのが相対リスク (RR: Relative Risk)。例えば、被爆線量がゼロの集団における疾患Aの発生リスクが10万人あたり2人で、被爆線量が1 Svの集団では10万人あたり5人だった場合、単純なRRは、5 ÷ 2 = 2.5 倍 となる。しかし、放射線被ばくのリスク評価では、放射線を浴びることによって単位線量当たりどのくらい過剰にリスクが上昇したのかをみる過剰相対リスク (ERR: Excess Relative Risk / 1Sv) という指標が用いられることが多い。上記の例では、単純なERRは、 (5 － 2) ÷ 2 ＝ 1.5 となり、被爆線量1Sv 浴びると1.5倍過剰に発生リスクが上昇することを意味する。
ほう。けど、なんでRRじゃダメなのか。もうちょっとググると、
Suissa S. Relative excess risk: an alternative measure of comparative risk. American journal of epidemiology. 1999;150(3):279-82. (PubMed, PDF)
この論文が出てきて、雑誌名から察するに疫学分野の有名な雑誌で、タイトルから察するにERRという指標を最初に提唱した論文かな。1999年ってわりと最近。で、中を見てみると、 
2種類の薬と高齢者ドライバーの自動車事故のリスクを考える。もともと、1000人の高齢者が1年間に起こす事故が10件だとする。もし、薬Aを飲むと8件余計に事故が起こり、薬Bを飲むと4件余計に事故が起こるとする。1年1000人あたり起こる事故の件数は、薬Aを飲む群で10 + 8 = 18件、 薬Bを飲む群では10 + 4 = 14件。なので、比較相対リスク (Comparative Relative Risk) は、 18 ÷ 14 = 1.3。一方で、 過剰相対リスク (Relative Excess Risk) は、(18 - 10) ÷ (14 -10)  = 2。
RRだとバックグラウンドのリスクに隠れてしまうので、暴露要因部分のリスクを比較する場合にはERRを使うのが適切だろう、という論文。そりゃそうですね。

たとえば、もともと10万人中1000人がかかる病気があって、条件Aではこれが1200件、条件Bでは1020件になるとしたら、比較相対リスクは1200 ÷ 1020 = 1.18 で、過剰比較リスクは(1200 - 1000) ÷ (1020 -1000)  = 10ですもんね。

で、ちょっと用語がすでに混乱してるけど、説明の中身からすると、比較相対リスク (Comparative Relative Risk) ＝ 相対リスク (Relative Risk) ですよね。あと、"Relative Excess Risk" と "Excess Relative Risk" は一緒で、過剰比較リスク（Comparative Excess Risk) はERRとそれぞれ同じもの？違うもの？日本語でも余剰と過剰が混在してるけど、一緒ですよね。違うのかな。整理すると、

疾病あり疾病なし
暴露XA人B人
暴露YC人D人
非暴露E人F人


として、
  • Xの危険度　R(X) = A / (A + B)
  • Yの危険度　R(Y) = C / (C + D)
  • バックグラウンド危険度　R(Bg) = E / (E + F)
とすると、
  • 相対危険度（Relative Risk） = R(X) / R(Y)
  • 寄与危険度(Attributable Risk）= R(X) - R(Y)
  • 過剰相対リスク (Excess Relative Risk) = {R(X) - R(Bg)} / {R(Y) - R(Bg)}
ですよね。

で、放射線影響協会というところの用語解説に、
死亡率（あるいは死亡数）や発生率（あるいは発生数）の観察値をO、期待値をEとすると、相対リスク（RR）、過剰相対リスク（ERR）はそれぞれ以下の式で示される。
RR ＝ O／E，　　ERR ＝ RR－1 ＝（O－E）／E
ってあるけど、ERR ＝ RR－1　になるかな…。式から判断して上の表で言い換えると、O = R(X), E = R(Y) ですよね。（O－E）ってつまりは寄与危険度ですね。それを R(X) で割ってERR...。ここでいうERRって、バックグラウンド危険度とか関係なくて（というか暴露Yがバックグラウンドで）、Suissa の「Relative excess risk」とは関係なくて、名前をつけるとすると、余剰寄与危険率比、ってことですね、うん。きっとそうだ。

単純に、危険度の上乗せ部分がもとの危険度の何倍か、の指標ですね。ERRが「2」というのは、暴露のせいで、もとのリスクの2倍が上乗せされて、「3」になってる、ってことですね。うん。そうだ。もともと10％の人がかかる病気があって、何かの暴露のERRが0.5の場合、0.1 + (0.1 * 0.5) = 0.15 、つまり15％になる、ってことですね。

ま、ERRは分かった（と思う）。なんだけど、広島や長崎のデータから、これ実際どうやって算出するのかな。またおいおい。

つづく。



Popular posts from this blog

エクセル1行をテキスト1ファイルに

ほんとによく働く。感心する。 調査票の自由記述、電子化されたデータを見るのと、手書きのものを読むの、疲れ方が全然違いますね...。もちろん後者の方が疲れるのですが、それは読みにくいからという理由だけではなくて、伝わってくるものが多いんですね。扱われている内容も内容ですし、忙しい中でわざわざ自由記述欄に書きこもうと思われるような言葉なので、まあ、感じること考えることたくさんあります。 さて、今の作業はだいぶテキストの量が多いので、 WeftQDA を使ってみようと企んでいます。初めての実戦での使用です。 データは.csvで手元にあるのですがこのままだと読み込めないので.txtにします。ずばっとコピペしてもいいのですが、「●●への言及が●●件あった。」なんて感じで数えたいので、1ケース1ファイルにしておくほうが便利そうです。500件ほどあるテキストをいちいちコピペして「名前をつけて保存」なんてやってられないのでググります。 エクセルに入力済みの連続したデータを1行単位で以下のような内容でテキストファイルを作成し、かつテキストファイル名を、セルに入力されている文字を使って保存するマクロを作りたいのですが、どのようにしたら良いのでしょうか？よろしくお願いします。（ 1行単位のデータをテキストファイルに保存するマクロ ） 良い質問です。そして良い回答です。自分用にメモ。 Sub ColumnOut2Text() Dim i As Long Dim j As Long Dim Fno As Integer Dim OutColumn As String 　 Const myPath As String = "C:\txt\" With Worksheets("Sheet1") For i = 2 To .Cells(Rows.Count, 1).End(xlUp).Row 　Fno = FreeFile() 　Open myPath & .Cells(i, 1).Value & ".txt" For Output As #Fno 　For j = 1 To .Cells(1, Columns.Count).End(xlToLeft).Column 　　
Keep reading

教師を授けられた

Image
教師を授けられた。 9月21日から30日、京都の本願寺西山別院で行われていた教師教習に参加してきた。感想というほどの感想もないのだけど、一応思ったことをメモしておこう。 初日、集合して簡単な説明があった後、居室に通されて俗服から法衣に着替え。そこに指導員の方も一緒についてきて着替えの様子をチェックするのだけど（驚くけど）、 「（真顔で）あれ、なんで下着が白じゃないの？」 「（半笑いで）え！？」 「（真顔で）案内に書いてあったよね？読んでない？」 「（ひきつった顔で）え、あぁ、書いてあったようにも...。」 「（真顔で）今日は電話を許可するので、電話して送ってもらうようにしてください。」 ま、確かに案内にそう書いてあった。だけど、わざわざ買うのももったいないし、別にどうでもいいことだろうと普段通りの下着で行ったわけだけど、まさかの展開...。公衆電話から家に電話すると父が出て、 「あ、お父さん？下着、白くないとあかんねんて。2～3枚買って送ってもらえるかな？」 とまぁ、なんというか、まぁ、なんというか...。 なぜ下着は白じゃないのいけないのか、理由の説明はないけれど、僧侶に華美な下着は不適切だから、みたいな感じだろうか。指導員の先生、赤いアルファロメオで出勤されてますけど...。ま、僧侶が所有する車の車種なんて、僧侶の下着の色と同じ程度にどうでもいいことですね。どうかな。 毎日5時半起床で、掃除、朝のお勤め、朝食、講義、昼食、講義、夕方のお勤め、夕食、夜のお勤め、翌日のお勤めの練習、課題、入浴、、、、と23時の消灯までみっちり予定がつまってた。ここ数年、人に叱られるという経験をほとんどしてなかったので、叱られっぱなしの生活はちょっと新鮮だった。 参加者は45名で、年齢層は、20代が半分、30代が2割、40代1割、50代1割、60代1割みたいな感じかな。男性が8割、女性が2割程度か。現役大学生、お寺でお坊さんをしてる人、普通の会社員をしてる人、定年退職後の人、いろんな人がいた。僧侶であることに熱い思いを持っている人もいれば、のっぴきならない事情でやむを得ず住職を引き受けることになってしまった人、なんとなくの人、参加の動機もいろいろ。 初日から「あと何日で帰れる」とカウントダウンを始める人が少なからずいるのだけど、若い人
Keep reading

サンプリングと本歌取り

Image
自由と 気持ちよさ 満ちた MUSIC 、Hiphopが好きなんですけども、好きな要素の一つがサンプリングです。ヒップホップって、もともとは楽器を使わずレコード2枚を2つのターンテーブルに乗せて、既存曲の一部をつないでトラックを作る、みたいなことをして作られた音楽なので、元ネタ（つながれてるレコードの曲、英語では sample とか sampling source とかいいます）があるんですね。貧しくて楽器が買えなくて、みたいな理由から生まれた手法ですが、これがとても素敵で、今でも基本的に全てのヒップホップの曲には元ネタがあるんです。 例えば分かりやすいところだと、 ↑この Isaac Hayes の1970年の曲、 Ike's Mood の1:45くらいから始まるピアノのフレーズが、 Mary J Blige の I Love You (1994) とか、 Foxy Brown の (Holy Matrimony) Letter from the Firm (1996) とか、 Cormega の Beautiful Mind (2004) とか、たぶん他にもいろいろありますが、サンプリングされてるんですね。 これは分かりやすいですが、1秒にも満たないサンプリングもあって、1曲に数十曲の元ネタが入ってることもあって、「わ、今の『イエー！』聞いたことある！」みたいなことになって面白いんです。というか気持ちいいんです。元ネタに気づくっていうのは、当然その曲を作ったミュージシャンと同じ曲を聞いていたことを意味するわけで、いきなり、なんていうか、体験を共有する仲間意識、みたいものを感じるんですね。これって和歌でいう本歌取りだと思うんです。またしても内田樹ですが、小説における本歌取りについて、 「先行作品を踏まえている」という事実それ自体が読者に対する「コールサイン」として機能するということである。書き手から読者への「コールサイン」はつねに「ダブル･ミーニング」として発信される。表層的に読んでもリーダブルである。でも、別の層をたどると「表層とは別の意味」が仕込んである。その層をみつけた読者は「書き手は私だけにひそかに目くばせをしている」という「幸福な錯覚」を感知することができる。 To the happy few. 自
Keep reading