Hi everyone, I am currently analyzing the association between several biomarkers (continuous variables) and a disease (a binary variable, 0 or 1) along with confounders using Cox proportional hazard models. Given the abundance of biomarkers, I should multiple testing with corrected p-values for each biomarker. I have obtained the estimate of the effective number of independent tests, which is 50. I understand that the significance threshold for p-values should be set at 0.05/50. Now, my question is: should I divide the calculated p-values by 50 as well? Also, if it should be is there any option to transform "<0.001" to detailed value? For example as code: proc phreg data=mydata;
class confounders;
model time*status(0)=confounders biomarker1/RL;
ods select ParameterEstimates;ods output ParameterEstimates=result;
run;
data result; set result;
corrected_p=ProbChiSq/50; /*This is my question*/
run;
/*repeat this as the number of biomarekrs */ Thanks in advance!
... View more