If I remember correctly, I noticed an effect that did give a p of slightly less than .05 was a hazard ratio of 3, which made me think of running that test, and then I think spower was the r function that I used to figure out what p they could get for a hazard ratio of 2 and 35 experimentals and 35 controls (or whatever the actual split was- I think it was slightly different?).
So you were using Hmisc::spower… I’m surprised that there was even such a function (however obtusely named) - why on earth isn’t it in the survival library?
I was going to try to replicate that estimate, but looking at the spower documentation, it’s pretty complex and I don’t think I could do it without the original paper (which is more work than I want to do).
If I remember correctly, I noticed an effect that did give a p of slightly less than .05 was a hazard ratio of 3, which made me think of running that test, and then I think spower was the r function that I used to figure out what p they could get for a hazard ratio of 2 and 35 experimentals and 35 controls (or whatever the actual split was- I think it was slightly different?).
So you were using
Hmisc::spower
… I’m surprised that there was even such a function (however obtusely named) - why on earth isn’t it in thesurvival
library?I was going to try to replicate that estimate, but looking at the spower documentation, it’s pretty complex and I don’t think I could do it without the original paper (which is more work than I want to do).