{"id":14,"date":"2024-07-27T11:55:41","date_gmt":"2024-07-27T02:55:41","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/how-to-calculate-median-survial-time-and-confidence-interval-of-iptw-km-curve-in-ezr\/"},"modified":"2024-08-25T20:21:39","modified_gmt":"2024-08-25T11:21:39","slug":"how-to-calculate-median-survial-time-and-confidence-interval-of-iptw-km-curve-in-ezr","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/how-to-calculate-median-survial-time-and-confidence-interval-of-iptw-km-curve-in-ezr\/","title":{"rendered":"EZR \u3067 IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc \u66f2\u7dda\u306e\u751f\u5b58\u671f\u9593\u4e2d\u592e\u5024\u3068 95 \uff05 \u4fe1\u983c\u533a\u9593\u3092\u8a08\u7b97\u3059\u308b\u65b9\u6cd5"},"content":{"rendered":"\n

EZR \u3067 IPTW\uff08\u9006\u78ba\u7387\u91cd\u307f\u4ed8\u3051\uff09\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306e\u751f\u5b58\u671f\u9593\u4e2d\u592e\u5024\u3068 95 \uff05 \u4fe1\u983c\u533a\u9593\u3092\u8a08\u7b97\u3059\u308b\u65b9\u6cd5\u3092\u7d39\u4ecb\u3059\u308b<\/p>\n\n\n\n\n\n\n\n

<\/div>\n\n\n\n\n\n\n\n

IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3068\u306f\u4f55\u304b\uff1f<\/h2>\n\n\n\n

IPTW\uff08\u9006\u78ba\u7387\u91cd\u307f\u4ed8\u3051\uff09 \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3068\u306f\u3001\u50be\u5411\u30b9\u30b3\u30a2\u306e\u9006\u6570\u3092\u7528\u3044\u305f IPTW \u3067\u4ea4\u7d61\u56e0\u5b50\u3092\u8abf\u6574\u3057\u305f\u3001\u7fa4\u9593\u6bd4\u8f03\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306e\u3053\u3068\u3067\u3042\u308b<\/p>\n\n\n\n

\u66f8\u304d\u65b9\u306e\u5177\u4f53\u7684\u306a\u65b9\u6cd5\u306f\u3001\u904e\u53bb\u8a18\u4e8b\u53c2\u7167<\/p>\n\n\n

\t\t\t
\n\t\t\t\t
\n\t\t\t\t\t\u3042\u308f\u305b\u3066\u8aad\u307f\u305f\u3044<\/span>\n\t\t\t\t\t
\"\"<\/figure><\/div>\t\t\t\t\t
\n\t\t\t\t\t\tR \u3067 IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u30b0\u30e9\u30d5\u306b Number at Risk \u3092\u66f8\u304d\u5165\u308c\u308b\u65b9\u6cd5<\/a>\n\t\t\t\t\t\tIPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306b\u304a\u3044\u3066\u3001\u4efb\u610f\u306e\u6642\u70b9\u306b\u304a\u3051\u308b Number at Risk \u3092\u66f8\u304d\u5165\u308c\u308b\u65b9\u6cd5 R \u3067\u4f5c\u6210\u3059\u308b\u65b9\u6cd5 IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda IPTW \u306f\u3001\u9006\u78ba\u7387\u91cd\u307f\u4ed8\u3051\u3068\u3082…<\/span>\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t<\/div>\n\t\t<\/div>\n\n\n\n\n\n
\t\t\t
\n\t\t\t\t
\n\t\t\t\t\t\u3042\u308f\u305b\u3066\u8aad\u307f\u305f\u3044<\/span>\n\t\t\t\t\t
\"\"<\/figure><\/div>\t\t\t\t\t
\n\t\t\t\t\t\tR \u3068 EZR \u3067 IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3092\u66f8\u304f\u65b9\u6cd5<\/a>\n\t\t\t\t\t\tIPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306e\u66f8\u304d\u65b9\u3002 R\u306e\u30d1\u30c3\u30b1\u30fc\u30b8\u3068EZR\u3092\u4f7f\u3046\u3068\u66f8\u3051\u308b\u3002 IPTW\u3068\u306f\uff1f IPTW \u3068\u306f\u3001Inverse Probability of Treatment Weights \u306e\u982d\u6587\u5b57\u8a9e\u3002 \u9006\u78ba\u7387\u91cd\u307f…<\/span>\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t<\/div>\n\t\t<\/div>\n\n\n\n\n\n\n

IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306e\u751f\u5b58\u671f\u9593\u4e2d\u592e\u5024\u3068 95 \uff05\u4fe1\u983c\u533a\u9593\u3092\u8a08\u7b97\u3059\u308b\u65b9\u6cd5<\/h2>\n\n\n\n

EZR \u3067\u3001IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306e\u751f\u5b58\u671f\u9593\u4e2d\u592e\u5024\u3068 95 \uff05\u4fe1\u983c\u533a\u9593\u3092\u8a08\u7b97\u3059\u308b\u65b9\u6cd5\u3092\u89e3\u8aac\u3059\u308b<\/p>\n\n\n\n

\u65b9\u91dd\u3068\u3057\u3066\u306f\u3001\u4e00\u90e8 EZR \u306e\u7d71\u8a08\u89e3\u6790\u30e1\u30cb\u30e5\u30fc\u3067\u51fa\u529b\u3055\u308c\u308b R \u30b9\u30af\u30ea\u30d7\u30c8\u3092\u501f\u308a\u308b<\/p>\n\n\n\n

\u4e00\u90e8\u52a0\u7b46\u3059\u308b\u3053\u3068\u3067\u8a08\u7b97\u3059\u308b<\/p>\n\n\n\n

\u307e\u305a\u3001\u901a\u5e38\u306e\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3092\u66f8\u304f<\/p>\n\n\n\n

\u7d71\u8a08\u89e3\u6790 \u2192 \u751f\u5b58\u671f\u9593\u306e\u89e3\u6790 \u2192 \u751f\u5b58\u66f2\u7dda\u306e\u8a18\u8ff0\u3068\u7fa4\u9593\u306e\u6bd4\u8f03\uff08Logrank \u691c\u5b9a\uff09\u3092\u9078\u629e\u3059\u308b<\/p>\n\n\n\n

\u6642\u9593\u5909\u6570\u3001\u30a4\u30d9\u30f3\u30c8\u5909\u6570\u3001\u7fa4\u5225\u5909\u6570\u3092\u9078\u629e\u3057\u3066\u3001OK \u3092\u30af\u30ea\u30c3\u30af\u3059\u308b\u3068\u3001\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u66f8\u3051\u308b<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

<\/span><\/p>\n\n\n\n

\u540c\u6642\u306b\u3001\u751f\u5b58\u671f\u9593\u4e2d\u592e\u5024\u3068 95\uff05 \u4fe1\u983c\u533a\u9593\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u8a08\u7b97\u3055\u308c\u308b<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

<\/span><\/p>\n\n\n\n\n\n\n\n

IPTW \u3092\u6307\u5b9a\u3057\u3066\u518d\u5ea6\u5b9f\u884c\u3059\u308b<\/h3>\n\n\n\n

IPTW \u3092 survfit \u306b\u66f8\u304d\u5165\u308c\u3066\u3001\u518d\u5ea6\u5b9f\u884c\u3059\u308b<\/p>\n\n\n\n

\u5177\u4f53\u7684\u306b\u306f\u3001weights \u3092\u8ffd\u8a18\u3059\u308b\uff08\u524d\u63d0\u3068\u3057\u3066\u3001IPTW \u304c\u8a08\u7b97\u3055\u308c\u3066\u3044\u3066\u3001weight.ATE.GLM.1 \u3068\u3044\u3046\u540d\u524d\u3067\u3042\u308b\uff09<\/p>\n\n\n\n

R \u30b9\u30af\u30ea\u30d7\u30c8\u67a0\u5185\u306b\u4ee5\u4e0b\u306e\u30b9\u30af\u30ea\u30d7\u30c8\u3092\u898b\u3064\u3051\u3066\u3001\u6700\u5f8c\u306e\u90e8\u5206\u306b\u3001weights = weight.ATE.GLM.1\uff08IPTW \u306e\u5909\u6570\u540d\u306b\u9069\u5b9c\u5909\u66f4\uff09\u52a0\u7b46\u3059\u308b<\/p>\n\n\n\n

km <-<\/span> survfit<\/span>(<\/span>Surv<\/span>((<\/span>time\/<\/span>1<\/span>),<\/span>status01==<\/span>1<\/span>)<\/span>~<\/span>ph.ecog01,<\/span> data=<\/span>lung_rev,<\/span> na.action =<\/span> na.omit,<\/span>\nconf.type=<\/span>\"log-log\"<\/span>,<\/span> weights=<\/span>weight.ATE.GLM.1)<\/span>\n<\/pre>\n\n\n\n
\u4e0a\u8a18\u3092\u9078\u629e\u3057\u3066\u5b9f\u884c\u5f8c\u3001\u4ee5\u4e0b\u306e\u90e8\u5206\u3092\u9078\u629e\u3057\u518d\u5ea6\u5b9f\u884c\u3059\u308b\u3068\u3001\u305d\u306e\u4e0b\u306e\u3088\u3046\u306a\u30b0\u30e9\u30d5\u304c\u66f8\u3051\u308b<\/p>\n\n\n\n
plot<\/span>(<\/span>km,<\/span> bty=<\/span>\"l\"<\/span>,<\/span> col=<\/span>1<\/span>:<\/span>32<\/span>,<\/span> lty=<\/span>1<\/span>,<\/span> lwd=<\/span>1<\/span>,<\/span> conf.int=<\/span>FALSE<\/span>,<\/span> mark.time=<\/span>TRUE<\/span>,<\/span> xlab=<\/span>\"time\"<\/span>,<\/span> ylab=<\/span>\"Probability\"<\/span>)<\/span>\nxticks <-<\/span> axTicks<\/span>(<\/span>1<\/span>)<\/span>\nn.atrisk <-<\/span> nrisk<\/span>(<\/span>km,<\/span> xticks)<\/span>\nfor<\/span> (<\/span>i in<\/span> 1<\/span>:<\/span>length<\/span>(<\/span>km$<\/span>strata)){<\/span>axis<\/span>(<\/span>1<\/span>,<\/span> at =<\/span> xticks,<\/span> labels =<\/span> n.atrisk[<\/span>i,],<\/span> line=<\/span>3<\/span>+<\/span>i,<\/span> tick =<\/span> FALSE<\/span>)}<\/span>\nfor<\/span> (<\/span>i in<\/span> 1<\/span>:<\/span>length<\/span>(<\/span>km$<\/span>strata)){<\/span>mtext<\/span>(<\/span>legend[<\/span>i],<\/span> at=-<\/span>(<\/span>xticks[<\/span>2<\/span>]<\/span>-<\/span>xticks[<\/span>1<\/span>])<\/span>\/<\/span>2<\/span>,<\/span> side=<\/span>1<\/span>,<\/span> line=<\/span>4<\/span>+<\/span>i,<\/span> cex=<\/span>1<\/span>)}<\/span>\ntitle<\/span>(<\/span>xlab =<\/span> \"Number at risk\"<\/span>,<\/span> line =<\/span> 3.5<\/span>,<\/span> adj =<\/span> 0<\/span>)<\/span>\nlegend <\/span>(<\/span>\"topright\"<\/span>,<\/span> legend,<\/span> col=<\/span>1<\/span>:<\/span>32<\/span>,<\/span> lty=<\/span>1<\/span>,<\/span> lwd=<\/span>1<\/span>,<\/span>  box.lty=<\/span>0<\/span>,<\/span> title=<\/span>\"ph.ecog01\"<\/span>)<\/span>\n<\/pre>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
Number at risk \u304c\u6574\u6570\u3067\u306f\u306a\u304f\u3001\u5c0f\u6570\u70b9\u4ee5\u4e0b\u304c\u3042\u308b\u305f\u3081\u3001\u75c7\u4f8b\u4e00\u4eba\u4e00\u4eba\u306b\u91cd\u307f\u304c\u304b\u304b\u3063\u305f\u3053\u3068\u304c\u308f\u304b\u308b<\/p>\n\n\n\n
\u751f\u5b58\u671f\u9593\u4e2d\u592e\u5024\u3068 95 \uff05\u4fe1\u983c\u533a\u9593\u306f\u3001\u4ee5\u4e0b\u306e\u7b87\u6240\u3092\u9078\u629e\u3057\u5b9f\u884c\u3059\u308b\u3068\u3001\u8a08\u7b97\u3055\u308c\u308b<\/p>\n\n\n\n
km.summary.table <-<\/span> summary.km<\/span>(<\/span>survfit=<\/span>km,<\/span> survdiff=<\/span>res)<\/span>\nkm.summary.table\n<\/pre>\n\n\n\n
\u751f\u5b58\u671f\u9593\u4e2d\u592e\u5024\u3068 95 \uff05\u4fe1\u983c\u533a\u9593\u306f\u3001\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u5c11\u3057\u3060\u3051\u5909\u308f\u3063\u305f<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u3061\u306a\u307f\u306b\u3001P \u5024\uff08\u30ed\u30b0\u30e9\u30f3\u30af\u691c\u5b9a\uff09\u306f\u3001IPTW \u524d\u306e\u307e\u307e\u306a\u306e\u3067\u3001\u4f7f\u308f\u306a\u3044\u3088\u3046\u306b\u6ce8\u610f<\/p>\n\n\n\n
IPTW \u30ed\u30b0\u30e9\u30f3\u30af\u691c\u5b9a\u306f\u3001\u3053\u3061\u3089\u3092\u53c2\u7167\u306e\u3053\u3068<\/p>\n\n\n
\t\t\t
\n\t\t\t\t
\n\t\t\t\t\t\u3042\u308f\u305b\u3066\u8aad\u307f\u305f\u3044<\/span>\n\t\t\t\t\t
\"\"<\/figure><\/div>\t\t\t\t\t
\n\t\t\t\t\t\tEZR \u3068 R \u3067 IPTW \u30ed\u30b0\u30e9\u30f3\u30af\u691c\u5b9a\u3092\u5b9f\u884c\u3059\u308b\u65b9\u6cd5<\/a>\n\t\t\t\t\t\tIPTW \u30ed\u30b0\u30e9\u30f3\u30af\u691c\u5b9a\u3092EZR\u3068R\u3092\u4f7f\u3063\u3066\u884c\u3046\u65b9\u6cd5\u3002 IPTW\u3068\u306f\uff1f IPTW \u3068\u306f\u3001Inverse Probability of Treatment Weights \u306e\u982d\u6587\u5b57\u8a9e\u3002 \u65e5\u672c\u8a9e\u3067\u306f\u3001\u9006\u78ba\u7387\u91cd\u307f\u3065\u3051\u3068\u8a00\u3046\u3002 …<\/span>\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t<\/div>\n\t\t<\/div>\n\n\n\n\n\n\n
\uff1e\uff1e\u3082\u3046\u7d71\u8a08\u3067\u60a9\u3080\u306e\u306f\u7d42\u308f\u308a\u306b\u3057\u307e\u305b\u3093\u304b\uff1f\u00a0<\/a><\/strong><\/span><\/p>\r\n\"\"<\/a>\r\n
\u21911\u4e07\u4eba\u4ee5\u4e0a\u306e\u533b\u7642\u5f93\u4e8b\u8005\u304c\u8cfc\u8aad\u4e2d<\/span><\/strong><\/span><\/p><\/div>
\u307e\u3068\u3081<\/h2>\n\n\n\n
EZR \u3067 IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306e\u751f\u5b58\u671f\u9593\u4e2d\u592e\u5024\u3068 95 \uff05\u4fe1\u983c\u533a\u9593\u3092\u8a08\u7b97\u3059\u308b\u65b9\u6cd5\u3092\u89e3\u8aac\u3057\u305f<\/p>\n\n\n\n
\u6700\u521d\u306b\u3001\u901a\u5e38\u306e\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3092\u66f8\u304d\u3001\u305d\u306e\u6642\u306b\u751f\u6210\u3055\u308c\u305f R \u30b9\u30af\u30ea\u30d7\u30c8\u306b IPTW \u3092\u6307\u5b9a\u3059\u308b\u3088\u3046\u306b\u52a0\u7b46\u3059\u308b\u3060\u3051\u3067\u3001\u7c21\u5358\u306b\u8a08\u7b97\u3067\u304d\u308b<\/p>\n\n\n\n
survfit \u3068\u3044\u3046\u95a2\u6570\u306e\u5f15\u6570\u3068\u3057\u3066\u3001weights \u3092\u6307\u5b9a\u3059\u308b\u3060\u3051<\/p>\n\n\n\n
\u53c2\u8003\u306b\u306a\u308c\u3070<\/p>\n","protected":false},"excerpt":{"rendered":"
EZR \u3067 IPTW\uff08\u9006\u78ba\u7387\u91cd\u307f\u4ed8\u3051\uff09\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306e\u751f\u5b58\u671f\u9593\u4e2d\u592e\u5024\u3068 95 \uff05 \u4fe1\u983c\u533a\u9593\u3092\u8a08\u7b97\u3059\u308b\u65b9\u6cd5\u3092\u7d39\u4ecb\u3059\u308b<\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"swell_btn_cv_data":"","_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_publicize_message":"","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":false,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2}},"categories":[4,7,8,9],"tags":[],"class_list":["post-14","post","type-post","status-publish","format-standard","hentry","category-ezr","category-iptw-","category-8","category-9"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/14","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/comments?post=14"}],"version-history":[{"count":2,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/14\/revisions"}],"predecessor-version":[{"id":708,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/14\/revisions\/708"}],"wp:attachment":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/media?parent=14"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/categories?post=14"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/tags?post=14"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}