タミフルと異常行動の因果性(再掲)

　今日は私の59歳の誕生日ですが,それとは全く関係なく,丁度インフルエンザの季節で新薬などのニュースもありました。

　この新薬は試験段階でまだ使用不可とのことなので,200７年の入院前日(3月23日)に書いた関連の特効薬タミフルの記事(「タミフルと異常行動の因果性,」)と言っても内容は統計学の仮説検定の話が中心ですが,それを再掲して手抜き記事でお茶を濁します。

　(※ Happy Birthday to me 。。Happy Birthday to me 。。

　Happy Birthday dear TOSHI 。。

  Happy Birthday to me 。。

　ちょっと淋しいね。何かテレビ朝日と同じ誕生日らしい。。。。

　生年月日が全く同じ芸能人なら山本譲二。。

　１歳年下なら中村雅俊がいるけど。。

　いい男ばかりだな。。。

　ちょっと離れた歳だと,女性だけど磯野貴理(子)くらいかな。。※)

※以下は再掲記事全文です。最後で入院する話を書いてますが,ほぼ２年前なので勘違いしないでください。

　今日は,今話題となっているインフルエンザの特効薬タミフルと,飛び降りや飛び込みなどの異常行動の間に,有意な因果関係があるかどうか？を判断する１つの統計的手法:仮説検定について述べてみます。

　最近のタミフル服用者の異常行動は,一見では自殺行動と判別できないように感じています。そこで因果関係の有無を調べる手段として,次のような仮説検定が有効ではないかと考えます。

　"タミフル服用者全員という総母数ｎ人に対し,これが原因で自殺行動をしたのではないかと予想される人の人数がｒ人である場合の比率ｐ＝ｒ/ｎが,日本人全員の平均自殺率:ｐ₀と比較して有意な差であると認められるかどうか"という仮説検定です。

　そのため,まず,未知の出現確率ｐを持つある確率事象がｎ回の独立試行中にｒ回起こったものとして,この未知確率ｐが確率が既知のｐ＝ｐ₀の事象と比較して差がないという仮説「Ｈ₀:ｐ＝ｐ₀」を検定する一般的問題を設定します。

　これの対立仮説としては「Ｈ₁：ｐ≠ｐ₀」を取ります。

(※この事象は有意な差があるだろう,とかこの薬は他と違って有効だろうとかいう意味のある予想(作業仮説＝対立仮説:Ｈ₁)を否定する(無に帰する)仮説が帰無仮説:Ｈ₀です。※)　

　ｎが小さいときには,こうした分布は二項分布で与えられますが,これはｎが大きいときには,中心極限定理により正規分布で近似できます。

　すなわち,母数ｎが大きい場合には,注目している事象が起これば１,起こらなければ 0 とする確率変数をＸ_i(ｉ＝1,2,...,ｎ)として,その平均を＜Ｘ＞≡ΣＸ_i/ｎとすると,仮説「Ｈ₀：ｐ＝ｐ₀」の下で,確率変数:Ｔ≡(＜Ｘ＞－ｐ₀)/[ｐ₀(1－ｐ₀)/ｎ]^1/2は,ほぼ正規分布Ｎ[0,1]に従うはずです。(大数の法則)

　そして,仮説Ｈ₀を真とみなしたとき,Ｔ∈Ｒが成り立つ確率:Ｐr(Ｔ∈Ｒ|Ｈ₀)がαに等しい,すなわち,Ｐr(Ｔ∈Ｒ|Ｈ₀)＝αとなるＴの変域Ｒを定めます。

　変域Ｒは,仮説Ｈ₀が真の場合の分布ならＲに入らない確率が(1－α)(α＝0.05なら95％の範囲)を除く領域に取ります。

　このＲを有意水準α,または危険率αの棄却域と呼びます。

　そして,実際に検定するデータでのＴがＴ^*のとき,Ｔ^*∈Ｒの場合は,仮説Ｈ₀を棄却し,そうでない場合にはＨ₀を採択するというわけです。

　今の場合では,例えばＴ^*∈Ｒのときは,通常の自殺率の範囲内からずれている,ということで「タミフルと異常行動の間に因果関係がない」という仮説は棄却されるわけです。

　具体的な例として,１つのサイコロを500回振ったとき１の目が66回出たとき,この結果から「このサイコロには異常がない」という仮説を有意水準α＝0.05で検定してみましょう。

　「Ｈ₀:ｐ＝1/6」,「Ｈ₁:ｐ≠1/6」とすると,Ｔ≡[＜Ｘ＞－(1/6)]/[{(1/6)×(5/6)}/500]^1/2ですが,この変数Ｔはほぼ正規分布Ｎ[0,1]に従います。

　そして,有意水準α＝0.05に対応して,[1/(2π)^1/2]∫₀^xexp(－ｔ²/2)ｄｔ＝(1－α)/2＝0.475となるｘの値は,ｘ＝1.96で与えられます。そこで,棄却域:ＲはＲ＝(－∞,－1.96)∪(1.96,∞)で与えられます。

　　　

そして,「500回投げたとき１の目が66回出た」という事象は,Ｔ^*＝[(66/500)－(1/6)]/[{(1/6)×(5/6)}/500]^1/2＝－2.10となります。

そこで,Ｔ^*＝-2.10∈Ｒが満足されるので,この場合,このサイコロについては仮説"Ｈ₀:ｐ＝1/6"は棄却されます。

すなわち,この検定によれば,「このサイコロは異常なし,とは言えない」という結論となります。

本題に戻ると,日本人の平均自殺率は,2003年のデータによると,10万人当たりで27人です。すなわち,ｐ₀＝27/100000ですね。

一方,タミフルの服用者の総数ｎについては,私はよく知らないし,そのうちの自殺行動者の数ｒについても現在流動的でデータを入手してはいません。

もしも,正確な異常行動者の比率ｐ＝ｒ/ｎのデータが入手できれば,こうした検定は容易に可能です。

すなわち,仮説棄却の有意水準を0.05とした場合,ｐ＝ｒ/ｎを変数Ｔに代入したときのＴの値:Ｔ^*＝(ｐ－ｐ₀)/[ｐ₀(1－ｐ₀)/ｎ]^1/2が,棄却域Ｒ＝(－∞,－1.96)∪(1.96,∞)に入る場合には,「タミフルと異常行動の間に因果関係がない」という仮説は棄却される,つまり微妙な言い方ですが「因果関係がないとは言えない」ということになります。

一方,そうでないならば仮説は採択される,つまり,「タミフルと異常行動の間に因果関係はない」ということになります。

もっとも,平均値ｐ₀としては日本人全体の自殺者数だけではなく,自殺失敗者数も含める必要があるとも考えられ,そうするとｐ₀は倍以上になると予想されます。

また,タミフルを使用するのは,インフルエンザにかかった人だけですから,日本人全体を対象とするのでなく,インフルエンザの罹患者のみを対象とした自殺的異常行動者の比率と,タミフル関係の異常行動者の比率との差異を問題にすべきではないか？とも考えられます。

実際,インフルエンザ脳症などにより,薬の副作用ではなくインフルエンザ自体が原因で精神に異常をきたす例もあるらしいです。

明日,心筋梗塞(心不全)で板橋区加賀の帝京大学病院に入院なので,もしかしたらこれが最後の記事になるかもしれません。キーボードを叩くという作業も心臓(ハート)がおかしいと結構苦しいものですね。

参考文献；藤沢武久 著「確率と統計」(日本理工出版会)

(再掲記事終わり※)

ＰＳ:結局,日曜日もやってる近くの店３軒に押しかけて無理やり自分の誕生日を押し付けて,料金もサービスしてもらって初めて会った人にもHappy Birthdayを唄ってもらいました。世間にはつくづく厚かましいヤツがいるもんだと呆れ果てましたね。。。あ？俺のことかよ。。

http://folomy.jp/heart/「folomy 物理フォーラム」サブマネージャーです。 

人気blogランキングへ　← クリックして投票してください。(１クリック＝１投票です。１人１日１投票しかできません。クリックすると人気blogランキングに跳びます。）

← クリックして投票してください。(ブログ村科学ブログランキング投票です。１クリック＝１投票です。１人１日１投票しかできません。クリックするとブログ村の人気ランキング一覧のホ　－ムページームに跳びます。）

　Attension!!　広告宣伝です。

　http://www.rakuten.co.jp/trs-kenko-land/　 健康商品の店「TRS健康ランド」－ SCS(食品洗浄剤),黒ウコン,酒の帝王,鵜鶏王などの専売店  ←この店は私TOSHIが店長をしています。(楽天ショップです。TRSのTはTOSHIのTです。)

「TRS健康ランド」では2008年１月10日よりお徳用SCS500mlを新発売！！当店の専売です。

　そこのお酒のみの方，いろいろと飲食の機会の増えたあなた,悪酔いを防止すると言われているウコンがいいですよ！！　そして特に今回提供する沖縄原産の純粋な黒ウコンは当店が専売の新製品ですが古くから沖縄地方ではいわゆる男性の力に効果があると言われています。

　おやおや,そこの静電気バチバチの人、いいものありますよ。。。

　それから農薬を落とした後の皮がピカピカに光っているリンゴなど商品として販売する際の見栄えをよくするなどのために化学処理をした食品を安全に洗浄する新商品の洗浄液SCSはいかがですか。。。農薬ジクロルポスも食品専用の洗浄液SCSで落ちて安全になります。(厚労省試験済み)

http://www.mediator.co.jp/category/pages.php?id=115「中古パソコン！メディエーター巣鴨店」

Dell-個人のお客様ページ

(Dellの100円パソコン(Mini9)↓私も注文しました。)

ヤーマン　プラチナゲルマローラー　1日3分コロコロエステ！ローラー型プラチナ配合美顔器  

　

お売りください。ブックオフオンラインのインターネット買取 展開へ！ ▼コミック　尾田栄一郎 「ONE PIECE（52）」 ▼コミック 「ONE PIECE」をオトナ買い

三国志特集　▼コミック 横山光輝 「三国志全巻セット」 「三国志（文庫版）全巻セット」   「三国志（ワイド版）全巻セット」 　▼書籍 「三国志」/吉川英治   「三国志」/北方謙三   「三国志」/宮城谷昌光