+ Reply to Thread
Results 1 to 4 of 4

Thread: How to obtain an unbiased post-stratified estimator for survival probabilities?

  1. #1
    Points: 4,825, Level: 44
    Level completed: 38%, Points required for next Level: 125

    Posts
    40
    Thanks
    8
    Thanked 0 Times in 0 Posts

    How to obtain an unbiased post-stratified estimator for survival probabilities?




    Suppose I have estimated the survival probabilities $\hat{S(t)}_1, \hat{S(t)}_2, \ldots, \hat{S(t)}_H$ from $H$ strata. Now I want to get an unbiased estimate of population survival probability $S(t)$ using these stratum-specific estimates. Then what should be the estimator?

    Should it be $\hat{S(t)}_{pst} = \sum_{i=1}^h \dfrac{N_h}{N} \times \hat{S(t)}_h$?

    If I want to adjust for under or over representation in the sample, then what should be the estimator (suppose I know the population and sample proportions)?

    Thanks for your kind guidance.

  2. #2
    Points: 3,631, Level: 37
    Level completed: 88%, Points required for next Level: 19
    staassis's Avatar
    Location
    New York
    Posts
    226
    Thanks
    2
    Thanked 41 Times in 39 Posts

    Re: How to obtain an unbiased post-stratified estimator for survival probabilities?

    No easy way to adjust. At any time horizon T, the event of survival up to T is co-dependent with the stratum membership. Not only do conditional probabilities differ from the unconditional probabilities but they change with T.

    The easiest way out is simulating a representative sample of the whole population out of your strata multiple times. For each simulated sample, run survival analysis and get the survival curve. Use the survival curves from multiple simulations to draw the confidence band for the true survival curve.

  3. #3
    Points: 2, Level: 1
    Level completed: 3%, Points required for next Level: 48

    Posts
    1
    Thanks
    0
    Thanked 0 Times in 0 Posts

    Re: How to obtain an unbiased post-stratified estimator for survival probabilities?

    Hi
    i want do stratified cox model in R :
    library(survival)
    library(foreign)
    leuk.data<-read.table("C://Users/Ehsan/Desktop/sur2.txt",header=T)
    leuk.data
    cox1<-coxph(Surv(t,status)~Rx+logWBC+strata(sex),data=leuk.data,method="breslow")
    but it shows this error :

    Error in Surv(t, status) : Time variable is not numeric
    pleas help me how can numeric time and whats your idea ?
    TNX

  4. #4
    Omega Contributor
    Points: 38,374, Level: 100
    Level completed: 0%, Points required for next Level: 0
    hlsmith's Avatar
    Location
    Not Ames, IA
    Posts
    6,998
    Thanks
    398
    Thanked 1,186 Times in 1,147 Posts

    Re: How to obtain an unbiased post-stratified estimator for survival probabilities?


    Not a regular R user, but have you looked to see how time is formatted in the dataset?


    Code: 
     is.numeric(leuk.data$t)
    Stop cowardice, ban guns!

+ Reply to Thread

           




Tags for this Thread

Posting Permissions

  • You may not post new threads
  • You may not post replies
  • You may not post attachments
  • You may not edit your posts






Advertise on Talk Stats