| summary.afttest | R Documentation |
Summary Method for Semiparametric AFT Model Diagnostics
## S3 method for class 'afttest'
summary(object, ...)
object |
An object of class |
... |
Other arguments passed to methods. |
summary.afttest returns a detailed summary of an afttest fit.
library(survival)
library(aftgee)
library(afttest)
datgen <- function(n = 100) {
z1 <- rbinom(n, 1, 0.5)
z2 <- rnorm(n)
e <- rnorm(n)
tt <- exp(2 + z1 + z2 + 0.5 * z2^2 + e)
cen <- runif(n, 0, 100)
data.frame(Time = pmin(tt, cen), status = 1 * (tt < cen),
z1 = z1, z2 = z2, id = 1:n)
}
set.seed(1)
simdata <- datgen(300)
# --------------------------------------------------
# Method 1: Formula (Runs quickly for CRAN tests)
# --------------------------------------------------
result_form <- afttest(Surv(Time, status) ~ z1 + z2, data = simdata,
npath = 50, testType = "covForm", estMethod = "rr",
eqType = "ns", covTested = "z2", npathsave = 50,
linApprox = TRUE, seed = 1)
print(result_form)
plot(result_form, std = TRUE)
# --------------------------------------------------
# Method 2: Fitted aftsrr object (Induced Smoothing)
# --------------------------------------------------
fit_srr <- aftsrr(Surv(Time, status) ~ z1 + z2, data = simdata,
eqType = "is", rankWeights = "gehan")
result_srr <- afttest(fit_srr, data = simdata, npath = 100, testType = "covForm",
covTested = "z2", npathsave = 50,
linApprox = TRUE, seed = 1)
summary(result_srr)
plot(result_srr, std = FALSE)
# --------------------------------------------------
# Method 3: Fitted aftgee object (Least Squares)
# --------------------------------------------------
fit_gee <- aftgee(Surv(Time, status) ~ z1 + z2, data = simdata)
result_gee <- afttest(fit_gee, data = simdata, npath = 100, testType = "covForm",
covTested = "z2", npathsave = 50,
linApprox = TRUE, seed = 1)
print(result_gee)
# --------------------------------------------------
# Method 4: Standard Resampling (linApprox = FALSE)
# --------------------------------------------------
result_resamp <- afttest(Surv(Time, status) ~ z1 + z2, data = simdata,
npath = 100, testType = "covForm", estMethod = "rr",
eqType = "ns", covTested = "z2", npathsave = 50,
linApprox = FALSE, seed = 1)
summary(result_resamp)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.