Thanks Paige, using your feedback readjusted code to below:   DATA HNC; SET SOURCE; IF SUBSITE IN = (C00.3, C00.4, C00.5, C00.6, C00.8, C00.9, C02.0, C02.1, C02.2, C02.3, C02.8, C02.9 C03.0, C03.1, C03.9 C04.0, C04.1, C04.8, C04.9 C05.0 C05.1, C05.2, C05.8, C05.9 C06.0, C06.1, C06.2 C06.8, C06.9) THEN SUB = 'ORAL CAVITY CANCER'; IF SUBSITE IN = (C01.X, C02.4, C09.0, C09.1, C09.8, C09.9, C10.0, C10.1, C10.2, C10.3, C10.4, C10.8, C10.9 C14.0, C14.2) THEN SUB = 'OROPHARYNGEAL CANCER'; IF SUBSITE IN = (C32.0, C32.1, C32.2, C32.3, C32.8, C32.9) THEN SUB = 'LARYNGEAL CANCER'; IF SUBSITE IN = (C00.3, C00.4, C00.5, C00.6, C00.8, C00.9 C01.X, C02.0, C02.1, C02.2, C02.3, C02.4, C02.8, C02.9 C03.0, C03.1, C03.9 C04.0, C04.1, C04.8, C04.9 C05.0, C05.1, C05.2, C05.8, C05.9 C06.0, C06.1, C06.2, C06.8, C06.9 C09.0, C09.1, C09.8, C09.9 C10.0, C10.1, C10.2, C10.3, C10.4, C10.8, C10.9 C11.0, C11.1, C11.2, C11.3, C11.8, C11.9 C12.X, C13.0, C13.1, C13.2, C13.8, C13.9 C14.0, C14.2, C14.8 C32.0, C32.1, C32.2, C32.3, C32.8, C32.9) THEN SUB = 'HEAD AND NECK CANCER'; RUN;   This produced more errors though 😞    531 DATA HNC; 532 SET SOURCE; ERROR: File WORK.SOURCE.DATA does not exist. 533 IF SUBSITE IN = (C00.3, C00.4, C00.5, C00.6, C00.8, C00.9, - 390 76 ERROR 390-185: Expecting an relational or arithmetic operator. ERROR 76-322: Syntax error, statement will be ignored. 534 C02.0, C02.1, C02.2, C02.3, C02.8, C02.9 535 C03.0, C03.1, C03.9 536 C04.0, C04.1, C04.8, C04.9 537 C05.0 C05.1, C05.2, C05.8, C05.9 538 C06.0, C06.1, C06.2 C06.8, C06.9) 539 THEN SUB = 'ORAL CAVITY CANCER'; 540 541 IF SUBSITE IN = (C01.X, - 22 ----- 557 ERROR: The right-hand operand must be an array name or a constant value list. The specified name NAME, is not an array. ERROR: DATA STEP Component Object failure. Aborted during the COMPILATION phase. ERROR 22-322: Expecting a name. ERROR 557-185: Variable C01 is not an object. NOTE: The SAS System stopped processing this step because of errors. NOTE: DATA statement used (Total process time): real time 0.04 seconds cpu time 0.01 seconds 542 C02.4, 543 C09.0, C09.1, C09.8, C09.9, 544 C10.0, C10.1, C10.2, C10.3, C10.4, C10.8, C10.9 545 C14.0, C14.2) 546 THEN SUB = 'OROPHARYNGEAL CANCER'; 547 548 IF SUBSITE IN = (C32.0, C32.1, C32.2, C32.3, C32.8, C32.9) 549 THEN SUB = 'LARYNGEAL CANCER'; 550 551 IF SUBSITE IN = (C00.3, C00.4, C00.5, C00.6, C00.8, C00.9 552 C01.X, 553 C02.0, C02.1, C02.2, C02.3, C02.4, C02.8, C02.9 554 C03.0, C03.1, C03.9 555 C04.0, C04.1, C04.8, C04.9 556 C05.0, C05.1, C05.2, C05.8, C05.9 557 C06.0, C06.1, C06.2, C06.8, C06.9 558 C09.0, C09.1, C09.8, C09.9 559 C10.0, C10.1, C10.2, C10.3, C10.4, C10.8, C10.9 560 C11.0, C11.1, C11.2, C11.3, C11.8, C11.9 561 C12.X, 562 C13.0, C13.1, C13.2, C13.8, C13.9 563 C14.0, C14.2, C14.8 564 C32.0, C32.1, C32.2, C32.3, C32.8, C32.9) 565 THEN SUB = 'HEAD AND NECK CANCER'; 566 567 568 RUN; ... View more

Hey sorry for delayed reply.   Got dpt. statistician to load it o=in eventually once they got back from leave.   All good now.   BW,   Craig  ... View more

Thanks Quentin,   Hmmm I'm without any stats support in my department till the 20th.   Was hoping calibration curve would be relatively straightforward to do myself.   Shall keep code form forums and see if they can modify it.    BW,   Craig   ... View more

Hi Quentin,   Tried adapting code from https://blogs.sas.com/content/iml/2019/02/20/easier-calibration-plot-sas.html to my data / logistic model    /* NEW in SAS/STAT 15.1 (SAS 9.4M6): PLOTS=CALIBRATION option in PROC LOGISTIC */ title "Calibration Plot for a Quadratic Model"; title2 "Created by PROC LOGISTIC"; proc logistic data=ARCAGE plots=calibration(CLM ShowObs); model y(Event='1') = x x*x / GOF; /* New in 15.1: More goodness-of-fit statistics */ run;   Gave log error of    3547 /* NEW in SAS/STAT 15.1 (SAS 9.4M6): PLOTS=CALIBRATION 3547! option in PROC LOGISTIC */ 3548 title "Calibration Plot for a Quadratic Model"; 3549 title2 "Created by PROC LOGISTIC"; NOTE: PROCEDURE LOGISTIC used (Total process time): real time 0.00 seconds cpu time 0.00 seconds ERROR: Variable Y not found. NOTE: The SAS System stopped processing this step because of errors. 3550 proc logistic data=ARCAGE plots=calibration(CLM ShowObs); 3551 model y(Event='1') = x x*x / GOF; /* New in 15.1: 3551! More goodness-of-fit statistics */ 3552 run;   I'm essentially unsure how to adapt the code from that to mine (zero coding background) thought DATA = ARCAGE would fix it but apparently not.  ... View more

Hi Quentin,   Thanks for your answers - hope you're well. Have tried the calibration code but unsure how to make it work to my model - do you know what I'm doing wrong.    Best wishes,   Craig      /* NEW in SAS/STAT 15.1 (SAS 9.4M6): PLOTS=CALIBRATION option in PROC LOGISTIC */ title "Calibration Plot for a Quadratic Model"; title2 "Created by PROC LOGISTIC"; proc logistic data=LogiSim plots=calibration(CLM ShowObs); model y(Event='1') = x x*x / GOF; /* New in 15.1: More goodness-of-fit statistics */ run; PROC LOGISTIC DATA = ARCAGE; plots=calibration(CLM ShowObs); model y(Event='1') = x x*x / GOF; /* New in 15.1: More goodness-of-fit statistics */ run;   ... View more

Thanks Quentin,   Sorry had some tech issues. Shall give this a crack and report back.   KR,   Craig  ... View more

Brilliant Tom - this is what I meant. Sorry for the delayed reply - all the best. Craig ... View more

Hi Paige,   Thanks for your reply. Sorry if this seems counterintuitive but as someone of a non-coding background when I say "merge" - what I mean is to have all my data in one place, ready for analysis. Not split between 2 files. Let's take hypothetical person X - I have their sex, dietary, smoking, *insert here* info. In this separate file I also have the HPV serostatus(es) for person X. I don't wish to manipulate or analyse the data in either just yet - just to have all my potential variables in the one place for person X.    I don't wish to touch the variables or "merge" these. I'd simply like these variables to all be together in the one programme so I can mix and match for regression models e.g. below:   proc logistic data= ARCAGE; FORMAT SMOK SMOK. DRINK DRINK. SEX SEX.; CLASS SMOK ( REF = 'Never smoker') DRINK (REF = '0: Never') SEX ( REF = '2: Female' ) BMIC ( REF = '2' ); model CASE = AGE10 SEX prolev SMOK SMOKP DRINK alcoIC VEGETC FRUITC DENTIS BMIC /rocci RL SELECTION = F; run;   Thanks for your understanding.    Craig  ... View more

Hi Paige,   Thanks for your reply.    I see - I'm essentially trying to merge the datatsets so I have all the information in one place for study participants and analysis - i.e. combining the demographic, behavioural and clinical data I already have with the HPV serology in the excel file (both are from same study).    Not sure what step code is sorry (sorry am a total newbie at this). Nothing crazy coding wise mostly descriptive analysis and logistic regression modeling.    Many thanks,   Craig  ... View more

Thanks so much Paige - got it working!   Many thanks,   Craig  ... View more