Does Surv in psm handle interval2 data?  The argument list seems to indicate it does but I get an error.

Thanks,
Chris

# code
library('survival')
left <- c(1, 3, 5, NA)
right <-c(2, 3, NA, 4)
Surv(left, right, type='interval2')
survreg(Surv(left, right, type='interval2') ~ 1)

library('rms')
Surv(left, right, type='interval2') # error
args(Surv)
psm(Surv(left, right, type='interval2') ~ 1) # same error (of course)
psm(survival::Surv(left, right, type='interval2') ~ 1) # runs


# output

R version 2.15.2 (2012-10-26) -- "Trick or Treat"

Copyright (C) 2012 The R Foundation for Statistical Computing

ISBN 3-900051-07-0

Platform: x86_64-w64-mingw32/x64 (64-bit)



R is free software and comes with ABSOLUTELY NO WARRANTY.

You are welcome to redistribute it under certain conditions.

Type 'license()' or 'licence()' for distribution details.



  Natural language support but running in an English locale



R is a collaborative project with many contributors.

Type 'contributors()' for more information and

'citation()' on how to cite R or R packages in publications.



Type 'demo()' for some demos, 'help()' for on-line help, or

'help.start()' for an HTML browser interface to help.

Type 'q()' to quit R.



> library('survival')

Loading required package: splines

> left <- c(1, 3, 5, NA)

> right <-c(2, 3, NA, 4)

> Surv(left, right, type='interval2')

[1] [1, 2] 3      5+     4-

> survreg(Surv(left, right, type='interval2') ~ 1)

Call:

survreg(formula = Surv(left, right, type = "interval2") ~ 1)



Coefficients:

(Intercept)

   1.317943



Scale= 0.6098782



Loglik(model)= -5.3   Loglik(intercept only)= -5.3

n= 4

> library('rms')

Loading required package: Hmisc

Hmisc library by Frank E Harrell Jr



Type library(help='Hmisc'), ?Overview, or ?Hmisc.Overview')

to see overall documentation.



NOTE:Hmisc no longer redefines [.factor to drop unused levels when

subsetting.  To get the old behavior of Hmisc type dropUnusedLevels().





Attaching package: 'Hmisc'



The following object(s) are masked from 'package:survival':



    untangle.specials



The following object(s) are masked from 'package:base':



    format.pval, round.POSIXt, trunc.POSIXt, units





Attaching package: 'rms'



The following object(s) are masked from 'package:survival':



    Surv



> Surv(left, right, type='interval2') # error

Error in Surv(left, right, type = "interval2") :

  argument "event" is missing, with no default

> args(Surv)

function (time, time2, event, type = c("right", "left", "interval",

    "counting", "interval2"), origin = 0)

NULL

> psm(Surv(left, right, type='interval2') ~ 1) # same error (of course)

Error in Surv(left, right, type = "interval2") :

  argument "event" is missing, with no default

> psm(survival::Surv(left, right, type='interval2') ~ 1) # runs



Parametric Survival Model: Weibull Distribution



psm(formula = survival::Surv(left, right, type = "interval2") ~

    1)



                Model Likelihood    Discrimination

                   Ratio Test          Indexes

Obs        4    LR chi2     0.00    R2       0.000

Events     6    d.f.           0    g        0.000

sigma 0.6099                        gr       1.000



            Coef    S.E.   Wald Z Pr(>|Z|)

(Intercept)  1.3179 0.3598  3.66  0.0002

Log(scale)  -0.4945 0.5977 -0.83  0.4081



**********************************************************
Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues 

	[[alternative HTML version deleted]]

