{"id":212,"date":"2023-01-08T19:58:03","date_gmt":"2023-01-08T10:58:03","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/how-to-keep-some-variables-in-stepwise-selection-procedure-in-ezr\/"},"modified":"2024-09-02T20:14:16","modified_gmt":"2024-09-02T11:14:16","slug":"how-to-keep-some-variables-in-stepwise-selection-procedure-in-ezr","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/how-to-keep-some-variables-in-stepwise-selection-procedure-in-ezr\/","title":{"rendered":"EZR \u3067\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u5909\u6570\u9078\u629e\u3092\u884c\u3046\u3068\u304d\u306b\u3042\u308b\u5909\u6570\u3060\u3051\u306f\u5fc5\u305a\u9078\u629e\u3055\u308c\u308b\u3088\u3046\u306b\u3059\u308b\u65b9\u6cd5"},"content":{"rendered":"\n

\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u5909\u6570\u9078\u629e\u306e\u969b\u306b \u3042\u308b\u5909\u6570\u3060\u3051\u306f\u6b8b\u3063\u3066\u3044\u3066\u307b\u3057\u3044\u5834\u5408\u306e\u6307\u5b9a\u65b9\u6cd5<\/p>\n\n\n\n\n\n\n\n

\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u5909\u6570\u9078\u629e\u3092 EZR \u3067\u5b9f\u884c\u3059\u308b\u65b9\u6cd5<\/h2>\n\n\n\n

\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u5909\u6570\u9078\u629e\u3092 EZR \u3067\u5b9f\u884c\u3059\u308b\u306b\u306f\u3001\u30c7\u30fc\u30bf\u3092\u8aad\u307f\u8fbc\u3093\u3060\u5f8c\u3001\u4f8b\u3048\u3070\u91cd\u56de\u5e30\u306e\u30e1\u30cb\u30e5\u30fc\u304b\u3089\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u3092\u9078\u629e\u3059\u308b<\/p>\n\n\n\n

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

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

\u5b9f\u884c\u3059\u308b\u3068\u3001\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u5909\u6570\u9078\u629e\u306f\u3001\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u5b9f\u65bd\u3055\u308c\u308b<\/p>\n\n\n\n

><\/span> res <-<\/span> stepwise<\/span>(<\/span>RegModel.1,<\/span> direction=<\/span>\"backward\/forward\"<\/span>,<\/span> criterion=<\/span>\"AIC\"<\/span>)<\/span>\nDirection:<\/span>  backward\/<\/span>forward\nCriterion:<\/span>  AIC\nStart:<\/span>  AIC=-<\/span>147.32<\/span>\nLogCD34perBW ~<\/span> age +<\/span> gender +<\/span> Group +<\/span> WBC.entry +<\/span> Plt.entry +<\/span>\nApheresisVolume\nDf Sum of Sq    RSS     AIC\n-<\/span> Plt.entry        1<\/span>   0.00864<\/span> 7.3280<\/span> -<\/span>149.24<\/span>\n-<\/span> gender           1<\/span>   0.03892<\/span> 7.3583<\/span> -<\/span>148.95<\/span>\n-<\/span> age              1<\/span>   0.12298<\/span> 7.4423<\/span> -<\/span>148.14<\/span>\n<<\/span>none><\/span>                         7.3194<\/span> -<\/span>147.32<\/span>\n-<\/span> Group            1<\/span>   0.25607<\/span> 7.5754<\/span> -<\/span>146.88<\/span>\n-<\/span> ApheresisVolume  1<\/span>   0.49402<\/span> 7.8134<\/span> -<\/span>144.69<\/span>\n-<\/span> WBC.entry        1<\/span>   0.80924<\/span> 8.1286<\/span> -<\/span>141.88<\/span>\nStep:<\/span>  AIC=-<\/span>149.24<\/span>\nLogCD34perBW ~<\/span> age +<\/span> gender +<\/span> Group +<\/span> WBC.entry +<\/span> ApheresisVolume\nDf Sum of Sq    RSS     AIC\n-<\/span> gender           1<\/span>   0.03840<\/span> 7.3664<\/span> -<\/span>150.87<\/span>\n-<\/span> age              1<\/span>   0.11755<\/span> 7.4455<\/span> -<\/span>150.11<\/span>\n<<\/span>none><\/span>                         7.3280<\/span> -<\/span>149.24<\/span>\n-<\/span> Group            1<\/span>   0.25043<\/span> 7.5784<\/span> -<\/span>148.85<\/span>\n+<\/span> Plt.entry        1<\/span>   0.00864<\/span> 7.3194<\/span> -<\/span>147.32<\/span>\n-<\/span> ApheresisVolume  1<\/span>   0.48539<\/span> 7.8134<\/span> -<\/span>146.69<\/span>\n-<\/span> WBC.entry        1<\/span>   0.96837<\/span> 8.2964<\/span> -<\/span>142.43<\/span>\nStep:<\/span>  AIC=-<\/span>150.87<\/span>\nLogCD34perBW ~<\/span> age +<\/span> Group +<\/span> WBC.entry +<\/span> ApheresisVolume\nDf Sum of Sq    RSS     AIC\n-<\/span> age              1<\/span>   0.12085<\/span> 7.4872<\/span> -<\/span>151.71<\/span>\n<<\/span>none><\/span>                         7.3664<\/span> -<\/span>150.87<\/span>\n-<\/span> Group            1<\/span>   0.23941<\/span> 7.6058<\/span> -<\/span>150.60<\/span>\n+<\/span> gender           1<\/span>   0.03840<\/span> 7.3280<\/span> -<\/span>149.24<\/span>\n+<\/span> Plt.entry        1<\/span>   0.00811<\/span> 7.3583<\/span> -<\/span>148.95<\/span>\n-<\/span> ApheresisVolume  1<\/span>   0.66779<\/span> 8.0342<\/span> -<\/span>146.71<\/span>\n-<\/span> WBC.entry        1<\/span>   0.99045<\/span> 8.3568<\/span> -<\/span>143.91<\/span>\nStep:<\/span>  AIC=-<\/span>151.71<\/span>\nLogCD34perBW ~<\/span> Group +<\/span> WBC.entry +<\/span> ApheresisVolume\nDf Sum of Sq    RSS     AIC\n-<\/span> Group            1<\/span>   0.19479<\/span> 7.6820<\/span> -<\/span>151.89<\/span>\n<<\/span>none><\/span>                         7.4872<\/span> -<\/span>151.71<\/span>\n+<\/span> age              1<\/span>   0.12085<\/span> 7.3664<\/span> -<\/span>150.87<\/span>\n+<\/span> gender           1<\/span>   0.04170<\/span> 7.4455<\/span> -<\/span>150.11<\/span>\n+<\/span> Plt.entry        1<\/span>   0.00283<\/span> 7.4844<\/span> -<\/span>149.74<\/span>\n-<\/span> ApheresisVolume  1<\/span>   0.60062<\/span> 8.0879<\/span> -<\/span>148.23<\/span>\n-<\/span> WBC.entry        1<\/span>   1.08102<\/span> 8.5683<\/span> -<\/span>144.14<\/span>\nStep:<\/span>  AIC=-<\/span>151.89<\/span>\nLogCD34perBW ~<\/span> WBC.entry +<\/span> ApheresisVolume\nDf Sum of Sq    RSS     AIC\n<<\/span>none><\/span>                         7.6820<\/span> -<\/span>151.89<\/span>\n+<\/span> Group            1<\/span>   0.19479<\/span> 7.4872<\/span> -<\/span>151.71<\/span>\n+<\/span> age              1<\/span>   0.07623<\/span> 7.6058<\/span> -<\/span>150.60<\/span>\n+<\/span> gender           1<\/span>   0.03050<\/span> 7.6515<\/span> -<\/span>150.17<\/span>\n+<\/span> Plt.entry        1<\/span>   0.00072<\/span> 7.6813<\/span> -<\/span>149.90<\/span>\n-<\/span> ApheresisVolume  1<\/span>   0.60501<\/span> 8.2870<\/span> -<\/span>148.51<\/span>\n-<\/span> WBC.entry        1<\/span>   0.99257<\/span> 8.6746<\/span> -<\/span>145.26<\/span>\n><\/span> summary<\/span>(<\/span>res)<\/span>\nCall:<\/span>\nlm<\/span>(<\/span>formula =<\/span> LogCD34perBW ~<\/span> WBC.entry +<\/span> ApheresisVolume,<\/span> data =<\/span> TempDF)<\/span>\nResiduals:<\/span>\nMin       1Q   Median       3Q      Max\n-<\/span>0.76985<\/span> -<\/span>0.22692<\/span>  0.03478<\/span>  0.20694<\/span>  0.66902<\/span>\nCoefficients:<\/span>\nEstimate  Std. Error t value Pr<\/span>(<\/span>>|<\/span>t|<\/span>)<\/span>\n(<\/span>Intercept)<\/span>     -<\/span>0.90751397<\/span>  0.32874553<\/span>  -<\/span>2.761<\/span>  0.00741<\/span> **<\/span>\nWBC.entry        0.00009510<\/span>  0.00003208<\/span>   2.964<\/span>  0.00418<\/span> **<\/span>\nApheresisVolume  0.05761960<\/span>  0.02489841<\/span>   2.314<\/span>  0.02369<\/span> *<\/span>\n---<\/span>\nSignif. codes:<\/span>  0<\/span> '***'<\/span> 0.001<\/span> '**'<\/span> 0.01<\/span> '*'<\/span> 0.05<\/span> '.'<\/span> 0.1<\/span> ' '<\/span> 1<\/span>\nResidual standard error:<\/span> 0.3361<\/span> on 68<\/span> degrees of freedom\nMultiple R-<\/span>squared:<\/span>  0.1725<\/span>,<\/span> Adjusted R-<\/span>squared:<\/span>  0.1482<\/span>\nF<\/span>-<\/span>statistic:<\/span>  7.09<\/span> on 2<\/span> and 68<\/span> DF,<\/span>  p-<\/span>value:<\/span> 0.001597<\/span>\n<\/code><\/pre>\n\n\n\n
\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u5909\u6570\u9078\u629e\u3067 \u3042\u308b\u5909\u6570\u3060\u3051\u306f\u6b8b\u3057\u3066\u304a\u304d\u305f\u3044\u5834\u5408<\/h2>\n\n\n\n
\u3053\u306e\u3068\u304d\u3001Group\u306f\u8131\u843d\u3057\u3066\u3057\u307e\u3063\u3066\u3044\u308b\u304c\u3001\u672c\u5f53\u306f Group \u306f\u3069\u3046\u3057\u3066\u3082\u8aac\u660e\u5909\u6570\u3068\u3057\u3066\u6b8b\u3057\u3066\u304a\u304d\u305f\u304b\u3063\u305f\u3068\u3059\u308b<\/p>\n\n\n\n
\u305d\u306e\u6642\u306f\u3069\u3046\u3059\u308b\u304b\u3068\u3044\u3046\u3068\u3001step() \u3068\u3044\u3046\u95a2\u6570\u3092\u4f7f\u3046<\/p>\n\n\n\n
R \u30b9\u30af\u30ea\u30d7\u30c8\u7a93\u306b\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u66f8\u3044\u3066\u5b9f\u884c\u3059\u308b<\/p>\n\n\n\n
step<\/span>(<\/span>RegModel.1,<\/span> scope=<\/span>list<\/span>(<\/span>lower=~<\/span>Group))<\/span>\n<\/code><\/pre>\n\n\n\n
RegModel.1 \u306f\u3001EZR \u304c\u81ea\u52d5\u3067\u540d\u4ed8\u3051\u305f\u91cd\u56de\u5e30\u30e2\u30c7\u30eb\u30aa\u30d6\u30b8\u30a7\u30af\u30c8<\/p>\n\n\n\n
scope \u304c\u3001\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u306e\u9069\u5fdc\u7bc4\u56f2\u3067\u3001list() \u306e\u4e2d\u306b\u3001lower= \u3067\u3001\u30c1\u30eb\u30c0 ( ~ ) \u306e\u3042\u3068\u306b\u6b8b\u3057\u305f\u3044\u5909\u6570\u540d\u3092\u66f8\u304d\u5165\u308c\u308b<\/p>\n\n\n\n
\u3053\u306e\u3088\u3046\u306b\u3059\u308b\u3068\u3001\u4e0b\u8a18\u306e\u3088\u3046\u306b Group \u3092\u6700\u5f8c\u307e\u3067\u6b8b\u305b\u308b<\/p>\n\n\n\n
><\/span> step<\/span>(<\/span>RegModel.1,<\/span> scope=<\/span>list<\/span>(<\/span>lower=~<\/span>Group))<\/span>\nStart:<\/span>  AIC=-<\/span>147.32<\/span>\nLogCD34perBW ~<\/span> age +<\/span> gender +<\/span> Group +<\/span> WBC.entry +<\/span> Plt.entry +<\/span>\nApheresisVolume\nDf Sum of Sq    RSS     AIC\n-<\/span> Plt.entry        1<\/span>   0.00864<\/span> 7.3280<\/span> -<\/span>149.24<\/span>\n-<\/span> gender           1<\/span>   0.03892<\/span> 7.3583<\/span> -<\/span>148.95<\/span>\n-<\/span> age              1<\/span>   0.12298<\/span> 7.4423<\/span> -<\/span>148.14<\/span>\n<<\/span>none><\/span>                         7.3194<\/span> -<\/span>147.32<\/span>\n-<\/span> ApheresisVolume  1<\/span>   0.49402<\/span> 7.8134<\/span> -<\/span>144.69<\/span>\n-<\/span> WBC.entry        1<\/span>   0.80924<\/span> 8.1286<\/span> -<\/span>141.88<\/span>\nStep:<\/span>  AIC=-<\/span>149.24<\/span>\nLogCD34perBW ~<\/span> age +<\/span> gender +<\/span> Group +<\/span> WBC.entry +<\/span> ApheresisVolume\nDf Sum of Sq    RSS     AIC\n-<\/span> gender           1<\/span>   0.03840<\/span> 7.3664<\/span> -<\/span>150.87<\/span>\n-<\/span> age              1<\/span>   0.11755<\/span> 7.4455<\/span> -<\/span>150.11<\/span>\n<<\/span>none><\/span>                         7.3280<\/span> -<\/span>149.24<\/span>\n-<\/span> ApheresisVolume  1<\/span>   0.48539<\/span> 7.8134<\/span> -<\/span>146.69<\/span>\n-<\/span> WBC.entry        1<\/span>   0.96837<\/span> 8.2964<\/span> -<\/span>142.43<\/span>\nStep:<\/span>  AIC=-<\/span>150.87<\/span>\nLogCD34perBW ~<\/span> age +<\/span> Group +<\/span> WBC.entry +<\/span> ApheresisVolume\nDf Sum of Sq    RSS     AIC\n-<\/span> age              1<\/span>   0.12085<\/span> 7.4872<\/span> -<\/span>151.71<\/span>\n<<\/span>none><\/span>                         7.3664<\/span> -<\/span>150.87<\/span>\n-<\/span> ApheresisVolume  1<\/span>   0.66779<\/span> 8.0342<\/span> -<\/span>146.71<\/span>\n-<\/span> WBC.entry        1<\/span>   0.99045<\/span> 8.3568<\/span> -<\/span>143.91<\/span>\nStep:<\/span>  AIC=-<\/span>151.71<\/span>\nLogCD34perBW ~<\/span> Group +<\/span> WBC.entry +<\/span> ApheresisVolume\nDf Sum of Sq    RSS     AIC\n<<\/span>none><\/span>                         7.4872<\/span> -<\/span>151.71<\/span>\n-<\/span> ApheresisVolume  1<\/span>   0.60062<\/span> 8.0879<\/span> -<\/span>148.23<\/span>\n-<\/span> WBC.entry        1<\/span>   1.08102<\/span> 8.5683<\/span> -<\/span>144.14<\/span>\nCall:<\/span>\nlm<\/span>(<\/span>formula =<\/span> LogCD34perBW ~<\/span> Group +<\/span> WBC.entry +<\/span> ApheresisVolume,<\/span>\ndata =<\/span> TempDF)<\/span>\nCoefficients:<\/span>\n(<\/span>Intercept)<\/span>            Group        WBC.entry  ApheresisVolume\n-<\/span>1.09045844<\/span>       0.10552746<\/span>       0.00009989<\/span>       0.05741134<\/span>\n<\/code><\/pre>\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>
\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u5909\u6570\u9078\u629e\u3067\u8907\u6570\u306e\u5909\u6570\u3092\u6b8b\u3057\u305f\u3044\u5834\u5408<\/h2>\n\n\n\n
\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u5909\u6570\u9078\u629e\u3067\u8907\u6570\u306e\u5909\u6570\u3092\u6b8b\u3057\u305f\u3044\u5834\u5408\u306f\u3001\u30d7\u30e9\u30b9\u3067\u3064\u306a\u3044\u3067\u6307\u5b9a\u3059\u308b<\/p>\n\n\n\n
\u4f8b\u3048\u3070\u3001lower = ~Group + age + gender \u306a\u3069\u3068\u3059\u308b<\/p>\n\n\n\n
\u5b9f\u884c\u7d50\u679c\u306f\u4ee5\u4e0b\u306e\u901a\u308a<\/p>\n\n\n\n
><\/span> step<\/span>(<\/span>RegModel.1,<\/span> scope=<\/span>list<\/span>(<\/span>lower=~<\/span>Group+<\/span>age+<\/span>gender))<\/span>\nStart:<\/span>  AIC=-<\/span>147.32<\/span>\nLogCD34perBW ~<\/span> age +<\/span> gender +<\/span> Group +<\/span> WBC.entry +<\/span> Plt.entry +<\/span>\nApheresisVolume\nDf Sum of Sq    RSS     AIC\n-<\/span> Plt.entry        1<\/span>   0.00864<\/span> 7.3280<\/span> -<\/span>149.24<\/span>\n<<\/span>none><\/span>                         7.3194<\/span> -<\/span>147.32<\/span>\n-<\/span> ApheresisVolume  1<\/span>   0.49402<\/span> 7.8134<\/span> -<\/span>144.69<\/span>\n-<\/span> WBC.entry        1<\/span>   0.80924<\/span> 8.1286<\/span> -<\/span>141.88<\/span>\nStep:<\/span>  AIC=-<\/span>149.24<\/span>\nLogCD34perBW ~<\/span> age +<\/span> gender +<\/span> Group +<\/span> WBC.entry +<\/span> ApheresisVolume\nDf Sum of Sq    RSS     AIC\n<<\/span>none><\/span>                         7.3280<\/span> -<\/span>149.24<\/span>\n-<\/span> ApheresisVolume  1<\/span>   0.48539<\/span> 7.8134<\/span> -<\/span>146.69<\/span>\n-<\/span> WBC.entry        1<\/span>   0.96837<\/span> 8.2964<\/span> -<\/span>142.43<\/span>\nCall:<\/span>\nlm<\/span>(<\/span>formula =<\/span> LogCD34perBW ~<\/span> age +<\/span> gender +<\/span> Group +<\/span> WBC.entry +<\/span>\nApheresisVolume,<\/span> data =<\/span> GCSF_CD34)<\/span>\nCoefficients:<\/span>\n(<\/span>Intercept)<\/span>              age      gender[<\/span>T.M]<\/span>            Group        WBC.entry\n-<\/span>0.96183059<\/span>      -<\/span>0.00308517<\/span>       0.05025930<\/span>       0.12130618<\/span>       0.00009524<\/span>\nApheresisVolume\n0.05569276<\/span>\n<\/code><\/pre>\n\n\n\n
\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u5909\u6570\u9078\u629e\u3067\u306f\u306a\u304f\u3001\u6e1b\u5c11\u6cd5\u3084\u5897\u52a0\u6cd5\u3067\u884c\u3044\u305f\u3044\u5834\u5408<\/h2>\n\n\n\n
step() \u95a2\u6570\u306e\u5f15\u6570\u306b direction= \u3068\u3044\u3046\u306e\u304c\u3042\u308b\u306e\u3067\u3001\u3053\u308c\u3092\u4f7f\u3046<\/p>\n\n\n\n
\u6e1b\u5c11\u6cd5\u306e\u5834\u5408\u306f\u3001direction=”backward” \u3068\u3059\u308b<\/p>\n\n\n\n
\u5897\u52a0\u6cd5\u306e\u5834\u5408\u306f\u3001direction=”forward” \u3068\u3059\u308b<\/p>\n\n\n\n
\u307e\u3068\u3081<\/h2>\n\n\n\n
\u30b9\u30c6\u30c3\u30d7\u30ef\u30a4\u30ba\u5909\u6570\u9078\u629e\u306e\u3068\u304d\u306b\u3001\u3042\u308b\u5909\u6570\u3060\u3051\u306f\u3069\u3046\u3057\u3066\u3082\u6b8b\u3057\u305f\u3044\u3068\u3044\u3046\u5834\u5408\u306e\u6307\u5b9a\u306e\u4ed5\u65b9\u3092\u89e3\u8aac\u3057\u305f<\/p>\n\n\n\n
\u30e2\u30c7\u30eb\u30aa\u30d6\u30b8\u30a7\u30af\u30c8\u3092\u4f5c\u6210\u3057\u305f\u5f8c\u3001step() \u306e\u5f15\u6570\u3068\u3057\u3066\u3001scope=list(lower=~) \u306e\u30c1\u30eb\u30c0\u306e\u53f3\u306b\u6b8b\u3057\u305f\u3044\u5909\u6570\u3092\u66f8\u304d\u5165\u308c\u3066 step() \u3092\u5b9f\u884c\u3059\u308b<\/p>\n\n\n\n
\u8907\u6570\u306e\u5909\u6570\u3092\u6b8b\u3057\u305f\u3044\u5834\u5408\u306f\u30d7\u30e9\u30b9\u8a18\u53f7\u3067\u4ed8\u3051\u8db3\u3059<\/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
\u53c2\u8003\u30b5\u30a4\u30c8<\/h2>\n\n\n\n
[R] for parameter ‘keep’ in ‘step’<\/a><\/p>\n\n\n\n
r – Force step() to keep a certain valuable – Stack Overflow<\/a><\/p>\n\n\n\n