黒木玄 Gen Kuroki
- いいね数 389,756/311,170
- フォロー 995 フォロワー 14,556 ツイート 293,980
- 現在地 (^-^)/
- Web https://genkuroki.github.io/documents/
- 自己紹介 私については https://twilog.org/genkuroki と https://genkuroki.github.io と https://github.com/genkuroki と https://github.com/genkuroki/public を見て下さい。
2016年07月03日(日)
次々と繰り出される面白メカ（ギヤ）！たまらん♪
/ 回らないはずの3つの密着した歯車、3Dプリンターを使うことで回転させることに成功 : ギズモード・ジャパン www.gizmodo.jp/2016/06/3_1120... @gizmodojapanさんから
タグ：
posted at 23:05:20
まあ、いずれにしてもこういう解析を行うこともできる、という例であって、これを以て何か確定的な結論を導こうとか、そういうものではない。それに、データの解析のプライオリティはあくまで名古屋市、名古屋市大にある。
(この話、ここまで）
タグ：
posted at 22:41:05
結果はこんな感じ。年齢を共変量としたときの結果と大体同じ、という感じ。なお、本当は正確確率を求めるアルゴリズムでやるべきなのだけど、それやると天国的時間の長さになるので、Efronの近似アルゴリズムを採用。 pic.twitter.com/pPly4SZOTT
タグ：
posted at 22:40:33
年齢を共変量としても、層別化してもサンプルサイズの偏りは解消できないので、こういう場合は条件付きロジスティック回帰を行うのが方法としては最も妥当と思われる(勿論、層別化条件によってデータに変な偏りがない、という前提ではあるけど)。
タグ：
posted at 22:40:11
それに、接種が勧奨されていることによる偏りとして接種しないのは特段の健康的事情がある人に限られていたという仮定と同時に、逆に接種勧奨が中止された後だって、今度は接種した群がその後の健康異常に注意深くなってより「症状あり」としやすくなるという可能性も排除できない。
タグ：
posted at 22:39:48
ただ、気を付けなければいけないのは、このように層別化した場合に、層別でのサンプルサイズが違う、ということ。これは信頼区間の幅に影響するので、これを同列に見るのは統計的に妥当とは言えない。サイズの違うデータのp値を比較することがナンセンスなのと同じこと。
タグ：
posted at 22:39:25
ということで、接種群と非接種群の比率が逆転する16歳以下か17歳以上かでまずは層別に見てみることにした。すると、症状17～24で結構層別化することによってオッズ比に違いがあるっぽいということが。
タグ：
posted at 22:38:41
ただ、牧野さんのこの指摘は全くの的外れということはないとは考えられる。勿論、ここに書かれているような方法で対応するのはあり得ないと思われるけど。接種の積極的勧奨が中止されたことによるバイアスだって排除できないわけで。
twitter.com/jun_makino/sta...
タグ：
posted at 22:37:58
で、コレスポンデンス分析を各「症状あり」に対してやってみてプロットするとこんな感じ。明らかに各症状の発生は独立ではない。従って、データがこのような挙動を示すのはそれぞれの症状が強い共起関係にあるということでかなり説明可能と思われ。 pic.twitter.com/a2h3xP2E6l
タグ：
posted at 22:37:31
まず、大量にオッズ比が＞1から＜1になるようなものが発生する原因について、これは「偶然そうなったのか」どうかは検証しておく必要があって、それぞれの症状が独立ならこんなに＜1側に傾くことは確率的にはまれそうではあるけど、そもそもこの症状、独立とは思えない性質のものが多い。
タグ：
posted at 22:36:45
ちなみに、分割表を層別解析すると全く逆の結果になってしまうことがある、というのは「シンプソンのパラドックス」という極めて有名な話で，統計屋にとっては常識の部類に入るもので、それ自体は「謎な結果」でも何でもない。
jtsutsui.hatenablog.com/entry/20100125...
タグ：
posted at 22:36:18
まあ、それだけではこういった人の結果に文句だけつけるような人がいるので、もうちょい解析をしてみる。まあ、自称研究者なら自分で解析しろ、の一言でもよいような気もするけど。
twitter.com/jun_makino/sta...
twitter.com/jun_makino/sta...
タグ：
posted at 22:35:42
まずは症状1～24に関して、HPVVの接種有無を説明変数として症状ありかなしかについてロジスティック回帰、次に年齢を共変量としてロジスティック回帰を行った。傾向としては消された元の結果と大体同じである様に見える。 pic.twitter.com/nGK9ayOJ5P
タグ：
posted at 22:35:12
ということで、こんな感じでまずデータを読み込んで、解析対象の変数を抽出、加工した。ちなみにクソコードなのであまり良い子は真似をしないように。解析はRを使用。 pic.twitter.com/EUZN9cDdDW
タグ：
posted at 22:33:44
それで出した結論が「ニセ科学批判」批判なわけかｗ この期に及んで、すげー、貧弱な政治的決断じゃな。
twitter.com/hirakawah/stat...
タグ：
posted at 22:33:08
ただ、この資料のp22に記載されている条件を加えたけどそれでも一致しなかったので、とりあえずは判断保留。いずれにしてもきちんとしたデータのQCが必要だし、これはPDFでデータが出ている限りは我々には難しいかと思われ。
タグ：
posted at 22:32:22
ちなみに、先生が匿名の方からいただいたというデータと先生が作ったデータでも整合性は確認した(ここには示さないけど)。ただ、名古屋市の発表しているデータと少し数字がずれている。どっちが正しいのかはわからない。
www.city.nagoya.jp/kenkofukushi/c...
タグ：
posted at 22:31:56
で、素直にread.csv関数を使った。この程度のサイズなら読み込みには時間はかからなかった。症状20のクロス集計表はこんな感じになり、数値は一致。
[1] "symptom 20"
1 2
1 9107 20571
2 58 198
タグ：
posted at 22:31:35
理由は、ファイルの読み込みに時間がかかることが予想されたので、data.tableパッケージのfread関数を使ったこと。これ使うと途中で列がずれたりしておかしな結果になり、数値が先生が提示したものと合わなくなる、という現象が発生。
oku.edu.mie-u.ac.jp/~okumura/stat/...
タグ：
posted at 22:30:32
名古屋のHPVVの調査結果のデータ、奥村@h_okumura 先生が作成したcsvファイルを用いてざっと解析を行った。ここに記載されている2つのファイルで結果を比較することからやろうとしたが、
oku.edu.mie-u.ac.jp/~okumura/stat/...
のっけからドツボにはまった。
タグ：
posted at 22:30:03
ええ機会なんで、@hirakawahさんが返答くださらなかったままになってるものを列挙しておきます。
〇エートスについて
twitter.com/shanghai_ii/st...
〇rather thanについて
twitter.com/shanghai_ii/st...
（他にもあります）
タグ：
posted at 21:57:41
件の南京事件の紛糾時には、コロラドさんにはこういうコメントを頂いたのをお忘れのようだね。
twitter.com/BB45_Colorado/...
ぼくはとても危うい意見だと思ったので、次のように諌めました。
twitter.com/shanghai_ii/st...
彼には不思議な発言が多い。
タグ：
posted at 21:31:46
山田太郎 ⋈（参議院議員・全国比例） @yamadataro43
本日、秋葉原でボランティアの方が金属製の扇子のようなもので暴力行為を受けたということがネット上にあがっています。この件については、警察の方で話をしていますので、詳細分かり次第また、ご連絡します。私（山田太郎）は特にケガなどありません。元気です。
タグ：
posted at 21:01:58
「問答無用の政権を黙認する議員をこれ以上増やせば、私たちは自由を奪われる。自分の声を国会に届けてくれる使者が必要だ。」までは同意なんだが、それで推すのが山本太郎に三宅洋平というのは……知性ってなんだったっけ、という気分になる。 twitter.com/SdaMhiko/statu...
タグ：
posted at 21:00:34
三宅洋平を批判するにあたっては「右」とか「左」とかというスタンスは関係ないですねえ。彼は「下」です。 twitter.com/0556saitoh/sta...
タグ：
posted at 19:22:14
非公開
タグ：
posted at xx:xx:xx
昨日の錯視立体、もう一個作ってみた。そのまま真似ても芸がないから、180度回しても左が円筒に見えるようにしてみた pic.twitter.com/hb44QcHE0g
タグ：
posted at 17:25:45
非公開
タグ：
posted at xx:xx:xx
三宅洋平や山本太郎が #自分らしくあれる社会 の実現のために社会保証を厚くしろ、富裕層から金を取れとか言ってるけど、それならまず先にスピリチュアル系やマルチ商法系の人脈を切ってから言えよ。スピ系・マルチ系こそが末端の金のない層から金を吸い上げるシステムなんだからよ。
タグ： 自分らしくあれる社会
posted at 13:06:36
OokuboTact 大久保中二病中年 @OokuboTact
@temmusu_n 今、書店に置いてある月刊誌『新しい算数研究 ７月号』に掲載の「かけ算九九」の指導の根本的な見直し／東京フェスム支部 34」が興味深いです。 www.toyokan.co.jp/book/b239545.h...
タグ：
posted at 12:08:35
山本弘 『BIS ビブリオバトル部』 @hirorin0015
なんか改憲派も護憲派も、憲法九条の効果をとてつもなく過大視してないか？ 「九条があるから日本人は襲われない」とか「改正しなければ日本人は守れない」とか。いや、そもそもテロリストは、日本国憲法なんてぜんぜん気にしてませんから。
タグ：
posted at 11:20:03
山本弘 『BIS ビブリオバトル部』 @hirorin0015
案の定、今回の事件を日本国憲法と強引に絡めてる声が多いなあ。「憲法九条は日本人を守れませんでした」とか。いや、九条なくても今回の事件は防げてないだろ。
タグ：
posted at 11:15:55
twitter.com/shangdai1225/s...
文字を利用して自分の気持ちを表現できている時点で教育が役に立つことを自ら証明しちゃっているんですけどね。私がツイッターでたくさんの人から情報を受け取ることができているのは公教育のおかげ。あなたの勉強は私にとって役に立っている。
タグ：
posted at 09:21:57
暑い季節になったので恒例ツイート。
妊婦は冷やしちゃだめってさんざん言われると思いますが、熱中症の方が悪いので、腹巻にカイロとかレッグウォーマーとかしなくていいです。自分が快適と思う服装がベスト。冷やしたから不妊や逆子や難産になるというのは迷信です。指導している医療者は勉強を。
タグ：
posted at 09:21:55
以前、「母親が四大卒の子どもにアトピーや食物アレルギーが多い」とやってくれたのも、この手のイシキタカイ系穀潰しメディアだった。 twitter.com/a5gt/status/74...
タグ：
posted at 07:49:25
#掛算 朗報もありまず。東京書籍はサイトにこんなアホなことを書いていましたが、削除した模様です。東書の本心はともかく、おかしな意見が引っ込むことは望ましいこと。
8254.teacup.com/kakezannojunjo...
タグ： 掛算
posted at 00:12:43
@sekibunnteisuu
> 包含除の場合はａ÷ｂ＝ｃあまりｄのたしかめは，ｂ×ｃ＋ｄ＝ａ
> 等分除の場合のたしかめは，ｃ×ｂ＋ｄ＝ａ
> 計算結果のみをたしかめるときには，場合分けをしないで，除数×商＋あまり＝被乗数
相変わらず狂っているようです。 #掛算
タグ： 掛算
posted at 00:08:06
