topic Re: Format and ref value in SAS Programming

Format and ref value

wajs — Wed, 13 Mar 2019 02:23:17 GMT

Hi, I am working on a dataset with microorganisms. For Proc format, my values are CMV_IgG_rev 1.'Present', 2.'Absent'. In Proc logistic, I put (ref='Absent') or (ref='2'), log turns out invalid reference value. Please guide me.

Thanks

Re: Format and ref value

PGStats — Wed, 13 Mar 2019 02:46:50 GMT

Please post the code, or better, the log of the format and logistic proc calls.

Re: Format and ref value

wajs — Wed, 13 Mar 2019 02:51:57 GMT

proc logistic data= pe. matched_21_plus_co_inf_updated;
strata group;
class CMV_IgG_rev (ref='2');
model PEstatus (event='1')= CMV_IgG_rev;
run;

316 proc logistic data= pe. matched_21_plus_co_inf_updated;
NOTE: Data file PE.MATCHED_21_PLUS_CO_INF_UPDATED.DATA is in a format that is native to another
host, or the file encoding does not match the session encoding. Cross Environment Data Access
will be used, which might require additional CPU resources and might reduce performance.
317 strata group;
318 class CMV_IgG_rev (ref='2');
319 model PEstatus (event='1')= CMV_IgG_rev;
320 run;

ERROR: Invalid reference value for CMV_IgG_rev.
NOTE: The SAS System stopped processing this step because of errors.
NOTE: There were 78 observations read from the data set PE.MATCHED_21_PLUS_CO_INF_UPDATED.
NOTE: PROCEDURE LOGISTIC used (Total process time):
real time 0.03 seconds
cpu time 0.03 seconds

Re: Format and ref value

Tom — Wed, 13 Mar 2019 03:07:52 GMT

If you think the variable is numeric with single digit values then try using a format statement to make sure.

proc logistic data= pe. matched_21_plus_co_inf_updated;
  strata group;
  class CMV_IgG_rev (ref='2');
  format CMV_IgG_rev 1.;
  model PEstatus (event='1')= CMV_IgG_rev;
run;

Re: Format and ref value

wajs — Wed, 13 Mar 2019 04:52:38 GMT

990 proc logistic data= pe. matched_21_plus_co_inf_updated;
NOTE: Data file PE.MATCHED_21_PLUS_CO_INF_UPDATED.DATA is in a format that is native to another
host, or the file encoding does not match the session encoding. Cross Environment Data Access
will be used, which might require additional CPU resources and might reduce performance.
991 strata group;
992 class CMV_IgG_rev (ref='2');
993 format CMV_IgG_rev 1.;
994 model PEstatus (event='1')= CMV_IgG_rev;
995 run;

Re: Format and ref value

Kurt_Bremser — Wed, 13 Mar 2019 08:34:51 GMT

Maxim 3: Know Your Data.

Run a proc contents on your dataset (to see the type of CMV_IgG_rev) and a proc freq for variable CMV_IgG_rev.

Re: Format and ref value

Tom — Wed, 13 Mar 2019 13:55:30 GMT

Look at what is in your dataset.

data test;
  set PE.MATCHED_21_PLUS_CO_INF_UPDATED;
  formatted=quote(trim(vvalue( CMV_IgG_rev)));
  raw =  CMV_IgG_rev  ;
run;
proc freq ;
  tables raw*formatted / list missing;
run;

Re: Format and ref value

wajs — Wed, 13 Mar 2019 15:48:49 GMT

Hi,

Thanks for the guidance! There was issue with my values in throriginal dataset which I could find after running proc freq.

wajs

Re: Format and ref value

wajs — Wed, 13 Mar 2019 15:51:56 GMT

Hi Tom,

Thanks for this advice!

Here is the output; what should do wth "."?

The SAS System

The FREQ Procedure

raw formatted Frequency Percent CumulativeFrequency CumulativePercent. " ."0 " 0"1 " 1"

5	6.41	5	6.41
40	51.28	45	57.69
33	42.31	78	100.00