[R] interval censoring case 1 or current status data

[Chao Zhu]

Dear all,

Is there a function in R dealing with the NPMLE for current status data or 
interval censored case 1 data?

Thanks,

Jimmy
[[alternative HTML version deleted]]

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] Cox model qustion in R

[Chao Zhu]

- Original Message - 
From: Thomas Lumley [EMAIL PROTECTED]
To: Chao Zhu [EMAIL PROTECTED]
Cc: r-help@stat.math.ethz.ch
Sent: Thursday, March 31, 2005 11:09 AM
Subject: Re: [R] Cox model qustion in R


 On Wed, 30 Mar 2005, Chao Zhu wrote:

 But it only gives me the curve with combined treatment. How do I get a 
 plot with two curves for different treatments. I don't want to do 
 new-data.frame(treat=c(0,1)) and plot(survfit(fit,newdata=new))
 cause  it gives you predicted curves for two treatments.

 
 I think you are going to have to be clearer about what you want.  If you 
 want the Kaplan--Meier estimator for each treatment separately then
   plot(survfit(Surv(t)~treat))
 will do it. Otherwise I don't know what you're trying to do.
 
 By the way, do you really have no censoring in the time variable?
 
  -thomas
 


__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] Cox model qustion in R

[Chao Zhu]

Dear all,

I have a beginner's R question for coxph() and survfit().
I have a data set with survival time (no censoring) and treatment 
indicator(0/1). I did something like

fit-coxph( Surv(t)~treat, method=breslow,data=data)
 fit
Call:
coxph(formula = Surv(t) ~ treat, data = data, method = breslow)


   coef exp(coef) se(coef) z   p
treat -1.96 0.1400.243 -8.07 6.7e-16

Likelihood ratio test=80.3  on 1 df, p=0  n= 144 

 plot(survfit(fit))

But it only gives me the curve with combined treatment. How do I get a plot 
with two curves for different treatments.
I don't want to do 
new-data.frame(treat=c(0,1))
and plot(survfit(fit,newdata=new)) cause it gives you predicted curves for two 
treatments.

Thanks you so much for your help.

Jimmy Chu
[[alternative HTML version deleted]]

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] a cox model question

[Chao Zhu]

Hi,

Is there anyone can tell me how to estimate baseline cumulative hazard function 
and baseline survival function 
by using Breslow's method in R?
Thanks,

Jimmy
[[alternative HTML version deleted]]

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] batch jobs question

[Chao Zhu]

Hi,there

I'm doing some R batch jobs in Unix. 
Something like
R prog1 output1 --save 
R prog2 output2 --save 

prog1 and prog2 are running at the same time and they are essentially same 
except it contains different parameter values. 
I was wondering if two processes will affect each other? Hopefully they are two 
independent jobs.

It's a beginner's question and thanks for your help.

Chao

[[alternative HTML version deleted]]

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html