I do not know what is wrong with my code: Below is the distribution of deposi: deposi Frequency Percent Cumulative Frequency Cumulative Percent 0 498 59.50 498 59.50 1 239 28.55 737 88.05 2 100 11.95 837 100.00   I fitted the SAS procedure below and had the error proc genmod data=stroke; class gender(ref="0") race_ethnic(ref="0") deposi(ref="0") pay_source(ref="0") cmg_tier(ref="a") ; model mrs_discharge= age gender race_ethnic deposi pay_source cmg_tier report_cmg weight_cmg expected_los/dist=gamma link=log type3; run;   ERROR: Invalid reference value for deposi. ERROR: No valid observations due to invalid or missing values in the response, explanatory, offset, frequency, or weight variable. NOTE: The SAS System stopped processing this step because of errors. NOTE: PROCEDURE GENMOD used (Total process time): real time 0.01 seconds cpu time 0.01 seconds   ... View more

I wish to compute the difference between the following two timepoints, that is Time2-Time2. But the time are characters, that is Char $16. I was expecting it to be Datetime. Anydtdtm40 but they are not. Please can you show me how to convert but Time1 and Time2 from Char $16 to Datetime. Anydtdtm40 and perform the subtraction.      Time1 Time2   1/3/18 1:30 PM .   . .   1/15/18 10:05 AM 1/16/18 1:31 AM   2/7/18 3:15 PM .   2/15/18 10:30 PM 2/16/18 11:50 AM   3/5/18 8:20 AM 3/6/18 2:02 AM   3/19/18 9:30 AM 3/19/18 1:55 PM             ... View more

Hello -,   I have this longitudinal data with population as outcome, I wish to ask the best statistical procedure. Treatment, sex, Pre/post, time are covariate and population is the outcome.     Treatment Sex Pre/Post time Population id Yes Female Pre 1 196 1 Yes Female Pre 3 132 1 Yes Female Pre 5 216 1 Yes Female Pre 10 520 1 Yes Female Post 1 2021 2 Yes Female Post 3 2551 2 Yes Female Post 5 3050 2 Yes Female Post 10 4811 2 Yes Male Post 1 6777 3 Yes Male Post 3 1165 3 Yes Male Post 5 2142 3 Yes Male Post 10 3135 3 No Female Pre 1 5376 4 No Female Pre 3 12736 4 No Female Pre 5 15278 4 No Female Pre 10 16896 4 No Female Post 1 68548 5 No Female Post 3 74211 5 No Female Post 5 101223 5 No Female Post 10 30569 5 No Male Pre 1 12141 6 No Male Pre 3 13280 6 No Male Pre 5 13869 6 No Male Pre 10 41440 6 ... View more

Hello,   I am trying to transform this data from wide to long but kept having issues: data long; set mental ; array Avg_hrs_sleep_month(0:2) Avg_hrs_sleep_month0-Avg_hrs_sleep_month2; array hrs_of_sleep_yest(0:2) hrs_of_sleep_yest0-hrs_of_sleep_yest2; *Insomnia0 Insomnia1 Insomnia2 see_yourself_6month0 see_yourself_6month1 see_yourself_6month2; do time=0 to 2; Avg_hrs_sleep_month = Avg_hrs_sleep_month(time); hrs_of_sleep_yest= hrs_of_sleep_yest(time); output; end; run;     215 data long; 216 set mental ; 217 array Avg_hrs_sleep_month(0:2) Avg_hrs_sleep_month0-Avg_hrs_sleep_month2; 218 array hrs_of_sleep_yest(0:2) hrs_of_sleep_yest0-hrs_of_sleep_yest2; 219 220 *Insomnia0 Insomnia1 Insomnia2 see_yourself_6month0 see_yourself_6month1 see_yourself_6month2 220! ; 221 do time=0 to 2; 222 Avg_hrs_sleep_month = Avg_hrs_sleep_month(time); ERROR: Illegal reference to the array Avg_hrs_sleep_month. 223 hrs_of_sleep_yest= hrs_of_sleep_yest(time); ERROR: Illegal reference to the array hrs_of_sleep_yest. 224 output; 225 end; 226 run; NOTE: The SAS System stopped processing this step because of errors. WARNING: The data set WORK.LONG may be incomplete. When this step was stopped there were 0 observations and 46 variables. WARNING: Data set WORK.LONG was not replaced because this step was stopped. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds   ... View more

I have the DateTime variable in a longitudinal data and i want to create two variable t1 and t2 as below:   id                     Datetime                                 t1                                                          t2        1               8/13/2020 10:47                            0 1               8/13/2020 15:04           8/13/2020 15:04 - 8/13/2020 10:47                8/13/2020 15:04 - 8/13/2020 10:47 1               8/13/2020  19:08          8/13/2020 15:04 - 8/13/2020  19:08              8/13/2020 10:47 - 8/13/2020  19:08            NB:  DateTime Num 8 DATETIME. ANYDTDTM40 ... View more

I have no issues running these codes below . But it seems I didn’t create the “antibiotics_flag" variable in my data, I just have the label, so every time I will have to re-run the whole codes (I have 20 codes blocks like this). Is there a way to save the “antibiotics_flag" as a permanent variable in my datasets?    data kid2012; set Kid_2012_core; array diagvar $ dx1-dx25; label antibiotics_flag = "Diagnosis code for antibiotics ADE"; antibiotics_flag = 0; do over diagvar;       if diagvar ="76074" or "9600 "=<diagvar<= "9609 " or diagvar="E856 " or "E9300 "=<diagvar<= "E9306 " or "E9308 "=<diagvar<= "E9309 "       or diagvar= "E9310 " or diagvar= "E9313 " then antibiotics_flag = 1;       end; PROC SURVEYFREQ DATA=kid2012; WEIGHT DISCWT; CLUSTER HOSP_KID ; STRATA KID_STRATUM ; table antibiotics_flag; run;     ... View more

I have this data with three classes that I wish to plot a forest plot on the same axis with three colors to distinguish between the three classes. The data is as follow;   Diagnosis       est.    LCL      UCL      Type Acne             2        1            3          Class I head        Acne             3        2.5         4.5       Class II head        Acne             3.2     2.5         4.3       Class III head      Asthma          4        1            6          Class I head        Asthma           5       3.5        7.5       Class II head        Asthma          7.2     4.5         8.3       Class III head  ..... The forest plot is for the diagnosis with point estimate (est.), lower confidence limit (LCL) and upper confidence limit (UCL) for each class (distinguish by a color for class I, a color for class II and a color for class III). Thank you for your help. ... View more