I'm running proc report on a proc freq output table and the resulting table looks good except the rows under the "year" headers are offset (see data example below).   The code is running without errors -- any ideas what's going wrong?   data test; input year biomarkername $ biomarkerstatus $ count pct_row; cards; 2017 ALK N 2716 90.8969 2017 ALK P 87 2.9116 2017 ALK U 185 6.1914 2017 BRA N 1786 92.1569 2017 BRA P 91 4.6956 2017 BRA U 61 3.1476 2018 ALK N 2839 92.3552 2018 ALK P 72 2.3422 2018 ALK U 163 5.3025 2018 BRA N 2364 92.0202 2018 BRA P 136 5.2939 2018 BRA U 69 2.6859 ; run; options missing=' '; proc report data=test missing; column biomarkername year,(biomarkerstatus count pct_row); define biomarkername / group ; define year / across ; run;   ... View more

I'm outputting the frequency counts and %s from a crosstab table to a data set (out=test_freq in example code below).   The conditional percentages by x values are correct in the SAS Results tab, however the %s in the test_freq output dataset are not correct.   Is my syntax incorrect?   data test; input x $ y $; cards; Alice blue Alice blue Alice blue Alice green Alice green Alice red Bob blue Bob green Bob green Bob green Bob red Bob red Bob red Bob red ; run; proc freq data=test; table x*y / nopercent nocol out=test_freq; run; ... View more

Using the following dataset for illustration:   data test; input id $ x y; cards; a1 1 2 a1 1 4 a1 2 8 a1 2 16 b1 2 2 b1 2 4 b1 2 8 b1 3 16 b1 4 32 ; run; I'd like to retain all records sharing the lowest values of x by subject id: id x y a1 1 2 a1 1 4 b1 2 2 b1 2 4 b1 2 8 There's a simple solution I'm missing -- tried the following code but it didn't work: data test2; set test; by id x; if first.x; run; ... View more

I'm running a repeated measures multiple regression for a healthcare dataset. Each member has two records, one pre-intervention and one post-intervention outcome value (# of skilled nursing facility visits).   My code looks something like this: ods output GEEEmpPEst=parms; proc genmod data=health_data;    class member_id;    y = time*class time class x1 x2 x3 / dist=negbin ;    repeated subject=member_id / type=ind; run; ods output close;   When I run the procedure the regression coefficients are calculated but the stnd errors, confidence intervals, and p-values are blank.   Any ideas what's going wrong?   Algorithm converged. GEE Model Information Correlation Structure Independent Subject Effect hce_member_id (1199 levels) Number of Clusters 1199 Correlation Matrix Dimension 8 Maximum Cluster Size 8 Minimum Cluster Size 2 ERROR: Error in estimation routine. GEE Fit Criteria QIC 0.0000 QICu 0.0000 Analysis Of GEE Parameter Estimates Empirical Standard Error Estimates Standard 95% Confidence Parameter Estimate Error Limits Z Pr > |Z| Intercept 223254.9 . . . . . time*class -48.2170 . . . . . time -130.091 . . . . . class 150.1242 . . . . . x1 -120.939 . . . . . x2 -11.0064 . . . . . x3 -5.7767 . . . . .   ... View more

I'd like to edit an ICD-10 diagnosis category format file by adding new categories.   I'm running code written by a former coworker and was unable to find the original proc format code that created the format catalog.   However, running the following code:     libname fmt "D:\Project_Files\395_Matched_Cohort_Study\engaged_performance"; options obs=max compress=yes fmtsearch=(fmt); proc format library=fmt cntlout=FormatDef; run; outputted a copy of the format dataset FormatDef into my work folder:   Can I edit this file by adding new lines to the dataset, then save as fmt.FormatDef (or fmt.FormatDef_new)?      ... View more