ウルヴァリン WOLVERINE ADDISON BOOT TWO TONE WINGTIP W05947 BLACK × BROWNウルバリン アディソン ブーツブラック × ブラウン ウィングチップ1000マイル ツートーン:SHOWA CL-50 レバーハンドル錠 HL15A型 NXシリンダー キー3本付属 ドアノブ 交換 取替え【ショウワ CLシリーズ】

BR-ケブラーパーカー ライトブルー(女性用):mobus(モーブス) MSW6010J/MSW6011P ユニクロスジャケット・パンツ上下セット サッカー フットサル トレーニングウェア ピンク

LACOSTE ラコステ スニーカー FAIRLEAD FRX WT/DB (MAB022-QW4)【bp】:ブーツ メンズフェイクスウェード サイド ジップ ブーツ(靴 メンズ靴 ブーツ エンジニア メンズimproves imp) インプローブス 秋物 メンズ 夏物 夏服 夏男

ELPA(エルパ) 増設用ワイヤレス防犯カメラ CMS-C71 1818700:木製ウエディングビッグスプーン 60cm【名入れなし】

防犯カメラ 監視カメラ 4台 録画 セット 屋外高画質 1TB 家庭用防犯カメラ スターターセット 【防犯カメラ】【録画機】【防水】:【国内正規品】NEW BALANCE ニューバランス M1400 BKS black ブラック メンズ レディース スニーカー 【ST】

VM10S-B/3 太陽電池式LED回転灯 Ф105 3点留め バッテリー式 ニコソーラー 日恵製作所 G17071 goleador シグナルスリング >

« 読了:Zeeelenberg & Pieters (2004) 顧客不満の裏にある感情の効果 | メイン | 読了:Hu, et al. (2001) 固定電話対象のRDD調査と携帯電話対象のRDD調査を両方やってひとつにまとめました »

ウォーキング シューズ YONEX ( ヨネックス ) ( MC30 ) パワークッション メンズ ウォーキング シューズ 送料無料 予約受付中 【送料無料100215 smtb-k w1 sybp w1】:svolme【スボルメ】ジャージロングパンツ〈フットサル・サッカー・ジャージ〉121-58550

パトライト LED流動表示灯■型式:KLE-24S-R■定格電圧:AC/DC24V■色:R(赤)■電流:5.8W■閃光数:50~200回/分■質量:1.2kg■ボディ仕様:ステンレスボディ:Tenvis ネットワークカメラ ベビーモニター 100万画素数 日本語対応 動体検知機能付き 遠隔操作【iPhone 6 6Plus/スマホ/Wi-Fi対応】◇ALW-IPROBOT3

クリスマスツリー スレンダーツリー150cm:goleador リュック [ドットバックパック] ゴレアドール

防犯カメラ SDカード録画 屋外 監視カメラ セット 防水 小型 暗視 【 SDカード 録画 完全防水】赤外線 録画 防水 リモコン付属 【ロングセラーモデル 初心者でも扱いが簡単】 05P06Aug16:品番:PB4526≪2014春夏≫【ペナルティ】ベンチレーションパック≪PENALTY2014SS≫

交換用シリンダー(鍵) FBロック MIWA LA・MA用 シルバー【防犯対策】【取替用】:【即納 】ラブライブ!スクフェス SR高坂穂乃果 コスプレ衣装 覚醒後 忍者

VM10S-B/3 太陽電池式LED回転灯 Ф105 3点留め バッテリー式 ニコソーラー 日恵製作所 G17071 goleador シグナルスリング

2015年6月 8日 (月)

 仕事の都合で、左側切断と右側打ち切りがある生存時間データについて考えていた。恥ずかしながら、両方ある奴を扱うのははじめてだ。わかりやすく書いてしまうといろいろ差し障りがあるので書けないけど、マーケティング・データ解析ではそんなに珍しい問題状況じゃないかもしれない、と反省。
 意外なことに、この状況についてきちんと説明した参考書がなかなか見当たらない。社会科学の文脈で生存分析を扱った名著Singer & Willett(2003, 最近訳書が出た)には切断の話が出てこないように思う。Rのsurvivalパッケージの中の人Therneauさんによる教科書Therneau & Grambsch (2010) にもほとんど出てこない(これはラッキーともいえる。私にはこの本はかなりわかりにくい)。いまのところ見つけたのはクライン&メシュベルガー(原著2003)なんだけど、先生、文系の私にはこれでさえハードルが高いんです...

Cain, K.C., Harlow, S.D., Little, R.J., Nan, B., Yosef, M., Taffe, J.R., & Elliott, J.R. (2011) Bias Due to Left Truncation and Left Censoring in Longitudinal Studies of Developmental and Disease Processes.

日本製 本革 コインローファー [SARABANDE サラバンド]8608 [送料無料]ローファー カジュアル ビジネス 革靴 天然皮革 アイビー ローファー メンズ スリッポン【RCP】02P06Aug16:【ポイント15倍】スピングルムーブ SPINGLE MOVE スニーカー SPM-101 メンズ レディース (ホワイト)(120810)

リーボック Reebok インスタ ポンプフューリー スニーカー INSTAPUMP FURY TECH V63499 メンズ レディース 靴 ネイビー:通販/正規品 おすすめ ブーツ エミュー エミュ モカシンシューズ レディース ファー ボア Moccasin Moonah ムーナ モカシン 最終!EMU

200円OFFクーポン発行中!MIZUNO ミズノ フットサルシューズ ウエーブグレヴィスタPRO3 12KF310:ニュースター ドアクローザードアチェックP182 P-182 パラレル段付アーム段付ブラケットタイプシルバー色 ストップ付☆☆NS NEWSTAR ニュースター☆P182AK☆NS NEWSTAR ニュースター☆☆ 玄関 ドア 扉 修理 交換 部品 パーツ

[送料無料/交換再送料無料]RegettaCanoe/ライトフット/ストライプクロスベルトサンダル/CJLF2116/日本製/リゲッタ:【公式】adidas アディダス 2016-17 マンチェスターユナイテッドFC プレゼンテーション パンツ メンズ AP3861【mail0806】

よさこい衣装 よさこいコスチューム:DalPonte(ダウポンチ) ピステ上下セットアップ(裏地メッシュ仕様)(NVY)【フットサルウェア/ピステ/ピステパンツ/上下セット】【150926】【DalPonte15fw】【evi】

American Journal of Epidemiology, 173(9), 1078-1084.
 というわけで、出張の新幹線でメモを取りながら読んだ。

 ええと...
 縦断研究で、関心がある最初のマイルストーン(発症とか)から最後のマイルストーン(死亡とか)までを観察しているのをインシデント・コホート・デザインという。これに対して、なんらかのマイルストーンを観察できなかった人を除外しちゃうのをプリバレント・コホート・デザインという。後者のデザインにはバイアス源が3つあることになる:

生存分析による右側打ち切りの対処は広く活用されているが、左側打ち切り・左側切断はそうでもない。
 左側切断によるバイアスは、切断に関わっているマイルストーンがひとつならまだわかりやすい。ややこしいのは複数のマイルストーンと関連しているときだ。たとえば、女性の更年期についての3つのコホート研究を比較すると、後期更年期の開始時点の分布が全然違う。これは研究によっては参加時にもう後期更年期が始まっている人がいたり(左側打ち切り)、更年期が終わっちゃった人が研究に参加できなかったり(その後のマイルストーンによる左側切断)するからだ。[←ああ、そうか。この例ならすぐ気が付くだろうけど、こうやって対象者条件が複雑になればなるほど、思わぬ落とし穴も増えるなあ...]

 終末イベント時の年齢をY、中間イベント時の年齢をXとする(たとえば、更年期の終了と後期更年期の開始)。D=Y-Xとする。Y, X, Dの分布を推定するというのが目標である。
 研究参加時の年齢をA, 退出時の年齢をBとする[←ここではすべてをカレンダー時間ではなく年齢で定義している点に注意。いやあ、途中で混乱したぜ。Figure 1をぼーっと見ていると誤解するぞ]

景品 二次会 ディズニーランドorシーペア パスポート チケット 景品3点セット 毛がに ハーゲンダッツ&ブランブリュン ( パネル付き!) 【送料無料・】 【 二次会 景品 目録 ビンゴ ゴルフコンペ 】:ニューバランス M996NAV NEWBALANCE M996 NAV D:width ネイビー NAVY メンズ スニーカー  『送料無料』

adidas(アディダス) エックス15.1 CG AF4806 コアBLK/クリスタルW 29.0cm:屋外用ダミーカメラ 【RD-3947】 ダミー防犯カメラ/ダミー監視カメラ/ダミー カメラ

★送料無料★コレガ(corega) PoE 暗視 屋外 防塵 防水 耐衝撃 対応 Full HD バンダルドーム型 監視・防犯 有線ネットワークカメラ(05P06Aug16) CG-NCVD031A:ソリッドカメラ 屋外用防犯カメラ(バリフォーカル対応) ANC-201WPIRV

ディズニー コスチューム 大人 男性 女性 用 ドナルド モコモコ 仮装 在庫限り:【スペシャルプライス】[ティンバーランド] Timberland TIMB A118U メンズ カジュアルシューズ レースアップシューズ シンプル スタイリッシュ オシャレ ブランド 人気 イエロー

YKK アルミサッシ 引戸錠 戸先 KH-65 キー3本付属 引き戸 鍵 交換 取替え【ウエスト WEST】【KH-65 KH65】【引戸 引違戸】:レイコップふとん専用ダニクリーナー(ホワイト) 景品3点セットB(景品 二次会 コンペ 新年会 忘年会 結婚式二次会 ゴルフ イベント 参加賞 抽選会 ビンゴ 粗品 賞品 選べる)

。可能な順序は6つある: AXYB、AXBY、ABXY、XABY、XAYB, AYAB。
 ここでは次のデザインに焦点を当てる。

 

最高品質VM10S-B/3 太陽電池式LED回転灯 Ф105 3点留め バッテリー式 ニコソーラー 日恵製作所 G17071 goleador シグナルスリング

デザインIについて。[←混乱するので自分のために例を挙げておくと、分析対象期間前に退会した会員はデータから抜け落ちているけど、現会員についての入会日はわかる、会員の在籍期間(入会日を原点にとってX)の平均を推定せよ、といわれているような状況だ。ありそうありそう]
 Xの左側切断を無視してYなりXなりDなりの分布を推定したとして、仮にAが固定されていたなら、その推定は X > Aという条件のもとでの条件付き分布となる。実際にはAはたいてい変動するけど。
 このデザインでは、Aを組み込んでXの周辺分布を不偏推定するノンパラなproduct-limit推定量が存在する[←なんのことかと思ったが、カプラン・マイヤー推定量のことをproduct-limit推定量というらしい。へー]。Stata, SAS, Rに載っている。ただし、Aの最小値がXの最小値より早いことが必要。さらに、Aの最小値とXの最小値が近いと不安定になる(リスク集合が小さすぎるから)。こういうときは、Aの最小値の右側にA0という点を定義し、X < A0 の下での条件付き分布を推定するという手もある [←へぇー!]。もちろん周辺分布に対しては不偏でなくなる。
 不幸にして、Xの多くがA (ないしA0)より左だったとしよう。Xの周辺分布を不偏推定するためには、もはやXのパラメトリックな分布を知るよりほかに手がない。それがわかれば、尤度を Pr(X < A)で割ればいい[←理解するまで1分ほど固まりました... いやぁ、専門家の奴らは頭がいい]。

正規代理店 日本国内送料・代引き手数料無料 REDWING (レッドウィング) 9875 6inch CLASSIC MOC TOE ブーツ ゴールドラセットセコイア 犬タグ【10P06Aug16】:★送料無料★ マジカルスパイク#900 オレンジ スパイクピン付作業長靴

よさこいポリエステル長半纏 =よさこい衣装=:【H-GANG】デジタルドアロック PRISMA M701◆電子鍵 デジタルキー 保安キー セキュリティー 【生活用品】【デジタル雑貨】

【ディアドラ】 DD-NA 中綿入り フィールドジャケット・パンツ セット(ブラック)FW5150 9998 99 サッカー◆W368 ウインドブレーカー 上下セット:USJペアチケットをメインにふく刺し・ステーキと参加者大盛り上がり!景品10点セット 【ユニバーサル・スタジオ・ジャパン】【景品目録】二次会・結婚式・ゴルフコンペ・イベント(ゴルフ コンペ 2次会 ペアギフト)

【全品送料無料】[スボルメ]SVOLME モビライトスタンドジップ 153-59601【フットサル サッカー】:8ホールレースアップブーツ/メンズ ブーツ レースアップブーツ エンジニア スタンダード レザー 本革 ライダース プレーン ロック ミリタリー レザーブーツ パンク モード

景品 二次会 【当日出荷OK&送料無料】 全てパネル&目録!ルンバと豪華人気景品5点セット ビンゴ 忘年会 新年会 パーティー・イベント用品・販促品 景品:UMBRO(アンブロ) アクセレイター KL Dブルー×ホワイト×ネイビー USS7602BWN F 28.0cm

。StataとRのサンプルコードを示すから読むがよい。
 ところで、すべての潜在的参加者について、Aの前にXがあったかどうかだけはわかっているとしたらどうか。この場合は、Xで左側切断されているのではなく、Xで左側打ち切りしていることになる。これもコードを示すので読め。
 さて、ここまでの話はXが正しくわかる場合の話だ。実際には往々にしてXはあいまいである。更年期の例だと、生理周期が60日を超えた初回を後期更年期の開始と定義するんだけど、ある参加者の生理周期が「はじめて」60日を超えたと参加者が報告したとして、実は観察期間の前にもそれは生じていたかもしれない。以下、参加者の報告を信じる方法をナイーブ法と呼ぶ。

 デザインIIについて。おさらいすると、イベントYで切断、Xで打ち切りが生じてます。[←分析対象期間前に退会した会員はデータから抜け落ちており、入会日(カレンダー時間でX)については新規会員についてはわかるが既存会員についてはもうわからん、会員の在籍期間(退会日をカレンダー時間でYとしてY-X)の平均を推定せよ、でも新規会員のみの平均じゃだめよといわれているようなものか。悪夢ではあるがありそうな話でもある]
 このときは、まずXとYの同時分布のモデルが必要になる。XとYが独立でないかぎり、Yの切断がXの分布の推定にバイアスをもたらすからだ。二変量正規分布を仮定したモデルのコード例を示すから読め。
 なお Jiang et al.(2005 Biometrics)は、XとYの関係についてはパラメトリックに記述しそれぞれの周辺分布についてはノンパラに推定するというモデルを提案している。[←なんでそんな変なことを考えたのかと思ったら、Xは糖尿病性網膜症、Yは糖尿病による死亡だそうだ。なるほど、観察開始前の死亡者はデータから抜け落ちる、網膜症を発症したかどうかはわかるが時期は観察期間でない限りわからない、発症からの平均余命を推定したい、発症年齢の確率分布はわからないけど発症から死亡までの期間についてはなんらかの仮定ができそうだ]

 シミュレーション。 
 TREMINというコホート研究をつかう。参加者は女性で、ずっと生理を記録している。40歳の時点で参加していた人の40歳以降の記録を分析対象とする。後期更年期の開始年齢をX, 閉経年齢をYとする。XとYの周辺分布は正規分布に従っているが、Dは歪んでいる(定義上 X > Yにはならないから当然だ)。でも二変量正規性を仮定して分析するぞ。詳細は付録をみよ。[←すんません、いまちょっとその時間がないっす]
 TREMINには左側打ち切りも左側切断もない。このデータで推定したX, Y, Dの分布を、このデータから以下の手順で抜き出したデータで推定した結果と比較する。(1)ある期間(たとえば40歳~50歳)に参加した人のデータだけ切りだす。(2)ブートストラップ抽出。(3)個々の対象者に一様乱数Aを与え、Aから10年分を切り出す。(4)X < Aの対象者を除外(デザインI)。ないし、Y < Aの対象者を除外(デザインII)。
 結果[きちんと読んでないけど...]。左側切断・左打ち切りをちゃんとモデル化しないとバイアスを受ける。正しく扱っても切断されている割合が50%あたりを超えると不安定になる。云々。

 考察。関心あるイベントが切断を受けるイベントである場合は簡単で、RやSASやStataで扱える。しかし関心があるイベントと切断のあるイベントがずれている場合はややこしくて、なんらかのパラメトリックな仮定が必要になる。
 実際の切断はそんなにシャープじゃなくて、むしろ確率的に生じることも多い。その調整方法は今後の課題。[←ヘックマン・モデルみたいな話か...]
 研究デザインを変えて切断をやり過ごせることもある。たとえば自然流産までのイベント時間分析だったら、すべての潜在的参加者について最初の評価の際に過去の自然流産の有無を調べておくとよい。このクロスセクショナルなデータがあれば、縦断データのほうがY > Aな対象者だけであっても、左側切断が左側打ち切りに変わる。

 正直なところ本文よりappendixのコード例のほうに惹かれて読んだんだけど、いやー、これはホントに勉強になった。ありがたや、ありがたや。

NSC-AHD921AHD SeriesAHDミニカメラボードレンズ1280×720解像度AHDカメラマイクも搭載!音声付きで記録を残せます:スピングルムーブ スニーカー SPM-110 ダークブルー SPINGLE MOVE SPM-110 SPM110 DARK BLUEスピングルムーヴ

サッカー フットサルゴール 2台セット スチール製 (固定用ペグ、ハンマー、専用バッグ付き) フィールドフォース 送料無料 サッカーゴール セール SALE:アシックス フットサル専用シューズ 屋内用TST216 DESTAQUE 6 J デスタッキ フラッシュイエロー×オニキス(0799)

★送料無料★ 【フォーマル靴】 紳士靴  本革メンズシューズ 【フォーマルにOK】2E(EE狭め) オシャレなエナメル素材!(メンズ シューズ 白) (メンズ シューズ ホワイ) 紳士靴 小さいサイズ24.0~27.5まで!YN318-1 【円高還元YDKG-kd】:2/19再入荷 【 こだわりの made in JAPAN 】 CONVERSE ONE STAR コンバース ワンスター J V-3 ベルクロ レザー 限定 白 黒

penetrar【ペネトラール】ジャージセットアップ〈フットサル サッカー〉243-06509:オプテックス LEDセンサライト シルバー LA-11LED S [LA11LEDS]【05P06Aug16】

優勝カップ【文字彫刻無料】ガラス台「桜色」のカップ GA221-Bサイズ●高さ195mm:キーン(KEEN) Yogui Arts ヨギ アーツ メンズ Graphite 1002036 【アウトドア スポーツサンダル】【SNDL15】

論文:データ解析 - 読了:Cain, et al. (2011) 左側打ち切りと左側切断が生存分析に与えるバイアス

{yahoojp} {elsur.jpn.org} nt07-myp01-gui-304