search for: flexsurv

Dear all, I am fitting a parametric regression model to survival data using the flexsurvreg function from the flexsurv package. I am using a Gompertz distribution (a 2-parameter distribution) to describe the hazard function and I want to compare two groups. The model is very simple: flexsurvreg(formula = Surv(day, censored) ~ group, dist = "gompertz") and the output is also...

Dear R I am trying to understand and use the flexible parametric survival model suggested by Royston and Parmar. However I am stuck trying to plot the adjusted survival curves for different covariates in the following code: library(flexsurv) library(graphics) spl <- flexsurvspline(Surv(futime, fustat) ~ rx+ecog.ps+resid.ds+age, data = ovarian, k=2, scale="odds") spl the code "plot(spl)" produces a red spline which is difficult for me to interpret. I would like to plot survival based on for example age or rx....

...ives or general internet. Is it possible to implement the gompertz distribution as survreg.distribution to use with survreg of the survival library? I haven't found anything and recent attempts from my side weren't succefull so far. I know that other packages like 'eha' and 'flexsurv' offer functions similar to survreg with gompertz support. However, due to the run-time environment were this needs to be running in the end, I can't use these packages :( Same questions for the gompertz-makeham distribution. Many thanks! Matthias

Dear all, I am using R 4.4.0 via RStudio (2024.04.0) on a Windows PC. The code below worked on the previous version of the flexsurv but is not working since my recent update (version 2.3). The code that was working is: library(survival) rfs <- pmax(rotterdam$recur, rotterdam$death) rfstime <- with(rotterdam, ifelse(recur==1, rtime, dtime)) fit1 <- flexsurvspline(Surv(rfstime, rfs) ~ age + meno + size + er, k=2, da...

I'm working on a survival analysis project and, based on previous work of my colleagues on similar projects, am trying to estimate an AFT model with a generalized gamma distribution. When I execute the following code: >flexsurvreg(timesurv ~ cycles + rcycle + cyctime, dist="gengamma") I get the following error message: Error in optim(optpars, minusloglik.flexsurv, t = Y[, "time"], dead = Y[, : non-finite finite-difference value [4] I get an almost identical error (a 3 in the brackets instead of a...

...vctrs cli fansi hms modeltools promises rstudioapi viridisLite clipr farver htmltools mstate ps rvest webshot coda.base flexmix htmlwidgets muhaz purrr scales whisker colorspace flexsurv httr munsell qpdf selectr withr crayon flextable isoband mvtnorm quadprog spatial xfun crosstalk forcats jsonlite nlme R6 stringi xml2 curl foreign kableExtra officer RColorB...