Hi All,
I use the stepAIC to select the variables of Cox ph model. But I got one
error as the following:
Start: AIC= 239.32
Surv(SURVP, STATP) ~ AGEBMT_med + SEX + DXCODE + STAGEGRP + BSYMDX +
BMINVDX + PS + IPIGROUP + INITXGRP + INICHEMO + PRICHEMO +
PRIRESP + DXBMTMON + CON1_3
Error in dropterm.default(fit, scope$drop, scale = scale, trace = max(0,
:
number of rows in use has changed: remove missing values?
Could you please help me with it?
Here is my code:
> whole<-read.table("E://my file/BC project/Project II/UVA and
MVA/whole.txt",header=TRUE,sep="\t")
>
full1<-coxph(Surv(SURVP,STATP)~AGEBMT_med+SEX+DXCODE+STAGEGRP+BSYMDX+BMI
NVDX+PS+IPIGROUP+INITXGRP+INICHEMO+PRICHEMO+PRIRESP+DXBMTMON+CON1_3,data
=whole)
> full<-stepAIC(full1)
> full
Thank you!
P.S.: The attached is the data.
UPN STATP SURVP STATNRMP STATREP STATEFP SURVEFP AGEBMT_med
AGEBMT_dec SEX DXCODE DX1 DXBMT? STAGEGRP BSYMDX BMINVDX
PS IPIGROUP INITXGRP PRIXRT INICHEMO PRIPURAN
PRIRITU PRICHOP PRICHEMO PRIRESP DXBMTMON CON1_3 CON1_2 STATUS
318 1 41 1 0 1 41 0 2 F
2 2 COMP ADV 0 0 0 1 2 1
1 0 0 1 2 P 1 3 0 1REL/S
492 1 2078 1 0 1 2078 0 2 M
2 2 DIS ADV 1 1 0 2 2 1
1 0 0 1 1 P 0 2 1 1PR
625 1 128 1 0 1 128 0 2 M
1 1 TRAN-H ADV 1 1 0 1 1 0
1 0 0 1 3 P 0 2 1 1REL/S
631 1 1548 0 1 1 656 0 1 F
2 2 DIS ADV 0 1 0 2 1 0
1 0 0 1 3 P 0 1 1 1REL/S
637 1 16 1 0 1 16 0 1 F
2 2 COMP ADV 0 1 0 2 1 0
1 0 0 1 1 C 1 2 1 1REL/S
652 1 17 1 0 1 17 1 2 M
1 1 TRAN-H ADV 0 1 0 2 2 1
1 0 0 1 2 C 0 2 1 1REL/S
759 1 156 0 1 1 66 0 2 M
1 1 TRAN-H ADV 0 0 0 1 2 1
1 0 0 1 3 C 0 1 1 1REL/S
826 1 52 1 0 1 52 0 1 F
1 1 TRAN-H ADV 0 1 0 1 2 1
1 0 0 1 2 P 0 1 1 1REL/S
918 1 289 0 1 1 121 1 3 M
2 2 COMP ADV 1 0 0 1 2 1
1 0 0 1 4 P 1 2 1 1REL/S
985 1 166 0 1 1 152 1 2 M
1 1 TRAN-C ADV 0 0 0 1 1 0
1 0 0 1 2 P 0 1 1 IF
1033 1 103 0 1 1 85 0 1 M
1 2 TRAN-C ADV 0 1 0 2 1 0
3 1 0 0 2 N 0 3 0 IF
1222 1 967 0 1 1 799 0 1 F
1 1 TRAN-C ADV 0 1 0 1 2 1
1 0 0 1 4 P 0 3 0 2REL/S
1225 1 639 0 1 1 182 0 1 M
1 1 TRAN-H ADV 0 1 0 2 1 1
3 1 0 0 1 P 0 1 1 1REL/U
1228 0 2575 0 0 0 2575 1 2 F
2 2 COMP ADV 0 0 0 1 2 1
1 0 0 1 2 C 1 2 1 1REL/U
1319 1 106 1 0 1 106 1 3 M
1 1 TRAN-H LIM 0 0 0 1 2 1
1 0 0 1 2 P 0 1 1 2CR
1321 0 2622 0 1 1 962 0 0 M
2 2 COMP ADV 1 0 0 1 2 1
1 0 0 1 2 P 0 2 1 1PR
1370 1 68 0 1 1 61 1 3 F
2 2 DIS ADV 1 1 0 2 1 0
1 0 0 1 1 P 0 2 1 IF
1468 0 2212 0 0 0 2212 1 2 F
2 2 DIS ADV 0 1 0 1 1 0
1 0 0 1 1 C 0 1 1 1PR
1528 1 456 0 1 1 245 0 0 F
1 2 TRAN-H ADV 0 0 0 2 2 1
1 0 0 1 3 C 0 1 1 2REL/R
1634 1 2074 0 1 1 1324 0 2 F
1 1 TRAN-C ADV 1 0 1 2 1 1
1 0 0 1 2 P 0 1 1 1REL/R
1687 1 1461 0 1 1 769 0 1 M
1 1 TRAN-H ADV 1 1 0 2 1 0
1 1 0 1 3 P 0 1 1 1REL/S
1810 1 173 0 1 1 104 1 2 M
1 1 TRAN-H ADV 0 0 0 1 1 0
1 0 0 1 4 C 0 1 1 1REL/S
1824 1 133 0 1 1 92 0 0 M
1 2 TRAN-C ADV 0 1 0 2 1 0
1 1 0 1 3 P 0 1 1 IF
1867 1 285 1 0 1 285 1 2 M
2 2 DIS ADV 1 1 1 3 1 0
1 0 0 1 2 P 1 1 1 2CR
1889 0 1525 0 0 0 1525 0 1 M
1 1 TRAN-H ADV 0 0 0 2 1 0
3 1 1 0 4 P 1 1 1 1REL/R
1893 0 1529 0 0 0 1529 0 1 M
2 2 DIS ADV 0 1 0 2 2 1
1 1 0 1 3 C 1 1 1 1REL/S
1895 0 1009 0 0 0 1009 0 1 M
1 1 TRAN-H LIM 0 1 0 1 1 0
3 1 1 0 4 P 0 1 1 1REL
1917 1 334 1 0 1 334 1 2 M
1 1 TRAN-C ADV 0 1 0 2 1 0
1 1 0 1 3 P 0 1 1 3PR
1958 1 38 1 0 1 38 0 1 F
2 2 COMP ADV 0 1 0 2 1 0
1 0 0 1 2 C 1 1 1 1REL/S
1972 1 197 1 0 1 197 1 2 M
1 1 TRAN-H ADV 0 1 0 2 2 0
1 1 0 1 3 G 1 1 1 1REL/S
1979 0 1275 0 0 0 1275 0 2 M
1 1 TRAN-C ADV 0 1 0 1 1 0
3 1 0 0 2 P 0 1 1 1REL
2009 1 223 0 1 1 82 1 3 F
2 2 COMP ADV 0 1 0 1 1 0
1 0 0 1 3 C 1 3 0 1REL/S
2039 1 124 0 1 1 1 1 3 M
1 1 TRAN-C ADV 1 0 0 2 1 0
1 1 0 1 3 P 0 3 0 1REL/S
2061 1 91 1 0 1 91 1 3 M
1 2 TRAN-H ADV 1 1 0 1 1 0
1 1 1 1 3 N 0 1 1 1PR
2114 0 754 0 0 0 754 1 3 M
1 1 TRAN-H ADV 1 1 0 2 2 1
1 0 1 1 3 P 0 1 1 1PR
2115 0 332 0 0 0 332 1 2 F
2 2 COMP ADV 0 1 0 2 1 0
1 0 1 1 3 P 1 1 1 1REL/S
2117 1 98 0 1 1 85 0 0 M
2 2 DIS ADV 0 1 0 2 1 0
2 0 0 0 1 P 0 1 1 1PR
2119 1 41 1 0 1 41 1 3 M
2 2 COMP ADV 0 1 0 1 2 1
1 0 1 1 2 C 1 1 1 1REL
2135 1 512 0 1 1 174 1 3 F
2 2 COMP LIM 0 0 0 1 1 0
1 0 1 1 2 G 1 3 0 2CR
2141 1 679 1 0 1 679 1 2 M
1 1 TRAN-H ADV 0 0 0 1 1 0
1 0 1 1 4 G 0 1 1 1PR
2159 0 722 0 0 0 722 1 3 F
1 1 TRAN-H LIM 0 0 0 1 2 1
1 0 1 1 6 C 1 3 0 3CR
2167 0 733 0 0 0 733 0 2 F
2 2 DIS ADV 1 1 0 2 1 0
1 1 1 1 4 C 1 1 1 2REL/S
2173 0 643 0 0 0 643 0 1 F
2 2 COMP ADV 1 1 0 2 1 0
1 0 1 1 4 P 1 1 1 1REL/S
2179 0 699 0 0 0 699 0 2 F
2 2 COMP ADV 0 0 0 1 2 1
1 1 0 1 3 C 1 1 1 3CR
2236 0 563 0 0 0 563 1 3 F
1 1 TRAN-C ADV 0 1 0 1 1 0
1 1 0 1 4 P 1 1 1 3REL/S
2250 0 386 0 0 0 386 1 2 M
1 1 TRAN-C ADV 1 1 0 1 1 0
1 0 1 1 3 P 0 2 1 3REL/S
2271 1 25 1 0 1 25 1 3 M
1 1 TRAN-H ADV 0 1 0 1 1 0
1 1 1 1 3 C 0 1 1 1PR
2295 0 396 0 0 0 396 1 2 M
1 1 TRAN-H LIM 0 0 0 1 2 1
1 1 1 1 5 N 0 3 0 1REL/S
2304 1 295 1 0 1 295 1 2 M
2 2 COMP ADV 1 1 0 1 1 0
3 1 1 0 2 P 0 1 1 1PR
2307 0 348 0 0 0 348 0 0 M
2 2 COMP ADV 1 0 1 2 1 0
1 0 1 1 2 C 1 3 0 1REL/S
2311 1 107 0 1 1 77 0 2 M
2 2 COMP LIM 0 0 0 1 2 1
1 0 1 1 2 C 1 1 1 1REL/S
2340 0 317 0 1 1 174 1 3 M
1 1 TRAN-H ADV 1 0 0 2 1 0
1 0 0 1 2 P 0 2 1 1REL/S
2360 0 37 0 0 0 37 0 2 M
1 2 TRAN-H LIM 0 0 0 1 2 1
1 0 0 1 3 C 1 3 0 2REL/S
2361 0 263 0 0 0 263 1 3 M
2 2 COMP ADV 0 0 1 1 0
1 1 1 1 4 G 1 3 0 IF
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
