名古屋市の子宮頸癌ワクチンアンケート調査データの解析に挑む
- parasite2006
- 70641
- 664
- 38
- 136
@atsushimiyahara 「年齢補正の誤り」とか言うてますけどねえyakugai.gr.jp/topics/topic.p… 実態は接種群と非接種群には年齢分布の大差がtwitter.com/tonkyo_Vc/stat… (2013年=平成25年6月14日の接種勧奨中止が原因)
2016-07-06 11:53:03名古屋市のHPVVの調査データ、接種群(赤)と非接種群(青)の年齢分布をグラフにしてみたけど、これは明らかに年齢分布が異なっているので、年齢調整したオッズ比出さなければいかんですね。 city.nagoya.jp/kenkofukushi/c… pic.twitter.com/cGHHzmTue8
2016-06-27 23:29:05@atsushimiyahara 年齢補正前なら接種群により高率で出る症状があるがtwitter.com/keigomi29/stat… 年齢補正後は接種群に非接種群より有意に高いオッズ比で出た症状はないとの@tonkyo_Vc さんの再解析twilog.org/tonkyo_Vc/date…
2016-07-06 12:05:29子宮頸がん予防接種、名古屋市の調査から単純集計で症状ありの割合を比較。年齢調整は回答データを読み込まないとできないが、単純集計は集計結果のP22に出てたので取り急ぎ。「接種あり」のほうが多いのはちょうど半分の12項目。 pic.twitter.com/g8oYB8PJMV
2016-06-27 14:05:42@parasite2006 @atsushimiyahara これは<意見の概要>の2がまずおかしいですね。誰も因果関係など推論していないのに「因果関係を推論するには適さない」などと話をずらしています(続)。
2016-07-06 23:04:30@parasite2006 @atsushimiyahara (承前)3もデータの構造上「HPVワクチン接種後の症状に着目した実態分析」を行うことは無理(それをやるにはRCTが必要)なのに、それを無視して「行っていない」と書いています。このデータでその解析を行う方が不適切です。
2016-07-06 23:07:57@parasite2006 @atsushimiyahara (承前)4も、実際にどのように年齢調整を行ったかが速報に記載されていないことと、奥村先生が起こしてくださったデータと数値で若干のずれがあるので何とも言えないですが、私が先日ツイートしたように(続く)
2016-07-06 23:12:43@parasite2006 @atsushimiyahara (承前)条件付ロジスティック回帰を行えば、年齢によるポピュレーションの偏りも同時に調整され、その条件下でも年齢を共変量として調整した場合と同じ結果であったこと、コレスポンデンス分析で各症状の発生に相関がある(続く)
2016-07-06 23:14:47@parasite2006 @atsushimiyahara (承前)ことが示されているので別に不思議なことではないと考えられます。加えて「ワクチン接種の影響がないのであれば差は現れないはず」は論理的におかしいですね。同じであっても5%の確率で有意な差を検出するわけですから。
2016-07-06 23:17:11@parasite2006 @atsushimiyahara (承前)もしこのように批判をするなら、より適切な解析プロトコールを示さなければ「様々な角度から徹底した分析をおこなうこと」などと書いても画餅にすらならないな、という感じです。
2016-07-06 23:20:42大崎市の開業医の反子宮頸がんワクチンブログ記事を見つけてびっくり
これは新手だなあ。この人名古屋の子宮頸がんワクチンアンケート調査の速報結果web.archive.org/web/2015122319… 見たのはいいけど、接種群と比較せずに非接種群の有症状率だけ取り上げて何が言いたいの?satouclk.jp/2016/07/post-7…
2016-07-06 13:26:18名古屋の子宮頸がんワクチンアンケート調査データを見る時は、接種群と非接種群の人数(接種なし9,245、接種あり21,034)twitter.com/tonkyo_Vc/stat… と年齢分布twitter.com/tonkyo_Vc/stat… に差があることを頭において読まないと足をすくわれる
2016-07-06 13:32:22数字は嘘をつかないが、嘘つきは数字を使う、とか言うよね。 いや、このドクターがそうだってことじゃなく、一般論としてだが。。(爆 twitter.com/parasite2006/s…
2016-07-06 13:34:47ああ「俺、戦後生まれ。今までワクチン打ってないけど、死んでないぜ。だからワクチンいらね」って言ってる、某内科循環器科の人じゃないですか...。思春期診療してない事と、自分に都合の悪い事を「ねつ造」っていう典型的な「煽り」ですね twitter.com/parasite2006/s…
2016-07-06 17:13:06@io302 なるほど、有名人でしたか。それにしてもこの人、陰謀論だーいすきのようですねsatouclk.jp/cat16/
2016-07-06 22:34:49@io302 「小児科は特に標榜しておりませんが、最近まで休日当番で小児を診ていました」satouclk.jp/cat2/ とありますが、休日に急病でお医者さんに連れてこられるこどもといえば小学校入学以前の可能性が高そう。
2016-07-06 22:52:10『小児科は特に標榜しておりませんが、最近まで休日当番で小児を診ていました』 こういうの一番危ないですね。つい最近も「小児も可」という他科医が、嘔吐で尿ケトン陽性の児に「絶対安静!」って言って1ヶ月登園禁止されてた例ありました。 twitter.com/parasite2006/s…
2016-07-07 07:59:12@parasite2006 陰謀論ですね...多変量データ(アンケートデータ)を偽造するってのは,本当はすごく難しいだよね。
2016-07-06 15:17:16(↑この医師のブログ記事には
「表4より、基準となる非接種群の有症状率を大きくし、接種群とのオッズ比が大きくならないように"操作した形跡"がわかるのである。そうして、以下のようなウソの結論をデッチ上げたのである。(中略)それでは具体的な捏造の手口はどのようなものか。アンケートは無記名である。非接種群の回答として、上の症状にチェックを入れた捏造アンケートを混ぜ込めば良い。それほど多数は必要ない。」
とあります)
@anonymous124816 非接種群の年齢分布twitter.com/tonkyo_Vc/stat… を見れば高年齢ほど人数が少ないのですから、症状が出たと回答した人数が同じだったとしても21歳(500人そこそこ)と15歳(目分量で3700人ぐらい)では有症状率が7倍違って当り前。
2016-07-06 22:42:49.@parasite2006 @anonymous124816 このブログもかなり豪快にデータを誤読していますね(ちょっと速報の結果の示し方にも問題があるとは思いますが)。これは発症がλ=年数×発症率のポアソン分布に従うと近似できるので、年齢が違えば累積の発症率が上昇(続く)
2016-07-06 23:30:00.@parasite2006 @anonymous124816 (承前)するというだけの話ですよね。試しに年齢を説明変数としてポアソン回帰もやってみたのですが、年齢の係数はいずれもほぼ0、つまり、exp(0)=1で年齢に依存せず一定確率で症状が発生しているだけでした。
2016-07-06 23:32:48