{"id":251,"date":"2022-09-15T22:28:48","date_gmt":"2022-09-15T13:28:48","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/how-to-draw-iptw-km-curve\/"},"modified":"2024-09-01T20:01:57","modified_gmt":"2024-09-01T11:01:57","slug":"how-to-draw-iptw-km-curve","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/how-to-draw-iptw-km-curve\/","title":{"rendered":"R \u3068 EZR \u3067 IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3092\u66f8\u304f\u65b9\u6cd5"},"content":{"rendered":"\n

IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306e\u66f8\u304d\u65b9\u3002<\/p>\n\n\n\n

R\u306e\u30d1\u30c3\u30b1\u30fc\u30b8\u3068EZR\u3092\u4f7f\u3046\u3068\u66f8\u3051\u308b\u3002<\/p>\n\n\n\n\n\n\n\n

IPTW\u3068\u306f\uff1f<\/h2>\n\n\n\n

IPTW \u3068\u306f\u3001Inverse Probability of Treatment Weights \u306e\u982d\u6587\u5b57\u8a9e\u3002<\/p>\n\n\n\n

\u9006\u78ba\u7387\u91cd\u307f\u4ed8\u3051\u306e\u610f\u5473\u3002<\/p>\n\n\n\n

\u50be\u5411\u30b9\u30b3\u30a2 (Propensity Score) \u3092\u4f7f\u3063\u305f\u3001\u4ea4\u7d61\u8981\u56e0\u8abf\u6574\u306e\u4e00\u7a2e\u3002<\/p>\n\n\n\n

\u51e6\u65b9\u3055\u308c\u306b\u304f\u3044\u5206\u753b\u306b\u3001\u51e6\u65b9\u610f\u5411\u78ba\u7387 (Propensity Score\u3068\u547c\u3070\u308c\u308b) \u306e\u9006\u6570\u3067\u91cd\u307f\u3092\u4ed8\u3051\u3066\u3001\u51e6\u65b9\u3055\u308c\u3084\u3059\u3044\u5206\u753b\u3068\u30d0\u30e9\u30f3\u30b9\u3092\u3068\u3063\u3066\u3001\u6bd4\u8f03\u53ef\u80fd\u6027\u3092\u4e0a\u3052\u308b\u65b9\u6cd5\u3002<\/p>\n\n\n\n

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

IPTW\u3067\u30d0\u30e9\u30f3\u30b9\u3092\u3063\u305f\u4e0a\u3067\u306e\u7fa4\u9593\u6bd4\u8f03\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306e\u3053\u3068\u3002<\/p>\n\n\n\n

\u3053\u308c\u3092\u66f8\u304f\u306b\u306f\u3069\u3046\u3059\u308c\u3070\u3088\u3044\u304b\uff1f<\/p>\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>

R\u3067IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3092\u66f8\u304f\u65b9\u6cd5<\/h2>\n\n\n\n

RISCA\u30d1\u30c3\u30b1\u30fc\u30b8\u306e\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb<\/h3>\n\n\n\n

RISCA\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3057\u3066\u3001ipw.survival() \u304c\u4f7f\u3048\u308b\u3088\u3046\u306b\u3059\u308b\u3002<\/p>\n\n\n\n

install.packages<\/span>(<\/span>\"RISCA\"<\/span>)<\/span>\nlibrary<\/span>(<\/span>RISCA)<\/span>\n<\/code><\/pre>\n\n\n\n
\u3044\u304f\u3064\u304b\u30d1\u30c3\u30b1\u30fc\u30b8\u306e\u66f4\u65b0\u304c\u5fc5\u8981\u306b\u306a\u308b\u5834\u5408\u304c\u3042\u308b\u306e\u3067\u3001\u6839\u6c17\u3088\u304f\u4e00\u3064\u305a\u3064\u66f4\u65b0\u3059\u308b\u3002<\/p>\n\n\n\n
R\u30b3\u30f3\u30bd\u30fc\u30eb\u30e1\u30cb\u30e5\u30fc\u306e\u300c\u30d1\u30c3\u30b1\u30fc\u30b8\u300d\u2192\u300c\u30d1\u30c3\u30b1\u30fc\u30b8\u306e\u66f4\u65b0\u300d\u3092\u4f7f\u3046\u3068\u66f4\u65b0\u3067\u304d\u308b\u3002<\/p>\n\n\n\n
\u4ee5\u4e0b\u306e\u30ea\u30f3\u30af\u5148\u306eR\u30b9\u30af\u30ea\u30d7\u30c8\u304c\u5fc5\u8981\u306a ipw.survivial() \u95a2\u6570\u306a\u306e\u3067\u3001\u3053\u308c\u3092\u30b3\u30f3\u30bd\u30fc\u30eb\u306b\u30b3\u30d4\u30da\u3057\u3066\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u306e\u4ee3\u308f\u308a\u306b\u3057\u3066\u3082\u3088\u3044\u3002<\/p>\n\n\n\n
EZR\u3067IPTW\u3092\u4f5c\u6210\u3059\u308b<\/h3>\n\n\n\n
EZR\u306e\u30ed\u30b8\u30b9\u30c6\u30a3\u30c3\u30af\u56de\u5e30\u5206\u6790\u306e\u30e1\u30cb\u30e5\u30fc\u3092\u501f\u308a\u3066\u3001IPTW\u306e\u5909\u6570\u3092\u4f5c\u6210\u3059\u308b\u3002<\/p>\n\n\n\n
\u51e6\u65b9\u306e\u6709\u7121\uff08\u80cc\u666f\u3092\u306a\u308b\u3079\u304f\u63c3\u3048\u305f\u3044\u56e0\u5b50\uff09\u3092\u76ee\u7684\u5909\u6570\u306b\u3057\u3066\u3001\u305d\u306e\u51e6\u65b9\u306b\u5f71\u97ff\u3092\u4e0e\u3048\u308b\u56e0\u5b50\u3067\u8aac\u660e\u3059\u308b\u30ed\u30b8\u30b9\u30c6\u30a3\u30c3\u30af\u56de\u5e30\u30e2\u30c7\u30eb\u3092\u4f5c\u6210\u3057\u3001IPTW\u3092\u8a08\u7b97\u3059\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u8d64\u67a0\u306e\u30c1\u30a7\u30c3\u30af\u304cIPTW\u4f5c\u6210\u306e\u305f\u3081\u306e\u6307\u5b9a\u7b87\u6240\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
OK\u3092\u62bc\u3057\u3066\u8a08\u7b97\u3059\u308b\u3068\u3001weight.ATE.GLM.1\u3068\u3044\u3046\u5909\u6570\u304c\u4f5c\u6210\u3055\u308c\u308b\u3002<\/p>\n\n\n\n
\u3053\u308c\u304cIPTW\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u3053\u3053\u307e\u3067\u6e96\u5099\u3057\u3066\u3001ipw.survival() \u95a2\u6570\u306b\u623b\u308b\u3002<\/p>\n\n\n\n
ipw.survival\u95a2\u6570\u3067IPTW \u751f\u5b58\u78ba\u7387\u3092\u8a08\u7b97\u3059\u308b<\/h3>\n\n\n\n
R\u30b3\u30f3\u30bd\u30fc\u30eb\u3067\u3082EZR\u306eR\u30b9\u30af\u30ea\u30d7\u30c8\u7a93\u3067\u3082\u3088\u3044\u306e\u3067\u3001\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u66f8\u3044\u3066\u5b9f\u884c\u3059\u308b\u3068\u3001IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u751f\u5b58\u66f2\u7dda\u306e\u305f\u3081\u306e\u6570\u5024\u304c\u8a08\u7b97\u3067\u304d\u308b\u3002<\/p>\n\n\n\n
\u3053\u306e\u95a2\u6570\u3067\u306f IPTW \u3092 IPW \u3068\u8868\u73fe\u3057\u3066\u3044\u308b\u306e\u3067\u6ce8\u610f\u3002<\/p>\n\n\n\n
IPTW \u3068 IPW \u306f\u540c\u3058\u3082\u306e\u3092\u6307\u3057\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
res <-<\/span> with<\/span>(<\/span>ElderlyAML,<\/span> ipw.survival<\/span>(<\/span>times =<\/span> Days,<\/span>\nfailures =<\/span> Survival,<\/span> variable =<\/span> Anthracyclines,<\/span>\nweights =<\/span> weight.ATE.GLM.1))<\/span>\n<\/code><\/pre>\n\n\n\n
\u6b20\u640d\u304c\u3042\u308b\u3068\u8a08\u7b97\u3055\u308c\u306a\u3044\u306e\u3067\u3001\u3042\u3089\u304b\u3058\u3081\u6b20\u640d\u304c\u3042\u308b\u30b1\u30fc\u30b9\u306f\u524a\u9664\u3059\u308b\u51e6\u7406\u3092\u884c\u3063\u3066\u304a\u304f\u3002<\/p>\n\n\n\n
\u89e3\u6790\u7d50\u679c\u3092\u5c11\u3057\u898b\u3066\u307f\u308b\u3068\u3001\u4ee5\u4e0b\u306e\u3088\u3046\u306a\u69cb\u9020\u306b\u306a\u3063\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
times \u306f\u751f\u5b58\u6642\u9593\u3001n.risk \u306f number at risk\u3001n.event \u306f number of events\u3001survival \u306f IPTW \u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u751f\u5b58\u78ba\u7387\u3001variable \u306f\u7fa4\u5206\u3051\u5909\u6570\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
n.risk \u3068 n.event \u304c\u5c0f\u6570\u307e\u3067\u3042\u308a\u3001IPTW\u304c\u304b\u304b\u3063\u3066\u3044\u308b\u3053\u3068\u304c\u308f\u304b\u308b\u3002<\/p>\n\n\n\n
\u3042\u3068\u306f\u3053\u306e\u30c7\u30fc\u30bf\u3092\u30b0\u30e9\u30d5\u306b\u66f8\u304f\u306e\u307f\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
IPTW\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3092\u66f8\u304f<\/h3>\n\n\n\n
\u307e\u305a\u3001\u7a7a\u767d\u306e\u30b0\u30e9\u30d5\u7bc4\u56f2\u3092\u4f5c\u6210\u3059\u308b\u3002<\/p>\n\n\n\n
plot<\/span>(<\/span>NULL<\/span>,<\/span> xlim=<\/span>c<\/span>(<\/span>0<\/span>,<\/span>2100<\/span>),<\/span> ylim=<\/span>c<\/span>(<\/span>0<\/span>,<\/span>1<\/span>),<\/span> ylab=<\/span>\"Survival\"<\/span>,<\/span>\nxlab=<\/span>\"Time\"<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
X\u8ef8\u306e\u7bc4\u56f2\u306f\u3042\u3089\u304b\u3058\u3081\u8abf\u3079\u3066\u304a\u304f\u3068\u3088\u3044\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u4e00\u672c\u76ee\u306f\u51e6\u65b9\u3042\u308a variable == 1 \u306e\u7dda\u3092\u66f8\u304f\u3002<\/p>\n\n\n\n
lines<\/span>(<\/span>x =<\/span> res$<\/span>table.surv$<\/span>times[<\/span>res$<\/span>table.surv$<\/span>variable==<\/span>1<\/span>],<\/span>\ny =<\/span> res$<\/span>table.surv$<\/span>survival[<\/span>res$<\/span>table.surv$<\/span>variable==<\/span>1<\/span>],<\/span>\ntype =<\/span> \"s\"<\/span>,<\/span> col=<\/span>1<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u4e8c\u672c\u76ee\u306f\u3001\u51e6\u65b9\u306a\u3057\u3067\u751f\u5b58\u78ba\u7387\u304c\u901f\u304f\u4e0b\u304c\u3063\u3066\u3057\u307e\u3046 variable == 0 \u3092\u66f8\u304f\u3002<\/p>\n\n\n\n
lines<\/span>(<\/span>x =<\/span> res$<\/span>table.surv$<\/span>times[<\/span>res$<\/span>table.surv$<\/span>variable==<\/span>0<\/span>],<\/span>\ny =<\/span> res$<\/span>table.surv$<\/span>survival[<\/span>res$<\/span>table.surv$<\/span>variable==<\/span>0<\/span>],<\/span>\ntype =<\/span> \"s\"<\/span>,<\/span> col=<\/span>2<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u6700\u5f8c\u306b\u3001\u53f3\u4e0a\u306e\u307b\u3046\u306b Legend \u3092\u633f\u5165\u3057\u3066\u5b8c\u6210\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
legend<\/span>(<\/span>\"topright\"<\/span>,<\/span> legend=<\/span>c<\/span>(<\/span>\"Anthracyclines Treated\"<\/span>,<\/span>\n\"Anthracyclines Untreated\"<\/span>),<\/span> col=<\/span>c<\/span>(<\/span>1<\/span>,<\/span>2<\/span>),<\/span> lty=<\/span>c<\/span>(<\/span>1<\/span>,<\/span>1<\/span>))<\/span>\n<\/code><\/pre>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
IPTW\u3067\u91cd\u307f\u4ed8\u3051\u3092\u3057\u3066\u3044\u306a\u3044\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3092\u66f8\u3044\u3066\u307f\u308b\u3068\u4ee5\u4e0b\u306e\u3068\u304a\u308a\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u9ed2\u306e\u7dda\u306f\u3042\u307e\u308a\u5909\u308f\u3063\u3066\u3044\u306a\u3044\u304c\u3001\u8d64\u306e\u7dda\u304c\u82e5\u5e72\u5909\u308f\u3063\u3066\u3044\u308b\u306e\u304c\u308f\u304b\u308b\u3002<\/p>\n\n\n\n
\u751f\u5b58\u78ba\u73870.4\u3042\u305f\u308a\u306e\u30a4\u30d9\u30f3\u30c8\u306e\u8d77\u304d\u65b9\u304c\u7570\u306a\u3063\u3066\u3044\u308b\u3088\u3046\u306b\u898b\u3048\u308b\u3002<\/p>\n\n\n\n
\u3053\u306e\u3088\u3046\u306bIPTW\u3067\u91cd\u307f\u4ed8\u3051\u3092\u3057\u305f\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306f\u3001\u91cd\u307f\u4ed8\u3051\u306a\u3057\u306e\u66f2\u7dda\u3068\u306f\u82e5\u5e72\u7570\u306a\u308b\u66f2\u7dda\u306b\u306a\u308b\u308f\u3051\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\u5b9f\u306f\u3001\u3082\u3063\u3068\u7c21\u5358\u306b\u66f8\u3051\u308b\u3002<\/p>\n\n\n\n
\u3042\u3068\u306b\u306a\u3063\u3066\u6c17\u3065\u3044\u305f\u3002<\/p>\n\n\n\n
### plot.survrisca<\/span>\nplot<\/span>(<\/span>res,<\/span> col=<\/span>2<\/span>:<\/span>1<\/span>,<\/span>\nxlab=<\/span>'Time'<\/span>,<\/span> ylab=<\/span>'Survival'<\/span>)<\/span>\nlegend<\/span>(<\/span>\"topright\"<\/span>,<\/span> legend=<\/span>c<\/span>(<\/span>\"Anthracyclines Treated\"<\/span>,<\/span>\n\"Anthracyclines Untreated\"<\/span>),<\/span>\ncol=<\/span>c<\/span>(<\/span>1<\/span>,<\/span>2<\/span>),<\/span> lty=<\/span>c<\/span>(<\/span>1<\/span>,<\/span>1<\/span>))<\/span>\n<\/code><\/pre>\n\n\n\n
\u3053\u306e\u3088\u3046\u306b\u30b9\u30af\u30ea\u30d7\u30c8\u3092\u66f8\u304f\u3060\u3051\u3067\u3001\u4ee5\u4e0b\u306e\u30b0\u30e9\u30d5\u304c\u63cf\u3051\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u6253\u3061\u5207\u308a\u30de\u30fc\u30af\u3092\u5165\u308c\u305f IPTW KM \u66f2\u7dda\u306e\u66f8\u304d\u65b9<\/h3>\n\n\n\n
\u6253\u3061\u5207\u308a\u304c\u767a\u751f\u3057\u305f\u3068\u3053\u308d\u306b\u7e26\u7dda\u3092\u4ed8\u3051\u305f\u3044\u5834\u5408\u306f\u3001\u9055\u3046\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u4f7f\u3046\u3068\u3067\u304d\u308b\u3002<\/p>\n\n\n\n
adjustedCurves \u3068\u3044\u3046\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u4f7f\u3046\u3002<\/p>\n\n\n\n
\u95a2\u6570\u306f\u3001adjustedsurv() \u3067\u3042\u308b\u3002<\/p>\n\n\n\n
library<\/span>(<\/span>adjustedCurves)<\/span>\nlibrary<\/span>(<\/span>survival)<\/span>\nlibrary<\/span>(<\/span>ggplot2)<\/span>\n### \u7fa4\u5225\u306e\u5909\u6570\u3092\u56e0\u5b50\u306b<\/span>\nElderlyAML$<\/span>Anthracyclines.Factor <-<\/span>\nas.factor<\/span>(<\/span>ElderlyAML$<\/span>Anthracyclines)<\/span>\n### \u50be\u5411\u30b9\u30b3\u30a2\u8a08\u7b97\u306e\u305f\u3081\u306e\u30ed\u30b8\u30b9\u30c6\u30a3\u30c3\u30af\u56de\u5e30<\/span>\nGLM.1 <-<\/span> glm<\/span>(<\/span>Anthracyclines ~<\/span>\nSex +<\/span> Age +<\/span> PS +<\/span> WBC +<\/span> BM_blast,<\/span>\nfamily=<\/span>binomial<\/span>(<\/span>logit),<\/span> data=<\/span>ElderlyAML)<\/span>\n### iptw_km \u3092\u6307\u5b9a\u3057\u3066\u3001\u50be\u5411\u30b9\u30b3\u30a2\u8a08\u7b97\u30e2\u30c7\u30eb\u3082\u6307\u5b9a\u3059\u308b<\/span>\nadjsurv <-<\/span> adjustedsurv<\/span>(<\/span>data=<\/span>ElderlyAML,<\/span>\nvariable=<\/span>'Anthracyclines.Factor'<\/span>,<\/span>\nev_time=<\/span>'Days'<\/span>,<\/span> event=<\/span>'Survival'<\/span>,<\/span>\nmethod=<\/span>'iptw_km'<\/span>,<\/span>\ntreatment_model=<\/span>GLM.1)<\/span>\n### censoring_ind='lines' \u3067\u6253\u3061\u5207\u308a\u4f8b\u767a\u751f\u6642\u70b9\u306b\u7dda\u304c\u5f15\u3051\u308b<\/span>\nplot<\/span>(<\/span>adjsurv,<\/span> censoring_ind=<\/span>'lines'<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
\u4ee5\u4e0b\u306e\u3088\u3046\u306a\u30b0\u30e9\u30d5\u306b\u306a\u308b<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
EZR \u306a\u3069\u3067\u5225\u9014\u4f5c\u6210\u3057\u305f IPTW \u3092\u4f7f\u7528\u3059\u308b\u3053\u3068\u3082\u3067\u304d\u308b<\/p>\n\n\n\n
weight.ATE.GLM.1 \u3068\u3044\u3046\u5909\u6570\u540d\u306e IPTW \u3092 treatment_model \u3067\u6307\u5b9a\u3059\u308b<\/p>\n\n\n\n
### iptw_km \u3092\u6307\u5b9a\u3057\u3066\u3001\u8a08\u7b97\u6e08\u307f\u306e IPTW \u3092\u6307\u5b9a\u3059\u308b<\/span>\nadjsurv <-<\/span> adjustedsurv<\/span>(<\/span>data=<\/span>ElderlyAML,<\/span>\nvariable=<\/span>'Anthracyclines.Factor'<\/span>,<\/span>\nev_time=<\/span>'Days'<\/span>,<\/span> event=<\/span>'Survival'<\/span>,<\/span>\nmethod=<\/span>'iptw_km'<\/span>,<\/span>\ntreatment_model=<\/span>ElderlyAML$<\/span>weight.ATE.GLM.1)<\/span>\n### censoring_ind='lines' \u3067\u6253\u3061\u5207\u308a\u4f8b\u767a\u751f\u6642\u70b9\u306b\u7dda\u304c\u5f15\u3051\u308b<\/span>\nplot<\/span>(<\/span>adjsurv,<\/span> censoring_ind=<\/span>'lines'<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
\u305d\u3046\u3059\u308b\u3068\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u30b0\u30e9\u30d5\u304c\u66f8\u3051\u308b<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
IPTW \u7d2f\u7a4d\u767a\u751f\u7387\u66f2\u7dda\u306e\u66f8\u304d\u65b9<\/h2>\n\n\n\n
\u7d2f\u7a4d\u767a\u751f\u7387\u306f\u3001\u751f\u5b58\u7387\u3092 1 \u304b\u3089\u5f15\u3051\u3070\u3088\u3044\u306e\u3067\u3001y = 1 – res$ … \u3068\u3059\u308c\u3070\u3088\u3044\u3002<\/p>\n\n\n\n
plot<\/span>(<\/span>NULL<\/span>,<\/span> xlim=<\/span>c<\/span>(<\/span>0<\/span>,<\/span>2100<\/span>),<\/span> ylim=<\/span>c<\/span>(<\/span>0<\/span>,<\/span>1<\/span>),<\/span>\nylab=<\/span>\"Cumulative Incidence\"<\/span>,<\/span> xlab=<\/span>\"Time\"<\/span>)<\/span>\nlines<\/span>(<\/span>x =<\/span> res$<\/span>table.surv$<\/span>times[<\/span>res$<\/span>table.surv$<\/span>variable==<\/span>1<\/span>],<\/span>\ny =<\/span> 1<\/span>-<\/span>res$<\/span>table.surv$<\/span>survival[<\/span>res$<\/span>table.surv$<\/span>variable==<\/span>1<\/span>],<\/span>\ntype =<\/span> \"s\"<\/span>,<\/span> col=<\/span>1<\/span>)<\/span>\nlines<\/span>(<\/span>x =<\/span> res$<\/span>table.surv$<\/span>times[<\/span>res$<\/span>table.surv$<\/span>variable==<\/span>0<\/span>],<\/span>\ny =<\/span> 1<\/span>-<\/span>res$<\/span>table.surv$<\/span>survival[<\/span>res$<\/span>table.surv$<\/span>variable==<\/span>0<\/span>],<\/span>\ntype =<\/span> \"s\"<\/span>,<\/span> col=<\/span>2<\/span>)<\/span>\nlegend<\/span>(<\/span>\"bottomright\"<\/span>,<\/span> legend=<\/span>c<\/span>(<\/span>\"Anthracyclines Treated\"<\/span>,<\/span>\n\"Anthracyclines Untreated\"<\/span>),<\/span> col=<\/span>c<\/span>(<\/span>1<\/span>,<\/span>2<\/span>),<\/span> lty=<\/span>c<\/span>(<\/span>1<\/span>,<\/span>1<\/span>))<\/span>\n<\/code><\/pre>\n\n\n\n
\u3059\u308b\u3068\u3001\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u751f\u5b58\u7387\u3092\u9006\u3055\u306b\u3057\u305f\u30b0\u30e9\u30d5\u304c\u63cf\u3051\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u901a\u5e38\u306e\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3082\u3001plot(km, …) \u3068\u3042\u308b\u30b9\u30af\u30ea\u30d7\u30c8\u306b\u3001fun = ‘event’ \u3068\u52a0\u7b46\u3059\u308b\u3060\u3051\u3067\u3001\u7d2f\u7a4d\u767a\u751f\u7387\u306e\u66f2\u7dda\u304c\u66f8\u3051\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u307e\u3068\u3081<\/h2>\n\n\n\n
\u7d71\u8a08\u30bd\u30d5\u30c8R\u3067\u3001RISCA \u30d1\u30c3\u30b1\u30fc\u30b8\u306e ipw.survival() \u95a2\u6570\u3092\u7528\u3044\u305f\u3001IPTW\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3092\u66f8\u3044\u3066\u307f\u305f\u3002<\/p>\n\n\n\n
EZR\u306e\u30ed\u30b8\u30b9\u30c6\u30a3\u30c3\u30af\u56de\u5e30\u5206\u6790\u306e\u30e1\u30cb\u30e5\u30fc\u3092\u501f\u308a\u308b\u3068\u7c21\u5358\u306bIPTW\u304c\u8a08\u7b97\u3067\u304d\u308b\u306e\u3067\u3001EZR\u306e\u529b\u3092\u501f\u308a\u308b\u306e\u304c\u304a\u3059\u3059\u3081\u3002<\/p>\n\n\n\n
\u3042\u3068\u306f\u3001\u5c11\u3057\u3060\u3051R\u30b9\u30af\u30ea\u30d7\u30c8\u3092\u66f8\u304b\u306a\u3044\u3068\u3044\u3051\u306a\u3044\u304c\u3001\u4e0a\u8a18\u3092\u53c2\u8003\u306b\u3001\u5909\u6570\u540d\u4ee5\u5916\u306f\u30b3\u30d4\u30da\u3067OK\u306a\u306e\u3067\u3001\u305c\u3072\u8a66\u3057\u3066\u307f\u3066\u307b\u3057\u3044\u3002<\/p>\n\n\n\n
\u53c2\u8003\u30b5\u30a4\u30c8<\/h2>\n\n\n\n
ipw.survival: Adjusted Survival Curves by Using IPW. in RISCA: Causal Inference and Prediction in Cohort-Based Analyses<\/a><\/p>\n\n\n\n
IPWsurvival \u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u7528\u3044\u305f\u7d71\u8a08\u7684\u56e0\u679c\u63a8\u8ad6 x \u751f\u5b58\u5206\u6790 – YuRAN-HIKO<\/a><\/p>\n\n\n\n
\u304a\u3059\u3059\u3081\u66f8\u7c4d<\/h2>\n\n\n\n
EZR\u516c\u5f0f\u30de\u30cb\u30e5\u30a2\u30eb<\/p>\n\n\n\n
\"EZR\u3067\u3084\u3055\u3057\u304f\u5b66\u3076\u7d71\u8a08\u5b66<\/a>\n
\n
EZR\u3067\u3084\u3055\u3057\u304f\u5b66\u3076\u7d71\u8a08\u5b66 \u6539\u8a023\u7248 \u301cEBM\u306e\u5b9f\u8df5\u304b\u3089\u81e8\u5e8a\u7814\u7a76\u307e\u3067\u301c<\/a><\/p>\n
RISCA source: R\/ipw.survival.R<\/a><\/p>\n\n\n\n