<div dir="ltr"><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763">Saudações amigos do R, </div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763">Estou as voltas de estimar taxas de eventos e estou batendo cabeça. Antes de escrever uma função eu mesmo para fazer uma tabela com alguns valores gostaria de uma luz dos amigos de R.</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763">O banco aml está no pacote survival então bastaria carregar o pacote pra reproduzir o exemplo. Eu gostaria de ver em formato de tabela alguns momentos especÃficos da tabela de sobrevivência. Só que o summary.survfit só retorna a sobrevivência e não a taxa cumulativa. </div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763"><br></div><span class="gmail_default" style="color:rgb(7,55,99)"><font face="monospace" style="" size="1">library("survival")</font></span><div><span class="gmail_default" style="color:rgb(7,55,99)"><font face="monospace" style="" size="1"></font></span> <font face="monospace" size="1">fit <span class="gmail_default" style="color:rgb(7,55,99)"></span><- survfit(Surv(time, status) ~ x, data = aml) <br>> fit<br>Call: survfit(formula = Surv(time, status) ~ x, data = aml)<br><br>         n events median 0.95LCL 0.95UCL<br>x=Maintained   11    7   31    18    NA<br>x=Nonmaintained 12   11   23    8    NA<br><span class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:rgb(7,55,99)"># Summary com todos os momentos de eventos</span><br></font></div><div><font face="monospace" size="1">> summary(fit)<br>Call: survfit(formula = Surv(time, status) ~ x, data = aml)<br><br>        x=Maintained <br> time n.risk n.event survival std.err lower 95% CI upper 95% CI<br>  9   11    1   0.909  0.0867    0.7541     1.000<br>  13   10    1   0.818  0.1163    0.6192     1.000<br>  18    8    1   0.716  0.1397    0.4884     1.000<br>  23    7    1   0.614  0.1526    0.3769     0.999<br>  31    5    1   0.491  0.1642    0.2549     0.946<br>  34    4    1   0.368  0.1627    0.1549     0.875<br>  48    2    1   0.184  0.1535    0.0359     0.944<br><br>        x=Nonmaintained <br> time n.risk n.event survival std.err lower 95% CI upper 95% CI<br>  5   12    2  0.8333  0.1076    0.6470     1.000<br>  8   10    2  0.6667  0.1361    0.4468     0.995<br>  12    8    1  0.5833  0.1423    0.3616     0.941<br>  23    6    1  0.4861  0.1481    0.2675     0.883<br>  27    5    1  0.3889  0.1470    0.1854     0.816<br>  30    4    1  0.2917  0.1387    0.1148     0.741<br>  33    3    1  0.1944  0.1219    0.0569     0.664<br>  43    2    1  0.0972  0.0919    0.0153     0.620<br>  45    1    1  0.0000   NaN      NA      NA<br><br><span class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:rgb(7,55,99)"># Summary com os momentos desejados</span><br>> <span class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:rgb(7,55,99)"></span>summary(fit, times = c(14,28,35))<br>Call: survfit(formula = Surv(time, status) ~ x, data = aml)<br><br>        x=Maintained <br> time n.risk n.event survival std.err lower 95% CI upper 95% CI<br>  14    8    2   0.818  0.116     0.619     1.000<br>  28    6    2   0.614  0.153     0.377     0.999<br>  35    3    2   0.368  0.163     0.155     0.875<br><br>        x=Nonmaintained <br> time n.risk n.event survival std.err lower 95% CI upper 95% CI<br>  14    7    5   0.583  0.142    0.3616     0.941<br>  28    4    2   0.389  0.147    0.1854     0.816<br>  35    2    2   0.194  0.122    0.0569     0.664<br><br></font><span style="font-family:monospace;font-size:x-small">> plot(fit)</span><font face="monospace" size="1"><br></font></div><div><font face="monospace" size="1">> plot(fit, cumhaz = T) </font><br><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763"><span style="font-family:Arial,Helvetica,sans-serif;font-size:small;color:rgb(34,34,34)">></span> </div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763">Reparem que há uma opção para o gráfico cumhaz = T, o que significa que o cumhaz está depositado no objeto, inclusive dentro do summary também. Tipo summary(fit)$cumhaz. Só que não há uma opção summary(fit, cumhaz = T) que retorne o cumhaz ao invés da sobrevivência. Alguém tem algum bizu pra fazer isso organizado, extrair a mesma tabela só que com o cumhaz, como <span class="gmail_default"></span><span style="color:rgb(34,34,34);font-family:monospace;font-size:x-small">summary(fit, times = c(14,28,35)) </span>sem muito trabalho?  </div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:#073763">Abraço forte,</div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><font face="'courier new', monospace" style="background-color:rgb(255,255,255)" color="#000066" size="4">Pedro Brasil</font></div></div></div></div></div></div></div></div></div></div>