That would be better to demonstrate your question by posing some data and the picture you want to see.   data have; call streaminit(123); array x{*} opt1-opt6; do id=1 to 100; do i=1 to dim(x); x{i}=rand('bern',i/10); end; output; end; drop i; run; proc transpose data=have out=want; by id; var opt:; run; proc sgplot data=want pctlevel=group; vbar _name_/group=col1 stat=percent ; run;   ... View more

I think you could use ESTIMATE to get RR and its CI ,after that polt it to get the graph you want. http://support.sas.com/kb/24/188.html ... View more

@art297 Arthur.T, Long time no see ! Did you see Perter.E in recently years ? ... View more

Did you try : ("MedDRA System Organ Class /* Preferred Term" (" " aesoc)) --> ( "09"x ("MedDRA System Organ Class /* Preferred Term" ("09"x aesoc)) ) ... View more

Yes. You can. Just list all these variable in sub-query of PROC SQL.   proc sql; create table want as select id,South ,North ,East,West, (select mean(South) from ( select South from have union all select North from have union all select East from have union all select West from have )) as mean , (select std(South) from ( select South from have union all select North from have union all select East from have union all select West from have )) as std from have; quit; ... View more

You could drop these variables firstly, after that adding these variables in this table . Here is an example:   data have; set sashelp.class; run; proc sql; alter table have drop name,age add name char(20),age num; quit; ... View more

data have; input id South North east west ; cards; 1 2 4 12 9 2 3 . 11 10 3 2 4 12 10 4 6 4 13 12 ; proc sql; create table want as select id,South ,North , (select mean(South) from (select South from have union all select North from have)) as mean_south_north , (select std(South) from (select South from have union all select North from have)) as SD_south_north, east,west, (select mean(east) from (select east from have union all select west from have)) as mean_east_west , (select std(east) from (select east from have union all select west from have)) as SD_east_west from have; quit; ... View more

Check @Rick_SAS blogs: https://blogs.sas.com/content/iml/2024/03/27/likelihood-ratio-test.html ... View more

Check @Rick_SAS blogs: https://blogs.sas.com/content/iml/2018/10/10/fit-growth-curve-sas.html ... View more

Sorry. I have no idea about it. If you want symmetric CI, you could try option " fisher(BIASADJ= no) ". But that is still not get exact symmetric CI. proc reg data=sashelp.heart(obs=100) ; weight ageatstart; model weight= height Diastolic Systolic; output out=want predicted=p; quit; ods output FisherPearsonCorr=FisherPearsonCorr; proc corr data=want outp=outp fisher(BIASADJ= no); var weight; with p; weight ageatstart; run; data FisherPearsonCorr; set FisherPearsonCorr; RSquare=Corr**2; lcl_rsquare=lcl**2; ucl_rsquare=ucl**2; _lcl2=lcl_rsquare-RSquare; _ucl2=ucl_rsquare-RSquare; run; proc print noobs;run;   Or you could try Bootstrap Method to get RSquare CI. Maybe @Rick_SAS @StatDave knew something you want. ... View more

Yeah. That is what supposed to do. Did you try @DanH_sas code ? ... View more

Did you try LSMESTIMATE ,not ESTIMATE ? ... View more