Dear All,   I would like to draw vlines with mean and standard deviation, instead of standard errors. The below code gives me the vline with StdErr. Could anybody help me with this? Thanks in advance!   * Plot ; data temp; set sashelp.cars; if origin in ("Asia" , "Europe") then group= 1; else group= 2; if type= "Hybrid" then delete; run; proc sgplot data= temp; vline Type / response= EngineSize stat= mean markers limitstat= stderr group= group; ; run; * Summary Stat ; proc sort data= temp; by type group ; run; proc means data= temp n sum mean std; class type group ; var EngineSize; run; ... View more

I would like to estimate regression with three sets of fixed effects and clustered standard errors. PROC SURVEYREG works, but it has memory issues when there are many fixed effects and/or clusters. So, I decided to use RPOC STDIZE to de-mean my variables. But, I don't think my code is proper. Could anybody help me to point out what is wrong in my code?   My code includes two PROC SURVEYREG. 1) The original data with explicit fixed effects. 2) The de-meaned data (i.e., implicit fixed effects).   I expected the results to be the same, but they are not.     proc freq data= test order=freq; table f1; run; proc freq data= test order=freq; table f2; run; proc freq data= test order=freq; table f3; run; * Centering data by f1, f2, and f3; data test_centered; set test; run; proc sort data= test_centered; by f1; run; proc stdize data= test_centered out= test_centered method= mean oprefix= old sprefix; by f1; var dv iv; run; proc sort data= test_centered; by f2; run; proc stdize data= test_centered out= test_centered method= mean oprefix= old sprefix; by f2; var dv iv; run; proc sort data= test_centered; by f3; run; proc stdize data= test_centered out= test_centered method= mean oprefix= old sprefix; by f3; var dv iv; run; proc means data= test; var dv iv; run; proc means data= test_centered; var dv iv; run; proc surveyreg data= test; class f1 f2 f3; model dv= iv f1 f2 f3/ solution; run; proc surveyreg data= test_centered; model dv= iv / solution; run; proc print data= test; run;     My Original DataSet:    Obs f1 f2 f3 dv iv 1 20202 37 2016 14.1690 4 2 20202 37 2017 14.3602 6 3 20202 37 2018 14.4787 5 4 293110377 38 2016 12.6248 3 5 293110377 38 2017 13.0541 3 6 293110377 38 2018 12.9945 3 7 2020245 28 2016 11.2252 3 8 2020245 28 2017 11.3145 4 9 2020245 28 2018 11.2848 4 10 40023 35 2016 15.6811 2 11 40023 35 2017 15.7506 2 12 40023 35 2018 15.8072 2 13 265160055 36 2016 13.5722 6 14 265160055 36 2017 13.5261 7 15 265160055 36 2018 13.4228 6 16 151340376 28 2016 13.0350 3 17 151340376 28 2017 12.9871 3 18 151340376 28 2018 13.2472 4 19 293110377 38 2016 13.5696 3 20 293110377 38 2017 13.6979 3 21 272200375 38 2016 11.6168 1 22 272200375 38 2017 11.6243 2 23 272200375 38 2018 11.6396 2 24 130376 36 2016 12.4459 1 25 130376 36 2017 12.8828 1 26 130376 36 2018 12.8970 1 27 20342 28 2016 12.9657 4 28 20342 28 2017 13.4190 3 29 20342 28 2018 13.1806 4 30 10421 38 2016 15.4309 2 31 10421 38 2017 15.2439 2 32 10421 38 2018 15.1993 2 33 265160089 33 2016 11.8636 2 34 265160089 33 2017 12.0782 2 35 265160089 33 2018 12.0895 2 36 30111 38 2016 15.1774 5 37 30111 38 2017 15.2131 7 38 30111 38 2018 15.2249 6 39 60435 42 2016 12.8841 2 40 60435 42 2017 12.8980 1 41 60435 42 2018 13.0379 2 42 20023 73 2016 15.8567 5 43 20023 73 2017 15.8304 5 44 20023 73 2018 15.9903 5 45 30245 38 2016 13.9173 8 46 30245 38 2017 15.0267 9 47 30245 38 2018 15.1380 11 48 82560245 73 2016 12.4837 1 49 82560245 73 2017 12.7665 2 50 82560245 73 2018 12.9239 1 51 30198 56 2016 14.5976 2 52 30198 56 2017 14.6309 2 53 30198 56 2018 14.6419 2 54 30322 73 2016 14.6623 1 55 30322 73 2017 14.6135 2 56 30322 73 2018 14.6001 2 57 18510373 38 2016 11.9250 4 58 18510373 38 2017 11.9576 4 59 18510373 38 2018 11.9184 4 60 20111 38 2016 14.6730 8 61 20111 38 2017 14.6226 9 62 20111 38 2018 14.7825 10 63 158380359 73 2017 11.2266 1 64 158380359 73 2018 11.0898 1 ... View more

Is there a way to check whether there is a new person in my grouped data?   HospitalID PersonID Year 1 1 2005 1 2 2005 1 3 2005 1 1 2006 1 2 2006 1 3 2006 1 4 2006 2 5 2005 2 6 2005 2 7 2005 2 5 2006 2 6 2006 2 8 2006 3 9 2005  3 10 2005 3 9 2006  3 10 2006  As shown above, there were 1, 2, and 3 in hospital A in year 2005. There is a new person, 4, in year 2006. In this case, I would like to give a value of 1 to observations of hospital 1 and year 2006 (i.e., 1 for the observations of groups where there is a new person). In hospital 2, there are 5, 6, and 7 in year 2005. But next year, 7 left while 8 joined. In this case, I would like to assign missing to the observations of this group (i.e., hospital 2 - year 2006). In hospital 3, there are 9 and 10 in years 2005 and 2006. For the observations of this group, I would like to assign 0 to the observations of this group.   The resulting dataset would be either at hospital-year level or at hospital-person-yearl level, whichever is OK for me. Is there a simple way to make this dataset? Thanks in advance!   ... View more

Dear All,   I have two sets, one is at an aggregated-level while the other is at a disaggregated-level. I would like to merge summary stat (e.g., mean, sum) from the second set with the first set by using proc SQL. I did this as shown below. I think it could be shorter, but I don't know how to do so. Please keep in mind that the aggregated-level has another variable that is not included in the the other data set. So I would like to use proc SQL to get summary information from the disaggregated set.   To illustrate, I generated the two datasets. - ID_set: Set at an aggregated-level. - Data_set: Set at a disaggregated-level.   *Generating ID_Set; data ID_set; set sashelp.cars; keep Make Type; run; proc sort data= ID_set nodupkey; by Make; run; *N= 38; *Generating Data_Set; data Data_set; set sashelp.cars; keep Make EngineSize; run; *N= 428;   Using the two datasets, I would like to calculate sum and average of EngineSize and attach them to the first data set. I did PROC SQL and PROC SORT to keep the unique observations at an aggregated-level. I feel like they can be done at once in PROC SQL (without PROC SORT), but I don't know how to do so. Any help would be appreciated!      proc sql; create table Merged as select a.*, sum(b.EngineSize) as Sum_Engine, Mean(b.EngineSize) as Mean_Engine from ID_set a left join Data_set b on a.Make= b.Make group by a.Make; quit; proc sort data= Merged nodupkey; by Make; run;     ... View more

Dear All,   I was wondering if dm statement can set up placement and size of "Results Viewer - SAS output". I use a keyboard shortcut to clear my log, output, and odsresults windows from time to time. After this, Results Viewer - SAS output gets always maximized. Is there a way for me to set up the placement of Results Viwer by revising the below dm statement?     log; clear; output; clear; odsresults; clear; ... View more