情報更新
ツイートの記録を停止しています

 

ツイート検索

 

@tonkyo_Vc
サイトメニュー
Twilogユーザー検索
新規ツイートの取得を再開しました!取得再開にはログインが必要です。

Twilog

ツイートの記録を停止しています

このアカウントはTwitter APIの仕様変更の影響でツイートの記録を停止しています。
記録を再開するには、Twilogにログインしてください。

 

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

Stats Twitter歴
4,763日(2011/03/15より)
ツイート数
222,024(46.6件/日)

ツイートの並び順 :

表示するツイート :

2016年07月03日(日)20 tweetssource

7月3日

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

まあ、いずれにしてもこういう解析を行うこともできる、という例であって、これを以て何か確定的な結論を導こうとか、そういうものではない。それに、データの解析のプライオリティはあくまで名古屋市、名古屋市大にある。
(この話、ここまで)

posted at 22:41:05

7月3日

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

結果はこんな感じ。年齢を共変量としたときの結果と大体同じ、という感じ。なお、本当は正確確率を求めるアルゴリズムでやるべきなのだけど、それやると天国的時間の長さになるので、Efronの近似アルゴリズムを採用。 pic.twitter.com/pPly4SZOTT

posted at 22:40:33

7月3日

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

年齢を共変量としても、層別化してもサンプルサイズの偏りは解消できないので、こういう場合は条件付きロジスティック回帰を行うのが方法としては最も妥当と思われる(勿論、層別化条件によってデータに変な偏りがない、という前提ではあるけど)。

posted at 22:40:11

7月3日

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

それに、接種が勧奨されていることによる偏りとして接種しないのは特段の健康的事情がある人に限られていたという仮定と同時に、逆に接種勧奨が中止された後だって、今度は接種した群がその後の健康異常に注意深くなってより「症状あり」としやすくなるという可能性も排除できない。

posted at 22:39:48

7月3日

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

ただ、気を付けなければいけないのは、このように層別化した場合に、層別でのサンプルサイズが違う、ということ。これは信頼区間の幅に影響するので、これを同列に見るのは統計的に妥当とは言えない。サイズの違うデータのp値を比較することがナンセンスなのと同じこと。

posted at 22:39:25

7月3日

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

ということで、接種群と非接種群の比率が逆転する16歳以下か17歳以上かでまずは層別に見てみることにした。すると、症状17~24で結構層別化することによってオッズ比に違いがあるっぽいということが。

posted at 22:38:41

7月3日

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

で、コレスポンデンス分析を各「症状あり」に対してやってみてプロットするとこんな感じ。明らかに各症状の発生は独立ではない。従って、データがこのような挙動を示すのはそれぞれの症状が強い共起関係にあるということでかなり説明可能と思われ。 pic.twitter.com/a2h3xP2E6l

posted at 22:37:31

7月3日

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

まず、大量にオッズ比が>1から<1になるようなものが発生する原因について、これは「偶然そうなったのか」どうかは検証しておく必要があって、それぞれの症状が独立ならこんなに<1側に傾くことは確率的にはまれそうではあるけど、そもそもこの症状、独立とは思えない性質のものが多い。

posted at 22:36:45

7月3日

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

まずは症状1~24に関して、HPVVの接種有無を説明変数として症状ありかなしかについてロジスティック回帰、次に年齢を共変量としてロジスティック回帰を行った。傾向としては消された元の結果と大体同じである様に見える。 pic.twitter.com/nGK9ayOJ5P

posted at 22:35:12

7月3日

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

ただ、この資料のp22に記載されている条件を加えたけどそれでも一致しなかったので、とりあえずは判断保留。いずれにしてもきちんとしたデータのQCが必要だし、これはPDFでデータが出ている限りは我々には難しいかと思われ。

posted at 22:32:22

7月3日

@tonkyo_Vc

地下楽師@Ph.D@tonkyo_Vc

で、素直にread.csv関数を使った。この程度のサイズなら読み込みには時間はかからなかった。症状20のクロス集計表はこんな感じになり、数値は一致。
[1] "symptom 20"
1 2
1 9107 20571
2 58 198

posted at 22:31:35

このページの先頭へ

×