1. reasonsforreferral versus locationoffistula
## 1 2 3 p.value
## 1 12 22 7 0.01
## 2 0 23 12 0.02
## 3 10 28 12 0.32
## 4 0 11 6 0.23
## 5 2 19 5 0.36
## 6 1 7 6 0.49
## 7 3 14 8 0.95
## 8 5 7 5 0.12
## 9 3 22 13 0.47
2. The distribution of number of lesions in each fistulogram
## 1 2 3 p.value
## 0 0 4 4 0.39
## 1 19 93 49 0.40
## 2 11 44 21 0.96
## 3 6 10 0 0.00
## 4 0 2 0 1.00
4. Angioplasty survival for each lesion:
## 25 50 75
## lesion.group=l12l13l14 91 191 NA
## lesion.group=l2l4 182 NA NA
## lesion.group=l5l6l7l8 94 252 NA
## P-value of three groups comparison is 0.00088
## Call: survfit(formula = Surv(surv.days, as.numeric(event)) ~ lesion.group,
## data = subset(dat.surv.lesion, !is.na(lesion.group)))
##
## lesion.group=l12l13l14
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 5 38 1 0.974 0.0260 0.924 1.000
## 41 37 2 0.921 0.0437 0.839 1.000
## 54 35 1 0.895 0.0498 0.802 0.998
## 66 34 1 0.868 0.0548 0.767 0.983
## 75 33 1 0.842 0.0592 0.734 0.966
## 84 32 1 0.816 0.0629 0.701 0.949
## 89 31 2 0.763 0.0690 0.639 0.911
## 91 29 1 0.737 0.0714 0.609 0.891
## 98 28 1 0.711 0.0736 0.580 0.870
## 100 27 1 0.684 0.0754 0.551 0.849
## 105 26 1 0.658 0.0770 0.523 0.827
## 111 25 1 0.632 0.0783 0.495 0.805
## 112 24 2 0.579 0.0801 0.441 0.759
## 138 22 1 0.553 0.0807 0.415 0.736
## 182 21 1 0.526 0.0810 0.389 0.712
## 187 20 1 0.500 0.0811 0.364 0.687
## 195 19 1 0.474 0.0810 0.339 0.662
## 205 18 1 0.447 0.0807 0.314 0.637
## 215 16 1 0.419 0.0803 0.288 0.610
## 235 13 1 0.387 0.0804 0.258 0.582
## 272 11 1 0.352 0.0804 0.225 0.551
## 436 6 1 0.293 0.0858 0.165 0.520
##
## lesion.group=l2l4
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 8 60 1 0.983 0.0165 0.951 1.000
## 15 59 1 0.967 0.0232 0.922 1.000
## 72 58 2 0.933 0.0322 0.872 0.999
## 75 56 1 0.917 0.0357 0.849 0.989
## 91 55 1 0.900 0.0387 0.827 0.979
## 93 54 1 0.883 0.0414 0.806 0.968
## 94 53 1 0.867 0.0439 0.785 0.957
## 98 51 1 0.850 0.0462 0.764 0.945
## 99 50 1 0.833 0.0483 0.743 0.933
## 105 49 1 0.816 0.0502 0.723 0.920
## 112 48 1 0.799 0.0520 0.703 0.907
## 154 47 1 0.782 0.0536 0.683 0.894
## 178 44 1 0.764 0.0552 0.663 0.880
## 182 43 1 0.746 0.0567 0.643 0.866
## 213 38 1 0.727 0.0585 0.620 0.851
## 255 35 1 0.706 0.0604 0.597 0.835
##
## lesion.group=l5l6l7l8
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 2 67 1 0.985 0.0148 0.956 1.000
## 15 66 1 0.970 0.0208 0.930 1.000
## 28 65 1 0.955 0.0253 0.907 1.000
## 36 64 1 0.940 0.0289 0.885 0.999
## 41 63 2 0.910 0.0349 0.845 0.981
## 43 61 2 0.881 0.0396 0.806 0.962
## 54 59 2 0.851 0.0435 0.770 0.940
## 70 57 2 0.821 0.0468 0.734 0.918
## 75 55 1 0.806 0.0483 0.717 0.906
## 76 54 1 0.791 0.0497 0.699 0.895
## 85 53 1 0.776 0.0509 0.682 0.883
## 89 52 1 0.761 0.0521 0.666 0.870
## 94 51 1 0.746 0.0532 0.649 0.858
## 98 50 2 0.716 0.0551 0.616 0.833
## 100 48 1 0.701 0.0559 0.600 0.820
## 105 47 1 0.687 0.0567 0.584 0.807
## 111 46 1 0.672 0.0574 0.568 0.794
## 112 45 1 0.657 0.0580 0.552 0.781
## 114 44 1 0.642 0.0586 0.537 0.768
## 123 43 1 0.627 0.0591 0.521 0.754
## 154 41 1 0.612 0.0596 0.505 0.740
## 168 40 1 0.596 0.0600 0.490 0.726
## 170 37 1 0.580 0.0605 0.473 0.712
## 192 35 1 0.564 0.0610 0.456 0.697
## 197 34 1 0.547 0.0614 0.439 0.682
## 201 33 1 0.530 0.0618 0.422 0.666
## 216 30 1 0.513 0.0622 0.404 0.650
## 252 26 2 0.473 0.0634 0.364 0.615
## 267 21 1 0.451 0.0642 0.341 0.596
## 322 20 1 0.428 0.0648 0.318 0.576
## P-value of l5+l6+l7+l8 vs. l12+l13+l14 is 0.40984
## P-value of l5+l6+l7+l8 vs. l2+l4 is 0.00351
## P-value of l2+l4 vs. l12+l13+l14 is 0.00067
## 25 50 75
## 100 300 NA
## Call: survfit(formula = Surv(surv.days, as.numeric(event)) ~ 1, data = dat.surv.lesion)
##
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 2 185 1 0.995 0.00539 0.984 1.000
## 5 184 1 0.989 0.00760 0.974 1.000
## 8 183 1 0.984 0.00929 0.966 1.000
## 12 182 1 0.978 0.01069 0.958 1.000
## 15 181 2 0.968 0.01302 0.942 0.993
## 18 179 1 0.962 0.01403 0.935 0.990
## 28 177 2 0.951 0.01584 0.921 0.983
## 36 175 1 0.946 0.01665 0.914 0.979
## 41 174 4 0.924 0.01950 0.887 0.963
## 42 170 1 0.919 0.02013 0.880 0.959
## 43 169 2 0.908 0.02131 0.867 0.951
## 54 167 3 0.891 0.02291 0.848 0.938
## 66 164 1 0.886 0.02341 0.841 0.933
## 70 163 2 0.875 0.02435 0.829 0.924
## 72 161 3 0.859 0.02565 0.810 0.911
## 75 158 3 0.843 0.02684 0.792 0.897
## 76 155 1 0.837 0.02721 0.785 0.892
## 84 154 1 0.832 0.02757 0.779 0.888
## 85 153 1 0.826 0.02792 0.773 0.883
## 89 152 3 0.810 0.02891 0.755 0.869
## 91 149 2 0.799 0.02953 0.743 0.859
## 93 147 1 0.794 0.02982 0.737 0.854
## 94 146 2 0.783 0.03039 0.725 0.845
## 98 143 4 0.761 0.03145 0.702 0.825
## 99 139 1 0.755 0.03170 0.696 0.820
## 100 138 2 0.744 0.03217 0.684 0.810
## 105 136 3 0.728 0.03283 0.666 0.795
## 109 133 1 0.723 0.03303 0.661 0.790
## 111 132 2 0.712 0.03343 0.649 0.780
## 112 130 6 0.679 0.03447 0.614 0.750
## 114 124 1 0.673 0.03462 0.609 0.745
## 123 123 1 0.668 0.03477 0.603 0.740
## 138 121 1 0.662 0.03492 0.597 0.734
## 144 120 1 0.657 0.03506 0.592 0.729
## 154 119 2 0.646 0.03533 0.580 0.719
## 168 117 1 0.640 0.03546 0.574 0.714
## 170 114 1 0.635 0.03559 0.569 0.708
## 178 110 1 0.629 0.03573 0.563 0.703
## 182 109 2 0.617 0.03599 0.551 0.692
## 185 107 1 0.612 0.03612 0.545 0.687
## 187 105 1 0.606 0.03624 0.539 0.681
## 192 104 1 0.600 0.03636 0.533 0.676
## 195 103 1 0.594 0.03647 0.527 0.670
## 197 101 1 0.588 0.03658 0.521 0.664
## 201 100 2 0.576 0.03678 0.509 0.653
## 205 97 1 0.570 0.03688 0.503 0.648
## 213 92 1 0.564 0.03699 0.496 0.642
## 215 91 1 0.558 0.03710 0.490 0.636
## 216 90 1 0.552 0.03720 0.484 0.630
## 226 88 1 0.546 0.03731 0.477 0.624
## 235 85 1 0.539 0.03742 0.471 0.618
## 252 79 2 0.526 0.03769 0.457 0.605
## 255 77 1 0.519 0.03782 0.450 0.598
## 267 73 1 0.512 0.03796 0.442 0.592
## 272 72 1 0.504 0.03809 0.435 0.585
## 300 68 1 0.497 0.03825 0.427 0.578
## 322 63 1 0.489 0.03845 0.419 0.571
## 401 36 1 0.476 0.03971 0.404 0.560
## 436 30 1 0.460 0.04143 0.385 0.549
## P-value of l2 vs. l4 is 0.40984