微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

删除缺失的数据值

如何解决删除缺失的数据值

删除一个原始帖子,以便能够发布更大版本的数据集。实际上总共有 418 行。

这是我正在进行的生存分析的数据。第一列是 ID 号,其他列标记为 V2 - V20。有很多缺失的数据用“.”表示。

我使用 coxph() 函数获得以下信息:

#Saves survival time in vector "time".
time = surv.df[,"V2"]

#Saves information about censoring status in vector "status" in Cox's regression; death code should always be set to 1.
state = 1 - surv.df[,"V3"]

#Conduct Cox's regression analysis.
library(survival)
surv.cox = coxph(Surv(time,state == 1) ~ V4 + V5 + V6 + V7 + V8 + V9 + V10 + V11 + V12 + V13 + V14 + V15 + V16 + V17 + V18 + V19 + V20,data = surv.df)

Error in fitter(X,Y,istrat,offset,init,control,weights = weights,: 
  NA/NaN/Inf in foreign function call (arg 6)
In addition: Warning messages:
1: In fitter(X,:
  Ran out of iterations and did not converge
2: In fitter(X,:
  one or more coefficients may be infinite

我假设问题在于缺失值。我尝试使用以下方法非常粗略地清理数据:

#Allows us to use the select() function.
library(tidyverse)

#Find the missing value in each column.
surv.df[! is.na(surv.df)]

#Select the data column needed in an automatic data frame.
s1 = select(surv.df,V2)
s1[! is.na(s1)]

s2 = select(surv.df,V3)
s2[! is.na(s2)]

s3 = select(surv.df,V4)
s3[! is.na(s3)]
s3 = s3[-c(313:418),]

等等。然而,我收效甚微。我想我的第一个问题将是,如何分析具有这么多缺失值的数据(由于字符限制,其中很多都在我无法发布的行中)?如何删除所有缺失值?一旦缺失值‘.’,数据是否可分析?被删除了?

  1  400 2 1 21464 1 1 1 1 1.0 14.5  261 2.60 156  1718.0 137.95 172 190 12.2 4
  2 4500 0 1 20617 1 0 1 1 0.0  1.1  302 4.14  54  7394.8 113.52  88 221 10.6 3
  3 1012 2 1 25594 0 0 0 0 0.5  1.4  176 3.48 210   516.0  96.10  55 151 12.0 4
  4 1925 2 1 19994 1 0 1 1 0.5  1.8  244 2.54  64  6121.8  60.63  92 183 10.3 4
  5 1504 1 2 13918 1 0 1 1 0.0  3.4  279 3.53 143   671.0 113.15  72 136 10.9 3
  6 2503 2 2 24201 1 0 1 0 0.0  0.8  248 3.98  50   944.0  93.00  63   . 11.0 3
  7 1832 0 2 20284 1 0 1 0 0.0  1.0  322 4.09  52   824.0  60.45 213 204  9.7 3
  8 2466 2 2 19379 1 0 0 0 0.0  0.3  280 4.00  52  4651.2  28.38 189 373 11.0 3
  9 2400 2 1 15526 1 0 0 1 0.0  3.2  562 3.08  79  2276.0 144.15  88 251 11.0 2
 10   51 2 2 25772 1 1 0 1 1.0 12.6  200 2.74 140   918.0 147.25 143 302 11.5 4
 11 3762 2 2 19619 1 0 1 1 0.0  1.4  259 4.16  46  1104.0  79.05  79 258 12.0 4
 12  304 2 2 21600 1 0 0 1 0.0  3.6  236 3.52  94   591.0  82.15  95  71 13.6 4
 13 3577 0 2 16688 1 0 0 0 0.0  0.7  281 3.85  40  1181.0  88.35 130 244 10.6 3
 14 1217 2 2 20535 0 1 1 0 1.0  0.8    . 2.27  43   728.0  71.00   . 156 11.0 4
 15 3584 2 1 23612 1 0 0 0 0.0  0.8  231 3.87 173  9009.8 127.71  96 295 11.0 3
 16 3672 0 2 14772 1 0 0 0 0.0  0.7  204 3.66  28   685.0  72.85  58 198 10.8 3
 17  769 2 2 19060 1 0 1 0 0.0  2.7  274 3.15 159  1533.0 117.80 128 224 10.5 4
 18  131 2 1 19698 1 0 1 1 1.0 11.4  178 2.80 588   961.0 280.55 200 283 12.4 4
 19 4232 0 1 18102 1 0 1 0 0.5  0.7  235 3.56  39  1881.0  93.00 123 209 11.0 3
 20 1356 2 2 21898 1 0 1 0 0.0  5.1  374 3.51 140  1919.0 122.45 135 322 13.0 4
 21 3445 0 2 23445 0 0 1 1 0.0  0.6  252 3.83  41   843.0  65.10  83 336 11.4 4
 22  673 2 1 20555 1 0 0 1 0.0  3.4  271 3.63 464  1376.0 120.90  55 173 11.6 4
 23  264 2 2 20442 1 1 1 1 1.0 17.4  395 2.94 558  6064.8 227.04 191 214 11.7 4
 24 4079 2 1 16261 0 0 1 0 0.0  2.1  456 4.00 124  5719.0 221.88 230  70  9.9 2
 25 4127 0 2 16463 1 0 0 0 0.0  0.7  298 4.10  40   661.0 106.95  66 324 11.3 2
 26 1444 2 2 19002 1 0 1 1 0.0  5.2 1128 3.68  53  3228.0 165.85 166 421  9.9 3
 27   77 2 2 19884 1 1 1 1 0.5 21.6  175 3.31 221  3697.4 101.91 168  80 12.0 4
 28  549 2 2 16417 1 1 1 1 1.0 17.2  222 3.23 209  1975.0 189.10 195 144 13.0 4
 29 4509 0 2 23331 1 0 0 0 0.0  0.7  370 3.78  24  5833.0  73.53  86 390 10.6 2
 30  321 2 2 15116 1 0 1 1 0.0  3.6  260 2.54 172  7277.0 121.26 158 124 11.0 4
 31 3839 2 2 15177 1 0 1 0 0.0  4.7  296 3.44 114  9933.2 206.40 101 195 10.3 2
 32 4523 0 2 19722 1 0 1 0 0.0  1.8  262 3.34 101  7277.0  82.56 158 286 10.6 4
 33 3170 2 2 18731 1 0 0 0 0.0  0.8  210 3.19  82  1592.0 218.55 113 180 12.0 3
 34 3933 0 1 19015 1 0 0 0 0.0  0.8  364 3.70  37  1840.0 170.50  64 273 10.5 2
 35 2847 2 2 17758 1 0 0 0 0.0  1.2  314 3.20 201 12258.8  72.24 151 431 10.6 3
 36 3611 0 2 20604 1 0 0 0 0.0  0.3  172 3.39  18   558.0  71.30  96 311 10.6 2
 37  223 2 1 22546 1 1 1 0 1.0  7.1  334 3.01 150  6931.2 180.60 118 102 12.0 4
 38 3244 2 2 13378 1 0 1 1 0.0  3.3  383 3.53 102  1234.0 137.95  87 234 11.0 4
 39 2297 2 1 20232 1 0 1 0 0.0  0.7  282 3.00  52  9066.8  72.24 111 563 10.6 4
 40 4467 0 1 17046 1 0 0 0 0.0  1.3    . 3.34 105 11046.6 104.49   . 358 11.0 4
 41 1350 2 1 12285 1 0 1 0 0.0  6.8    . 3.26  96  1215.0 151.90   . 226 11.7 4
 42 4453 0 2 12307 1 0 1 1 0.0  2.1    . 3.54 122  8778.0  56.76   . 344 11.0 4
 43 4556 0 1 17850 1 0 0 0 0.0  1.1  361 3.64  36  5430.2  67.08  89 203 10.6 2
 44 3428 2 2 13727 1 0 1 1 1.0  3.3  299 3.55 131  1029.0 119.35  50 199 11.7 3
 45 4025 0 2 15265 1 0 0 0 0.0  0.6    . 3.93  19  1826.0  71.30   . 474 10.9 2
 46 2256 2 1 16728 1 0 1 0 0.0  5.7  482 2.84 161 11552.0 136.74 165 518 12.7 3
 47 2576 0 2 17323 1 0 0 0 0.0  0.5  316 3.65  68  1716.0 187.55  71 356  9.8 3
 48 4427 0 2 17947 0 0 0 0 0.0  1.9  259 3.70 281 10396.8 188.34 178 214 11.0 3
 49  708 2 2 22336 1 0 1 0 0.0  0.8    . 3.82  58   678.0  97.65   . 233 11.0 4
 50 2598 2 1 19544 1 0 1 0 0.0  1.1  257 3.36  43  1080.0 106.95  73 128 10.6 4
 51 3853 2 2 19025 1 0 0 0 0.0  0.8  276 3.60  54  4332.0  99.33 143 273 10.6 2
 52 2386 2 1 18460 0 0 0 0 0.0  6.0  614 3.70 158  5084.4 206.40  93 362 10.6 1
 53 1000 2 1 24621 1 0 1 0 0.0  2.6    . 3.10  94  6456.2  56.76   . 214 11.0 4
 54 1434 2 1 14317 1 1 1 1 1.0  1.3  288 3.40 262  5487.2  73.53 125 254 11.0 4
 55 1360 2 1 24020 0 0 0 0 0.0  1.8  416 3.94 121 10165.0  79.98 219 213 11.0 3
 56 1847 2 2 12279 1 0 1 1 0.0  1.1  498 3.80  88 13862.4  95.46 319 365 10.6 2
 57 3282 2 1 19567 1 0 1 0 0.5  2.3  260 3.18 231 11320.2 105.78  94 216 12.4 3
 58 4459 0 1 16279 0 0 0 0 0.0  0.7  242 4.08  73  5890.0  56.76 118   . 10.6 1
 59 2224 2 1 14754 1 0 1 1 0.0  0.8  329 3.50  49  7622.8 126.42 124 321 10.6 3
 60 4365 0 1 21324 1 0 0 0 0.0  0.9  604 3.40  82   876.0  71.30  58 228 10.3 3
 61 4256 0 2 16034 0 0 0 0 0.0  0.6  216 3.94  28   601.0  60.45 188 211 13.0 1
 62 3090 2 2 22173 1 1 0 0 0.0  1.3  302 2.75  58  1523.0  43.40 112 329 13.2 4
 63  859 2 2 17031 1 0 0 1 1.0 22.5  932 3.12  95  5396.0 244.90 133 165 11.6 3
 64 1487 2 2 22977 1 0 1 0 0.0  2.1  373 3.50  52  1009.0 150.35 188 178 11.0 3
 65 3992 0 1 14684 1 0 0 0 0.0  1.2  256 3.60  74   724.0 141.05 108 430 10.0 1
 66 4191 2 1 16967 0 0 1 0 0.0  1.4  427 3.70 105  1909.0 182.90 171 123 11.0 3
 67 2769 2 2 18733 1 0 0 0 0.0  1.1  466 3.91  84  1787.0 328.60 185 261 10.0 3
 68 4039 0 1 11912 1 0 0 0 0.0  0.7  174 4.09  58   642.0  71.30  46 203 10.6 3
 69 1170 2 1 18021 1 0 1 1 0.5 20.0  652 3.46 159  3292.0 215.45 184 227 12.4 3
 70 3458 0 1 20600 1 0 0 0 0.0  0.6    . 4.64  20   666.0  54.25   . 265 10.6 2
 71 4196 0 2 17841 1 0 1 0 0.0  1.2  258 3.57  79  2201.0 120.90  76 410 11.5 4
 72 4184 0 2 11868 1 0 0 0 0.0  0.5  320 3.54  51  1243.0 122.45  80 225 10.0 3
 73 4190 0 2 14060 1 0 0 0 0.0  0.7  132 3.60  17   423.0  49.60  56 265 11.0 1
 74 1827 2 1 18964 1 0 1 1 0.0  8.4  558 3.99 280   967.0  89.90 309 278 11.0 4
 75 1191 2 1 15895 1 1 1 1 0.5 17.1  674 2.53 207  2078.0 182.90 598 268 11.5 4
 76   71 2 1 18972 1 0 1 1 0.5 12.2  394 3.08 111  2132.0 155.00 243 165 11.6 4
 77  326 2 2 18199 1 0 1 1 0.5  6.6  244 3.41 199  1819.0 170.50  91 132 12.1 3
 78 1690 2 1 17512 1 0 1 0 0.0  6.3  436 3.02  75  2176.0 170.50 104 236 10.6 4
 79 3707 0 1 16990 1 0 1 0 0.0  0.8  315 4.24  13  1637.0 170.50  70 426 10.9 3
 80  890 2 2 24622 0 0 1 0 0.0  7.2  247 3.72 269  1303.0 176.70  91 360 11.2 4
 81 2540 2 1 23107 1 0 1 1 0.0 14.4  448 3.65  34  1218.0  60.45 318 385 11.7 4
 82 3574 2 1 24585 1 0 0 0 0.0  4.5  472 4.09 154  1580.0 117.80 272 412 11.1 3
 83 4050 0 1 20459 1 0 1 0 0.5  1.3  250 3.50  48  1138.0  71.30 100  81 12.9 4
 84 4032 0 2 20392 1 0 0 0 0.0  0.4  263 3.76  29  1345.0 137.95  74 181 11.2 3
 85 3358 2 2 17246 1 0 1 0 0.0  2.1  262 3.48  58  2045.0  89.90  84 225 11.5 4
 86 1657 2 1 19270 1 0 1 1 0.0  5.0 1600 3.21  75  2656.0  82.15 174 181 10.9 3
 87  198 2 1 13616 1 0 0 0 0.0  1.1  345 4.40  75  1860.0 218.55  72 447 10.7 3
 88 2452 0 2 15119 1 0 0 0 0.5  0.6  296 4.06  37  1032.0  80.60  83 442 12.0 3
 89 1741 2 1 19155 1 0 1 0 0.0  2.0  408 3.65  50  1083.0 110.05  98 200 11.4 2
 90 2689 2 1 12227 0 0 0 0 0.0  1.6  660 4.22  94  1857.0 151.90 155 337 11.0 2
 91  460 2 2 16658 1 0 1 1 0.5  5.0  325 3.47 110  2460.0 246.45  56 430 11.9 4
 92  388 2 1 28018 1 1 0 0 1.0  1.4  206 3.13  36  1626.0  86.80  70 145 12.2 4
 93 3913 0 1 13344 1 0 0 0 0.0  1.3  353 3.67  73  2039.0 232.50  68 380 11.1 2
 94  750 2 1 19693 1 0 1 1 0.0  3.2  201 3.11 178  1212.0 159.65  69 188 11.8 4
 95  130 2 2 16944 1 1 1 1 1.0 17.4    . 2.64 182   559.0 119.35   . 401 11.7 2
 96 3850 0 1 17841 1 0 0 0 0.0  1.0    . 3.70  33  1258.0  99.20   . 338 10.4 3
 97  611 2 2 26259 0 0 1 0 0.5  2.0  420 3.26  62  3196.0  77.50  91 344 11.4 3
 98 3823 0 1 10550 1 0 0 0 0.0  1.0  239 3.77  77  1877.0  97.65 101 312 10.2 1
 99 3820 0 2 17703 0 0 0 0 0.0  1.8  460 3.35 148  1472.0 108.50 118 172 10.2 2
100  552 2 2 18799 0 0 1 0 0.0  2.3  178 3.00 145   746.0 178.25 122 119 12.0 4
101 3581 0 2 16418 1 0 0 0 0.0  0.9  400 3.60  31  1689.0 164.30 166 327 10.4 3
102 3099 0 1 20662 1 0 0 0 0.0  0.9  248 3.97 172   646.0  62.00  84 128 10.1 1
103  110 2 2 17884 1 1 1 1 1.0  2.5  188 3.67  57  1273.0 119.35 102 110 11.1 4
104 3086 2 1 15712 1 0 0 0 0.0  1.1  303 3.64  20  2108.0 128.65  53 349 11.1 2
105 3092 1 2 12433 1 0 1 0 0.0  1.1  464 4.20  38  1644.0 151.90 102 348 10.3 3
106 3222 2 1 25023 1 1 1 0 0.0  2.1    . 3.90  50  1087.0 103.85   . 137 10.6 2
107 3388 0 2 22836 1 0 0 0 0.0  0.6  212 4.03  10   648.0  71.30  77 316 17.1 1
108 2583 2 1 18393 1 0 0 0 0.0  0.4  127 3.50  14  1062.0  49.60  84 334 10.3 2
109 2504 0 2 16094 1 0 0 0 0.0  0.5  120 3.61  53   804.0 110.05  52 271 10.6 3
110 2105 2 1 14212 1 0 1 1 0.0  1.9  486 3.54  74  1052.0 108.50 109 141 10.9 3
111 2350 1 1 15031 1 0 0 0 0.0  5.5  528 4.18  77  2404.0 172.05  78 467 10.7 3
112 3445 2 2 20256 1 0 1 1 0.0  2.0  267 3.67  89   754.0 196.85  90 136 11.8 4
113  980 2 1 18713 1 0 1 1 0.0  6.7  374 3.74 103   979.0 128.65 100 266 11.1 4
114 3395 2 2 19295 0 0 0 0 0.0  3.2  259 4.30 208  1040.0 110.05  78 268 11.7 3
115 3422 0 2 15574 1 0 0 1 0.0  0.7  303 4.19  81  1584.0 111.60 156 307 10.3 3
116 3336 0 1 22306 1 0 0 1 0.5  3.0  458 3.63  74  1588.0 106.95 382 438  9.9 3
117 1083 2 1 18137 1 0 1 1 0.0  6.5  950 3.11 111  2374.0 170.50 149 354 11.0 4
118 2288 2 1 17844 1 0 1 0 0.0  3.5  390 3.30  67   878.0 137.95  93 207 10.2 3
119  515 2 1 19817 1 0 0 1 0.0  0.6  636 3.83 129   944.0  97.65 114 306  9.5 3
120 2033 1 1 12839 0 0 0 0 0.0  3.5  325 3.98 444   766.0 130.20 210 344 10.6 3
121  191 2 2 24803 0 1 1 0 1.0  1.3  151 3.08  73  1112.0  46.50  49 213 13.2 4
122 3297 0 1 20248 1 0 0 0 0.0  0.6  298 4.13  29   758.0  65.10  85 256 10.7 3
123  971 2 1 16736 1 0 1 1 1.0  5.1    . 3.23  18   790.0 179.80   . 104 13.0 4
124 3069 0 1 19318 0 0 1 0 0.0  0.6  251 3.90  25   681.0  57.35 107 182 10.8 4
125 2468 1 2 17233 1 0 1 0 0.0  1.3  316 3.51  75  1162.0 147.25 137 238 10.0 4
126  824 2 1 19577 1 1 1 1 0.0  1.2  269 3.12   .  1441.0 165.85  68 166 11.1 4
127 3255 0 2 16109 1 0 0 0 0.0  0.5  268 4.08   9  1174.0  86.80  95 453 10.0 2
128 1037 2 1 15322 1 0 1 1 0.0 16.2    . 2.89  42  1828.0 299.15   . 123 12.6 4
129 3239 0 1 23235 1 0 1 0 0.0  0.9  420 3.87  30  1009.0  57.35 232   .  9.7 3
130 1413 2 2 16154 1 0 1 1 0.0 17.4 1775 3.43 205  2065.0 165.85  97 418 11.5 3
131  850 2 2 22646 1 0 1 1 0.0  2.8  242 3.80  74   614.0 136.40 104 121 13.2 4
132 2944 0 1 14812 1 0 0 0 0.0  1.9  448 3.83  60  1052.0 127.10 175 181  9.8 3
133 2796 2 2 22881 0 0 0 0 0.0  1.5  331 3.95  13   577.0 128.65  99 165 10.1 4
134 3149 0 2 15463 1 0 0 0 0.0  0.7  578 3.67  35  1353.0 127.10 105 427 10.7 2
135 3150 0 1 15694 1 0 0 0 0.0  0.4  263 3.57 123   836.0  74.40 121 445 11.0 2
136 3098 0 1 20440 1 0 0 0 0.0  0.8  263 3.35  27  1636.0 116.25  69 206  9.8 2
137 2990 0 1 22960 1 0 0 0 0.0  1.1  399 3.60  79  3472.0 155.00 152 344 10.1 2
138 1297 2 1 18719 0 0 1 0 0.0  7.3  426 3.93 262  2424.0 145.70 218 252 10.5 3
139 2106 0 2 17080 1 0 1 0 0.0  1.1  328 3.31 159  1260.0  94.55 134 142 11.6 4
140 3059 0 1 19751 1 0 1 0 0.0  1.1  290 4.09  38  2120.0 186.00 146 318 10.0 3
141 3050 0 1 17180 1 0 0 0 0.0  0.9  346 3.77  59   794.0 125.55  56 336 10.6 2
142 2419 2 2 20354 1 0 1 0 0.0  1.0  364 3.48  20   720.0 134.85  88 283  9.9 2
143  786 2 2 16839 1 0 1 0 0.0  2.9  332 3.60  86  1492.0 134.85 103 277 11.0 4
144  943 2 2 19098 1 0 1 0 0.5 28.0  556 3.26 152  3896.0 198.40 171 335 10.0 3
145 2976 0 2 18701 1 0 0 1 0.0  0.7  309 3.84  96   858.0  41.85 106 253 11.4 3
146 2615 0 2 12369 1 0 0 0 0.5  1.2    . 3.89  58  1284.0 173.60   . 239  9.4 3
147 2995 0 1 27398 1 0 0 0 0.5  1.2  288 3.37  32   791.0  57.35 114 213 10.7 2
148 1427 2 2 11273 1 0 1 0 0.0  7.2 1015 3.26 247  3836.0 198.40 280 330  9.8 3
149  762 2 1 22574 0 0 1 1 0.5  3.0  257 3.79 290  1664.0 102.30 112 140  9.9 4
150 2891 0 2 12779 1 0 0 1 0.0  1.0    . 3.63  57  1536.0 134.85   . 233 10.0 1
151 2870 0 1 20104 1 0 0 0 0.0  0.9  460 3.03  57   721.0  85.25 174 301  9.4 2
152 1152 2 1 25546 0 0 1 0 0.0  2.3  586 3.01 243  2276.0 114.70 126 339 10.9 3
153 2863 0 1 18118 1 0 0 0 0.0  0.5  217 3.85  68   453.0  54.25  68 270 11.1 1
154  140 2 1 25340 0 0 0 1 1.0  2.4  168 2.56 225  1056.0 120.90  75 108 14.1 3
155 2666 0 2 15909 1 0 1 1 0.5  0.6  220 3.35  57  1620.0 153.45  80 311 11.2 4
156  853 2 2 21699 1 0 1 0 0.0 25.5  358 3.52 219  2468.0 201.50 205 151 11.5 2
157 2835 0 2 17809 1 0 0 0 0.0  0.6  286 3.42  34  1868.0  77.50 206 487 10.0 2
158 2475 1 1 13329 1 0 0 0 0.0  3.4  450 3.37  32  1408.0 116.25 118 313 11.2 2
159 1536 2 2 16714 0 0 0 0 0.0  2.5  317 3.46 217   714.0 130.20 140 207 10.1 3
160 2772 0 2 20955 1 0 0 0 0.0  0.6  217 3.62  13   414.0  75.95 119 224 10.5 3
161 2797 0 2 15612 1 0 0 0 0.0  2.3  502 3.56   4   964.0 120.90 180 269  9.6 2
162  186 2 2 21483 1 0 1 1 0.0  3.2  260 3.19  91   815.0 127.10 101 160 12.0 4
163 2055 2 1 19540 1 0 0 0 0.0  0.3  233 4.08  20   622.0  66.65  68 358  9.9 3
164  264 2 2 15857 1 0 1 1 0.5  8.5    . 3.34 161  1428.0 181.35   .  88 13.3 4
165 1077 2 1 19470 0 0 1 0 0.0  4.0  196 3.45  80  2496.0 133.30 142 212 11.3 4
166 2721 0 2 15105 1 0 1 0 0.0  5.7 1480 3.26  84  1960.0 457.25 108 213  9.5 2
167 1682 2 1 22265 0 0 1 0 0.0  0.9  376 3.86 200  1015.0  83.70 154 238 10.3 4
168 2713 0 2 17442 1 0 1 0 0.0  0.4  257 3.80  44   842.0  97.65 110   .  9.2 2
169 1212 2 2 12963 1 0 0 0 0.0  1.3  408 4.22  67  1387.0 142.60 137 295 10.1 3
170 2692 0 1 17774 1 0 0 0 0.0  1.2  390 3.61  32  1509.0  88.35  52 263  9.0 3
171 2574 0 1 19237 1 0 0 0 0.0  0.5    . 4.52  31   784.0  74.40   . 361 10.1 3
172 2301 0 2 18215 1 0 0 1 0.0  1.3  205 3.34  65  1031.0  91.45 126 217  9.8 3
173 2657 0 1 11058 1 0 1 1 0.0  3.0  236 3.42  76  1403.0  89.90  86 493  9.8 2
174 2644 0 1 20296 1 0 0 0 0.0  0.5    . 3.85  63   663.0  79.05   . 311  9.7 1
175 2624 0 2 19049 1 0 0 0 0.0  0.8  283 3.80 152   718.0 108.50 168 340 10.1 3
176 1492 2 1 15198 1 0 1 1 0.0  3.2    . 3.56  77  1790.0 139.50   . 149 10.1 4
177 2609 0 2 20254 1 0 0 0 0.0  0.9  258 4.01  49   559.0  43.40 133 277 10.4 2
178 2580 0 1 25569 1 0 0 0 0.0  0.6    . 4.08  51   665.0  74.40   . 325 10.2 4
179 2573 0 2 16050 1 0 1 0 0.0  1.8  396 3.83  39  2148.0 102.30 133 278  9.9 4
180 2563 0 2 15548 1 0 0 0 0.0  4.7  478 4.38  44  1629.0 237.15  76 175 10.4 3
181 2556 0 1 16279 1 0 1 1 0.0  1.4  248 3.58  63   554.0  75.95 106  79 10.3 4
182 2555 0 1 20799 1 0 1 0 0.0  0.6    . 3.69 161   674.0  26.35   . 539  9.9 2
183 2241 1 2 14705 1 0 0 0 0.0  0.5  201 3.73  44  1345.0  54.25 145 445 10.1 2
184  974 2 2 13736 1 0 1 0 0.0 11.0  674 3.55 358  2412.0 167.40 140 471  9.8 3
185 2527 0 1 17664 1 0 0 0 0.0  0.8  256 3.54  42  1132.0  74.40  94 192 10.5 3
186 1576 2 1 25873 1 0 0 1 0.5  2.0  225 3.53  51   933.0  69.75  62 200 12.7 3
187  733 2 2 13073 1 0 1 0 0.0 14.0  808 3.43 251  2870.0 153.45 137 268 11.5 3
188 2332 0 1 22873 1 0 1 0 0.0  0.7  187 3.48  41   654.0 120.90  98 164 11.0 4
189 2456 0 2 18499 1 0 1 0 0.0  1.3  360 3.63  52  1812.0  97.65 164 256  9.9 3
190 2504 0 1 19916 1 0 0 1 0.0  2.3    . 3.93  24  1828.0 133.30   . 327 10.2 2
191  216 2 2 19246 1 1 1 1 0.0 24.5 1092 3.35 233  3740.0 147.25 432 399 15.2 4
192 2443 0 1 19256 1 0 1 0 0.0  0.9  308 3.69  67   696.0  51.15 101 344  9.8 4
193  797 2 2 20736 1 0 0 0 0.0 10.8  932 3.19 267  2184.0 161.20 157 382 10.4 4
194 2449 0 1 16216 1 0 0 0 0.0  1.5  293 4.30  50   975.0 125.55  56 336  9.1 2
195 2330 0 1 10795 1 0 1 0 0.0  3.7  347 3.90  76  2544.0 221.65  90 129 11.5 4
196 2363 0 1 20834 1 0 1 1 0.0  1.4  226 3.36  13   810.0  72.85  62 117 11.6 4
197 2365 0 1 16300 1 0 0 0 0.0  0.6  266 3.97  25  1164.0 102.30 102 201 10.1 2
198 2357 0 2 13075 1 0 0 1 0.0  0.7  286 2.90  38  1692.0 141.05  90 381  9.6 2
199 1592 0 1 14872 1 0 0 0 0.0  2.1  392 3.43  52  1395.0 184.45 194 328 10.2 3
200 2318 0 2 11773 1 0 0 1 0.0  4.7  236 3.55 112  1391.0 137.95 114 332  9.9 3
201 2294 0 2 15009 1 0 1 0 0.0  0.6  235 3.20  26  1758.0 106.95  67 228 10.8 4
202 2272 0 1 22514 1 0 0 0 0.0  0.5  223 3.80  15  1044.0  80.60  89 514 10.0 2
203 2221 0 2 13535 1 0 1 0 0.0  0.5  149 4.04 227   598.0  52.70  57 166  9.9 2
204 2090 2 2 22857 1 0 0 0 0.0  0.7  255 3.74  23  1024.0  77.50  58 281 10.2 3
205 2081 2 1 17889 1 1 0 0 0.0  2.5  382 3.55 108  1516.0 238.70   . 126 10.3 3
206 2255 0 1 22642 1 0 0 0 0.0  0.6  213 4.07  12  5300.0  57.35  68 240 11.0 1
207 2171 0 1 26580 1 0 0 0 0.5  0.6    . 3.33  14   733.0  85.25   . 259 10.1 4
208  904 2 1 22388 1 0 1 0 0.0  3.9  396 3.20  58  1440.0 153.45 131 156 10.0 4
209 2216 0 2 19221 1 0 1 1 0.0  0.7  252 4.01  11  1210.0  72.85  58 309  9.5 2
210 2224 0 2 18176 0 0 1 0 0.0  0.9  346 3.37  81  1098.0 122.45  90 298 10.0 2
211 2195 0 2 19327 1 0 0 0 0.0  1.3    . 3.76  27  1282.0 100.75   . 114 10.3 3
212 2176 0 2 17263 1 0 0 0 0.0  1.2  232 3.98  11  1074.0 100.75  99 223  9.9 3
213 2178 0 1 18337 1 0 0 1 0.0  0.5  400 3.40   9  1134.0  96.10  55 356 10.2 3
214 1786 2 2 25329 1 0 1 0 0.0  0.9  404 3.43  34  1866.0  79.05 224 236  9.9 3
215 1080 2 2 15037 1 0 0 0 0.0  5.9 1276 3.85 141  1204.0 203.05 157 216 10.7 3
216 2168 0 1 21610 1 0 0 0 0.0  0.5    . 3.68  20   856.0  55.80   . 146 10.4 3
217  790 2 2 13178 1 0 1 0 0.0 11.4  608 3.31  65  1790.0 151.90 210 298 10.8 4
218 2170 0 1 12636 1 0 0 0 0.0  0.5    . 3.89  29   897.0  66.65   . 423 10.1 1
219 2157 0 2 15601 1 0 0 0 0.0  1.6  215 4.17  67   936.0 134.85  85 176  9.6 3
220 1235 2 1 23241 1 0 0 1 0.0  3.8  426 3.22  96  2716.0 210.80 113 228 10.6 2
221 2050 0 2 20684 1 0 1 0 0.0  0.9  360 3.65  72  3186.0  94.55 154 269  9.7 4
222  597 2 2 16898 1 0 1 0 0.0  4.5  372 3.38 227  2310.0 167.40 135 240 12.4 3
223  334 2 1 22369 1 1 1 0 1.0 14.1  448 2.43 123  1833.0 134.00 155 210 11.0 4
224 1945 0 1 14106 1 0 0 0 0.0  1.0  309 3.66  67  1214.0 158.10 101 309  9.7 3
225 2022 0 1 14161 1 0 0 0 0.0  0.7  274 3.66 108  1065.0  88.35 135 251 10.1 2
226 1978 0 2 20708 1 0 1 0 0.0  0.5  223 3.70  39   884.0  75.95 104 231  9.6 3
227  999 2 1 21532 0 0 0 0 0.0  2.3  316 3.35 172  1601.0 179.80  63 394  9.7 2
228 1967 0 2 13486 1 0 0 0 0.0  0.7  215 3.35  41   645.0  93.00  74 165  9.6 3
229  348 2 1 22797 1 1 1 0 0.5  4.5  191 3.05 200  1020.0 175.15 118 139 11.4 4
230 1979 0 2 12641 1 0 1 1 0.0  3.3  302 3.41  51   310.0  83.70  44  95 11.5 4
231 1165 2 2 21307 1 0 1 1 0.0  3.4  518 1.96 115  2250.0 203.05  90 190 10.7 4
232 1951 0 1 18329 1 0 1 0 0.0  0.4  267 3.02  47  1001.0 133.30  87 265 10.6 3
233 1932 0 1 15591 1 0 1 1 0.0  0.9  514 3.06 412  2622.0 105.40  87 284  9.8 4
234 1776 0 2 12557 1 0 0 0 0.0  0.9  578 3.35  78   976.0 116.25 177 322 11.2 2
235 1882 0 2 12120 1 0 1 0 0.0 13.0 1336 4.16  71  3510.0 209.25 111 338 11.9 3
236 1908 0 1 14019 1 0 1 1 0.0  1.5  253 3.79  67  1006.0 139.50 106 341  9.7 3
237 1882 0 1 21828 1 0 1 0 0.0  1.6  442 2.95 105   820.0  85.25 108 181 10.1 3
238 1874 0 2 24257 1 0 0 0 0.5  0.6  280 3.35   .  1093.0 128.65  81 295  9.8 2
239  694 2 1 17090 1 0 1 1 0.0  0.8  300 2.94 231  1794.0 130.20  99 319 11.2 4
240 1831 0 1 20483 1 0 0 0 0.0  0.4  232 3.72  24   369.0  51.15 139 326 10.1 3
241  837 1 2 15112 1 0 1 1 0.0  4.4  316 3.62 308  1119.0 114.70 322 282  9.8 4
242 1810 0 1 23585 1 0 1 0 0.0  1.9  354 2.97  86  1553.0 196.85 152 277  9.9 3
243  930 2 2 24650 1 0 1 0 0.0  8.0  468 2.81 139  2009.0 198.40 139 233 10.0 4
244 1690 2 1 16374 1 0 0 1 0.0  3.9  350 3.22 121  1268.0 272.80 231 270  9.6 3
245 1790 0 2 16718 1 0 1 0 0.0  0.6  273 3.65  48   794.0  52.70 214 305  9.6 3
246 1435 1 1 12035 1 0 1 0 0.0  2.1  387 3.77  63  1613.0 150.35  33 185 10.1 4
247  732 1 1 15056 1 0 1 0 0.0  6.1 1712 2.83  89  3681.0 158.10 139 297 10.0 3
248 1785 0 2 20241 1 0 1 0 0.0  0.8  324 3.51  39  1237.0  66.65 146 371 10.0 3
249 1783 0 1 17525 1 0 0 1 0.0  1.3  242 3.20  35  1556.0 175.15  71 195 10.6 4
250 1769 0 2 14899 1 0 1 0 0.0  0.6  299 3.36  23  2769.0 220.10  85 303 10.9 4
251 1457 0 1 20810 1 0 0 0 0.0  0.5  227 3.61  40   676.0  83.00 120 249  9.9 2
252 1770 0 1 25006 1 0 1 1 0.0  1.1  246 3.35 116   924.0 113.15  90 317 10.0 4
253 1765 0 1 28650 0 1 1 1 0.0  7.1  243 3.03 380   983.0 158.10 154  97 11.2 4
254  737 1 1 14558 1 0 1 1 0.0  3.1  227 3.75 121  1136.0 110.00  91 264 10.0 3
255 1735 0 2 12897 1 0 1 1 0.0  0.7  193 3.85  35   466.0  53.00 118 156 10.3 3
256 1701 0 1 11485 1 0 0 0 0.0  1.1  336 3.74  48   823.0  84.00 108 242  9.7 3
257 1614 0 1 21281 1 0 0 0 0.0  0.5  280 4.23  36   377.0  56.00 146 227 10.6 2
258 1702 0 1 18806 1 0 0 0 0.0  1.1  414 3.44  80  1003.0  99.00  55 271  9.6 1
259 1615 0 2 21904 1 0 1 0 0.0  3.1  277 2.97  42  1110.0 125.00 126 221  9.8 3
260 1656 0 2 27220 0 0 1 0 0.0  5.6  232 3.59 188  1120.0  98.00 128 248 10.9 4
261 1677 0 2 19126 1 0 1 1 0.0  3.2  375 3.14 129   857.0  89.00   . 375  9.5 3
262 1666 0 2 15628 1 0 1 0 0.0  2.8  322 3.06  65  2562.0  91.00 209 231  9.5 3
263 1301 1 2 12738 1 0 1 1 0.5  1.1  432 3.57  45  1406.0 190.00  77 248 11.4 4
264 1542 1 2 16122 1 0 1 1 0.0  3.4  356 3.12 188  1911.0  92.00 130 318 11.2 3
265 1084 1 2 16941 1 0 1 0 0.0  3.5  348 3.20 121   938.0 120.00 146 296 10.0 4
266 1614 0 1 20567 1 0 0 0 0.0  0.5  318 3.32  52   613.0  70.00 260 279 10.2 3
267  179 2 1 25899 1 1 1 1 1.0  6.6  222 2.33 138   620.0 106.00  91 195 12.1 4
268 1191 2 1 20233 1 1 1 0 0.5  6.4  344 2.75  16   834.0  82.00 179 149 11.0 4
269 1363 0 2 16467 1 0 0 0 0.0  3.6  374 3.50 143  1428.0 188.00  44 151 10.1 2
270 1568 0 1  9598 1 0 1 1 0.0  1.0  448 3.74 102  1128.0  71.00 117 228 10.2 3
271 1569 0 2 18435 1 0 1 0 0.0  1.0  321 3.50  94   955.0 111.00 177 289  9.7 3
272 1525 0 1 14025 1 0 0 0 0.0  0.5  226 2.93  22   674.0  58.00  85 153  9.8 1
273 1558 0 2 17320 1 0 0 1 0.0  2.2  328 3.46  75  1677.0  87.00 116 202  9.6 3
274 1447 1 1 17525 1 0 0 0 0.0  1.6    . 3.07 136  1995.0 128.00   . 372  9.6 4
275 1349 0 1 13995 1 0 0 0 0.0  2.2  572 3.77  77  2520.0  92.00 114 309  9.5 4
276 1481 0 1 18302 1 0 0 0 0.0  1.0  219 3.85  67   640.0 145.00 108  95 10.7 2
277 1434 0 2 12816 1 0 0 0 0.5  1.0  317 3.56  44  1636.0  84.00 111 394  9.8 3
278 1420 0 2 11872 1 0 0 0 0.0  5.6  338 3.70 130  2139.0 185.00 193 215  9.9 4
279 1433 0 2 20510 1 0 0 0 0.0  0.5  198 3.77  38   911.0  57.00  56 280  9.8 2
280 1412 0 1 16858 1 0 0 0 0.0  1.6  325 3.69  69  2583.0 142.00 140 284  9.6 3
281   41 2 1 24064 1 1 0 0 1.0 17.9  175 2.10 220   705.0 338.00 229  62 12.9 4
282 1455 0 2 12398 1 0 1 0 0.0  1.3  304 3.52  97  1622.0  71.00 169 255  9.5 4
283 1030 0 2 22960 1 0 0 0 0.0  1.1  412 3.99 103  1293.0  91.00 113 422  9.6 4
284 1418 0 2 17738 1 0 0 0 0.0  1.3  291 3.44  75  1082.0  85.00 195 251  9.5 3
285 1401 0 1 16929 1 0 0 0 0.0  0.8  253 3.48  65   688.0  57.00  80 252 10.0 1
286 1408 0 1 14191 1 0 1 1 0.0  2.0  310 3.36  70  1257.0 122.00 118 143  9.8 3
287 1234 0 1 21421 1 0 0 1 0.0  6.4  373 3.46 155  1768.0 120.00 151 258 10.1 4
288 1067 1 2 17874 1 0 1 0 0.5  8.7  310 3.89 107   637.0 117.00 242 298  9.6 2
289  799 2 1 24681 0 0 1 0 0.5  4.0  416 3.99 177   960.0  86.00 242 269  9.8 2
290 1363 0 1 24101 1 0 0 0 0.0  1.4  294 3.57  33   722.0  93.00  69 283  9.8 3
291  901 1 1 14939 1 0 0 0 0.0  3.2  339 3.18 123  3336.0 205.00  84 304  9.9 4
292 1329 0 2 18352 0 0 1 0 0.0  8.6  546 3.73  84  1070.0 127.00 153 291 11.2 3
293 1320 0 2 20891 1 0 1 1 1.0  8.5  194 2.98 196   815.0 163.00  78 122 12.3 4
294 1302 0 1 22111 0 0 1 0 0.0  6.6 1000 3.07  88  3150.0 193.00 133 299 10.9 4
295  877 1 1 12912 0 0 0 0 0.0  2.4  646 3.83 102   855.0 127.00 194 306 10.3 3
296 1321 0 2 11462 1 0 0 0 0.0  0.8  328 3.31  62  1105.0 137.00  95 293 10.9 4
297  533 1 1 20449 0 0 1 0 0.0  1.2  275 3.43 100  1142.0  75.00  91 217 11.3 4
298 1300 0 2 19258 1 0 1 0 0.0  1.1  340 3.37  73   289.0  97.00  93 243 10.2 3
299 1293 0 1 13913 1 0 0 0 0.0  2.4  342 3.76  90  1653.0 150.00 127 213 10.8 3
300  207 2 2 21247 1 0 1 0 0.0  5.2    . 2.23 234   601.0 135.00   . 206 12.3 4
301 1295 0 2 16513 1 0 0 0 0.0  1.0  393 3.57  50  1307.0  74.00 103 295 10.5 4
302 1271 0 1 13806 1 0 0 0 0.0  0.7  335 3.95  43   657.0  52.00 104 268 10.6 2
303 1250 0 2 22156 1 0 1 1 0.0  1.0  372 3.25 108  1190.0 140.00  55 248 10.6 4
304 1230 0 1 12979 1 0 0 0 0.0  0.5  219 3.93  22   663.0  45.00  75 246 10.8 3
305 1216 0 2 15730 1 0 1 1 0.0  2.9  426 3.61  73  5184.0 288.00 144 275 10.6 3
306 1216 0 2 20597 1 0 1 0 0.0  0.6  239 3.45  31  1072.0  55.00  64 227 10.7 2
307 1149 0 2 11167 1 0 0 0 0.0  0.8  273 3.56  52  1282.0 130.00  59 344 10.5 2
308 1153 0 1 22347 1 0 1 0 0.0  0.4  246 3.58  24   797.0  91.00 113 288 10.4 2
309  994 0 2 21294 1 0 0 0 0.0  0.4  260 2.75  41  1166.0  70.00  82 231 10.8 2
310  939 0 1 22767 1 0 0 0 0.0  1.7  434 3.35  39  1713.0 171.00 100 234 10.2 2
311  839 0 1 13879 1 0 0 0 0.0  2.0  247 3.16  69  1050.0 117.00  88 335 10.5 2
312  788 0 2 12109 1 0 0 1 0.0  6.4  576 3.79 186  2115.0 136.00 149 200 10.8 2
313 4062 0 . 21915 1 . . . 0.0  0.7    . 3.65   .      .     .     . 378 11.0 .
314 3561 2 . 23741 1 . . . 0.5  1.4    . 3.04   .      .     .     . 331 12.1 4
315 2844 0 . 19724 1 . . . 0.0  0.7    . 4.03   .      .     .     . 226  9.8 4
316 2071 2 . 27394 1 . . . 0.5  0.7    . 3.96   .      .     .     .   . 11.3 4
317 3030 0 . 22646 1 . . . 0.0  0.8    . 2.48   .      .     .     . 273 10.0 .
318 1680 0 . 15706 1 . . . 0.0  0.7    . 3.68   .      .     .     . 306  9.5 2

解决方法

任何以“.”开头的列在过滤操作将保留在它从文本文件读取时拥有的字符或因子类之前。您需要在缩减的数据集上运行它:

surv.df[] <- lapply( surv.df,function(x) as.numeric(as.character(x)))

它将任何因子或字符分类的列转换为数字类。您可以在 na.strings="." 操作中指定 read 参数以避免对字符进行隐式和静默强制,或者您可以指定 colClasses="numeric"

这是经过测试的代码:

names(surv.df)[ 2:3] <- c("times","state") # better to stay in df

surv.cox = coxph(Surv(times,state == 1) ~ V4 + V5 + V6 + V7 + V8 + V9 + V10 + V11 + V12 + V13 + V14 + V15 + V16 + V17 + V18 + V19 + V20,data = surv.df)

#------------------------------------------


surv.cox
Call:
coxph(formula = Surv(time,state == 1) ~ V4 + V5 + V6 + V7 + 
    V8 + V9 + V10 + V11 + V12 + V13 + V14 + V15 + V16 + V17 + 
    V18 + V19 + V20,data = surv.df)

          coef  exp(coef)   se(coef)      z        p
V4  -5.421e-01  5.815e-01  5.616e-01 -0.965 0.334359
V5  -3.486e-04  9.997e-01  9.199e-05 -3.790 0.000151
V6  -1.151e+00  3.162e-01  9.316e-01 -1.236 0.216510
V7  -1.918e+01  4.658e-09  5.592e+03 -0.003 0.997263
V8   9.737e-01  2.648e+00  6.580e-01  1.480 0.138930
V9  -8.052e-01  4.470e-01  6.920e-01 -1.164 0.244590
V10  1.596e+00  4.935e+00  1.470e+00  1.086 0.277434
V11  5.874e-02  1.061e+00  1.027e-01  0.572 0.567217
V12  1.202e-03  1.001e+00  1.290e-03  0.931 0.351631
V13 -1.289e+00  2.757e-01  9.781e-01 -1.317 0.187686
V14  1.310e-03  1.001e+00  2.874e-03  0.456 0.648608
V15 -7.055e-04  9.993e-01  4.298e-04 -1.641 0.100732
V16  6.424e-06  1.000e+00  5.255e-03  0.001 0.999025
V17  4.114e-03  1.004e+00  4.627e-03  0.889 0.373941
V18  5.841e-03  1.006e+00  3.255e-03  1.795 0.072718
V19 -4.326e-01  6.488e-01  4.629e-01 -0.935 0.350031
V20  1.125e+00  3.081e+00  4.656e-01  2.416 0.015672

Likelihood ratio test=39.63  on 17 df,p=0.001458
n= 276,number of events= 18 
   (42 observations deleted due to missingness)

其中一些列可能真的应该转换为因子。这可以在读取文件时使用更复杂的 colClasses 规范完成,也可以在之后使用

完成
fac_vector <- # a character vector with names of columns to be coerced 
surv.df[ fac_vector ] <- lapply (surv.df[fac_vector],factor)

我还注意到您的“状态”向量至少有 3 个级别,因此这可能意味着您需要做进一步的工作来定义事件标记真正记录的关于现实的内容。

,

试试 complete.cases() 函数。来自描述:“返回一个逻辑向量,指示哪些案例是完整的,即没有缺失值。

编辑:我没有注意到缺失值由“.”表示。为了使其正常工作,您需要用正确的 NA 值替换句点,如 IRTMF 的帖子中所示。

surv.df[] <- lapply( surv.df,function(x) as.numeric(as.character(x)))

那么 complete.cases() 函数应该会按预期工作。

surv.df[complete.cases(surv.df),]

library(tidyverse)
surv.df %>%
    filter(complete.cases(.))

.括号中是指传入的数据集。

,

你可以这样做:

surv.df <- surv.df %>%
   dplyr::filter_all(dplyr::all_vars(. != '.'))

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。