Почему выживание дает две кривые возврата?
Я пытаюсь построить кривую выживаемости, и происходит что-то необычное, чего не было для других наборов данных. Когда я пытаюсь построить кривую, я получаю две кривые с CI каждая. Обычно у меня есть следующая сводка данных:
> summary(fit1)
Call: survfit(formula = Surv(Day, Status) ~ 1, data = d1)
time n.risk n.event survival std.err lower 95% CI upper 95% CI
1 1940 16 0.992 0.00205 0.988 0.996
2 1924 5 0.989 0.00235 0.985 0.994
3 1919 5 0.987 0.00261 0.981 0.992
4 1914 3 0.985 0.00276 0.980 0.990
5 1911 1 0.985 0.00280 0.979 0.990
6 1910 2 0.984 0.00289 0.978 0.989
8 1908 4 0.981 0.00306 0.975 0.987
11 1904 3 0.980 0.00319 0.974 0.986
14 1901 1 0.979 0.00323 0.973 0.986
18 1900 1 0.979 0.00327 0.972 0.985
Для другого набора данных я получаю две кривые, и сводка выглядит совсем иначе. Как я могу выбрать для построения только одну из кривых? Я использую автоплот для их построения. Вот краткое изложение:
> summary(fit10)
Call: survfit(formula = Surv(Day, Status) ~ Symptomatic, data = d2)
2 observations deleted due to missingness
Symptomatic=N
time n.risk n.event P((s0)) P(1)
1 985 2 0.998 0.00203
2 983 2 0.996 0.00406
3 981 3 0.993 0.00711
4 978 2 0.991 0.00914
6 976 1 0.990 0.01015
8 975 3 0.987 0.01320
11 972 2 0.985 0.01523
Symptomatic=Y
time n.risk n.event P((s0)) P(1)
1 955 14 0.985 0.0147
2 941 3 0.982 0.0178
3 938 2 0.980 0.0199
4 936 1 0.979 0.0209
5 935 1 0.978 0.0220
6 933 1 0.977 0.0230
8 932 1 0.976 0.0241
11 931 1 0.975 0.0251
14 930 1 0.974 0.0262
18 928 1 0.973 0.0272
Почему я не могу использовать ggsurvplot для второго набора данных (d2)? Вот сообщение об ошибке:
> ggsurvplot(fit10)
Error in data.frame(..., check.names = FALSE) :
arguments imply differing number of rows: 26, 0, 52