<div dir="ltr"><div><font face="arial, sans-serif">Isto não é suficiente?</font></div><div><font face="arial, sans-serif"><br></font></div><font face="monospace">> fit$cumhaz<br> [1] 0.09090909 0.19090909 0.31590909 0.45876623 0.45876623 0.65876623<br> [7] 0.90876623 0.90876623 1.40876623 1.40876623 0.16666667 0.36666667<br>[13] 0.49166667 0.49166667 0.65833333 0.85833333 1.10833333 1.44166667<br>[19] 1.94166667 2.94166667<br>></font><br><div><font face="arial, sans-serif">HTH</font></div><div><font face="arial, sans-serif">--</font></div><div><font face="arial, sans-serif">Cesar Rabak</font></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Dec 8, 2021 at 5:12 PM Pedro Emmanuel Alvarenga Americano do Brasil por (R-br) <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:rgb(7,55,99)">Saudações amigos do R, </div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:rgb(7,55,99)"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:rgb(7,55,99)">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:rgb(7,55,99)"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:rgb(7,55,99)">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:rgb(7,55,99)"><br></div><span class="gmail_default" style="color:rgb(7,55,99)"><font face="monospace" size="1">library("survival")</font></span><div><span class="gmail_default" style="color:rgb(7,55,99)"><font face="monospace" 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:rgb(7,55,99)"><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:rgb(7,55,99)">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:rgb(7,55,99)"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:large;color:rgb(7,55,99)">Abraço forte,</div><div><div dir="ltr"><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>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div>