BookmarkSubscribeRSS Feed
Maartje
Calcite | Level 5

Hello all,

 

For my school project I'm developing a risk prediction model for overall mortality. The overall equation to do so is: 1 - 10yearsurival probability^(exp(regression coefficient x1.. x2.. x3.. x3). Where the 10 year survival probability has to be calculated with proc phreg.

The syntax for my model is:

   proc phreg data=a.predict;

model time*died15(0)= age cur former mod med crp ;

run;

 

I'm not sure how to calculate the 10 year surival probability.. I know I should use the BASELINE statement and that I have to calculate it for my covariates set to zero?

 

I hope you can help me out :)!

Maartje

 

 

 

5 REPLIES 5
Reeza
Super User

I don’t know if you necessarily need a BASELINE statement, does the output dataset not have the survival estimates at various points in time? You can definitely use BASELINE to specify the values you’re interested in as well. 

 

Its not 0 though, it’s usually the mean for continuous variables and the reference level for categorical variables. 

Maartje
Calcite | Level 5

Ohh thanks for your response. I didnt specify the reference levels for my categorical variables and thats why I didnt get the survival probability.

 

I still have a small question about filling in the equation for the 10-year risk prediction for mortality. Now I've calculated the survival probability with the baseline statement where the mean for continuous variables and the reference level for categorical variables was used. If I want to calculate the 10 year risk for mortality for a person of age 70 and the mean age used to calculate the baseline survival probability was 60, should i use the difference in age between the mean used for baseline survival calculations and the actual age (in this case age = 10? Or should i just fill in the actual age of this person?

 

Thanks again 🙂

Reeza
Super User

The model uses the actual age so you would use that. 

Maartje
Calcite | Level 5

Hmm, thanks and sorry to bother u again.. But my predicted 10 year mortality risks are really high (also compared to the number of events) and does not vary much in between participants and I can't seem to find out why..

Should I calculate the baseline survival probability for ALL variables available in the data set? Now Ive calculated the baseline survival probability with the predictors which are in the end model.

 

Thanks!

 

 

Reeza
Super User

@Maartje wrote:

Hmm, thanks and sorry to bother u again.. But my predicted 10 year mortality risks are really high (also compared to the number of events) and does not vary much in between participants and I can't seem to find out why..

Should I calculate the baseline survival probability for ALL variables available in the data set? Now Ive calculated the baseline survival probability with the predictors which are in the end model.

 

Thanks!

 

 


No, you first build a model with the variables that you think are relevant. Once your model is built you then would look at the 10 year survival. If you N is small, especially for the subgroup analysis you're results will be all over the place. 

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

Click image to register for webinarClick image to register for webinar

Classroom Training Available!

Select SAS Training centers are offering in-person courses. View upcoming courses for:

View all other training opportunities.

Discussion stats
  • 5 replies
  • 1849 views
  • 1 like
  • 2 in conversation