黒木玄 Gen Kuroki(@genkuroki)/2017年12月/Page 3

Cornell Alliance for Science「反ワクチン活動家に立ち向かった日本の医師が国際的な賞を受賞」 allianceforscience.cornell.edu/blog/japanese-... 村中璃子先生がジョン・マドックス賞を受賞した件。断片的に聞いていたよりえぐい話だったので、訳してみた。長文注意。

タグ：

posted at 01:26:34

「日本の反ワクチン活動に異議を唱えた同国の医師が、2017年『科学のため立ち上がった人物のためのジョン・マドックス賞』を受賞した。

タグ：

posted at 01:27:24

「京都大学医師で文筆家の村中璃子医師は、高名な医師による反ワクチン的主張を『捏造』と呼んだのち、訴訟の標的となり、嫌がらせを受け、メディアから排除され、著書の出版は撤回された。」

タグ：

posted at 01:27:43

「村中氏は『受賞は個人的にも名誉だが、それより、これが政府やメディアが態度を変える助けになることを望む。科学のために立ち上がるのは簡単ではないが、代わりの多くの友人を得た。メディアの中にすら。』」

タグ：

posted at 01:28:14

「『科学のために立ち上がった人物のためのジョン・マドックス賞』は6年目。世界のどこであれ、困難や敵意に直面しながら、公共の利益のために健全な科学や根拠を推進する人物を表彰する。世界をリードする科学誌『サイエンス』等が年1回1～2名を選出する。」

タグ：

posted at 01:29:26

「村中氏は、日本における激しい反HPVワクチンキャンペーンに挑戦した努力を認められた。反ワクチン活動家は、若い女性達がワクチン接種後に痛みや車椅子生活など謎の副作用を被ったと称するビデオを流布し、接種率は2013年の70％から現在の1％未満まで激減した。」

タグ：

posted at 01:30:11

「WHOは世界のHPV関連がん発生は年間63万件、死亡は26万6千名と推計している。同機関はHPVワクチンを手に入りうる最も安全なワクチンのひとつと認識し、各国政府に使用の義務化を推奨している。」

タグ：

posted at 01:30:59

蒼き夜に @aokiyoruni

肝心の受賞理由を書かないのはなぜでしょうね。続報を期待しましょう。村中さんに政治性がないとは思いませんが、だからこそこつこつと事実を積み重ねバイアスに留意した情報発信が大事と言うことでしょう。真っ当なジャーナリストかどうかの境目。
www.sankei.com/life/news/1712...

タグ：

posted at 01:31:21

「しかし、日本でHPVパニックが巻き起こったため、同国厚労省は『被害者への配慮で』HPVワクチンの優先的推奨を停止した。医療専門家によると、これは誰も確かめていない戯言である。 www.washingtonpost.com/news/worldview... 」

タグ：

posted at 01:31:30

「こうして日本は、反ワクチン活動家が議論に効果的に勝利するという、常軌を逸した事例を提供した。HPVにまつわる恐怖の利用は世界に広がり、医療コミュニティは、MMRワクチンと自閉症の捏造のように、これが公衆衛生を害するのではないかと恐れている」

タグ：

posted at 01:32:02

「村中氏は数か月を費やして、高名な医師である池田修一氏の主張を精査した。池田氏は人気のテレビニュース番組で、ワクチンを注射されたマウスの脳に異常な自己抗体が沈着したと称するスライドを示していた。」

タグ：

posted at 01:32:31

「スライドの根拠を追求した結果、実験に加わった研究者は、実験には各種ワクチンを注射された一匹のマウスしか含まれておらず、スライドに示されたマウスの脳はHPVワクチンを注射されたものですらないと彼女に告げた。」

タグ：

posted at 01:32:58

「また、実験のマウスは加齢と共に自然に自己抗体を産生するよう遺伝的に改変されていた。これらマウスから採取された、抗体を大量に含む血清が正常マウスの脳に噴霧され、HPVが原因と思われる『脳障害』を示す写真とされた」

タグ：

posted at 01:34:00

注）実験の経緯はWedgeの記事が詳しい。 wedge.ismedia.jp/articles/-/7080

タグ：

posted at 01:34:26

「しかし、村中氏がこの無根拠な『捏造』を暴露したのち、池田氏は東京地裁に名誉棄損の訴訟を起こした、池田氏は厚労省のHPVワクチン副作用調査班長であり、信州大学脳神経内科学の前教授で、現在は大学病院でHPVワクチン被害者と思われる患者のための外来を運営している。」

タグ：

posted at 01:35:07

「名誉棄損訴訟の結果、村中氏の書いたコラムはdropされ(訳語がわからない)、メディアから排除され、書きあがっていたHPVワクチン議論の著書は出版を拒否された。またワクチン『被害者』グループは彼女の記事を掲載するあらゆるメディアに猛烈な苦情を殺到させた。」

タグ：

posted at 01:35:41

SYSTEM Y すごすぎる！ twitter.com/yonemotominoru...

タグ：

posted at 01:35:59

「『池田医師による訴訟はメディアを怯えさせ、彼らは私を使うのをやめました』と村中氏。『被害者団体はとてもプロフェッショナルです。彼らの苦情は必ずしも編集部に行かず、出版社と利害関係にある企業の経営者や、それら企業のビジネスに影響力のある政治家に行きます。』」

タグ：

posted at 01:36:23

「『東京都知事の娘、NHKのディレクター、被害者団体の親しい友人が、私がどこに住んでいるか、誰が私の家族か、熱心に探し出しました。家族や私は脅迫を受けました。』」

タグ：

posted at 01:36:54

「『コラムを元にした著書の出版は中止されました。刊行日が発表され、著者近影が撮影され、表紙のデザインも決まっていたのに。他の8社に話しましたが、『よく書けているし非常に面白いが、今は出版できない』と同じことを言われました。』」

タグ：

posted at 01:37:50

「村中氏は、池田氏のスライドが放送された僅か2週間後に、『被害者』団体による国とワクチン製造企業への集団訴訟が起こされたという興味深い偶然の一致を指摘する。」

タグ：

posted at 01:38:12

「『日本では子宮頸がんで毎年3000名が亡くなり、1万件の子宮摘出が行われています。集団訴訟は10年かかると言われており、その間はどの政治家もHPVワクチン再開を決断しそうにありません。』」

タグ：

posted at 01:38:42

「『その10年の間に、日本の産婦人科医はいくつの子宮を摘出しなければいけないでしょうか。答えは10万です。摘出される女性、そこから生まれ母親を失う子供、そこから生まれるはずだった子供のことを想像してみてください』」(了)

タグ：

posted at 01:39:02

同じサイトに、ジョン・マドックス賞候補となった他の方々の紹介。投獄されたり、家族ともども襲撃されたり、検察から嫌がらせされたり。この中から選ばれたってことは、日本の反HPVワクチン活動は相当ヤバいと危惧されているんだと思う。 allianceforscience.cornell.edu/blog/science-u...

タグ：

posted at 01:46:55

訂正：『サイエンス』→『ネイチャー』

タグ：

posted at 01:49:00

ツリーにうまく反映されないので、訂正を再掲。ジョン・マドックス賞は、ネイチャー、the Kohn Foundation、the charity Sense about Scienceによる賞です。 twitter.com/ymori117/statu...

タグ：

posted at 02:19:03

#統計訂正。「自分だったらWBICとWAICの比較でベイズ自由エネルギーの近似計算用の逆温度β^*をどのように推定するか？」をWBIC原論文を見ながら自分なりに考えたら数分でそれらしき公式が得られました。
twitter.com/genkuroki/stat...

タグ：統計

posted at 05:16:42

#統計自分なりに得た公式は
twitter.com/genkuroki/stat...
で引用した別の人達の論文のdiscussion欄で発表されている公式と本質的に同じものに見えたので、その公式でいいのかなと判断して、
nbviewer.jupyter.org/gist/genkuroki...
の第1.6節で実装しました。

タグ：統計

posted at 05:20:58

#統計私が実装したケースは、最も簡単な確率モデルである正規分布のケースなので、WBICやi-WBICの推定先であるベイズ自由エネルギーおよびWBICとi-WBIC自身を表す公式があって、MCMCによる誤差を心配せずにそれらを比較できます。

ただし、私の計算がミスっているとバグが混入してしまう(笑)。

タグ：統計

posted at 05:25:47

#統計 nbviewer.jupyter.org/gist/genkuroki...
の第1.9, 1.10節では実際にその比較を実行しています。

すると、サンプルサイズが小さくて、事前分布の分布範囲が狭いときに、WBICとWAICのサンプルごとの揺らぎのせいで、β^*の推定値が負になる場合が出て来たので、その場合にはβ^*=0として計算しました。

タグ：統計

posted at 05:31:09

#統計 β^*が負になったりしない場合も比較をかなり行なっています。

私のコードがバグっている可能性はあると思う。しかし、私が計算した公式による計算とモンテカルロ法による別の方法による計算結果との一致も確認しているので、私の公式とコードが全く信用できないというわけではないと思う。

タグ：統計

posted at 05:34:53

#統計比較結果をみると、私の公式で計算したベイズ自由エネルギーの正確な値とWBICによるモデル選択はほぼ一致しているのですが、私の解釈によるi-WBICは補正に失敗しているように見える。これはわたしに原因があるのかどうかを知りたいので、iWBICについて詳細を知りたくなったのです。

タグ：統計

posted at 05:42:04

#統計私見では、WBICとWAICの差の揺らぎが結構大きなことが問題。WBICとWAICの差を取ると共通の "nL_n(w_0)" の成分はキャンセルします。しかし、WAICの中の <√tξ_n(t)>とV(ξ_n)の差の分の揺らぎは残ります。その揺らぎはβ^*の推定の誤差を増やします。

タグ：統計

posted at 05:47:38

#統計ベイズ統計でも

自由エネルギー=-(1/β)log分配函数

という定義は同じです。以上はこれの近似計算の話。ベイズ自由エネルギーはモデルへのインプットであるサンプルの違いで揺らぐのですが、その期待値はある種のKullback-Leibler情報量になります。続く

タグ：統計

posted at 05:52:49

#統計そのKL情報量は、モデルとサンプルを生成している未知の確率分布の違いを表す量になっています。ざくっと言えば、モデルがどれだけ誤っているかを表しています。

だから、もしも、ベイズ自由エネルギーの期待値がわかれば、その値の大小によって、真の確率分布に近いモデルを選択できる。

タグ：統計

posted at 05:56:31

#統計しかし、現実の実験・観察では、未知の確率分布に好きなだけサンプルを生成させることはできないので、ベイズ自由エネルギーの(サンプルを動かすときの)期待値は計算できません。

ベイズ自由エネルギーの期待値の推定値で代用することになります。続く

タグ：統計

posted at 05:59:15

#統計ベイズ自由エネルギー(周辺対数尤度×(-2))の期待値の推定値の候補の代表は以下の通り:

* 各サンプルごとのベイズ自由エネルギー

* 各サンプルごとのベイズ自由エネルギーの推定値
- BIC
- WBIC
- improved WBIC = i-WBIC
- sBIC

タグ：統計

posted at 06:03:56

#統計 20世紀の統計学には結構黒歴史があると思います。赤池弘次さんが1980年に述べているように、頻度主義 vs. ベイズ主義の論争は不毛でした。他にも赤池さんのAICに対するBICによる批判も不当で不毛な攻撃に過ぎなかったと言えると思います。続く

タグ：統計

posted at 06:07:13

#統計 AICとBICでは基礎となるKullback-Leibler情報量の種類が異なります。使い分ければ問題なし。

数学的道具を目的に合わせて使い分ければいいだけの話なのに、「AICよりもBICがよい」などと言う困った人達が出て来たせいで、赤池さんの貴重な時間が無駄に費やされた。不毛。

タグ：統計

posted at 06:10:18

#統計 AICやWAICの基礎は、モデルによる予測分布p^*(x)と未知の真の分布q(x)に関するKL情報量

D(q||p^*)=∫q(x)log(q(x)/p^*(x))dx

です。これが、モデルと真の分布の違いではなく、サンプルごとに作られる予測分布の予測誤差を意味しています。続く

タグ：統計

posted at 06:14:39

#統計 AICやWAICは、モデルと真の分布の違いの大きさを測るための道具ではなく、サンプルごとに作られる予測分布の誤差の大きさを測るための指標なので、真の分布に近いモデルを選択したい場合には適さない道具です。こういう自明な事柄を理由にAICが不当に攻撃された歴史がある。

タグ：統計

posted at 06:19:41

以前アメリカで「退行催眠」によって子供の頃の性虐待を「思い出した」人たちが、親を訴えるという件が相次ぎ、社会問題になったが、実はその中に、虚偽の記憶を植え付けられ実際には起きなかったことを「思い出した」ケースが少なからず存在した。それを研究したロフタス氏もマドックス賞の受賞者。

タグ：

posted at 06:22:46

#統計 BICやWBICの基礎は、ベイズ統計におけるサンプルx_1,…,x_nに関する逆温度1での分配函数(周辺尤度、これは確率密度函数になる)

Z_n=Z(x_1,…,x_n)=∫p(x_1|w)…p(x_n|w)φ(w)dw

と真の分布q(x)の直積

q_n=q(x_1)…q(x_n)

に関するKL情報量です。続く

タグ：統計

posted at 06:24:29

#統計ここで、p(x|w)は推定に用いられるパラメーターwを持つxに関する確率密度函数で、φ(w)はパラメーターに関する確率分布(所謂事前分布=テスト函数)です。例えば正規分布モデルでは、wは平均と分散の組み合わせになり、p(x|w)は正規分布の確率密度函数になります。

タグ：統計

posted at 06:27:42

そのロフタス氏を「大物政治家や権力者を救う役割を担う」と揶揄するツイートを子宮頸がんワクチン被害者連絡会の事務局長が「拡散希望」としている twitter.com/toshi2133/stat... 。科学的に間違っているからではなく、「政治家や権力者に有利だから」が判断基準。

タグ：

posted at 06:29:21

#統計分配函数と呼ばれる確率密度函数Z_nは、p(x|w)とφ(w)で記述されたモデル世界におけるサンプルx_1,…,x_nのランダム生成を記述しています。それとq_nに関するKL情報量D(q_n||Z_n)はモデル世界の確率分布Z_nによるサンプル生成による現実世界の確率分布q_nによるサンプル生成の～続く

タグ：統計

posted at 06:33:05

名取宏（なとろむ） @NATROM

@yusuke_tsugawa 以前に考察したことがあります。「九州や北海道などの原発事故の影響をほとんど受けていない小児集団を対照にする目的で甲状腺がん検診を行うべきではない理由」[ h.hatena.ne.jp/NATROM/2281746... ]。

タグ：

posted at 06:34:07

なるほど、「性的虐待の加害者」とされる人を弁護し、「被害者」とその支援者を「嘘つき」として糾弾するように見える内容を「科学的に正しいから」と主張するのは勇気の必要なことだろう。そのロフタス氏と、村中氏は同様のことをしたとみなされているわけだ。

タグ：

posted at 06:34:37

#統計～シミュレーションの誤差を表しています。ザクッと言えば、p(x|w)とφ(w)で記述されたモデル世界がどれだけ現実世界をうまく近似しているかを表しています(D(q_n|Z_n)が小さいほどモデル世界は現実世界をよく近似している)。続く

タグ：統計

posted at 06:36:00

#統計だから、BICやWBICによるモデル選択は、現実世界の未知の分布に近いモデルを選択するために適した道具。

実際の応用で、モデル間の相違が大きいときには、AIC、WAIC達とBIC、WBIC達のような精密な道具を使わなくても、尤度に巨大な違いが出るのでどちらのモデルが最もらしいか判明します。

タグ：統計

posted at 06:40:45

#統計 AIC、WAICのそれぞれと、BIC、WBICのn大での漸近挙動の最初の主要項は一致しています。それらと予測分布の対数尤度の-2n倍の漸近挙動の最初の主要項は一致します。

AIC、WAICやBIC、WBIC以前に、尤度によるモデル選択について知っておかないと、色々誤解してしまう可能性がある。

タグ：統計

posted at 06:47:33

#統計訂正：「最もらしい」→「もっともらしい」

タグ：統計

posted at 06:48:34

#統計以前も強調したことですが、赤池弘次さんも言っているように、AIC以前の問題として、尤度がどうしてモデルのもっともらしさの指標になっているかを理解するためには、Kullback-Leibler情報量に関するSanovの定理を理解しておく必要があります。

タグ：統計

posted at 06:51:05

#統計 KL情報量が予測誤差を意味すること(Sanovの定理)を理解していない人には以上の雑談は全て何も理解できなかったはずです。

KL情報量に関するSanovの定理は、大数の法則や中心極限定理に匹敵するくらい統計学を学ぶときの必須の教養だと思われます。

タグ：統計

posted at 06:53:48

#統計余談：Sanovの定理などに触れずに、いきなり「Kullback-Leibler情報量(もしくはShannon情報量)を最小化する分布」から出発して物理的なカノニカル分布を説明したつもりになる議論はひどい誤りです。計算が簡単になって印象的なので騙される可能性がある。続く

タグ：統計

posted at 06:58:51

#統計続き。exp(-βE)型のボルツマン因子の導出は、i.i.d.という特殊な状況における大偏差原理(Sanovの定理)に出て来るKL情報量の特別な形には依存しておらず、「大偏差原理が成立している」という一般的な物理的仮定だけから得られる結果です。普通の統計力学の教科書には～

タグ：統計

posted at 07:03:07

#統計～、「大偏差原理が成立している」という一般的な物理的仮定だけに依存するボルツマン因子の導出が書いてあります。KL情報量やShannon情報量最小化による簡単な計算によるボルツマン因子の導出は、統計力学の教科書的導出よりも適用可能範囲が狭いのです。計算の簡明さ騙されちゃダメ。

タグ：統計

posted at 07:06:47

WideExit @wideexit

水ぼうそう患者が激減　無料の予防ワクチン効果か www.asahi.com/articles/ASKCN... 明確なエビデンス出て大変結構。『特に１～４歳の子どもの患者が減った。原則無料で受けられるようになった予防ワクチンの効果』との事。子供の貧困状況を鑑みると予防接種の無料化は少なくとも幼児では重要。

タグ：

posted at 07:11:21

#統計赤池弘次さんによる尤度の概念の説明については以下のリンク先で紹介しておきました。1980年の論説なのでそろそろ常識になって良い頃だと思う。統計学教科書における尤度概念の説明では赤池さんの説明を引用するべき。

twitter.com/genkuroki/stat...

タグ：統計

posted at 07:15:35

#統計予測分布p^*(x)のサンプルX_1,…,X_nに関する対数尤度の定義は

Σ_{i=1}^n log p^*(X_i)

です。X_iの各々が確率密度函数q(x)の確率分布によって独立生成されているならば、対数尤度の1/n倍は大数の法則よりn→∞で

∫q(x) log p^*(x)dx

に収束します。続く

タグ：統計

posted at 07:25:02

#統計すなわち

(対数尤度) = n ∫q(x) log p^*(x) dx + o(n)

です。最初の主要項の係数の-1倍は

-∫q(x) log p^*(x) dx = D(q||p^*) + S(q)

と表せます。ここでD(q||p^*)はKL情報量でS(q)はShannon情報量です。続く

タグ：統計

posted at 07:28:15

#統計それらの定義は

D(q||p^*)=∫q(x) log p^*(x) dx,
S(q) = -∫q(x) log q(x) dx

です。KL情報量 D(q||p^*) は予測分布p^*(x)による未知の真の分布q(x)の予測誤差を意味することがSanovの定理からわかります。続く

タグ：統計

posted at 07:31:00

#統計 S(q)は予測分布p^*(x)によらない定数なので、予測分布の予測誤差の大小を比較する場合には影響しません。

対数尤度の漸近挙動の最初の主要項に予測誤差を表すKL情報量が出て来ます。これが「対数尤度が予測分布のもっともらしさを表す数学的理由」です。

タグ：統計

posted at 07:34:16

#統計実用的には、尤度(ゆうど、いぬどではない)は対数尤度の形式で使われるのですが、対数尤度が確率モデルのもっともらしさを表す理由は、

* 大数の法則と
* Kullback-Leibler情報量に関するSanovの定理

にあるわけです。

タグ：統計

posted at 07:37:24

砂___の___女 @vecchio_ciao

「教育」業界で偉くなるには、いかに根拠のないことでも、何となく感じたことを自信満々に言い切るかがポイント、っていうのがよくわかる。

タグ：

posted at 07:46:29

@yusuke_tsugawa @ROCKY_Eto @NATROM 「事実を明らかにする」ことは検診の目的として不適当でしょう。そういう理由で検診するのは倫理的に認められないと思います

タグ：

posted at 07:48:23

いよいよ明日です！特に学部生、院生の方歓迎！
田崎さん@Hal_Tasaki の東北大の一般講演です。
日時：１２月４日１７：３０～１９：００
会場：東北大学理学部物理講義棟３F第２講義室（MAPのH-24）
講師：田崎晴明(学習院大）
タイトル：Haldane現象と「その先」 pic.twitter.com/hFddwzxoz7

タグ：

posted at 08:10:01

明後日は英語ですが、時間の矢に関する講演も田崎さんにお願いしております。twitter.com/Hal_Tasaki/sta...

タグ：

posted at 08:11:49

佐々木俊尚 @sasakitoshinao

たいへん素晴らしい受賞。しかし新聞・テレビメディアではほとんど報道されてませんね。／海外の一流科学誌「ネイチャー」HPVワクチンの安全性を検証してきた医師・ジャーナリストの村中璃子さんを表彰 bzfd.it/2iAO3El

タグ：

posted at 08:19:00

koji hasegawa @myfavoritescene

@hottaqu @genkuroki @Hal_Tasaki 出張につき残念です。

タグ：

posted at 08:37:23

@myfavoritescene @genkuroki @Hal_Tasaki あらっ。残念です。

タグ：

posted at 08:37:57

@yusuke_tsugawa @ROCKY_Eto @NATROM 検診はあくまでも個人の利益が最優先でなくてはなりません。医学研究に関するヘルシンキ宣言にしたがうべきです。要するに「人体実験禁止」ということです

タグ：

posted at 08:48:54

ぶたやま@「ぶたやまかあさんのやり過ごし @Butayama3

＃超算数
近隣中学の数学期末試験問題

次の（　）にあてはまるものを書きいれなさい。

（①　　　　）が等しい２つの比と、等号を用いて表した式を（②　　　　）という。
また、②の性質として、a:b＝c:dならば方程式（③　　　　）が成り立つ。

タグ：超算数

posted at 08:53:14

Nature誌が選んだ賞ってことね。
Nature誌に載ったってことではないのね。
でもNatureに載るのがどれだけ難しいかわかると、その選択基準もわかるよね。 twitter.com/usg_ringo/stat...

タグ：

posted at 08:56:10

ジョン・マドックス賞、「それでも地球は回っている」と言う人に与えられる賞なのかな、という印象。

タグ：

posted at 08:58:46

@Butayama3 Natureが共催している元編集長の名前を冠したジョン・マドックス賞を村中さんが受賞して、Sienceがインタビュー「日本の医師が反ワクチン活動家との戦いで受賞」というタイトルの記事だしてるの。これは現時点では科学としては世界的にガチであることを意味するわけよ。

タグ：

posted at 09:02:31

非公開

タグ：

posted at xx:xx:xx

@Butayama3 ジョン・マドックス賞は「科学のために立ち上がるStanding up for Scienceことを推奨する賞」なのさ。

タグ：

posted at 09:06:12

「反ワクチン活動家」と名指された方々は本望かな？

タグ：

posted at 09:07:16

あおじるPPPP @kale_aojiru

新興宗教か？ twitter.com/Kageyama_hideo...

タグ：

posted at 09:10:06

@yusuke_tsugawa @ROCKY_Eto @NATROM 過剰診断のリスクは個人にとって大きな不利益であり、これを上回る公益など、今の問題では見当たりません

タグ：

posted at 09:18:41

これを報道できたらメディアを見直すんだけどな。 >RTs

タグ：

posted at 09:27:51

佐々木俊尚 @sasakitoshinao

反子宮頸がんワクチン運動と村中さんのジョンマドックス賞受賞から浮き彫りになった非常に重大な日本のオルトファクト問題。これを無視しておいて「フェイクニュースをなくせ」と呼びかけてる新聞は何言ってるんだか…自分の足元をもう一度見ましょう。

タグ：

posted at 09:35:27

#超算数 pic.twitter.com/HZzYHI2LUg

タグ：超算数

posted at 09:35:31

理由を知りたいので投げてみる。なぜ村中璃子氏のジョン・マドックス賞受賞をニュースとして扱わないのでしょうか。十分なニュースバリューがあると思うのですが。　@nhk @asahi @mainichi @asahi_apital @Yomiuri_Online @yomiDr @tbs_newsbird @NTV_shakaibu @tv_asahi_houdou @fujitv @tx_news

タグ：

posted at 09:41:18

#統計

問題：尤度を確率モデルのもっともらしさの指標のひとつとみなせる数学的理由を説明せよ。

この問題に答えられるようになってから、AICなどの(対数尤度より)精密な情報量基準について勉強するようにした方が良いと思う。上の問題に答えることができない人にはAICなどの勉強は早過ぎる。

タグ：統計

posted at 09:41:25

@nhk @asahi @mainichi @asahi_apital @Yomiuri_Online @yomiDr @tbs_newsbird @NTV_shakaibu @tv_asahi_houdou @fujitv @tx_news 反ワクチン側にとってもこの受賞はおおごとだから、どっちから見てもニュースバリューはあると思うのだよなあ。

タグ：

posted at 09:42:08

@nhk @asahi @mainichi @asahi_apital @Yomiuri_Online @yomiDr @tbs_newsbird @NTV_shakaibu @tv_asahi_houdou @fujitv @tx_news 今まで反ワクチンの姿勢で報道していたから無視、なんてのは理由にならんよねえ。

タグ：

posted at 09:43:49

#超算数　おまけ
これは、（x+60)5　がまずくて（x+60)×5　だったら良かった、ということだと思う。算数的な掛け算の順序としては間違っていないけどね^^ pic.twitter.com/aaFaSHVQJz

タグ：超算数

posted at 09:47:29

#統計赤池弘次さんは、(最尤法(尤度最大法)の開発普及者でもあるあの偉大な)フィッシャーさんの尤度概念への理解の足り無さを厳しく指摘しています。20世紀の数学の発展普及のおかげで、赤池さんの世代でやっと尤度概念をまともに理解できるようになったのだと思う。

ismrepo.ism.ac.jp/index.php?acti...

タグ：統計

posted at 09:47:46

@nhk @asahi @mainichi @asahi_apital @Yomiuri_Online @yomiDr @tbs_newsbird @NTV_shakaibu @tv_asahi_houdou @fujitv @tx_news 反ワクチン側アカウントを見ると、賞そのものやら選考方式やらジョン・マドックス氏(Natureの編集長を務めた)をくさしているんだが、でも腐ってもNatureですぜ。2大論文誌のひとつですぜ。賞の内実を別にしてもとそこの出す賞を日本人が受賞という事実にニュースバリューがないとは思えない。

タグ：

posted at 09:51:30

#超算数　そもそも私、「比例式」という言葉を知らなかったことに気づいた。y=3x　みたいなのを「比例式」と言うのかと思っていたけど違うらしい。　ja.wikipedia.org/wiki/%E6%AF%94...

a/b＝c/d　が比例式なら、2/3＝4/6も比例式で、√4/1＝2/1も比例式？
となると、√4＝2　も比例式？　　#どーでもいい

タグ：どーでもいい超算数

posted at 09:55:21

#統計用語解説【情報量】

(Shannonの意味での)情報量の最もシンプルな定義は

確率の対数の-1倍

です。KL情報量が

(確率の対数の-1倍)=n×(KL情報量)+o(n)

の形で現れるというのが、大偏差原理の一種のSanovの定理。nはサンプルサイズ。

補足：確率の対数そのものはエントロピー。

タグ：統計

posted at 09:58:53

#統計用語解説【大偏差原理】

大偏差原理が成立するとは、系のサイズを表すパラメーターλを大きくするとき、

(確率の対数) = -λ(何か) + o(λ)

の形の漸近挙動が成立していることである。「何か」の部分はよく「エントロピー」(もしくは「エントロピー密度」)と呼ばれる。続く

タグ：統計

posted at 10:03:15

禁煙ポンスケ＿心臓リハビリへ転職 @janeway_tlr2

#統計続き。「何か」の-1倍はよく「情報量」と呼ばれる。

i.i.d.の文脈では系のサイズλはサンプルサイズとすることが多い。

統計力学に文脈では系のサイズとして体積を採用したりする。

系のサイズを大きくしたときの確率の漸近挙動を見るのが大偏差原理。弱い意味でのラプラスの近似。

タグ：統計

posted at 10:06:59

村中璃子氏がジョン・マドックス賞とるまで、なぜ、今、ツイッター内で賞賛している医者たちは自信をもって推奨できなかったか？なんて、ツイートをみたが、産婦人科・小児科中心の医師たち、有識者の方々は、かなり推奨していましたよ。その患者を思う「必死のつぶやき」を見逃していただけでは？

タグ：

posted at 10:11:52

#統計用語解説【弱ラプラス近似もしくは弱ラプラスの方法】

あるパラメーターλを大きくしたとき

(ある量の対数) = λ×(何か) + o(λ)

が成立していれば、「(ある量の対数)/λ」を「何か」で近似できる。この近似を弱ラプラス近似と呼び、この近似法を弱ラプラスの方法と呼ぶ。

タグ：統計

posted at 10:12:41

⚡️ "森勇一(ymori117)さんが訳した『反ワクチン活動家に立ち向かった日本の医師が国際的な賞を受賞』"
twitter.com/i/moments/9371...

タグ：

posted at 10:13:49

#統計続き。弱ラプラス近似は

(ある量) = exp(λ×(何か) + o(λ))

を

exp(λ×(何か))

で荒く近似することと同等である。本来のラプラスの方法よりもずっと荒い近似なので「弱」という形容詞がついている。大偏差原理は確率が弱ラプラス近似を満たしていることを意味する。

タグ：統計

posted at 10:17:29

残り52件のツイートを見る

io302 @io302

hpvワクチンの議論は、ワクチン不信&反対者が、危険性や不安の煽り（根拠なし）を行い、それが次第にワクチン勧める人達の誹謗中傷に至り、私のようにIPアドレス探られたり、勤務先を探られる、脅迫行為に進展した事。

ワクチンの『事実』は必要なく、自分たちを支持する『印象』に飛び付いた事。

タグ：

posted at 10:23:55

#統計用語解説【スターリングの公式】階乗に関する

log n! = n log n - n + (1/2)log n + √(2π) + O(1/n)

という近似公式をスターリングの公式と呼ぶ。 n log nの項はShannonエントロピーの由来になっている項である。ラプラスの方法によるラプラス近似で容易に証明される。

タグ：統計

posted at 10:24:02

あらためて読み返すと、長い…。 .@tyurukichi_AA さんの「反HPVワクチン活動家に科学で立ち向かった日本医師への言論弾圧」をお気に入りにしました。 togetter.com/li/1177678

タグ：

posted at 10:26:36

#統計スターリングの公式の証明

n! = ∫_0^∞ x^n e^{-x} dx を部分積分で容易に示せる。

x^n e^{-x} dx に x=n+√n yを代入すると

n^n e^{-n} √n exp(-y^2/2 + o(1)) dy の形になる。

ゆえに、

n! = n^n e^{-n} √n ∫_{-√n}^∞ exp(-y^2/2 + o(1)) dy.

続く

タグ：統計

posted at 10:29:34

#統計続き。そしてn→∞で

∫_{-√n}^∞ exp(-y^2/2 + o(1)) dy
→ ∫_{-∞}^∞ exp(-y^2/2)dy = √(2π).

これで

n! = n^n e^{-n} √n √(2π) (1+o(1))

が示された。両辺の対数をとると、

log n! = n log n - n +(1/2)log n + log√(2π) + o(1).

q.e.d.

タグ：統計

posted at 10:33:25

#統計訂正：スターリングの公式の用語解説の式の √(2π) の前にあるべき log を書き忘れていた。正しい公式は

log n! = n log n - n +(1/2)log n + log√(2π) + 1/(12n) + O(1/n^2).

1/(12n)の項まで出すことは上の議論のちょっとした拡張で可能です。

タグ：統計

posted at 10:37:01

Yamanami @yamanamitakeshi

①は「知らん」②は「解らん」③は「a/b=c/d」だけど多分バツで、正解は「比の値」「比例式」「ad=bc」なんだろうなー、と思っていたら案の定だったが嬉しくない。
#超算数 twitter.com/sekibunnteisuu...

タグ：超算数

posted at 10:49:49

twitter.com/arapanman/stat...
私自身は小中の頃は、勉強させる目的は子どもに苦行を与えることだと割と素で思っていた。そうすることで働くようになったときの耐性をつけるためだと思っていた。

そう思ってしまったのは、私なりに合理的に解釈しようとした結果である。

タグ：

posted at 10:53:18

kei natsui @kei_natsui

田中康夫さんと浅田彰さんは、この件に関しては認識を改めてほしいし、間違った認識を拡散したことに対して何らかのアクションが欲しいところです。浅田さん、土人とか言って大衆のことを罵るんだから、自分の無知蒙昧もちゃんと反省してほしい。 twitter.com/genkuroki/stat...

タグ：

posted at 10:57:11

garmy @garmy

ちなみにBuzzfeed自身は子宮頸がんワクチンに関しては、反「反ワクチン」の記事のみでているが、これは石戸諭さんが2016/12に取り上げ始めてから時々掲載されている。石戸さんは毎日新聞出身だが、記事の中ではメディアが引き金と述べている www.google.co.jp/search?q=%E5%A...

タグ：

posted at 11:01:22

取りあえずxに色々当てはめる。-1より大きくて0より小さいまでは分かる。分数だと予想して何とか答えを出す。

この手の問題を色々やって、

ax+b＝cx+d　まで一般化する。

ここに至るまで相当の時間を費やしたが正解に行き着いた。

タグ：

posted at 11:11:21

以上はある生徒の例。

その後、学校でもこのあたりのことが教えられるが、余裕で理解できたようだ。そりゃそうだろう。

タグ：

posted at 11:12:10

この過程を経ないと、「なんだかよく分からないが手順を覚えてそれに従って答えを出せば正解になるらしい」となってしまって、下手したら「方程式とは、当てはまる未知数を求めること」ということすら認識しないままに、「正解になる答案」を書く術だけ身につけることになりかねない。

タグ：

posted at 11:14:50

the_spoiler @don_jardine

宮原篤 6th「小児科医ママとパパのやさ @atsushimiyahara

村中さんのケースは震災後強化されたマスコミとデマ拡散者の共同戦線に厚労省が加わった強大な相手と戦ったわけであり様々な妨害活動にも負けなかったその知力と胆力、何より理性の力は大いなる称賛に値する。これを評価したネイチャーも同様。。一方でこれを無視する言論人たちは全く信用に値しない。

タグ：

posted at 11:16:29

色々と意見はあると思いますが、読売新聞で森内教授のコラムが急に削除されたり、名古屋市のデータ解析が急に非公表になった背景がわかると思います。4年前の当初から..「「反HPVワクチン活動家に科学で立ち向かった日本医師への言論弾圧」」togetter.com/li/1177678#c44...

タグ：

posted at 11:20:35

MMR騒動の時にイギリスのメディアが果たした役割をHPVVでは日本のメディアがまるっと踏襲したのですよ。メディアは事実と向き合って下さい。あなたの文字が声が確実に誰かの将来を閉ざします。ひいては今不調に苦しんでいる少女たちが活動家に食い物にされることからも救います。

タグ：

posted at 11:49:43

我が家の娘2人HPVV接種済みです。下の娘と同級生達は接種の時期にいるのですが他のご家庭はだいたい、怖いからうちは接種させない！と言っています。
話をする機会があるお母さん達には、我が家は接種済みな事や息子にも接種する予定である事、なぜ必要と思い接種したのかをお話ししています。

タグ：

posted at 11:49:50

上のお嬢さんがHPVV接種していて片頭痛で苦しむのはそのせいかと心配するお母さんが下の子には接種させたくないと言っていました。その下のお子さん、接種していないけど片頭痛に苦しんでたりする。その状況でも上のお嬢さんの片頭痛の原因はHPVVじゃないのかって心配し続けているんだよ。

タグ：

posted at 11:54:24

体調が優れないと原因はHPVV接種じゃないかと思ってしまうような状況を作ったのは誰だ。

タグ：

posted at 11:56:33

今大学生の娘がHPVV接種する時と違い下の娘が接種を受けようとしたら凄い大変でした。すんなり受けられないの。自治体に問い合わせたり、意思確認したり、ワクチン取り寄せてもらったり。

タグ：

posted at 12:03:30

OokuboTact　大久保中二病中年 @OokuboTact

【「逆順」だけどマル、という答案】
拡散されないだけで、当然あるはず、という話になると考えます。
だって「心ある一部の教師は、順序強制には反対」とされているのですから。 #掛算 twitter.com/sekibunnteisuu...

タグ：掛算

posted at 12:40:25

＃超算数　　瀬山士郎の『読む数学記号』が角川ソフィア文庫になっている

タグ：超算数

posted at 13:02:08

#掛算　 #超算数　【#拡散希望】
「掛け算順序でバツ」の画像が拡散するのは晩秋の風物詩だけど、「逆順だけどマル」というケースもあるので、そういう例があれば画像をアップして欲しい。

異常採点だけじゃなく、まともな採点も、研究資料としては重要です。

タグ：掛算超算数

posted at 13:04:08

【最初からそれは無駄だと思う人もいればあった方がいいという人もいる】
それなのに、全員がいったん「仮組」を使うことを強制されているから、こうして批判が起きているのです。#掛算 twitter.com/Crossrib/statu...

タグ：掛算

posted at 13:49:51

@Crossrib 【みんな方法がバラバラだったらそれこそ批判を招くんじゃないですか】
体育の鉄棒で「順手でも逆手でも、好きなほうでやりなさい」とされる事が好ましいと私（達）は思うのですが、そうではない人もいるのですね。 #掛算

タグ：掛算

posted at 14:17:14

（あくまでも思考実験なので、リプは自分宛てのみ）
仮に「鉄棒で後々いろんな技を習得するには最初は順手でやった方が良い」などというエビデンスが出てきたりしたら、１年生も全員順手でやりなさいという指導になったりするのだろうか？
そしてそれは支持されるべきなんだろうか？　#掛算

タグ：掛算

posted at 14:34:49

Alsnova @alsnova

@ymori117 NHK勤務の猪瀬直樹の娘でしょうかね？

タグ：

posted at 14:46:03

結城浩 / Hiroshi Yuki @hyuki

このツイートに三つの写真が付いていますが、三枚目（楽しそうにやってるのは止めたくなっちゃう）に一番ゾッとしました。こんな親御さんも実際にいらっしゃるのですね。番組を盛り上げるために雇われたサクラだったらいいのですけれど… twitter.com/arapanman/stat...

タグ：

posted at 15:01:35

io302 @io302

ぐるぐるうづまき東V30b・2日目・C @guruguruuzumaki

受賞を聞いても『ジョン・マドックス賞ははじまって5年足らず』『トホホのお話』という反応で、相変わらず村中さんへの誹謗中傷を繰り返してますね。結局同じ事の繰り返しだな、この人 pic.twitter.com/BQIcOoYb4Q

タグ：

posted at 15:13:34

どこのメディアでもいいから、子宮頚がんワクチンに対する今までの報道は間違ってたってキチンと総括できる番組を作れば、少しは信頼回復できると思うんだけどなぁ。

タグ：

posted at 15:19:40

けるちゃ @kerutya

酔ってきたのでマドックス賞がマッドマックス賞に見えて脳内が一気にこうなった pic.twitter.com/MbY41ItglN

タグ：

posted at 16:02:52

訂正：原文を読み直したら単数形だったので、住所を探し出した人物は一人のようです。「NHKのディレクターで、被害者団体の親しい友人でもある、東京都知事の娘」に訂正します。

タグ：

posted at 16:30:16

訂正その2。 twitter.com/ymori117/statu...

タグ：

posted at 16:31:32

@hottaqu @myfavoritescene @Hal_Tasaki 私は行きます。

タグ：

posted at 16:41:57

@genkuroki @myfavoritescene @Hal_Tasaki 歓迎です。懇親会も是非。

タグ：

posted at 16:42:22

@yusuke_tsugawa @ROCKY_Eto @NATROM 検診をやめて、徹底的に説明する以外にないでしょうね。過剰診断リスクがなければ、検診を続けるという解もありますけど、検診そのものに害があるので

タグ：

posted at 17:04:00

BuzzFeed：「救えるはずの患者を救えない」子宮頸がんワクチン副作用「問題」はなぜ起きた？ www.buzzfeed.com/jp/satoruishid... @satoruishidoさんから改めて読み直すと、大事な論点がおおむね網羅されていて、本当に良い記事だと思う。

タグ：

posted at 17:33:40

照ZO(厨2助詞) @terzoterzo

伊藤隼也や内藤聡と言った連中、詐欺罪で立件出来んモンかねェ…。少しでもこう言う連中、法的に処罰せんと、次から次へと様々なモンを狙って喰い物にして行くやろし、歯止めも掛からん様に成ってくやろ…。→

タグ：

posted at 18:52:36

禁煙外来医師 @happylife0951

紙巻きタバコでも、IQOSでも、液体のニコチンを使ったPloomTECHでも、排出されたニコチンは空気中の亜硝酸と反応しニトロソアミンと言う、強力な発がん物質になりばら撒かれる。

#家庭内の大気汚染PM2.5を考える

タグ：家庭内の大気汚染PM2

posted at 18:54:35

@GreatDemon1701 @genkuroki ja.wikipedia.org/wiki/%E6%96%B9...　によれば、文字と等号を含んだ文字式は全部「方程式」となりそう。

そうするとx^2+y^2＝25　などを「円の方程式」と言うのも合点がいく。
「方程式」だからといって必ずしも未知数を求めることは要求されていなくて、未知数を求めるのは「方程式を解く」と言うみたい。

タグ：

posted at 19:24:49

Riko Muranaka/村中璃子 @rikomrnk

@GreatDemon1701 @genkuroki ＞恒等式とは異なり、方程式は変数の取り得るすべての値に対して等式が成り立つ必要はない

全てについて成り立っても構わないんでしょうね。
方程式のなかで、「方程式の解が、変数の取り得るすべての値」となるものが恒等式　となるのかな。

こんな言葉の用法など、あまり意識しない。#どーでもいい

タグ：どーでもいい

posted at 19:27:52

素晴らしい援護射撃のリード。岩永さん、ありがとう。

ネイチャーは日本の状況について「このワクチンの信頼性を貶める誤った情報キャンペーンが全国的に繰り広げられた」と厳しく批判しています。 twitter.com/nonbeepanda/st...

タグ：

posted at 19:28:50

Sudoku Smith @SudokuSmith

Riko Muranaka/村中璃子 @rikomrnk

日本でのHPVパニックって、明らかに人為的な扇動と似非科学の勝利でしかない。こんなことで人の命が失われることをなぜマスコミは問題視しない？ twitter.com/ymori117/statu...

タグ：

posted at 20:18:29

Robert Geller; ロバート・ @rjgeller

「オブザーバー」紙は受賞式から2日遅れの今日、私のマドックス賞授賞スピーチにも触れながら、SNSを用いた”度を超えた敵意に満ちた”反子宮頸がんワクチンキャンペーンついて記事を掲載しています。... fb.me/1rX3ne8Eb

タグ：

posted at 20:35:06

日本ではHPVワクチンの適応率は70％から1%に落ちた。原因は客観的根拠がない噂の広まり。@rikomrnk さんは正確な情報を社会に提供した努力が高く評価されマドックス賞を受賞した。おめでとうございます。適応率の改善を期待する。 www.theguardian.com/society/2017/d...

タグ：

posted at 21:05:37

Riko Muranaka/村中璃子 @rikomrnk

その始まって5年で他に誰がどんな理由で貰っているのか確認すれば良いのに。 twitter.com/io302/status/9...

タグ：

posted at 21:26:53

日本時間の明日（12月4日）夜10時半より、BBCワールドのニュース番組「インパクト」にスタジオから生出演し、子宮頸がんワクチン問題について話します。
www.bbc.co.uk/programmes/n13...

タグ：

posted at 21:30:18

ぶっちゃけジョン・マドックス賞で一番批判されてるのは国（厚労省）とメディアだ。だからって何のコメントも報道もしないのはもんのすごーーーーく異常な事態（官民一体の情報統制？）。このままなら、この国にはそういう行政とメディアばかりだってことは心に深くメモすることになる。

タグ：

posted at 21:35:27

杏子 @kyoko_shirota1

今まで子供には有料の予防接種も受けさせてきたが、HPVワクチンはテレビでの大量の薬害報道に怖気付き、受けさせなかった。
でも、もう受けさせよう、決めた。
www.sankei.com/life/news/1712...

タグ：

posted at 21:45:31

skon @skonn96

しっかし子宮頸がんワクチン(HPVワクチン)副反応の被害を訴える少女たちに“寄り添って”各方面に圧力をかけていた方々は「我々は反ワクチンじゃない」が決まり文句だったけど、結局ネイチャーからも反ワクチンのお墨付きを貰う羽目になりましたな。

タグ：

posted at 22:42:43