TheGofCens package include the following graphical tools and goodness-of-fit tests for right-censored data:
- Kolmogorov-Smirnov, Cramér-von Mises, and Anderson-Darling tests, which use the empirical distribution function for complete data and are extended for right-censored data.
- Generalized chi-squared-type test, which is based on the squared differences between observed and expected counts using random cells with right-censored data.
- A series of graphical tools such as probability or cumulative hazard plots to guide the decision about the most suitable parametric model for the data.
These functions share several features as they can handle both complete and right-censored data, and they provide parameter estimates for the distributions under study.
Brief Example
To conduct goodness-of-fit tests with right censored data we can use theKScens(),CvMcens(),ADcens() andchisqcens() functions. We illustrate this by means of thecolon dataset:
# Kolmogorov-Smirnovset.seed(123)KScens(Surv(time,status)~1,colon, distr="norm")# Cramér-von Misescolonsamp<-colon[sample(nrow(colon),300),]CvMcens(Surv(time,status)~1,colonsamp, distr="normal")# Anderson-DarlingADcens(Surv(time,status)~1,colonsamp, distr="normal")# Generalized chi-squared-type testchisqcens(Surv(time,status)~1,colonsamp, M=6, distr="normal")The graphical tools provide nice plots via the functionscumhazPlot(),kmPlot() andprobPlot(). See several examples using thenba data set:
License
- GPL (>= 2)