topic Re: how to use proc power to calculate power given alpha, hazard ratio and number of events in Statistical Procedures

how to use proc power to calculate power given alpha, hazard ratio and number of events

fengyuwuzu — Thu, 23 Jan 2020 06:52:45 GMT

In survival analysis, given alpha, hazard ratio and number of events, assuming no Interim Analysis, how to calculate the power?

Re: how to use proc power to calculate power given alpha, hazard ratio and number of events

Rick_SAS — Thu, 23 Jan 2020 14:00:36 GMT

I believe you can use the COXREG statement in PROC POWER to get the power.

proc power;
   coxreg
      alpha = 0.1
      hazardratio = 1.4
      stddev = 1.0
      ntotal = 80 to 120 by 10
      eventprob = 0.8
      power = .
   ;
run;

Re: how to use proc power to calculate power given alpha, hazard ratio and number of events

fengyuwuzu — Thu, 23 Jan 2020 15:31:02 GMT

Thank you very much Rick. I will look into it.

I just did it in R, using the following code

> z<-0.5*sqrt(n)*log(1/hr)-qnorm(1-alpha/2)
> power<-pnorm(z)

Re: how to use proc power to calculate power given alpha, hazard ratio and number of events

Rick_SAS — Fri, 24 Jan 2020 14:10:03 GMT

If you want to use an explicit formula, you can run the DATA step or PROC IML to get the same results as R:

data Power;
n = 20;
hr=0.3;
alpha=0.05;

z = 0.5*sqrt(n)*log(1/hr)-quantile("Normal", 1-alpha/2);
power = cdf("Normal", z);
run;

proc print; run;

Re: how to use proc power to calculate power given alpha, hazard ratio and number of events

fengyuwuzu — Fri, 24 Jan 2020 15:44:10 GMT

Thank you very much Rick for providing code for this option.

I have a question on stddev in proc power with coxreg: how do I decide what number for stddev to use?

stddev = 1.0

Re: how to use proc power to calculate power given alpha, hazard ratio and number of events

Rick_SAS — Fri, 24 Jan 2020 15:58:03 GMT

The parameters in the equation are usually based on pilot studies, previous research, or your best guess.

Re: how to use proc power to calculate power given alpha, hazard ratio and number of events

haoduonge — Fri, 23 Sep 2022 18:03:17 GMT

Hi all,

A follow-up question to the [proc power]: can it integrate correlation coefficient in multilevel data?

Thanks!

Hao

Re: how to use proc power to calculate power given alpha, hazard ratio and number of events

TomHsiung — Sun, 21 Apr 2024 06:18:36 GMT

Hello, Rick. I am not sure about the parameter stddev. Is it the standard error of a coefficient estimate from a PROC PHREG procedure, or something different? Thanks.

Re: how to use proc power to calculate power given alpha, hazard ratio and number of events

John_Castelloe — Mon, 22 Apr 2024 19:53:30 GMT

The STDDEV= option in the COXREG statement is the assumed standard deviation of the predictor being tested---that is, your best guess of the true standard deviation of this predictor's sampling distribution under the same conditions as expected in the study you're planning. This isn't the same as the standard error of a coefficient estimate, just a simple property of the predictor variable itself as it would be sampled (or perhaps taking on values that you control yourself) in your planned study.

I'd also like to point out that the original question in this thread didn't specify what kind of data analysis is being planned, and it could just as easily pertain to the TWOSAMPLESURVIVAL statement in PROC POWER if you're planning a logrank test rather than a Cox proportional hazards regression. These are two scenarios in which PROC POWER supports the calculation of power for a survival analysis given alpha, hazard ratio, and number of events.