Thanks for being so helpful already, hoping I can stretch it with two more questions on this topic, sorry I'm new to SAS and have a long way to go. Question 1 I'm trying to compare power of cox and logistic regression. I want to run a logistic regression on the same data however the intercept values are being included in the pvalues. I can't seem to use the 'by treatment' statement it says 'variable treatment not found'. /*Create dataset with ODS*/ ods output ParameterEstimates=logest; /*Run logistic regression specify descending to find P(event=1)*/ proc logistic data=subset2 descending; by sampleno; model event = treatment; run; /*Unsuppress ODS listing*/ ods exclude none; ods results; options notes; /*Create Power criteria*/ proc format; value signif low-0.05 = 'significant' 0.05<-high = 'not significant'; run; /*Display Power results*/ proc freq data=logest; format ProbChiSq signif.; tables ProbChiSq; run; /*calculate means of sampling dist for log regression */ proc means data=logest mean stddev; var estimate; run; Question 2 I am also trying to do the logistic regression with a glm to get a RR estimate so in place of logistic regression command I have written: /*Run log-binomial regression to look at RR*/ proc genmod data=subset descending; model event = treatment/dist=binomial link=log; run; however it seems to only run 1 sample simulation not the 10,100, or 1000 specified. AN ideas why it might be doing this? Thanks heaps, Ben
... View more