<!DOCTYPE html>
<html>
<head>
<title></title>
</head>
<body><div style="font-family:Arial;">Estou sem o R no momento. Talvez isto funcione:<br></div>
<div style="font-family:Arial;"> </div>
<div style="font-family:Arial;">[...]<br></div>
<div style="font-family:Arial;">f <- cph(S ~ rcs(age,4) + sex, x=TRUE, y=TRUE)<br></div>
<div style="font-family:Arial;">df <- data.frame(coef = coef(f), row.names = names(coef(f)))<br></div>
<div style="font-family:Arial;">write.csv(df, "meuarquivo.csv")<br></div>
<div style="font-family:Arial;"> </div>
<div style="font-family:Arial;">Espero ter ajudado!</div>
<div style="font-family:Arial;"> </div>
<div id="sig4487139"><div id="signature" class="signature"><a href="http://lattes.cnpq.br/9234772336296638" title="Currículo Lattes">Leonardo Ferreira Fontenelle</a><br></div>
</div>
<div> </div>
<div> </div>
<div>Em Qua 20 abr. 2016, às 17:36, Pedro Emmanuel Alvarenga Americano do Brasil escreveu:<br></div>
<blockquote type="cite"><div dir="ltr"><div style="font-family:'courier new', monospace;color:rgb(0, 0, 102);">Amigos de R,<br></div>
<div style="font-family:'courier new', monospace;color:rgb(0, 0, 102);"> </div>
<div style="font-family:'courier new', monospace;color:rgb(0, 0, 102);">Ajustei um modelo de sobrevivencia e gostaria de exportar o modelo uma planilha xlsx ou csv, de tal forma que os coeficientes fiquem numa coluna e na coluna ao lado fique pelo o que o coeficiente está multiplicando. Já quebrei a cabeça aqui e não consegui. A pegadinha é que o modelo possui algumas transformações por splines e por isso preciso dos valores dos nós pra cada variável contínua. Achei uma função chamada Function que faz isso. Mas depois disso eu não avancei. Eu gostaria de pegar da saída função Function o que está entre as chaves, quebrar nos "+" e nos "-" pra separar os termos e depois quebrar nos "*" pra separar as colunas. Organizar em formato de planilha e exportar.<br></div>
<div style="font-family:'courier new', monospace;color:rgb(0, 0, 102);"> </div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">library(rms)</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">n <- 1000</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">age <- 50 + 12*rnorm(n)</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">label(age) <- "Age"</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">sex <- factor(sample(c('Male','Female'), n, </span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'"> rep=TRUE, prob=c(.6, .4)))</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">cens <- 15*runif(n)</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">h <- .02*exp(.04*(age-50)+.8*(sex=='Female'))</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">dt <- -log(runif(n))/h</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">label(dt) <- 'Follow-up Time'</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">e <- ifelse(dt <= cens,1,0)</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">dt <- pmin(dt, cens)</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">units(dt) <- "Year"</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">dd <- datadist(age, sex)</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">options(datadist='dd')</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">S <- Surv(dt,e)</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">f <- cph(S ~ rcs(age,4) + sex, x=TRUE, y=TRUE)</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', ' monospace'">>Function(f)</span></span><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace"> </span></span><br></div>
<div><div style="font-family:Arial;"> </div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">>function(age = 48.800654,sex = "Male") {0.28612987-0.017321951* age+0.00012586581*pmax(age-29.600456,0)^3-0.00046285671*pmax(age-44.273198,0)^3+0.00042111763*pmax(age-53.769262,0)^3-8.4126731e-05*pmax(age-69.855597,0)^3-0.64448916*(sex=="Male") }</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace"><environment: 0x0000000015fa7378></span></span><br></div>
<div> </div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">Pra ficar mais ou menos assim</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">+0.28612987</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">-0.017321951 age</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">+0.00012586581 pmax(age-29.600456,0)^3</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">-0.00046285671 pmax(age-44.273198,0)^3</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">+0.00042111763 pmax(age-53.769262,0)^3</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">-8.4126731e-05 pmax(age-69.855597,0)^3</span></span><br></div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">-0.64448916 (sex=="Male") </span></span><br></div>
<div> </div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">Alguma luz?</span></span><br></div>
<div> </div>
<div style="font-family:Arial;"> </div>
</div>
<div><span class="colour" style="color:rgb(0, 0, 102)"><span class="font" style="font-family:'courier new', monospace">Abraço forte,</span></span><br></div>
<div style="font-family:'courier new', monospace;color:rgb(0, 0, 102);"> </div>
<div><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><span class="highlight" style="background-color:rgb(255, 255, 255)">Pedro Brasil</span><br></div>
</div>
</div>
</div>
</div>
</div>
</div>
<div><u>_______________________________________________</u><br></div>
<div>R-br mailing list<br></div>
<div><a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br></div>
<div><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br></div>
<div>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></div>
</blockquote><div style="font-family:Arial;"> </div>
</body>
</html>