similar to: Royston Parmar adjusted survival curves using flexsurv

Hi Dear R-users I have a database with 18000 observations and 20 variables. I am running cox regression on five variables and trying to use survfit to plot the survival based on a specific variable without success. Lets say I have the following coxph: >library(survival) >fit <- coxph(Surv(futime, fustat) ~ age + rx, data = ovarian) >fit what I am trying to do is plot a survival

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) ~

Hi I have been looking around for a couple of days without being able to find a solution. I am trying to plot the following table using the package ggplot2 Month Mean SEM cry3m 3 -0.05227273 0.1607422 cry6m 6 0.65826087 0.6638734 cry1y 12 0.38333333 0.3062907 cry2y 24 1.04811111 0.4921863 cry3y

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 quite

System: R 2.3.1 on a Windows XP computer. I am validating several cancer prognostic models that have been published with a large independent dataset. Some of the models report a probability of survival at a specified timepoint, usually at 5 and 10 years. Others report only the linear predictor of the Cox model. I have used Harrell's c index for censored data (rcorr.cens) as a measure of

Hi I am a new fan of R after getting mad with the graphical functional in SPSS. I have been able to create a nice looking Kaplan Meyer graph using Survplot function. However I have difficulties in turning the y axis to percent instead of the default 0-1 scale. Further I have tried the function yaxt="n" without any results. Any help in this matter will be appreciated. The code is

To whom it may concern, I'm Utsav Parmar, pursuing my B. Tech in Computer Engineering. I like to work on new technologies and am currently looking for open-source projects to contribute to. As it may turn out, I've got a college project in my curriculum this semester under “Software Development Practice”, and I'd like to work upon a project and/or a feature in pipeline spanning over

Dear R-list readers: Royston described the effect of sample size on the p-value obtained from the Shapiro-Francia test (Estimating departure from normality. Stat Med 1991;10:1283-93). He developed two indices from the Shapiro-Francia test (i) V' - an index of departure from normality and (ii) v' - a plot of the cumulative squared residuals. He mentioned the availability (in

Hi, I was reading a paper published in JCO "Prediction of risk of distant recurrence using 21-gene recurrence score in node-negative and node-positive postmenopausal patients with breast cancer treated with anastrozole or tamoxifen: a TransATAC study" (ICO 2010 28: 1829). The author uses a method to estimate the 9-year risk of distant recurrence as a function of continuous recurrence

Hi all, Sorry if this has been answered already, but I couldn't find it in the archives 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

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[,

Dear Sir, Thanks for all good effort . Recently we have installed Red Hat Linus 8.0 ,but when use samba in browser its not working.error unable to connect to Host:901. then we have try to find SWAT file and Dir but can not find any file .Please let us know how to add /install SWAT successfully. Samba ver=2.2.5 Rgds NP

hello everyone, i am getting problems in graph plotting. When i attach file after adding color attributes in my data set. i got problem of "GlobalEnv" and masked the followings. Like this >attach(machm) The following object(s) are masked _by_ '.GlobalEnv': coll, sp The following object(s) are masked from 'mach': angle, area, dis, plot, sp

Is there a way to tell glm() that rows in the data represent a certain number of observations other than one? Perhaps even fractional values? Using the weights argument has no effect on the standard errors. Compare the following; is there a way to get the first and last models to produce the same results? data(sleep) coef(summary(glm(extra ~ group, data=sleep))) coef(summary(glm(extra ~ group,

Hi, After loading R 4.0.0 on Ubuntu 20.04. R version 4.0.0 (2020-04-24) -- "Arbor Day" Copyright (C) 2020 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()'

In the next month or two I've got to upgrade a number of servers that are currently on an EOL'd version of 4-STABLE. I foresee that I'll have very limited time to do full OS upgrades on these systems in the coming several years, so I want to make sure I bring them onto an extended-life branch. Right now 4.11 has the furthest projected EOL date (Jan 31 2007), and the projected EOL

Howyd all? Seems that I have been routed. Possibly by a physical B&E, but who knows? Probably some of you do.... anyways, some politically sensitive email was deleted from a user account and the line low -tr & inserted into my .xinitrc . Duncan (Dhu) Campbell

Hi, Does anyone know if the D'agostino test is available with R ? Alex

Brokers Move On ERMX! EntreMetrix Inc. (ERMX) $0.18 Heavy trading today as ERMX announced its launch of digital support tools for its portfolio companies. Brokers are getting ahead of this steady climb as they grab up large blocks of shares for there clients. Look at the numbers and get on ERMX Friday morning! KNEEguru: Are there any contraindications for the procedure? if the laws of physics

Brokers Move On ERMX! EntreMetrix Inc. (ERMX) $0.18 Heavy trading today as ERMX announced its launch of digital support tools for its portfolio companies. Brokers are getting ahead of this steady climb as they grab up large blocks of shares for there clients. Look at the numbers and get on ERMX Friday morning! KNEEguru: Are there any contraindications for the procedure? if the laws of physics