{"id":384,"date":"2021-10-24T21:36:20","date_gmt":"2021-10-24T12:36:20","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/checking-proportional-hazards-assumption-in-ezr\/"},"modified":"2024-09-29T11:02:06","modified_gmt":"2024-09-29T02:02:06","slug":"checking-proportional-hazards-assumption-in-ezr","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/checking-proportional-hazards-assumption-in-ezr\/","title":{"rendered":"EZR \u3067\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u3092\u4e8c\u91cd\u5bfe\u6570\u30d7\u30ed\u30c3\u30c8\u3067\u78ba\u8a8d\u3059\u308b\u65b9\u6cd5"},"content":{"rendered":"\n

Cox \u56de\u5e30\u3067\u5fc5\u8981\u306a\u3001\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u3092\u78ba\u8a8d\u3059\u308b\u65b9\u6cd5\u306e\u89e3\u8aac<\/p>\n\n\n\n\n\n\n\n

\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u3092\u4e8c\u91cd\u5bfe\u6570\u30d7\u30ed\u30c3\u30c8\u3067\u78ba\u8a8d\u3059\u308b\u65b9\u6cd5<\/h2>\n\n\n\n

EZR\u3067\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u3092\u78ba\u8a8d\u3059\u308b\u65b9\u6cd5\u306e\u4e00\u3064\u76ee\u3068\u3057\u3066\u3001\u30b0\u30e9\u30d5\u3067\u78ba\u8a8d\u3059\u308b\u65b9\u6cd5\u304c\u3042\u308b\u3002<\/p>\n\n\n\n

\u4e8c\u91cd\u5bfe\u6570\u30d7\u30ed\u30c3\u30c8\uff08 Complementary log-log plot \uff09\u3068\u3044\u3046\u3082\u306e\u3092\u63cf\u304f\u3002<\/p>\n\n\n\n

\u901a\u5e38\u306e\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306eR\u30b9\u30af\u30ea\u30d7\u30c8\u3092\u5fdc\u7528\u3057\u3066\u7c21\u5358\u306b\u63cf\u3051\u308b\u3002<\/p>\n\n\n\n

\u307e\u305a\u89e3\u6790\u30c7\u30fc\u30bf\u3092\u7528\u610f\u3057\u305f\u3089\u3001\u300c\u7d71\u8a08\u89e3\u6790\u300d\u2192\u300c\u751f\u5b58\u671f\u9593\u306e\u89e3\u6790\u300d\u2192\u300cLogrank\u691c\u5b9a\u300d\u3092\u9078\u629e\u3059\u308b\u3002<\/p>\n\n\n\n

\"\"
\u300c\u7d71\u8a08\u89e3\u6790\u300d\u2192\u300c\u751f\u5b58\u671f\u9593\u306e\u89e3\u6790\u300d\u2192\u300cLogrank\u691c\u5b9a\u300d\u3092\u9078\u629e<\/figcaption><\/figure>\n\n\n\n

\u89b3\u5bdf\u671f\u9593\u5909\u6570\u3001\u30a4\u30d9\u30f3\u30c8\u5909\u6570\u3001\u7fa4\u5225\u5909\u6570\u3092\u9078\u629e\u3057OK\u3092\u30af\u30ea\u30c3\u30af\u3059\u308b\u3068\u3001\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u304c\u63cf\u304b\u308c\u308b\u3002<\/p>\n\n\n\n

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

\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u306e\u4f8b<\/p>\n\n\n\n

\u3053\u306e\u3068\u304dR\u30b9\u30af\u30ea\u30d7\u30c8\u306e\u67a0\u306b\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306a\u30b9\u30af\u30ea\u30d7\u30c8\u304c\u51fa\u529b\u3055\u308c\u3066\u3044\u308b\u3002<\/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>\nxlab=<\/span>\"Days\"<\/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>\nline=<\/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>\nside=<\/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>\ntitle=<\/span>\"Anthracyclines\"<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
\u3053\u3053\u3067\u6700\u521d\u306e\uff12\u884c\u306eplot()\u306e\u6700\u5f8c\u306b\u3001fun=\"cloglog\"<\/code>\u3092\u52a0\u7b46\u3057\u3066\u3001\u5168\u4f53\u3092\u9078\u629e\u3057\u3066\u5b9f\u884c\u3092\u30af\u30ea\u30c3\u30af\u3059\u308b\u3002<\/p>\n\n\n\n
\u52a0\u3048\u3066\u3001\u4e0b\u304b\u3089\uff12\u884c\u76ee\u306elegend()\u5185\u306etopright\u3092bottomright\u306b\u5909\u66f4\u3057\u3066\u3001legend (\u51e1\u4f8b) \u304c\u30b0\u30e9\u30d5\u3068\u88ab\u3089\u306a\u3044\u3088\u3046\u306b\u3059\u308b\u3068\u3088\u3044\u3002<\/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>\nxlab=<\/span>\"Days\"<\/span>,<\/span> ylab=<\/span>\"Probability\"<\/span>,<\/span> fun=<\/span>\"cloglog\"<\/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>\nline=<\/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>\nside=<\/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>\"bottomright\"<\/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>\ntitle=<\/span>\"Anthracyclines\"<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
\u305d\u3046\u3059\u308b\u3068\u4ee5\u4e0b\u306e\u3088\u3046\u306a\u30b0\u30e9\u30d5\u304c\u63cf\u304b\u308c\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
Complementary log-log plot\u306e\u4f8b<\/p>\n\n\n\n
\uff12\u672c\u306e\u968e\u6bb5\u72b6\u306e\u30b0\u30e9\u30d5\u304c\u3001\u304a\u304a\u3088\u305d\u5e73\u884c\u306b\u306a\u3063\u3066\u3044\u308c\u3070\u3001\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u304c\u6e80\u305f\u3055\u308c\u3066\u3044\u308b\u3068\u8a00\u3048\u308b\u3002<\/p>\n\n\n\n
\u3053\u306e\u4f8b\u306f\u3001\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u304c\u6e80\u305f\u3055\u308c\u3066\u3044\u308b\u3068\u8a00\u3048\u308b\u30ec\u30d9\u30eb\u3068\u601d\u3046\u3002<\/p>\n\n\n\n
\u3053\u306e\u3088\u3046\u306b\u3057\u3066\u3001\u30b0\u30e9\u30d5\u306b\u3088\u3063\u3066\u8996\u899a\u7684\u306b\u78ba\u8a8d\u3059\u308b\u65b9\u6cd5\u304c\u4e00\u3064\u76ee\u306e\u65b9\u6cd5\u3060\u3002<\/p>\n\n\n\n
\u4e8c\u91cd\u5bfe\u6570\u306e\u610f\u5473<\/h3>\n\n\n\n
Complementary log log \u306e\u610f\u5473\u306f\u30010 \u304b\u3089 1 \u306e\u78ba\u7387\u3092\u7121\u9650\u5c0f\u304b\u3089\u7121\u9650\u5927\u306b\u5909\u63db\u3059\u308b\u3068\u3044\u3046\u3053\u3068\u3067\u3042\u308b<\/p>\n\n\n\n
\u30b0\u30e9\u30d5\u3092\u63cf\u3044\u3066\u307f\u308b\u3068\u308f\u304b\u308b<\/p>\n\n\n\n
X \u8ef8\u306b\u78ba\u7387\uff080 \u304b\u3089 1\uff09\u3001Y \u8ef8\u306b complementary log log \u3092\u65bd\u3057\u305f \u78ba\u7387\uff08\uff0d\u221e \u304b\u3089 \u221e\uff09\u3092\u53d6\u3063\u3066\u30d7\u30ed\u30c3\u30c8\u3059\u308b\u3068\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u306a\u308b<\/p>\n\n\n\n
prob <- 1:999\/1000\ncloglog <- log(-log(prob))\nplot(cloglog ~ prob, xlab=\"Probability\", ylab=\"log(-log(Probability))\")<\/code><\/pre>\n\n\n\n
\"\"<\/figure>\n\n\n\n
\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u3092\u691c\u5b9a\u3067\u78ba\u8a8d\u3059\u308b\u65b9\u6cd5<\/h2>\n\n\n\n
EZR\u3067\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u3092\u78ba\u8a8d\u3059\u308b\u65b9\u6cd5\u306e\u4e8c\u3064\u76ee\u306f\u3001\u691c\u5b9a\u3092\u7528\u3044\u308b\u65b9\u6cd5\u3060\u3002<\/p>\n\n\n\n
\u300c\u7d71\u8a08\u89e3\u6790\u300d\u2192\u300c\u751f\u5b58\u671f\u9593\u306e\u89e3\u6790\u300d\u2192\u300cCox\u56de\u5e30\u300d\u3092\u9078\u629e\u3059\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
\u300c\u7d71\u8a08\u89e3\u6790\u300d\u2192\u300c\u751f\u5b58\u671f\u9593\u89e3\u6790\u300d\u2192\u300cCox\u56de\u5e30\u300d\u3092\u9078\u629e<\/p>\n\n\n\n
\u6642\u9593\u5909\u6570\u3001\u30a4\u30d9\u30f3\u30c8\u5909\u6570\u3001\u8aac\u660e\u5909\u6570\u3092\u9078\u629e\u6295\u5165\u3059\u308b\u3002<\/p>\n\n\n\n
\u30ab\u30d7\u30e9\u30f3\u30de\u30a4\u30e4\u30fc\u66f2\u7dda\u3092\u63cf\u3044\u305f\u6642\u306b\u7528\u3044\u305f\u3001\u89b3\u5bdf\u671f\u9593\u5909\u6570\u3001\u30a4\u30d9\u30f3\u30c8\u5909\u6570\u3001\u7fa4\u5225\u5909\u6570\u3068\u540c\u3058\u5909\u6570\u3092\u9078\u629e\u3059\u308b\u3002<\/p>\n\n\n\n
\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u306e\u5206\u6790\u3092\u884c\u3046\u306b\u30c1\u30a7\u30c3\u30af\u3092\u5165\u308c\u3066\u3001OK\u3092\u30af\u30ea\u30c3\u30af\u3059\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u306e\u5206\u6790\u3092\u884c\u3046\u306b\u30c1\u30a7\u30c3\u30af<\/p>\n\n\n\n
cox.zph() \u3068\u3044\u3046\u95a2\u6570\u3067\u3001\u691c\u5b9a\u3055\u308c\u308b\u3002<\/p>\n\n\n\n
\u5e30\u7121\u4eee\u8aac\u306f\u3001\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u304c\u6e80\u305f\u3055\u308c\u3066\u3044\u308b\u3001\u3068\u3044\u3046\u3082\u306e\u3002<\/p>\n\n\n\n
\u7d71\u8a08\u5b66\u7684\u6709\u610f\u3067\u306a\u3051\u308c\u3070\u3001\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u304c\u6e80\u305f\u3055\u308c\u3066\u3044\u306a\u3044\u3068\u306f\u8a00\u3048\u306a\u3044\u3001\u3064\u307e\u308a\u6d88\u6975\u7684\u306b\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u304c\u6e80\u305f\u3055\u308c\u3066\u3044\u308b\u3068\u8a00\u3063\u3066\u3082\u3044\u3044\u3068\u3044\u3046\u3053\u3068\u306b\u306a\u308b\u3002<\/p>\n\n\n\n
\u7d50\u679c\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u8868\u793a\u3055\u308c\u308b\u3002<\/p>\n\n\n\n
><\/span> print<\/span>(<\/span>cox.zph<\/span>(<\/span>CoxModel.1))<\/span>\nchisq df    p\nAnthracyclines 0.609<\/span>  1<\/span> 0.44<\/span>\nGLOBAL         0.609<\/span>  1<\/span> 0.44<\/span>\n<\/code><\/pre>\n\n\n\n
\u691c\u5b9a\u306e\u672c\u6765\u306e\u76ee\u7684\u306f\u3001\u5e30\u7121\u4eee\u8aac\u304c\u68c4\u5374\u3055\u308c\u308b\u3053\u3068\u306a\u306e\u306b\u3001\u5e30\u7121\u4eee\u8aac\u304c\u68c4\u5374\u3055\u308c\u306a\u3044\u3053\u3068\u304c\u76ee\u7684\u3067\u3042\u308b\u691c\u5b9a\u30b7\u30ea\u30fc\u30ba\uff08\u6b63\u898f\u6027\u306e\u691c\u5b9a\u3001\u7b49\u5206\u6563\u6027\u306e\u691c\u5b9a\u3001\u7403\u9762\u6027\u306e\u691c\u5b9a\u7b49\uff09\u306e\u4e00\u3064\u3067\u3001\u304a\u3059\u3059\u3081\u306e\u65b9\u6cd5\u3067\u306f\u306a\u3044<\/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>
\u307e\u3068\u3081<\/h2>\n\n\n\n
EZR\u3067\u751f\u5b58\u6642\u9593\u89e3\u6790\u3092\u5b9f\u65bd\u3059\u308b\u969b\u306b\u691c\u8a0e\u3057\u305f\u3044\u6bd4\u4f8b\u30cf\u30b6\u30fc\u30c9\u6027\u306e\u78ba\u8a8d\u65b9\u6cd5\u3092\uff12\u3064\u7d39\u4ecb\u3057\u305f\u3002<\/p>\n\n\n\n
\u500b\u4eba\u7684\u306b\u306f complementary log-log plot \u3067\u8996\u899a\u7684\u306b\u78ba\u8a8d\u3059\u308c\u3070\u5341\u5206\u3068\u601d\u3063\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
\u691c\u5b9a\u306f\u3001\u30b5\u30f3\u30d7\u30eb\u30b5\u30a4\u30ba\u304c\u5927\u304d\u3044\u5834\u5408\u3001\u307b\u3093\u306e\u5c11\u3057\u306e\u305a\u308c\u3067\u3082\u691c\u51fa\u3057\u3066\u3057\u307e\u3044\u3001\u9069\u5207\u3067\u306f\u306a\u3044\u304b\u3089\u3060\u3002<\/p>\n\n\n\n
\u6b63\u898f\u6027\u306e\u78ba\u8a8d\u3068\u540c\u69d8\u306b\u3001\u3069\u3093\u306a\u3068\u304d\u3082\u307e\u305a\u306f\u30b0\u30e9\u30d5\u3067\u8996\u899a\u7684\u306b\u78ba\u8a8d\u3059\u308b\u3053\u3068\u304c\u6700\u512a\u5148\u3068\u3044\u3046\u4e00\u4f8b\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\u4f55\u3089\u304b\u53c2\u8003\u306b\u306a\u308c\u3070\u3002<\/p>\n\n\n\n
\u304a\u3059\u3059\u3081\u66f8\u7c4d<\/h2>\n\n\n\n
\"EZR\u3067\u3084\u3055\u3057\u304f\u5b66\u3076\u7d71\u8a08\u5b66<\/a><\/figure>\n\n\n\n
<\/div>\n\n\n\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