Questions tagged [survival]

Survival analysis compares the statistical strength of predictors for events encoded as occuring or censored in concert with a time variable. It can be parametric or non-parametric.

445 questions
1
vote
0 answers

How to perform Prametric mixture cure modelling in R? Although semi-parametric modelling package is there called smcure

How to perform Prametric mixture cure modelling in R? Although semi-parametric modelling package is there called smcure. But i need to model my data by parametric mixture cure model. Or can anyone suggest other statistical package to do the same.…
0
votes
0 answers

How to get predicted and observed risk estimates from riskRegression::FGR() subdistribution competing risk model object?

I want to extract observed and predicted estimates from time to event riskRegression::FGR() competing risk model object. Below is some simulated…
0
votes
0 answers

How to combine ggsurvplot and ggplot in a same plot

I had drawn kaplan meier plot for stanford heart transplant…
0
votes
0 answers

Coxph ran out of iteration and did not Converge

I am trying to run a Cox proportional hazard model to determine the effects of covariates on the survival of neonates. Previously when I ran coxph with numeric Predictor Variables, I had a result. # Create a survival object using age and event…
0
votes
0 answers

Dealing with variables violating proportional hazard assumption, is it okay to use both step function and time-transform function simultaneously?

I tried to fit survival model with 4 variables violating proportional hazard assumption. I read following vignette (https://cran.r-project.org/web/packages/survival/vignettes/timedep.pdf) and seared online thoroughly, but could not find similar…
JY J
  • 21
  • 3
0
votes
1 answer

How to enhance the visualisation of my survplot?

I'm conducting an experiment (during 10 days of observation) to compare four different conditions named NTS, NTC, TS, and TC. I performed a logrank test, and the differences between the conditions seem to be statistically significant. I am aware…
lddc
  • 1
  • 2
0
votes
1 answer

Facet competing risks by outcome rather than group

I have a survival dataset with outcomes for multiple groups separated into causes of death, and would like to plot the competing risks of each outcome. This tutorial provides an excellent example of using the cpmrsk and survminer packages to plot…
Timmo83
  • 119
  • 5
0
votes
0 answers

Way to handle multiple variables that violate porportional hazard assumption. Step function and time transformation

I tried to fit survival model with several variables violating proportional hazard assumption. I read following vignette thoroughly, but cannot get satisfactory answers for my…
JY J
  • 21
  • 3
0
votes
1 answer

Rescale model linear predictor in rms::cph

I'm working with the rms package to derive a score for predicting response to a drug. I actually derived two scores: one including only pre-treatment variables (x and y), and a second including also variable that can be measured only after 6 months…
Devanto
  • 33
  • 5
0
votes
1 answer

Survival Curves in R

I am new to survival curves and trying to graph them in R. I have the following columns in my excel file: Sample (Experimental or Control), Age (Old or Young), Sex (Male or Female), Bacteria status (either positive or negative) and days 0-20 (21…
m.rodwell
  • 39
  • 5
0
votes
0 answers

Negative Numbers at Risk; Survfit doubling my subject sample size; Competing Risks Model

I'm evaluating competing risks of mortality of birds due to various causes (snakes, mammals, hawks, other) with staggered entry. I have 120 individuals; 60 entered at time 0 and 60 at time 63; I ran survcheck in the R survival package. The…
bkube
  • 11
  • 2
0
votes
0 answers

Plotting a KM plot when there's a strata pairing variable

I have survival data and I'm comparing two groups. The groups are optimal matched prior to analysis, so I need to include the pairings in the KM for the log rank test. I completed the KM with the code…
John Ryan
  • 343
  • 1
  • 9
0
votes
1 answer

Extracting Standard error from crq function in r

I am using censored quantile regression for survival data. For a specific qunatile (say 0.8), I get the summary result. I can extract the intercept and beta coefficient but I can not extract standard error of the coefficients from the summary…
UNI39
  • 21
  • 2
0
votes
1 answer

How to get survival table in R?

This question concerns survival data i.e. a triple that is (enter, exit, event). My data is structured such that for each individual id, I can have event = 1 only once. I want to generate a survival table which shows for any given age, how many…
Cloft X
  • 141
  • 7
0
votes
0 answers

ggsurvplot_facet p-values - error variable lengths differ

I am trying to create a faceted survival plot with confidence intervals and p-values. The plot should look essentially equivalent to the below but with p-values. However, I keep running into the error variable lengths differ. I have found several…
Pat Taggart
  • 321
  • 1
  • 9