{"id":40,"date":"2018-07-24T22:41:33","date_gmt":"2018-07-24T13:41:33","guid":{"rendered":"http:\/\/jazzman.pe.kr\/blog\/?p=40"},"modified":"2026-02-13T22:15:21","modified_gmt":"2026-02-13T13:15:21","slug":"r-pooled-cohort-equation-10-year-absolute-risk","status":"publish","type":"post","link":"https:\/\/jazzman.pe.kr\/blog\/index.php\/2018\/07\/24\/40\/","title":{"rendered":"[R] Pooled cohort equation (10-year absolute risk)"},"content":{"rendered":"\n<p>Pooled cohort equation \uc73c\ub85c 10\ub144 \uac04\uc758 \uc2ec\ub1cc\ud608\uad00\uc9c8\ud658\uc758 \ubc1c\ubcd1 \uc704\ud5d8\ub3c4\ub97c \uc608\uce21\ud558\ub294 \uacf5\uc2dd\uc744 R function \uc73c\ub85c \uad6c\ud604\ud558\uc600\uc2b5\ub2c8\ub2e4.<\/p>\n\n\n\n<div class=\"wp-block-kevinbatdorf-code-block-pro\" data-code-block-pro-font-family=\"Code-Pro-JetBrains-Mono\" style=\"font-size:clamp(14px, .875rem, 21px);font-family:Code-Pro-JetBrains-Mono,ui-monospace,SFMono-Regular,Menlo,Monaco,Consolas,monospace;line-height:clamp(20px, 1.25rem, 30px);--cbp-tab-width:2;tab-size:var(--cbp-tab-width, 2)\"><span style=\"display:block;padding:16px 0 0 16px;margin-bottom:-1px;width:100%;text-align:left;background-color:#2e3440ff\"><svg xmlns=\"http:\/\/www.w3.org\/2000\/svg\" width=\"54\" height=\"14\" viewBox=\"0 0 54 14\"><g fill=\"none\" fill-rule=\"evenodd\" transform=\"translate(1 1)\"><circle cx=\"6\" cy=\"6\" r=\"6\" fill=\"#FF5F56\" stroke=\"#E0443E\" stroke-width=\".5\"><\/circle><circle cx=\"26\" cy=\"6\" r=\"6\" fill=\"#FFBD2E\" stroke=\"#DEA123\" stroke-width=\".5\"><\/circle><circle cx=\"46\" cy=\"6\" r=\"6\" fill=\"#27C93F\" stroke=\"#1AAB29\" stroke-width=\".5\"><\/circle><\/g><\/svg><\/span><span role=\"button\" tabindex=\"0\" style=\"color:#d8dee9ff;display:none\" aria-label=\"Copy\" class=\"code-block-pro-copy-button\"><pre class=\"code-block-pro-copy-button-pre\" aria-hidden=\"true\"><textarea class=\"code-block-pro-copy-button-textarea\" tabindex=\"-1\" aria-hidden=\"true\" readonly>## Pooled cohort equation: implementing pooled cohort equation\n## AR10 (ten-year absolute risk) = 1 - S^exp(XB - meanXB)\n##     S = survival rate at 10 yrs\n## women and men, white\n## gender M:1, F:2\n## bpl = BP lowering treatment, untreated or treated: 0 or 1\n## Reference:\n## 2013 ACC\/AHA Guideline on the Assessment of Cardiovascular Risk\n## A Report of the American College of Cardiology\/American Heart\n## Association Task Force on Practice Guidelines\n## Circulation. 2014;129&#91;suppl 2&#93;:S49-S73\n\npce &lt;- function(gender, age, tc, hdl.c, sbp, bpl, cur.smoker, dm) {\n    b.surv &lt;- c(0.9144, 0.9665)  ## baseline survival in US cohort, for white men &amp; women\n    names(b.surv) &lt;- c(\"male\", \"female\")\n\n    coef.matrix &lt;- array(\n        dim = c(2, 11, 2),\n        dimnames = list(\n            c(\"bp.untreated\", \"bp.treated\"),\n            c(),\n            c(\"male\", \"female\")\n            )\n        )\n    coef.matrix&#91;\"bp.untreated\", , \"male\"&#93; &lt;-\n        c(12.344, 0, 11.853, -2.664, -7.990, 1.769, 1.764, 0, 7.837, -1.795, 0.658)\n    coef.matrix&#91;\"bp.treated\", , \"male\"&#93; &lt;-\n        c(12.344, 0, 11.853, -2.664, -7.990, 1.769, 1.797, 0, 7.837, -1.795, 0.658)\n    coef.matrix&#91;\"bp.untreated\", , \"female\"&#93; &lt;-\n        c(-29.799, 4.884, 13.54, -3.114, -13.578, 3.149, 1.957, 0, 7.574, -1.665, 0.661)\n    coef.matrix&#91;\"bp.treated\", , \"female\"&#93; &lt;-\n        c(-29.799, 4.884, 13.54, -3.114, -13.578, 3.149, 2.019, 0, 7.574, -1.665, 0.661)\n\n    d.xb &lt;- function(.df, c) {\n        .df$cvs &lt;- rowSums(data.frame(t(apply(.df, 1, function(x) {x * c}))))\n        return(.df$cvs - mean(.df$cvs, na.rm = TRUE))\n    }\n\n    ar10 &lt;- function(s, x) {\n        return(1 - s**exp(x))\n    }\n\n    df &lt;- data.frame(log(age), log(age)**2,\n                     log(tc), log(age) * log(tc),\n                     log(hdl.c), log(age) * log(hdl.c),\n                     log(sbp), log(age)*log(sbp),\n                     cur.smoker,\n                     log(age) * cur.smoker,\n                     dm,\n                     gender, bpl)\n\n    for (i in 1:2) {\n        for (j in 1:2) {\n            df&#91;gender == i &amp; bpl == j - 1, \"risk.ascvd.10\"&#93; &lt;-\n                ar10(b.surv&#91;i&#93;, d.xb(df&#91;gender == i &amp; bpl == j - 1, 1:11&#93;, coef.matrix&#91;j, , i&#93;))\n        }\n    }\n\n    return(df$risk.ascvd.10)\n}<\/textarea><\/pre><svg xmlns=\"http:\/\/www.w3.org\/2000\/svg\" style=\"width:24px;height:24px\" fill=\"none\" viewBox=\"0 0 24 24\" stroke=\"currentColor\" stroke-width=\"2\"><path class=\"with-check\" stroke-linecap=\"round\" stroke-linejoin=\"round\" d=\"M9 5H7a2 2 0 00-2 2v12a2 2 0 002 2h10a2 2 0 002-2V7a2 2 0 00-2-2h-2M9 5a2 2 0 002 2h2a2 2 0 002-2M9 5a2 2 0 012-2h2a2 2 0 012 2m-6 9l2 2 4-4\"><\/path><path class=\"without-check\" stroke-linecap=\"round\" stroke-linejoin=\"round\" d=\"M9 5H7a2 2 0 00-2 2v12a2 2 0 002 2h10a2 2 0 002-2V7a2 2 0 00-2-2h-2M9 5a2 2 0 002 2h2a2 2 0 002-2M9 5a2 2 0 012-2h2a2 2 0 012 2\"><\/path><\/svg><\/span><pre class=\"shiki nord\" style=\"background-color: #2e3440ff\" tabindex=\"0\"><code><span class=\"line\"><span style=\"color: #616E88\">## Pooled cohort equation: implementing pooled cohort equation<\/span><\/span>\n<span class=\"line\"><span style=\"color: #616E88\">## AR10 (ten-year absolute risk) = 1 - S^exp(XB - meanXB)<\/span><\/span>\n<span class=\"line\"><span style=\"color: #616E88\">##     S = survival rate at 10 yrs<\/span><\/span>\n<span class=\"line\"><span style=\"color: #616E88\">## women and men, white<\/span><\/span>\n<span class=\"line\"><span style=\"color: #616E88\">## gender M:1, F:2<\/span><\/span>\n<span class=\"line\"><span style=\"color: #616E88\">## bpl = BP lowering treatment, untreated or treated: 0 or 1<\/span><\/span>\n<span class=\"line\"><span style=\"color: #616E88\">## Reference:<\/span><\/span>\n<span class=\"line\"><span style=\"color: #616E88\">## 2013 ACC\/AHA Guideline on the Assessment of Cardiovascular Risk<\/span><\/span>\n<span class=\"line\"><span style=\"color: #616E88\">## A Report of the American College of Cardiology\/American Heart<\/span><\/span>\n<span class=\"line\"><span style=\"color: #616E88\">## Association Task Force on Practice Guidelines<\/span><\/span>\n<span class=\"line\"><span style=\"color: #616E88\">## Circulation. 2014;129&#91;suppl 2&#93;:S49-S73<\/span><\/span>\n<span class=\"line\"><\/span>\n<span class=\"line\"><span style=\"color: #88C0D0\">pce<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">function<\/span><span style=\"color: #ECEFF4\">(<\/span><span style=\"color: #D8DEE9\">gender<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #D8DEE9\">age<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #D8DEE9\">tc<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #D8DEE9\">hdl.c<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #D8DEE9\">sbp<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #D8DEE9\">bpl<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #D8DEE9\">cur.smoker<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #D8DEE9\">dm<\/span><span style=\"color: #ECEFF4\">)<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #ECEFF4\">{<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #D8DEE9\">b.surv<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #88C0D0\">c<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #B48EAD\">0.9144<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">0.9665<\/span><span style=\"color: #D8DEE9FF\">)  <\/span><span style=\"color: #616E88\">## baseline survival in US cohort, for white men &amp; women<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #88C0D0\">names<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">b.surv<\/span><span style=\"color: #D8DEE9FF\">) <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #88C0D0\">c<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">male<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">female<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #D8DEE9FF\">)<\/span><\/span>\n<span class=\"line\"><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #D8DEE9\">coef.matrix<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #88C0D0\">array<\/span><span style=\"color: #D8DEE9FF\">(<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        <\/span><span style=\"color: #D8DEE9\">dim<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">=<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #88C0D0\">c<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #B48EAD\">2<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">11<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">2<\/span><span style=\"color: #D8DEE9FF\">),<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        <\/span><span style=\"color: #D8DEE9\">dimnames<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">=<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">list<\/span><span style=\"color: #ECEFF4\">(<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">            <\/span><span style=\"color: #88C0D0\">c<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">bp.untreated<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">bp.treated<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #ECEFF4\">)<\/span><span style=\"color: #D8DEE9FF\">,<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">            <\/span><span style=\"color: #88C0D0\">c<\/span><span style=\"color: #D8DEE9FF\">(),<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">            <\/span><span style=\"color: #88C0D0\">c<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">male<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">female<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #D8DEE9FF\">)<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">            )<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        )<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #D8DEE9\">coef.matrix<\/span><span style=\"color: #ECEFF4\">&#91;<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">bp.untreated<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #D8DEE9FF\">, , <\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">male<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #ECEFF4\">&#93;<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        <\/span><span style=\"color: #88C0D0\">c<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #B48EAD\">12.344<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">0<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">11.853<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">2.664<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">7.990<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">1.769<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">1.764<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">0<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">7.837<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">1.795<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">0.658<\/span><span style=\"color: #D8DEE9FF\">)<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #D8DEE9\">coef.matrix<\/span><span style=\"color: #ECEFF4\">&#91;<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">bp.treated<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #D8DEE9FF\">, , <\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">male<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #ECEFF4\">&#93;<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        <\/span><span style=\"color: #88C0D0\">c<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #B48EAD\">12.344<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">0<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">11.853<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">2.664<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">7.990<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">1.769<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">1.797<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">0<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">7.837<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">1.795<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">0.658<\/span><span style=\"color: #D8DEE9FF\">)<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #D8DEE9\">coef.matrix<\/span><span style=\"color: #ECEFF4\">&#91;<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">bp.untreated<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #D8DEE9FF\">, , <\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">female<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #ECEFF4\">&#93;<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        <\/span><span style=\"color: #88C0D0\">c<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">29.799<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">4.884<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">13.54<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">3.114<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">13.578<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">3.149<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">1.957<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">0<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">7.574<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">1.665<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">0.661<\/span><span style=\"color: #D8DEE9FF\">)<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #D8DEE9\">coef.matrix<\/span><span style=\"color: #ECEFF4\">&#91;<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">bp.treated<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #D8DEE9FF\">, , <\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">female<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #ECEFF4\">&#93;<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        <\/span><span style=\"color: #88C0D0\">c<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">29.799<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">4.884<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">13.54<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">3.114<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">13.578<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">3.149<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">2.019<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">0<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">7.574<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #B48EAD\">1.665<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">0.661<\/span><span style=\"color: #D8DEE9FF\">)<\/span><\/span>\n<span class=\"line\"><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #88C0D0\">d.xb<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">function<\/span><span style=\"color: #ECEFF4\">(<\/span><span style=\"color: #D8DEE9FF\">.<\/span><span style=\"color: #D8DEE9\">df<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #D8DEE9\">c<\/span><span style=\"color: #ECEFF4\">)<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #ECEFF4\">{<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        .<\/span><span style=\"color: #D8DEE9\">df<\/span><span style=\"color: #81A1C1\">$<\/span><span style=\"color: #D8DEE9\">cvs<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #88C0D0\">rowSums<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #88C0D0\">data.frame<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #88C0D0\">t<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #88C0D0\">apply<\/span><span style=\"color: #D8DEE9FF\">(.<\/span><span style=\"color: #D8DEE9\">df<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">1<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #81A1C1\">function<\/span><span style=\"color: #ECEFF4\">(<\/span><span style=\"color: #D8DEE9\">x<\/span><span style=\"color: #ECEFF4\">)<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #ECEFF4\">{<\/span><span style=\"color: #D8DEE9\">x<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">*<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #D8DEE9\">c<\/span><span style=\"color: #ECEFF4\">}<\/span><span style=\"color: #D8DEE9FF\">))))<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        <\/span><span style=\"color: #81A1C1\">return<\/span><span style=\"color: #ECEFF4\">(<\/span><span style=\"color: #D8DEE9FF\">.<\/span><span style=\"color: #D8DEE9\">df<\/span><span style=\"color: #81A1C1\">$<\/span><span style=\"color: #D8DEE9\">cvs<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #88C0D0\">mean<\/span><span style=\"color: #D8DEE9FF\">(.<\/span><span style=\"color: #D8DEE9\">df<\/span><span style=\"color: #81A1C1\">$<\/span><span style=\"color: #D8DEE9\">cvs<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #D8DEE9\">na.rm<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">=<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">TRUE<\/span><span style=\"color: #ECEFF4\">)<\/span><span style=\"color: #D8DEE9FF\">)<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #ECEFF4\">}<\/span><\/span>\n<span class=\"line\"><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #88C0D0\">ar10<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">function<\/span><span style=\"color: #ECEFF4\">(<\/span><span style=\"color: #D8DEE9\">s<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #D8DEE9\">x<\/span><span style=\"color: #ECEFF4\">)<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #ECEFF4\">{<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        <\/span><span style=\"color: #81A1C1\">return<\/span><span style=\"color: #ECEFF4\">(<\/span><span style=\"color: #B48EAD\">1<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #D8DEE9\">s<\/span><span style=\"color: #81A1C1\">**<\/span><span style=\"color: #88C0D0\">exp<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">x<\/span><span style=\"color: #ECEFF4\">)<\/span><span style=\"color: #D8DEE9FF\">)<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #ECEFF4\">}<\/span><\/span>\n<span class=\"line\"><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #D8DEE9\">df<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #88C0D0\">data.frame<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">age<\/span><span style=\"color: #D8DEE9FF\">), <\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">age<\/span><span style=\"color: #D8DEE9FF\">)<\/span><span style=\"color: #81A1C1\">**<\/span><span style=\"color: #B48EAD\">2<\/span><span style=\"color: #D8DEE9FF\">,<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">                     <\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">tc<\/span><span style=\"color: #D8DEE9FF\">), <\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">age<\/span><span style=\"color: #D8DEE9FF\">) <\/span><span style=\"color: #81A1C1\">*<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">tc<\/span><span style=\"color: #D8DEE9FF\">),<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">                     <\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">hdl.c<\/span><span style=\"color: #D8DEE9FF\">), <\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">age<\/span><span style=\"color: #D8DEE9FF\">) <\/span><span style=\"color: #81A1C1\">*<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">hdl.c<\/span><span style=\"color: #D8DEE9FF\">),<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">                     <\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">sbp<\/span><span style=\"color: #D8DEE9FF\">), <\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">age<\/span><span style=\"color: #D8DEE9FF\">)<\/span><span style=\"color: #81A1C1\">*<\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">sbp<\/span><span style=\"color: #D8DEE9FF\">),<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">                     <\/span><span style=\"color: #D8DEE9\">cur.smoker<\/span><span style=\"color: #D8DEE9FF\">,<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">                     <\/span><span style=\"color: #88C0D0\">log<\/span><span style=\"color: #D8DEE9FF\">(<\/span><span style=\"color: #D8DEE9\">age<\/span><span style=\"color: #D8DEE9FF\">) <\/span><span style=\"color: #81A1C1\">*<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #D8DEE9\">cur.smoker<\/span><span style=\"color: #D8DEE9FF\">,<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">                     <\/span><span style=\"color: #D8DEE9\">dm<\/span><span style=\"color: #D8DEE9FF\">,<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">                     <\/span><span style=\"color: #D8DEE9\">gender<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #D8DEE9\">bpl<\/span><span style=\"color: #D8DEE9FF\">)<\/span><\/span>\n<span class=\"line\"><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #81A1C1\">for<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #ECEFF4\">(<\/span><span style=\"color: #D8DEE9\">i<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">in<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #B48EAD\">1<\/span><span style=\"color: #81A1C1\">:<\/span><span style=\"color: #B48EAD\">2<\/span><span style=\"color: #ECEFF4\">)<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #ECEFF4\">{<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        <\/span><span style=\"color: #81A1C1\">for<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #ECEFF4\">(<\/span><span style=\"color: #D8DEE9\">j<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">in<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #B48EAD\">1<\/span><span style=\"color: #81A1C1\">:<\/span><span style=\"color: #B48EAD\">2<\/span><span style=\"color: #ECEFF4\">)<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #ECEFF4\">{<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">            <\/span><span style=\"color: #D8DEE9\">df<\/span><span style=\"color: #ECEFF4\">&#91;<\/span><span style=\"color: #D8DEE9\">gender<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">==<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #D8DEE9\">i<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&amp;<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #D8DEE9\">bpl<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">==<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #D8DEE9\">j<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #B48EAD\">1<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #A3BE8C\">risk.ascvd.10<\/span><span style=\"color: #ECEFF4\">&quot;<\/span><span style=\"color: #ECEFF4\">&#93;<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&lt;-<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">                ar10<\/span><span style=\"color: #ECEFF4\">(<\/span><span style=\"color: #D8DEE9\">b.surv<\/span><span style=\"color: #ECEFF4\">&#91;<\/span><span style=\"color: #D8DEE9\">i<\/span><span style=\"color: #ECEFF4\">&#93;,<\/span><span style=\"color: #D8DEE9FF\"> d.xb<\/span><span style=\"color: #ECEFF4\">(<\/span><span style=\"color: #D8DEE9\">df<\/span><span style=\"color: #ECEFF4\">&#91;<\/span><span style=\"color: #D8DEE9\">gender<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">==<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #D8DEE9\">i<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">&amp;<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #D8DEE9\">bpl<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">==<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #D8DEE9\">j<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #81A1C1\">-<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #B48EAD\">1<\/span><span style=\"color: #D8DEE9FF\">, <\/span><span style=\"color: #B48EAD\">1<\/span><span style=\"color: #81A1C1\">:<\/span><span style=\"color: #B48EAD\">11<\/span><span style=\"color: #ECEFF4\">&#93;,<\/span><span style=\"color: #D8DEE9FF\"> <\/span><span style=\"color: #D8DEE9\">coef.matrix<\/span><span style=\"color: #ECEFF4\">&#91;<\/span><span style=\"color: #D8DEE9\">j<\/span><span style=\"color: #D8DEE9FF\">, , <\/span><span style=\"color: #D8DEE9\">i<\/span><span style=\"color: #ECEFF4\">&#93;))<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">        <\/span><span style=\"color: #ECEFF4\">}<\/span><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #ECEFF4\">}<\/span><\/span>\n<span class=\"line\"><\/span>\n<span class=\"line\"><span style=\"color: #D8DEE9FF\">    <\/span><span style=\"color: #81A1C1\">return<\/span><span style=\"color: #ECEFF4\">(<\/span><span style=\"color: #D8DEE9\">df<\/span><span style=\"color: #81A1C1\">$<\/span><span style=\"color: #D8DEE9\">risk.ascvd.10<\/span><span style=\"color: #ECEFF4\">)<\/span><\/span>\n<span class=\"line\"><span style=\"color: #ECEFF4\">}<\/span><\/span><\/code><\/pre><\/div>\n","protected":false},"excerpt":{"rendered":"<p>Pooled cohort equation \uc73c\ub85c 10\ub144 \uac04\uc758 \uc2ec\ub1cc\ud608\uad00\uc9c8\ud658\uc758 \ubc1c\ubcd1 \uc704\ud5d8\ub3c4\ub97c \uc608\uce21\ud558\ub294 \uacf5\uc2dd\uc744 R function \uc73c\ub85c \uad6c\ud604\ud558\uc600\uc2b5\ub2c8\ub2e4.<\/p>\n","protected":false},"author":1,"featured_media":1165,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[17,3],"tags":[14,13],"class_list":["post-40","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-biostatistics","category-research","tag-pooled-cohort-equation","tag-r"],"views":2393,"_links":{"self":[{"href":"https:\/\/jazzman.pe.kr\/blog\/index.php\/wp-json\/wp\/v2\/posts\/40","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/jazzman.pe.kr\/blog\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/jazzman.pe.kr\/blog\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/jazzman.pe.kr\/blog\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/jazzman.pe.kr\/blog\/index.php\/wp-json\/wp\/v2\/comments?post=40"}],"version-history":[{"count":4,"href":"https:\/\/jazzman.pe.kr\/blog\/index.php\/wp-json\/wp\/v2\/posts\/40\/revisions"}],"predecessor-version":[{"id":1169,"href":"https:\/\/jazzman.pe.kr\/blog\/index.php\/wp-json\/wp\/v2\/posts\/40\/revisions\/1169"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/jazzman.pe.kr\/blog\/index.php\/wp-json\/wp\/v2\/media\/1165"}],"wp:attachment":[{"href":"https:\/\/jazzman.pe.kr\/blog\/index.php\/wp-json\/wp\/v2\/media?parent=40"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/jazzman.pe.kr\/blog\/index.php\/wp-json\/wp\/v2\/categories?post=40"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/jazzman.pe.kr\/blog\/index.php\/wp-json\/wp\/v2\/tags?post=40"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}