CureRate-class {CR} | R Documentation |
"CureRate"
Class of objects like the output of function "curerate"
.
Objects can be created by calls of the form new("CureRate", ...)
.
cureobs
:Object of class "numeric"
medobs
:Object of class "numeric"
curerx
:Object of class "numeric"
medrx
:Object of class "numeric"
actime
:Object of class "vector"
futime
:Object of class "vector"
info
:Object of class "vector"
crits
:Object of class "vector"
alpha
:Object of class "numeric"
rho
:Object of class "numeric"
acrate
:Object of class "numeric"
probrx
:Object of class "numeric"
numreps
:Object of class "integer"
numobs
:Object of class "matrix"
timept
:Object of class "array"
deaths
:Object of class "array"
testname
:Object of class "character"
power
:Object of class "array"
beta
:Object of class "matrix"
indac
:Vector of objects of class "integer"
indfu
:Vector of objects of class "integer"
printflag
:Object of class "integer"
signature(object = "CureRate")
Emil A. Cornea, Bahjat F. Qaqish, and Joseph G. Ibrahim
Maintainer: Emil A. Cornea <ecornea@live.unc.edu>
curerate
, show-methods
, showcr
showClass("CureRate") ## Input Data ############### # rho - specifies the value of rho in the G-rho test (Harrington and Fleming, 1982). # rho = 0 gives the logrank test, and rho = 1 the Peto-Peto Wilcoxon test #(and rho = -1 the test discussed by Gray and Tsiatis, 1989). rho = 0 numreps = 500 ##cureobs - probability of cure on the observation arm of the study cureobs = .359 ##curerx - probability of cure on the treatment arm of the study curerx = .459 ##medobs - median survival among the non-cured on observation medobs = .747 ##medrx - median survival among the non-cured on treatment medrx = .859 ##acrate - rate at which patients accrue, combined over both ##observation and treatment arms of the study acrate = 232 ##probrx - probability of assignment to the treatment arm probrx = .5 ##actime - accrual duration actime = c(3.3, 3.5, 3.8); ##futime - followup duration futime = c(2.0, 2.5) ##info - vector of information times for the interim looks ##must be an increasing sequence in (0,1] info = c( .35, .61, .86, 1.0) ##crits - vector of critical values corresponding to the ##interim analysis time points specified in info crits = c(3.6128, 2.6506, 2.1894, 2.0536) ############################################################## ### Log-rank test: rho=0 (default) mycr<-curerate(rho ,numreps,cureobs,curerx,medobs,medrx,acrate, probrx,actime,futime,info,crits) mycr ### (*) show(mycr) ### same as above showcr(mycr) ### same as above unclass(mycr) showcr(mycr,full.results=TRUE) showcr(mycr,indac=3,indfu=1)