Hi all, I have a little problem. I make an weibull survival analysis using the survival package. It,s OK, them I have the functions. I plot this funcions with curve(). I want to make a plot with the real survival points (proportion of alive x time) and them add the curves to points. I have the time to dead, the censor data and my trataments. To analysis the model is: model1 <- survreg(Surv(time,censor)~trat) Exist any function in R to plot the real point (proportion of alive in time)? Thanks Ronaldo -- That's odd. That's very odd. Wouldn't you say that's very odd? -- | //|\\ [*****************************][*******************] || ( ? ? ) [Ronaldo Reis J?nior ][PentiumIII-600 ] | V [ESALQ/USP-Entomologia, CP-09 ][HD: 30 + 10 Gb ] || / l \ [13418-900 Piracicaba - SP ][RAM: 128 Mb ] | /(lin)\ [Fone: 19-429-4199 r.229 ][Video: SiS620-8Mb ] ||/(linux)\ [chrysopa at insecta.ufv.br ][Modem: Pctel-onboar] |/ (linux) \[ICQ#: 5692561 ][SO: CL 7.0 (2.2.19)] || ( x ) [*****************************][*******************] ||| _/ \_Powered by Conectiva Linux 7.0 D+:) | Lxuser#: 205366 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
On Thu, 2 May 2002, Ronaldo Reis Jr. wrote:> Hi all, > I have a little problem. > I make an weibull survival analysis using the survival package. It,s OK, them > I have the functions. I plot this funcions with curve(). I want to make a > plot with the real survival points (proportion of alive x time) and them add > the curves to points. I have the time to dead, the censor data and my > trataments. To analysis the model is: > model1 <- survreg(Surv(time,censor)~trat) > Exist any function in R to plot the real point (proportion of alive in time)?plot(survfit(Surv(time,censor)~trt)) or to superimpose on an existing plot lines(survfit(Surv(time,censor)~trt)) -thomas Thomas Lumley Asst. Professor, Biostatistics tlumley at u.washington.edu University of Washington, Seattle -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Em Qui 02 Mai 2002 14:42, Thomas Lumley escreveu:> > plot(survfit(Surv(time,censor)~trt)) > > or to superimpose on an existing plot > > lines(survfit(Surv(time,censor)~trt)) > > -thomas > > > Thomas Lumley Asst. Professor, Biostatistics > tlumley at u.washington.edu University of Washington, SeattleHi Thomas I try it, but the graphic make is not a point graphic, is a ?scale? plot, like this: |_ | |_ | | | |__ | |_ | |_ |_____________ I need the same graphic, but with point like this: | | * | * | * | * | * |____________ I try to change the type (plot(survfit(Surv(tempo,censor)~trat),type="p")), but dont work:(((. I try to use points(survfit(Surv(tempo,censor)~trat), but this only work over the plot graphics with ?scales lines? I try to plot using type="n" but dont work too. Whem I use type the error appear: Error in plot.default(tempx, tempy * yscale, type = "n", log = logax, : formal argument "type" matched by multiple actual arguments Exist a mode for save the data to be plotted, saving all points calculated by the function survfit(Surv(tempo,censor)~trat in a object? Thanks for all Inte mais Ronaldo -- Maternity pay? Now every Tom, Dick and Harry will get pregnant. -- Malcolm Smith -- | //|\\ [*****************************][*******************] || ( ? ? ) [Ronaldo Reis J?nior ][PentiumIII-600 ] | V [ESALQ/USP-Entomologia, CP-09 ][HD: 30 + 10 Gb ] || / l \ [13418-900 Piracicaba - SP ][RAM: 128 Mb ] | /(lin)\ [Fone: 19-429-4199 r.229 ][Video: SiS620-8Mb ] ||/(linux)\ [chrysopa at insecta.ufv.br ][Modem: Pctel-onboar] |/ (linux) \[ICQ#: 5692561 ][SO: CL 7.0 (2.2.19)] || ( x ) [*****************************][*******************] ||| _/ \_Powered by Conectiva Linux 7.0 D+:) | Lxuser#: 205366 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._