{"id":3703,"date":"2025-02-13T11:39:59","date_gmt":"2025-02-13T02:39:59","guid":{"rendered":"https:\/\/plaza.umin.ac.jp\/~OIO\/?p=3703"},"modified":"2025-02-13T15:44:27","modified_gmt":"2025-02-13T06:44:27","slug":"roc-%e3%82%ab%e3%83%bc%e3%83%96%e3%82%92r%e3%81%a7%e6%8f%8f%e3%81%84%e3%81%a6%e3%81%bf%e3%81%9f","status":"publish","type":"post","link":"https:\/\/plaza.umin.ac.jp\/~OIO\/?p=3703","title":{"rendered":"ROC \u30ab\u30fc\u30d6\u3092R\u3067\u63cf\u3044\u3066\u307f\u305f"},"content":{"rendered":"\n

pROC\u30d1\u30c3\u30b1\u30fc\u30b8\u306e\u4f7f\u3044\u65b9<\/h1>\n\n\n\n

\u307e\u305a\u3001pROC<\/strong>\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u4f7f\u7528\u3057\u3066ROC\u66f2\u7dda\u3092\u4f5c\u6210\u3057\u3001AUC\u3092\u8a08\u7b97\u3059\u308b\u4f8b\u3092\u793a\u3057\u307e\u3059\u3002<\/p>\n\n\n

\n# \u30d1\u30c3\u30b1\u30fc\u30b8\u306e\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3068\u8aad\u307f\u8fbc\u307f\ninstall.packages("pROC")\nlibrary(pROC)\n\n# \u30b5\u30f3\u30d7\u30eb\u30c7\u30fc\u30bf\u306e\u8aad\u307f\u8fbc\u307f\ndata(aSAH)\n\n# \u6563\u5e03\u56f3\u306e\u4f5c\u6210 \u30c7\u30fc\u30bf\u306e\u5916\u89b3\u3092\u78ba\u8a8d\u3059\u308b\nggplot(aSAH, aes(x=wfns, y=s100b, color=outcome)) +\n  geom_point() +\n  labs(title="aSAH Data: wfns vs s100b",\n       x="wfns value",\n       y="s100b value")+\n  scale_y_continuous(limits = c(0, 1))\n\n# ROC\u66f2\u7dda\u306e\u4f5c\u6210\u3000s100b\u8840\u6e05\u5024\u3067SAH\u306e\u4e88\u5f8c\u4e88\u6e2c\u3059\u308b\u8a2d\u5b9a\nroc_obj <- roc(aSAH$outcome, aSAH$s100b)\n\n# ROC\u66f2\u7dda\u306e\u30d7\u30ed\u30c3\u30c8\nplot(roc_obj, main="ROC Curve predicting outcome of SAH using s100b value")\n\n# AUC\u306e\u8a08\u7b97\nauc_value <- auc(roc_obj)\nprint(paste("AUC:", auc_value))\n\n# \u611f\u5ea6\u3001\u7279\u7570\u5ea6\u3001\u967d\u6027\u7684\u4e2d\u7387\u3001\u9670\u6027\u7684\u4e2d\u7387\u3092\u8a08\u7b97\ncoords_obj <- coords(roc_obj, x="best", ret=c("sensitivity", "specificity", "ppv", "npv"))\n\n# \u7d50\u679c\u3092\u51fa\u529b\nprint(paste("\u611f\u5ea6:", coords_obj["sensitivity"]))\nprint(paste("\u7279\u7570\u5ea6:", coords_obj["specificity"]))\nprint(paste("\u967d\u6027\u7684\u4e2d\u7387:", coords_obj["ppv"]))\nprint(paste("\u9670\u6027\u7684\u4e2d\u7387:", coords_obj["npv"]))\n\n# \u5c24\u5ea6\u6bd4\u3092\u7b97\u51fa\u3001\u305d\u306e\u524d\u306b\u6539\u3081\u3066\u611f\u5ea6\u3068\u7279\u7570\u5ea6\u3092\u53d6\u5f97\ncoords_obj <- coords(roc_obj, x="best", ret=c("sensitivity", "specificity"))\nsensitivity <- coords_obj["sensitivity"]\nspecificity <- coords_obj["specificity"]\n\n# \u967d\u6027\u5c24\u5ea6\u6bd4 (LR+) \u3068 \u9670\u6027\u5c24\u5ea6\u6bd4 (LR-) \u3092\u8a08\u7b97\npositive_likelihood_ratio <- sensitivity \/ (1 - specificity)\nnegative_likelihood_ratio <- (1 - sensitivity) \/ specificity\n\n# \u7d50\u679c\u3092\u51fa\u529b\nprint(paste("\u967d\u6027\u5c24\u5ea6\u6bd4 (LR+):", positive_likelihood_ratio))\nprint(paste("\u9670\u6027\u5c24\u5ea6\u6bd4 (LR-):", negative_likelihood_ratio))\n\n<\/pre><\/div>\n\n\n
aSAH\u30c7\u30fc\u30bf\uff08\u4e00\u90e8\uff09<\/h3>\n\n\n\n
\"\"<\/figure>\n\n\n\n
\u30c7\u30fc\u30bf\u306e\u6982\u89b3<\/h3>\n\n\n\n
\"\"<\/figure>\n\n\n\n
plot\u51fa\u529b<\/h3>\n\n\n\n
\"\"<\/figure>\n\n\n\n
AUC\u306e\u8a08\u7b97 <\/h3>\n\n\n\n
s100b\u3067outcome\u3092\u4e88\u6e2c\u3059\u308b\u4e88\u5f8c\u4e88\u6e2c\u306e\u6027\u80fd\u3068\u3057\u3066\u306e\u8a08\u7b97<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
\u611f\u5ea6\u30fb\u7279\u7570\u5ea6\u30fb\u967d\u6027\u7684\u4e2d\u7387\u30fb\u9670\u6027\u7684\u4e2d\u7387<\/h3>\n\n\n\n
\"\"<\/figure>\n\n\n\n
\u967d\u6027\u5c24\u5ea6\u6bd4\u30fb\u9670\u6027\u5c24\u5ea6\u6bd4<\/h3>\n\n\n\n
\"\"<\/figure>\n\n\n\n
\u5c0f\u6280<\/h3>\n\n\n\n
\u3053\u306e\u30d1\u30c3\u30b1\u30fc\u30b8\u3001\u30d1\u30a4\u30d7\u3092\u4f7f\u3046\u3068\u3001\u3060\u3044\u3076\u30b7\u30f3\u30d7\u30eb\u3067\u898b\u3084\u3059\u304f\u306a\u308a\u307e\u3059
\u3053\u3053\u3067\u306f\u4e88\u6e2c\u306b\u4e0a\u8ff0\u306es100b\u306e\u4ee3\u308f\u308a\u306bwfns\u3092\u4f7f\u7528\u3057\u3066\u307f\u307e\u3057\u305f<\/p>\n\n\n
\ndata(aSAH)\n\nlibrary(dplyr) #\u30d1\u30a4\u30d7\u3092\u4f7f\u7528\u3059\u308b\u305f\u3081\u306b\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u4f7f\u7528\n\nroc2.1 <- aSAH %>% roc(outcome, wfns)\nggroc(roc2.1)\n\n<\/pre><\/div>\n\n\n
\"\"<\/figure>\n\n\n\n
R\u306epROC\u30d1\u30c3\u30b1\u30fc\u30b8\u306e\u4e3b\u306a\u95a2\u6570<\/h2>\n\n\n\n
\u95a2\u6570\u540d<\/th>\u8aac\u660e<\/th>\u4f7f\u7528\u4f8b<\/th><\/tr><\/thead>
roc()<\/code><\/td>ROC\u66f2\u7dda\u3092\u4f5c\u6210<\/td>roc(response, predictor)<\/code><\/td><\/tr>
plot()<\/code><\/td>ROC\u66f2\u7dda\u3092\u63cf\u753b<\/td>plot(roc_object)<\/code><\/td><\/tr>
auc()<\/code><\/td>ROC\u66f2\u7dda\u306eAUC\u3092\u8a08\u7b97<\/td>auc(roc_object)<\/code><\/td><\/tr>
ci()<\/code><\/td>ROC\u66f2\u7dda\u306e\u4fe1\u983c\u533a\u9593\u3092\u8a08\u7b97<\/td>ci(roc_object)<\/code><\/td><\/tr>
coords()<\/code><\/td>ROC\u66f2\u7dda\u4e0a\u306e\u7279\u5b9a\u306e\u70b9\u306e\u5ea7\u6a19\u3092\u53d6\u5f97<\/td>coords(roc_object, \"best\")<\/code><\/td><\/tr>
ggroc()<\/code><\/td>ggplot2\u3092\u4f7f\u7528\u3057\u3066ROC\u66f2\u7dda\u3092\u63cf\u753b<\/td>ggroc(roc_object)<\/code><\/td><\/tr>
roc.test()<\/code><\/td>ROC\u66f2\u7dda\u306e\u7d71\u8a08\u7684\u306a\u6bd4\u8f03<\/td>roc.test(roc_object1, roc_object2)<\/code><\/td><\/tr>
smooth()<\/code><\/td>ROC\u66f2\u7dda\u3092\u5e73\u6ed1\u5316<\/td>smooth(roc_object)<\/code><\/td><\/tr>
var()<\/code><\/td>ROC\u66f2\u7dda\u306eAUC\u306e\u5206\u6563\u3092\u8a08\u7b97<\/td>var(roc_object)<\/code><\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n
<\/p>\n\n\n\n
ROCR\u30d1\u30c3\u30b1\u30fc\u30b8\u306e\u4f7f\u3044\u65b9<\/h1>\n\n\n\n
\u6b21\u306b\u3001ROCR<\/strong>\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u4f7f\u7528\u3057\u3066ROC\u66f2\u7dda\u3092\u4f5c\u6210\u3057\u3001\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u3092\u8a55\u4fa1\u3059\u308b\u4f8b\u3092\u793a\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
\u305d\u306e\u524d\u306biris \u30c7\u30fc\u30bf\u3092\u30d7\u30ed\u30c3\u30c8\u3002\u3053\u3053\u3067\u306fPetalLength\u304c5\u4ee5\u4e0a\u3067virginica\u3068\u5224\u65ad\u3059\u308b\u4e88\u6e2c\u65b9\u6cd5\u306e\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u3092\u8abf\u3079\u308b\u3088\u3046\u306b\u3057\u307e\u3059\u3000<\/p>\n\n\n
\n# \u57fa\u672c\u30d1\u30c3\u30b1\u30fc\u30b8\u306e\u8aad\u307f\u8fbc\u307f\nlibrary(ggplot2)\n\n# iris\u30c7\u30fc\u30bf\u30bb\u30c3\u30c8\u306e\u8aad\u307f\u8fbc\u307f\ndata(iris)\n\n# \u6563\u5e03\u56f3\u306e\u4f5c\u6210\nggplot(iris, aes(x=Petal.Length, y=Petal.Width, color=Species)) +\n  geom_point() +\n  labs(title="Iris Data: Petal Length vs Petal Width",\n       x="Petal Length",\n       y="Petal Width")\n\n<\/pre><\/div>\n\n\n
\"\"<\/figure>\n\n\n
\n# \u30d1\u30c3\u30b1\u30fc\u30b8\u306e\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3068\u8aad\u307f\u8fbc\u307f\n# install.packages("ROCR")\nlibrary(ROCR)\n\n# iris\u30c7\u30fc\u30bf\u30bb\u30c3\u30c8\u306e\u8aad\u307f\u8fbc\u307f\ndata(iris)\n\n# \u4e8c\u5024\u5206\u985e\u306e\u305f\u3081\u306e\u30c7\u30fc\u30bf\u6e96\u5099\niris_binary <- iris\niris_binary$Species <- ifelse(iris$Species == "virginica", 1, 0)\n\n# \u4e88\u6e2c\u5024\u3068\u3057\u3066Petal.Length\u3092\u4f7f\u7528\u3057\u3001\u95be\u5024\uff15\u3067\u4e8c\u5024\u5316\npredictions <- ifelse(iris_binary$Petal.Length >= 5, TRUE, FALSE)\nlabels <- iris_binary$Species\n\n# Prediction\u30aa\u30d6\u30b8\u30a7\u30af\u30c8\u306e\u4f5c\u6210\npred <- prediction(as.numeric(predictions), labels)\n\n# \u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u30aa\u30d6\u30b8\u30a7\u30af\u30c8\u306e\u4f5c\u6210\nperf <- performance(pred, measure = "tpr", x.measure = "fpr")\n\n# ROC\u66f2\u7dda\u306e\u30d7\u30ed\u30c3\u30c8\nplot(perf, main="ROC Curve with Petal Width Threshold")\n\n# AUC\u306e\u8a08\u7b97\nauc_value <- performance(pred, measure = "auc")\nauc_value <- auc_value@y.values[[1]]\nprint(paste("AUC:", auc_value))\n\n<\/pre><\/div>\n\n\n
\u7d50\u679c\u306e\u30d7\u30ed\u30c3\u30c8<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
\"\"<\/figure>\n\n\n\n
\u7d50\u679c\u304c\u5947\u9e97\u3059\u304e\uff08\u30c7\u30fc\u30bf\u306e\u30d7\u30ed\u30c3\u30c8\u898b\u305f\u3089\u308f\u304b\u308b\u3088\u3046\u306b\u3001\u4e00\u3064\u306e\u5909\u6570\u3067\u7a2e\u3092\u304d\u308c\u3044\u306b\u5206\u96e2\u3067\u304d\u307e\u3059\u304b\u3089\u306d\uff09<\/p>\n","protected":false},"excerpt":{"rendered":"
pROC\u30d1\u30c3\u30b1\u30fc\u30b8\u306e\u4f7f\u3044\u65b9 \u307e\u305a\u3001pROC\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u4f7f\u7528\u3057\u3066ROC\u66f2\u7dda\u3092\u4f5c\u6210\u3057\u3001AUC\u3092\u8a08\u7b97\u3059\u308b\u4f8b\u3092\u793a\u3057\u307e\u3059\u3002 aSAH\u30c7\u30fc\u30bf\uff08\u4e00\u90e8\uff09 \u30c7\u30fc\u30bf\u306e\u6982\u89b3 plot\u51fa\u529b AUC\u306e\u8a08\u7b97 s100b\u3067outcome\u3092\u4e88\u6e2c\u3059\u308b\u4e88\u5f8c…<\/p>\n","protected":false},"author":5,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[5],"tags":[],"class_list":["post-3703","post","type-post","status-publish","format-standard","hentry","category-miscellaneous"],"jetpack_featured_media_url":"","jetpack_shortlink":"https:\/\/wp.me\/p9b6zl-XJ","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/plaza.umin.ac.jp\/~OIO\/index.php?rest_route=\/wp\/v2\/posts\/3703","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/plaza.umin.ac.jp\/~OIO\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/plaza.umin.ac.jp\/~OIO\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/plaza.umin.ac.jp\/~OIO\/index.php?rest_route=\/wp\/v2\/users\/5"}],"replies":[{"embeddable":true,"href":"https:\/\/plaza.umin.ac.jp\/~OIO\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=3703"}],"version-history":[{"count":8,"href":"https:\/\/plaza.umin.ac.jp\/~OIO\/index.php?rest_route=\/wp\/v2\/posts\/3703\/revisions"}],"predecessor-version":[{"id":3726,"href":"https:\/\/plaza.umin.ac.jp\/~OIO\/index.php?rest_route=\/wp\/v2\/posts\/3703\/revisions\/3726"}],"wp:attachment":[{"href":"https:\/\/plaza.umin.ac.jp\/~OIO\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=3703"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/plaza.umin.ac.jp\/~OIO\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=3703"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/plaza.umin.ac.jp\/~OIO\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=3703"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}