Ch13-mysurvfit: My Survfit

In PwrGSD: Power in a Group Sequential Design

My Survfit

Description

Computes numbers at risk, numbers of events at each unique event time within levels of a blocking factor

Usage

  mysurvfit(formula = formula(data), data = parent.frame(), subset, na.action = na.fail)

Arguments

formula	Should be a formula of the form Surv(ti, ev) ~ block where block is the blocking factor. It need not be a factor per se but should have relatively few discrete levels. Sorry, no staggered entry allowed at present
data	a dataframe
subset	you can subset the analysis via logical expression in variables in the dataframe
na.action	pass a method for handling NA values in block such as na.omit, or na.fail

Value

A dataframe of 2*NLEV + 1 columns where NLEV is the number of levels of the factor block.

time	The sorted vector of unique event times from all blocks
nrisk1	The number at risk in block level 1 at each event time
nevent1	The number of events in block level 1 at each event time
...
nriskNLEV	The number at risk in block level NLEV at each event time
neventNLEV	The number of events in block level NLEV at each event time

Author(s)

Grant Izmirlian <[e-mail]>

Examples

  library(PwrGSD)
  data(lung)

  fit.msf <- mysurvfit(Surv(time, I(status==2)) ~ sex, data=lung)

  fit.msf
  ## Not run: 
  plot(fit.msf)
  
## End(Not run)

PwrGSD documentation built on Nov. 2, 2024, 1:11 a.m.