Hi I am trying to built a stacked bar chart that shows the frequency percentage of a variable, stacked to 100% and further grouped in categories This example might explain what I am trying to do http://www.cdc.gov/nchs/data/databriefs/db182_fig1.png So there are three levels Year x (population, discharge, days of care ) And each bar shows age categories stacked to 100% I can arrange the bars as a stack or group, but cant use both in the same chart Any suggestions please? ... View more

Hi I know this forum is mainly for help with the SAS code itself and I have received some fantastic answers here, I would really appreciate if anyone could please help me with this. I am doing a regression analysis, looking at the effect of 2 drugs (A or B) in a population where we stratified individuals by ethnic group: Estimates for European Population proc logistic data=Studydatabase; class sex drug study; model Efficacy (event='Good')=   drug age sex weight egfr; where population='European'; run;   Effect OddsRatioEst LowerCL UpperCL ProbChiSq Age 1.005 0.973 1.037 0.7729 Drug A vs B 1.862 1.229 2.821 0.0033 Sex Female vs Male 1.469 0.925 2.334 0.1033 Weight 1.015 1.001 1.028 0.0346 eGFR 1.006 0.995 1.016 0.2967 ________________________________________ Estimates for Asian Population proc logistic data=Studydatabase; class sex drug study; model Efficacy (event='Good')=   drug age sex weight egfr; where population='Asian'; run; Obs Effect OddsRatioEst LowerCL UpperCL ProbChiSq 1 Age 1.047 0.994 1.103 0.0839 2 Drug A vs B 0.719 0.217 2.382 0.5892 3 Sex Female vs Male 1.932 0.445 8.384 0.3795 4 Weight 0.987 0.954 1.020 0.4357 5 eGFR 1.014 0.991 1.036 0.2300 As you can notice, the influence of  drug is significant in the European, but not in the Asian population. One could interpret this as the effect of the drug being modified according to the study subpopulation, and hence expect the p value for interaction to be significant, but it is not at all! proc logistic data=Studydatabase; class sex drug study; model effectnum (event='Good')=  study drug age sex weight egfr population*drug; run; Interactin term P value ProbChiSq 0.1874 I know this forum is mainly for questions regarding the SAS code itself, but I would appreciate any help from all the SAS and statistics experts here Best wishes Am Any explanation please? Am I missing something ... View more

Hi all, I am using a proc glimmix to perform a poison regression, and using the "store" statement to store the output parameters proc glimmix data=Poisson_sim; class prog ; model num_awards = prog math / dist=poisson  s cl; store p1; run; Then using the following to estimate the IRR, and lower and upper CI ods output ParameterEstimates = est; proc plm source = p1;   show parameters; run; data est_exp (drop=stderr llci luci estimate );   set est;   IRR = exp(estimate);     llci=estimate - 1.96*stderr;   luci=estimate + 1.96*stderr;   UpperCI=exp (luci);   LowerCI= exp (llci); run; proc print data = est_exp; run; The problem is: I am not managing to get the p-values from the proc Glimmix to be displayed in the output Could anyone please advice me on how to do it? Best wishes ... View more