@di_niu0 wrote: Hi there. I am converting some Matlab code to SAS code. In Matlab, there is a built-in function that calculates the matrix exponential of a matrix called EXPM. There is no such built-in function for PROC IML. Does anyone know how to write the SAS code to produce the same result as EXPM? Thank you!  I think you are mistaken. The SAS IML language has supported the EXPMATRIX function, which computes the exponential of a matrix, for 20+ years. ... View more

E.X. proc format; value fmt 1='X' 2='Y' ; run; and the value of Trt_CC variable is 1,2 Therefore, 'X' ,'Y' is the FORMATTED values. 1,2 is the internal values. ... View more

Just to clarify something that was said further up, you can in check for the presence of an array value in another array (though it's correct that the way you specified 'array_1` was a syntax error.  For example, this works fine:   data test; array c {3} $5 _temporary_ ('dm', 'chf', 'ami'); length dx1-dx5 $5; array d {*} dx:; dx3='chf'; inC=(d[3] in c); run; ...the above produces a value of 1 for the variable 'inC', i.e., true.     ... View more

This worked beautifully!  Thank you! ... View more

How about a multlabel format.   data all_freqs; PARCAT3='HAEMATOLOGY';btoxgr='Grade0';TRT01A='LZA';PARAMCD='GRAN';cnt=12;output; PARCAT3='HAEMATOLOGY';btoxgr='Grade2';TRT01A='LZA';PARAMCD='GRAN';cnt=20;output; PARCAT3='HAEMATOLOGY';btoxgr='Grade1';TRT01A='LZA';PARAMCD='LYMP';cnt=19;output; PARCAT3='HAEMATOLOGY';btoxgr='Grade2';TRT01A='LZA';PARAMCD='KKSP';cnt=90;output; PARCAT3='LIVER AND KIDNEY';btoxgr='Grade0';TRT01A='LZA';PARAMCD='PLAT';cnt=12;output; PARCAT3='LIVER AND KIDNEY';btoxgr='Grade2';TRT01A='LZA';PARAMCD='WBC';cnt=22;output; PARCAT3='LIVER AND KIDNEY';btoxgr='Grade1';TRT01A='LZA';PARAMCD='WBC';cnt=24;output; PARCAT3='LIVER AND KIDNEY';btoxgr='Grade3';TRT01A='LZA';PARAMCD='AST';cnt=11;output; run; proc sort; by parcat3 paramcd trt01a; run; proc format; value $txg(notsorted multilabel) 'Grade0' = 'Grade 0' 'Grade1' = 'Grade 1' 'Grade2' = 'Grade 2' 'Grade3' = 'Grade 3' 'Grade4' = 'Grade 4' 'Grade0'-'Grade4' = 'Total' ; quit; proc summary data=all_freqs nway missing completetypes; by parcat3 paramcd trt01a; class btoxgr / mlf order=data preloadfmt; format btoxgr $txg.; freq cnt; output out=summary; run; proc print; run; ... View more

The response here have been of very high quality. Proc Studley thanks you all for your service! ... View more

data have; input ID age dx :$15. weight test :$15.; datalines; 1 29 covid 140 covid . . flu . blood_pressure . . diabetes . . 2 33 high_bp 165 covid . . gest_diab . blood_pressure . . preeclamp . . ; data temp; set have; retain _id _age; if not missing(id) then _id=id; if not missing(age) then _age=age; drop id age; run; proc sql noprint; select max(n) into :n from (select count(*) as n from temp group by _id,_age); quit; proc summary data=temp nway; class _id _age; output out=want idgroup(out[&n.] (dx weight test)=); run; ... View more

@E_W wrote: Thank you, Paige. Is there a way to do this without creating new column names? There are other programs that have data in the same column and I would like to maintain that. After you create a data set with the new column names, you can move the values into the old column names and then delete the new column names. ... View more

could you share results of this code with us: data test; set <your Data Set> (keep= <one variable that makes problems>); test1=<one variable that makes problems>; test2 = compress(test1,,'kw'); a=put(test1,$hex128.); b=put(test2,$hex128.); run; proc print data = test; run; Bart ... View more

Astounding, Don't be offended. I truely understand what you mean. We all provide some suggestion based on the data posted by OP . We don't know how it looks like or how it was in reality. So base on the data posted by OP and write code to get what OP want, That is OK. My real point: make my own code to be as strong/robust as it could be to take into acount of many accident scenarios . All the best. ... View more

An aleternative to using symget (a DATA step function  that will execute once per iteration of the loop), would be to use a macro quoting function, such as %superq, e.g.:   1 %let string = "This is sample string" & char(10) & "This is second line of sample string"; 2 %put &=string; STRING="This is sample string" & char(10) & "This is second line of sample string" 3 4 data x; 5 a="%superq(string)"; 6 put a= ; 7 run; a="This is sample string" & char(10) & "This is second line of sample string" NOTE: The data set WORK.X has 1 observations and 1 variables. ... View more

data have; input group $ order L1 L2 ; datalines; A1 1 45 36 A1 2 25 55 A2 1 88 99 A2 2 74 78 A3 1 4 . 25 A3 2 88 99 B4 1 78 65 B4 2 . . ; run; proc sql noprint; select max(n) into :n from (select count(*) as n from have group by group); quit; proc summary data=have nway; class group; output out=want idgroup(out[&n.] (L1 L2)=); run; ... View more

To get a variable with those missing values you will need an separate variable to keep the count.  So use the BY statement to find the groups. The SUM statement to increment the counter and an IF statement to assign the counter to create your wanted variable with the missing values otherwise. data want; set ds1; by id stop notsorted; if first.id then counter=0; counter + (first.stop and stop); if stop then sseq=counter; run; Result   ... View more

Thanks so much, it was helpful ... View more

The code snippet you've provided is creating a custom format in SAS, and I'll walk you through each step. Step-by-Step Explanation Creating the Dataset disp : data disp; set trcb.DispositionCode_Xref; length start $2 label $100; fmtname = '$disposition'; start = Disp_Code; label = Disp_Code_Desc; keep start label fmtname; run; data disp; : This starts a new dataset named disp . set trcb.DispositionCode_Xref; : This reads in data from an existing dataset called trcb.DispositionCode_Xref . length start $2 label $100; : This sets the length of two new variables: start (length 2 characters) label (length 100 characters) fmtname = '$disposition'; : This assigns the value '$disposition' to a new variable called fmtname . This variable will be used to define the name of the custom format. start = Disp_Code; : This assigns the value from the variable Disp_Code to the start variable. In the context of creating formats, start represents the values that will be formatted. label = Disp_Code_Desc; : This assigns the value from the variable Disp_Code_Desc to the label variable. The label represents the formatted output that corresponds to the start value. keep start label fmtname; : This keeps only the start , label , and fmtname variables in the final disp dataset, discarding any others from the original dataset. Summary of this step: You're creating a new dataset called disp that contains three variables: start (the value to format), label (the formatted output), and fmtname (the name of the format). Sorting the Dataset: proc sort data=disp nodupkey; by start; run;   proc sort data=disp nodupkey; : This sorts the disp dataset by the start variable. The nodupkey option removes any duplicate records based on the start variable, ensuring that each value to be formatted is unique. by start; : This specifies that the sorting should be done by the start variable. Summary of this step: You're ensuring that the dataset disp is sorted by the start values and that there are no duplicate start values. What Does This All Mean? This code is preparing a dataset that will be used to create a custom format in SAS. A format in SAS is a way to map values (like Disp_Code ) to corresponding labels (like Disp_Code_Desc ), so that when you apply this format to a variable, the value gets displayed as the label. In this specific case, a custom format called '$disposition' is being created, which will map each Disp_Code to its corresponding Disp_Code_Desc . Final Step (Not in the Code) Typically, after preparing this dataset ( disp ), you would use the PROC FORMAT procedure to actually create the format from this dataset:   proc format cntlin=disp; run;   proc format cntlin=disp; : This tells SAS to create a format based on the disp dataset you've created. Once this is done, you can use the '$disposition' format to display Disp_Code values as their corresponding Disp_Code_Desc labels in reports or datasets.       ... View more