I'm trying to do nominal and define the reference category. Thank you for your help! ... View more

Also, the documentation you referred me to states, "EVENT=’category’ | keyword specifies the event category for the binary response model. PROC LOGISTIC models the probability of the event category. The EVENT= option has no effect when there are more than two response categories." The response variable (mj_user) is NOT binary. It has 4 response levels. Should I be using a different PROC? ... View more

Here is the exact log:  1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 72 73 PROC LOGISTIC DATA=ucsd_ed.dataset_1 ORDER=FORMATTED; 74 CLASS med_or_rec (REF="For medical purposes only")/PARAM=REF; 75 MODEL mj_user (REF="Never") = med_or_rec /LINK=GLOGIT RSQ; 76 FORMAT mj_user mj_user_. med_or_rec med_or_rec_.; 77 TITLE "XXX"; 78 RUN; ERROR: Invalid reference value for mj_user. NOTE: The SAS System stopped processing this step because of errors. NOTE: There were 2577 observations read from the data set UCSD_ED.DATASET_1. NOTE: PROCEDURE LOGISTIC used (Total process time): real time 1.35 seconds cpu time 0.91 seconds 79 80 81 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 93 When I change REF to EVENT, the PROC LOGISTIC seems to run fine. Here is the log:  1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 72 73 PROC LOGISTIC DATA=ucsd_ed.dataset_1 ORDER=FORMATTED; 74 CLASS med_or_rec (REF="For medical purposes only")/PARAM=REF; 75 MODEL mj_user (EVENT="Never") = med_or_rec /LINK=GLOGIT RSQ; 76 FORMAT mj_user mj_user_. med_or_rec med_or_rec_.; 77 TITLE "XXX"; 78 RUN; NOTE: The EVENT= option for the response variable is ignored for LINK=GLOGIT. NOTE: PROC LOGISTIC is fitting the generalized logit model. The logits modeled contrast each response level against the reference level (mj_user='Current, regular'). Use the response variable option REF= if you want to change the reference level. NOTE: Convergence criterion (GCONV=1E-8) satisfied. NOTE: There were 2577 observations read from the data set UCSD_ED.DATASET_1. NOTE: PROCEDURE LOGISTIC used (Total process time): real time 4.99 seconds cpu time 3.40 seconds 79 80 81 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 93 "The logits modeled contrast each response level against the reference level (mj_user='Current, regular'). Use the response variable option REF= if you want to change the reference level."   I do want to change the reference level, but when I do...I get the original error message.   Yes, there is a lot of missing data in the variables that produce the error. These data are from a complex survey so MOST of the data are missing for these/certain variables. ... View more

This is great. Thank you! Another question: Is there a way to look at these responses across another variable (eg. gender, see attachment Desired output 2.0.pdf)? ... View more

Also...I was just trying to get some help on Step 1 of this process. My ultimate goal is to look at these responses across another variable (See Desired Output 2.0). Perhaps I should have posted that from the beginning...   I'm basically wanting to do a PROC FREQ on a multiple response variable where each response is it's own variable. I'm open to other ways to think about this.    Thank you.  ... View more

@ballardw Thanks for your help.   Apologies for the confusion.    <Your table of percentages shows a total of 100%. > The desired numerator is the total number of responses for that question. The desired denominator is the total number of responses overall - for all questions. As opposed to the total number of respondents.    Also, ideally each question would be output as a row in a table, like a single response variable in PROC FREQ output.    Your syntax is useful though and I may use it instead! Thank you.   Is there alternative syntax you might recommend for percent of responses?   🙏 ... View more

1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 72 73 DATA data_working_keep; 74 KEEP record_id provider_type___1-provider_type___14; 75 SET mchps.data_working; 76 77 array p provider_: ; ERROR: All variables in array list must be the same type, i.e., all numeric or character. ERROR: All variables in array list must be the same type, i.e., all numeric or character. 78 if sum(of p(*))=0 then call missing(of p(*)); 79 80 LABEL 81 provider_type___1= "Medical Doctor (i.e. MD)" 82 provider_type___2 = "Osteopathic Doctor (i.e. DO)" 83 provider_type___3 = "Naturopathic Doctor (i.e. ND)" 84 provider_type___4= "Physicians Assistant (i.e. PA)" 85 provider_type___5= "Nurse Practitioner (i.e. NP)" 86 provider_type___6= "Registered Nurse (i.e. RN)" 87 provider_type___7= "Pharmacist (i.e. Rph)" 88 provider_type___8= "Massage therapist (i.e. LMP)" 89 provider_type___9 = "Acupuncturist (i.e. LA.c)" 90 provider_type___10= "Chiropractor (i.e. DC)" 91 provider_type___11 = "Clinical social worker (i.e. LCSW)" 92 provider_type___12 = "Physical therapist (i.e. PT)" 93 provider_type___13= "Psychiatrist" 94 provider_type___14= "Other"; 95 RUN; NOTE: Character values have been converted to numeric values at the places given by: (Line):(Column). 78:15 NOTE: The SAS System stopped processing this step because of errors. WARNING: The data set WORK.DATA_WORKING_KEEP may be incomplete. When this step was stopped there were 0 observations and 15 variables. NOTE: DATA statement used (Total process time): real time 0.02 seconds cpu time 0.01 seconds ... View more