<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">Hola,</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">Mira el ejemplo de la ayuda:</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><div class="gmail_default">> library(survival)</div><div class="gmail_default">> fit2 <- coxph( Surv(stop, event) ~ size, data = bladder )</div><div class="gmail_default">> # single curve</div><div class="gmail_default">> ggadjustedcurves(fit2, data = bladder)</div><div class="gmail_default"><br></div><div class="gmail_default">Que produce este gráfico.</div><div class="gmail_default"><br></div><div class="gmail_default"><img src="cid:ii_161399419636bfba" alt="Imágenes integradas 1" width="562" height="246"><br></div><div class="gmail_default"><br></div><div class="gmail_default">Solamente tienes que pasar la variable donde guardas el modelo y los datos.</div><div class="gmail_default"><br></div><div class="gmail_default">Pero puedes valorar otras opciones, representar las curvas de supervivencia, estratificadas por otra variable (los grupos, etc, etc).</div><div class="gmail_default"><br></div><div class="gmail_default">Saludos,</div><div class="gmail_default">Carlos Ortega</div><div class="gmail_default"><a href="http://www.qualityexcellence.es">www.qualityexcellence.es</a></div><div class="gmail_default"><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">El 27 de enero de 2018, 15:50, Patricio Suárez Gil <span dir="ltr"><<a href="mailto:patricsg@gmail.com" target="_blank">patricsg@gmail.com</a>></span> escribió:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Tengo un modelo de regresión de Cox y quiero obtener el plot ajustado por una covariable (sexo) con la función ‘ggadjustedcurves’, pero me da el siguiente error:<br>
<br>
> cox2 <- coxph(os ~ imc_25 + sexo.1, data = datos)<br>
> cox2<br>
Call:<br>
coxph(formula = os ~ imc_25 + sexo.1, data = datos)<br>
<br>
coef exp(coef) se(coef) z p<br>
imc_25 -0.621 0.537 0.299 -2.08 0.038<br>
sexo.1M 0.714 2.042 0.387 1.84 0.065<br>
<br>
Likelihood ratio test=6.23 on 2 df, p=0.0444<br>
n= 76, number of events= 54<br>
> ggadjustedcurves(cox2, variable = datos$sexo.1, data = datos)<br>
Error in `[.data.frame`(data, , variable) : undefined columns selected<br>
<br>
Agradezco cualquier feedback.<br>
<br>
Saludos,<br>
<br>
<br>
Patricio<br>
<br>
<br>
Patricio Suárez Gil<br>
Unidad de Investigación Área V-Gijón<br>
Planta 5ª Impar<br>
Hospital Universitario de Cabueñes<br>
C/Prado, 395<br>
33394 Gijón (Asturias)<br>
Tfno: 985 185 000 (Ext. 85715)<br>
@uinvest_psg<br>
<a href="mailto:unidadinvestigacion.area5@sespa.es">unidadinvestigacion.area5@<wbr>sespa.es</a><br>
ESPAÑA<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
[[alternative HTML version deleted]]<br>
<br>
______________________________<wbr>_________________<br>
R-help-es mailing list<br>
<a href="mailto:R-help-es@r-project.org">R-help-es@r-project.org</a><br>
<a href="https://stat.ethz.ch/mailman/listinfo/r-help-es" rel="noreferrer" target="_blank">https://stat.ethz.ch/mailman/<wbr>listinfo/r-help-es</a></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><span style="font-family:verdana,sans-serif">Saludos,</span><br style="font-family:verdana,sans-serif">
<span style="font-family:verdana,sans-serif">Carlos Ortega</span><br style="font-family:verdana,sans-serif">
<span style="font-family:verdana,sans-serif"><a href="http://www.qualityexcellence.es" target="_blank">www.qualityexcellence.es</a></span></div>
</div>