2010-01-24 3 views
1

J'ai utilisé les données, aml dans le paquet de survie dans R et calculé une fonction de survie par "survfit". Étant donné que le résultat de survfit ne montre pas la valeur moyenne, j'ai utilisé le code suivant pour imprimer la moyenne:Je ne peux pas imprimer la valeur moyenne de survit (...)

> print(survfit(Surv(aml1$time,aml1$status)~1),show.rmean=T) 

(les données j'est aml1 <-aml[aml$x="Maintained",])

Le code ci-dessus travaillé dans mon ami PC, mais pas le mien. Donc, j'ai pensé à télécharger un paquet supplémentaire pour utiliser print(...show.rmean=T). Mais, "print" est basique, je n'ai donc besoin de rien pour imprimer. Alors pourquoi je ne peux pas obtenir la valeur moyenne?

+0

Mettez à jour votre survie de paquet ... cela a fonctionné pour moi! –

Répondre

3

print est une fonction générique, avec des options spécifiques à la chose à imprimer. Dans ce cas, les options sont fournies par print.survfit. Bien que ma version de survival n'ait pas de show.rmean. Peut-être que vous voulez print.rmean? Voir l'aide pour print.survfit pour plus d'informations.

> print(survfit(Surv(aml1$time,aml1$status)~1), print.rmean=T) 
Call: survfit(formula = Surv(aml1$time, aml1$status) ~ 1) 

    records  n.max n.start  events  *rmean *se(rmean)  median 
    23.00  23.00  23.00  18.00  36.36  9.85  27.00 
    0.95LCL 0.95UCL 
    18.00  45.00 
    * restricted mean with upper limit = 161 
Questions connexes