Skip to content

survival analysis

3 messages · p.b.pynsent@bham.ac.uk, Jun Yan

Original
p.b.pynsent@bham.ac.uk
# 
Has anybody written an actuarial (life) survival procedure, this does 
not appear to be an option in the survival package?
This approach is common in orthopaedic surgery to demonstrate the 
survival of prostheses.
I need to apply the "modified" lower conf.int because of the censoring 
over time.
I want a life table which I can then easily plot.
Many thanks

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Jun Yan
# 
I am not sure if this is what you want. The function lifetab in package
KMsurv creates a cohort life table as described in Klein and Moeschberger
(1997), Survival Analysis.
Jun
On Wed, 20 Nov 2002, Paul Pynsent wrote:

            

      
      
      
    

        
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

  
  


  


      

    

    

      
      

        


  

        

      p.b.pynsent@bham.ac.uk
    

    

      
      #
    

  


  

      
Jun,
Thank you for this information. It is half of what I wanted as there  
are no confidence limits calculated.
Survfit in the survival package gives very useful alternatives for  
calculating the confidence limits including adjustments to the lower  
limits, for data that is heavily censored (see Dorey and Korn as  
referenced in the survfit documentation). On the other hand,  
unfortunately survfit does not offer an actuarial method.
It looks as though I am going to have to write my own code for this.

        
On Wednesday, November 20, 2002, at 03:19  pm, Jun Yan wrote:

        

            

      
      
      
    

        
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._