This function returns a ggplot derived from an rpsftm object showing psi vs test statistic. In case of convergence should see a unique solution.

ggpsi(x)

Arguments

x

rpsftm object (rpsftm)

Examples

# use data included in rpsftm package
library(rpsftm)
immdef <- rpsftm::immdef
immdef$rx <- with(immdef, 1 - xoyrs/progyrs)

# fit the model
imm.fit <- rpsftm(Surv(progyrs, prog) ~ rand(imm, rx),
                  data=immdef,
                  censor_time=censyrs)

ggpsi(imm.fit)