Hello, I want to determine whether an individual completed any form (event_date) within 60 days of their program (end_date). I have a demo data below and am struggle with the code. Any help please.  data A; input ssn @4 start_date :date9. @14 end_date :date9. @24 event_date :date9. ; infile datalines firstobs=2 ; informat start_date end_date event_date date9.; format start_date end_date event_date date9.; datalines; id start_date end_date event_date 101 01DEC2021 28FEB2022 09NOV2020 101 16OCT2019 18NOV2019 03FEB2020 102 16OCT2019 18NOV2019 07MAY2020 103 03FEB2022 02MAY2022 08JUL2022 103 03FEB2022 12FEB2022 15JAN2023 ;; run;  For example, id 101 first end_date was 28FEB2022 but the event date was less than 60days there the first id will not be my desired output. BUT id 101 second end_date was 18NOV2019 but the event date is 60days more than the end_date therefore the second ID will be part of my desired output. I don't want to dedup any data by id so I can capture all multiple events within one person(id). Thanks in advance! ... View more

Hello, I want to import a text file data into SAS. The file contain some issues that you will see in my example.  Some of the data is missing  and the data is not well organized as you can see. Student_id and d_pin are character variable. I would like to learn how to use either infile or proc import or both to import the data CORRECTLY. Here is sample code. I have over 10,000 records. Student_ID|D_PIN|SERVICE_DATE |1624413098|11/08/2023 |1628543409|11/08/2023 001020084|1615413808|10/12/2023 001021665|1616220906|03/31/2023 002781314|1462647405|07/19/2023 |1331118504|11/02/2023 |1611325610|12/26/2023 080504901||11/08/2023   Here is my code: data STUDENT; *; infile 'c:\dog\dummydata.txt' truncover delimiter = '|' firstobs=2; format SERVICE_DATE mmddyy10.; input Student_ID $ 1-9 D_PIN $ 11-20 SERVICE_DATE :mmddyy10.; run; I am getting many of the notes below. When I checked it seems it's because the original text file is not well arranged.   NOTE: Invalid data for SERVICE_DATE in line 1756 21-22. RULE: ----+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8----+----9----+----0 1756 |1331118504|11/02/2023 22 SERVICE_DATE=. Student_ID=|13311185 D_PIN=4|11/02/20 _ERROR_=1 _N_=1755 NOTE: Invalid data for SERVICE_DATE in line 1757 21-22. 1757 |1611325610|12/26/2023 22 Thanks           ... View more

Hello, I would like to convert the following character values into numeric. My values are in character;  Interst Freeway Pri Art Min Art Maj Col Min Col Local First I want to convert to numeric then convert to 1,2,3 in the following order. 'Interst' = 1 'Freeway' = 2 'Pri Art' = 3 'Min Art' = 4 'Maj Col' = 5 'Min Col' = 6 'Local' = 7 I initially created a format and i ran an input statement in the data step but I didn't get the expected result. My results were blank for funtional_num. Please help!! proc format; value $functional 'Interst' = 1 'Freeway' = 2 'Pri Art' = 3 'Min Art' = 4 'Maj Col' = 5 'Min Col' = 6 'Local' = 7 ; run; data A; length functional $10.; input functional$ 1-10; datalines; Interst Freeway Pri Art Min Art Maj Col Min Col Local ; run; data B; set A; functional_num = input(functional, 8.); format functional_num functional.; run; Result obtained: Obs functional functional_num 1      Interst                          . 2     Freeway                      . 3    Pri Art                          . 4   Min Art                        . 5   Maj Col                        . 6   Min Col                        . 7   Local                            . ... View more

Hello,  Please can anyone help with correctly using the code for empirical (observed) logits for ordinal response data. SAS provided this code in this paper: 37944 - Plotting empirical (observed) logits for binary and ordinal response data (sas.com) I carefully copied the code into SAS on demand and only changed the data=?, y=? & x=? data= name of the input data set. Required. y= name of the response variable. Required. x= space-separated list of predictor variables. Required.  but my code is not working. Please I would really appreciate any help. I wanted to use the code for my publication report.   The empirical is a long code so I will not copy it here but here is my log when i run it. 1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; NOTE: ODS statements in the SAS Studio environment may disable some output features. 69 70 /* INSTRUCTIONS FOR USE 71 72 Version 1.1 or higher of the CtoN macro is required by the EmpiricalLogitPlot 73 macro and must be defined before using the EmpiricalLogitPlot macro. Download 74 it from this location: http://support.sas.com/kb/60678 and see the Downloads 75 tab at that link. The code defining the CtoN macro and the code defining the 76 macros below must all be submitted in your SAS session to make the macros 77 available for use in that SAS Session. Do not alter any of the macro code 78 to ensure successful macro execution. After running the macro code, check 79 your SAS log to verify that no errors were generated. 80 *******************************************************************************/ 81 82 %macro EmpiricalLogitPlot(version, data=stcp.cleanedpracticum2, y=issgrp, x=functional person_type gender dayofweek 82 ! drug_bac lnage, const=0.5, neighbors=50, 83 smooth=.3, contcutoff=20, options=DESCENDING PANEL) / minoperator; 84 /*-------------------------------------------------------------------------- 85 A macro extension of SAS Usage Note 37944 for both class and continuous 86 covariates. Calls the OneContPlot macro to plot the empirical logits defined 87 using the NEIGHBORS/2 observed covariate values preceding and following each 88 observation against the levels of a covariate with more than CONTCUTOFF 89 levels. Calls the OneClassPlot macro to plot the empirical logits computed 90 within each level of a covariate with fewer than CONTCUTOFF levels. 91 92 version Optional. Any text immediately following the open parenthesis and 93 followed by a comma will display the macro version. 94 data= name of the input data set. Required. 95 y= name of the response variable. Required. 96 x= space-separated list of predictor variables. Required. 97 const= small value between 0 and 1 to prevent zero counts. Default: 0.5. 98 neighbors= number of neighbors around each observation used to compute 99 logits for continuous effects. Default: Average population size but 100 at least 1% of the sample size with a minimum of 2. 101 smooth= smoothing parameter used for loess curve fit to logits for 102 continuous effects. Default: 0.3. 103 contcutoff= if an x= variable has more than the specified number of levels 104 then logits are computed within sets of neighboring observations; 105 otherwise computed within each level of the variable. Default: 20. 106 options= DESCENDING (high response levels are of interest), ASCENDING (low 107 response levels are of interest), DATAORDER (use for character 108 response after sorting data so that first occurrence of each level 109 is in logically ascending or descending order as desired with the 110 first levels are of interest). Specify only one of DESCENDING, 111 ASCENDING, or DATAORDER. PANEL (single panel of all predictor 112 plots), NOPANEL (separate plot for each predictor). Default: 113 DESCENDING PANEL. 114 --------------------------------------------------------------------------*/ 115 116 %macro existchk(data=, var=, dmsg=e, vmsg=e); 117 %global status; %let status=ok; 118 %if &dmsg=e %then %let dmsg=ERROR; 119 %else %if &dmsg=w %then %let dmsg=WARNING; 120 %else %let dmsg=NOTE; 121 %if &vmsg=e %then %let vmsg=ERROR; 122 %else %if &vmsg=w %then %let vmsg=WARNING; 123 %else %let vmsg=NOTE; 124 %if %quote(&data) ne %then %do; 125 %if %sysfunc(exist(&data)) ne 1 %then %do; 126 %put &dmsg: Data set %upcase(&data) not found.; 127 %let status=nodata; 128 %end; 129 %else %if &var ne %then %do; 130 %let dsid=%sysfunc(open(&data)); 131 %if &dsid %then %do; 132 %let i=1; 133 %do %while (%scan(&var,&i) ne %str() ); 134 %let var&i=%scan(&var,&i); 135 %if %sysfunc(varnum(&dsid,&&var&i))=0 %then %do; 136 %put &vmsg: Variable %upcase(&&var&i) not found in data %upcase(&data).; 137 %let status=novar; 138 %end; 139 %let i=%eval(&i+1); 140 %end; 141 %let rc=%sysfunc(close(&dsid)); 142 %end; 143 %else %put ERROR: Could not open data set &data.; 144 %end; 145 %end; 146 %else %do; 147 %put &dmsg: Data set not specified.; 148 %let status=nodata; 149 %end; 150 %mend; 151 152 %let elptime = %sysfunc(datetime()); 153 %let _version=1.2; 154 %if &version ne %then %put NOTE: &sysmacroname macro Version &_version; 155 %let _elpopts = %sysfunc(getoption(notes)); 156 %let version=%upcase(&version); 157 %if %index(&version,DEBUG) %then %do; 158 options notes mprint 159 %if %index(&version,DEBUG2) %then mlogic symbolgen; 160 ; 161 ods select all; 162 %put _user_; 163 %end; 164 %else %do; 165 options nonotes nomprint nomlogic nosymbolgen; 166 ods exclude all; 167 %end; 168 169 /* Check for newer version */ 170 %if %index(&version,NOCHECK)=0 %then %do; 171 %let _notfound=0; %let _newver=0; 172 filename _ver url 'http://ftp.sas.com/techsup/download/stat/versions.dat' 173 termstr=crlf; 174 data _null_; 175 infile _ver end=_eof; 176 input name:$18. ver; 177 if upcase(name)="&sysmacroname" then do; 178 call symput("_newver",ver); stop; 179 end; 180 if _eof then call symput("_notfound",1); 181 run; 182 options notes; 183 %if &syserr ne 0 or &_notfound=1 or &_newver=0 %then 184 %put NOTE: Unable to check for newer version of &sysmacroname macro.; 185 %else %if %sysevalf(&_newver > &_version) %then %do; 186 %put NOTE: A newer version of the &sysmacroname macro is available at; 187 %put NOTE- this location: http://support.sas.com/ ; 188 %end; 189 %if %index(%upcase(&version),DEBUG)=0 %then options nonotes;; 190 %end; 191 192 /* Input checks and process options */ 193 %if %index(&y,/) %then %do; 194 %put ERROR: Events/Trials syntax is not supported in Y=.; 195 %goto exit; 196 %end; 197 198 %existchk(data=&data, var=&y &x); 199 %if &status=nodata or &status=novar %then %goto exit; 200 201 %if %quote(&y)= %then %do; 202 %put ERROR: Y= is required.; 203 %goto exit; 204 %end; 205 206 %if &x= %then %do; 207 %put ERROR: One or more variables must be specified in X=.; 208 %goto exit; 209 %end; 210 211 %if %sysevalf(&neighbors ne) %then %do; 212 %if %sysevalf(%sysfunc(mod(&neighbors,1)) ne 0 or &neighbors<=0) %then %do; 213 %put ERROR: NEIGHBORS= must be an integer value greater than zero.; 214 %goto exit; 215 %end; 216 %end; 217 %else %let neighbors=50; 218 %let nbhd=&neighbors; 219 220 %if %sysevalf(&contcutoff ne) %then %do; 221 %if %sysevalf(%sysfunc(mod(&contcutoff,1)) ne 0 or &contcutoff<=0) %then %do; 222 %let contcutoff=%sysfunc(ceil(%sysfunc(abs(&contcutoff)))); 223 %if %index(&version,DEBUG)=0 %then options notes;; 224 %put NOTE: CONTCUTOFF= value must be a positive integer.; 225 %put NOTE- CONTCUTOFF= set to &contcutoff; 226 %if %index(&version,DEBUG)=0 %then options nonotes;; 227 %end; 228 %end; 229 %else %let contcutoff=20; 230 231 %if %quote(&const) ne %then %do; 232 %if (%sysfunc(verify(%sysfunc(catx( ,&const)),'0123456789.'))>0) + 233 (%sysfunc(count(&const,.))>1) %then %do; 234 %put ERROR: The CONST= value must be a positive value less than 1.; 235 %goto exit; 236 %end; 237 %else %if &const<0 or &const>=1 %then %do; 238 %put ERROR: The CONST= value must be a positive value less than 1.; 239 %goto exit; 240 %end; 241 %end; 242 %else %do; 243 %put ERROR: CONST= is required and must be a positive value less than 1.; 244 %goto exit; 245 %end; 246 247 %if %quote(&smooth) ne %then %do; 248 %if (%sysfunc(verify(%sysfunc(catx( ,&smooth)),'0123456789.'))>0) + 249 (%sysfunc(count(&smooth,.))>1) %then %do; 250 %put ERROR: The SMOOTH= value must be a positive value.; 251 %goto exit; 252 %end; 253 %else %if %sysevalf(&smooth<=0) %then %do; 254 %put ERROR: The SMOOTH= value must be greater than zero.; 255 %goto exit; 256 %end; 257 %end; 258 %else %do; 259 %put ERROR: SMOOTH= is required and must be positive value.; 260 %goto exit; 261 %end; 262 263 %let validopts=ASCENDING DESCENDING DATAORDER PANEL NOPANEL; 264 %let panel=1; %let plot=0; %let hievent=1; %let dataorder=0; 265 %let i=1; 266 %do %while (%scan(&options,&i) ne %str() ); 267 %let option&i=%upcase(%scan(&options,&i)); 268 %if &&option&i=DATAORDER %then %do; 269 %let dataorder=1; %let hievent=0; 270 %end; 271 %else %if &&option&i=ASCENDING %then %let hievent=0; 272 %else %if &&option&i=NOPANEL %then %do; %let panel=0; %let plot=1; %end; 273 %else %do; 274 %let chk=%eval(&&option&i in &validopts); 275 %if not &chk %then %do; 276 %put ERROR: Valid values of OPTIONS= are &validopts..; 277 %goto exit; 278 %end; 279 %end; 280 %let i=%eval(&i+1); 281 %end; 282 283 %if %index(&version,DEBUG)=0 %then options notes;; 284 %put NOTE: Logits are computed over lower Ordered Values of &y..; 285 %if %index(&version,DEBUG)=0 %then options nonotes;; 286 287 proc sql; 288 select count(*) into :n from &data where &y is not missing; 289 quit; 290 ods graphics / loessmaxobs=%sysfunc(max(5000,%eval(&n+1000))); 291 data _tempdata; 292 set &data; 293 run; 294 295 /* Order y= variable as requested */ 296 %if &dataorder=0 %then %do; 297 proc sort data=_tempdata; 298 by %if &hievent %then descending; &y; 299 run; 300 %end; 301 %let chary=0; %let yname=&y; 302 %let dsid=%sysfunc(open(_tempdata)); 303 %if &dsid %then %do; 304 %let varnum=%sysfunc(varnum(&dsid,&y)); 305 %if %sysfunc(vartype(&dsid,&varnum))=C %then %do; 306 %let rc=%sysfunc(close(&dsid)); 307 %CtoN(&version,data=_tempdata,var=&y,out=_tempdata,order=data, 308 options=noreplace noformat nodatanote nonewcheck) 309 data _tempdata; set _tempdata(rename=(&y=_chary)); 310 label _chary='00'x; 311 run; 312 %let y=&y._N; %let chary=1; 313 %end; 314 %let rc=%sysfunc(close(&dsid)); 315 %end; 316 proc freq nlevels data=_tempdata 317 %if &hievent or &dataorder %then order=data; 318 ; 319 where not missing(&y); 320 table &y %if &chary %then * _chary; / out=_ylevs noprint; 321 ods exclude nlevels; 322 ods output nlevels=_numy; 323 run; 324 data _null_; 325 set _numy; 326 call symputx('numy',nlevels); 327 run; 328 %if &chary=0 %then %do; 329 data _null_; 330 set _ylevs nobs=nlev; 331 call symputx(cats('l',_n_),&y); 332 run; 333 data _tempdata; 334 set _tempdata; 335 select (&y); 336 %do i=1 %to &numy; 337 when (&&l&i) &y = &i; 338 %end; 339 otherwise &y = .; 340 end; 341 run; 342 %end; 343 %if &numy=2 %then %let lgtlabl=Logit; 344 %else %let lgtlabl=CLogit; 345 346 /* Variables with > 50 levels will use neighborhood method to compute logits. 347 Variables with <=50 levels will have logits computed within each level. 348 */ 349 %let cont=; %let class=; 350 proc freq data=_tempdata nlevels; 351 table &x / noprint; 352 ods exclude nlevels; 353 ods output nlevels=_nx; 354 run; 355 data _null_; 356 set _nx; 357 call symputx(cats('nlev',_n_),nnonmisslevels); 358 run; 359 %let i=1; 360 %let dsid=%sysfunc(open(_tempdata)); 361 %if &dsid %then %do %while (%scan(&x,&i) ne %str() ); 362 %let varnum=%sysfunc(varnum(&dsid,%scan(%upcase(&x),&i))); 363 %if %sysfunc(vartype(&dsid,&varnum))=N and &&nlev&i>20 %then 364 %let cont=&cont %scan(&x,&i); 365 %else %let class=&class %scan(&x,&i); 366 %let i=%eval(&i+1); 367 %end; 368 %let rc=%sysfunc(close(&dsid)); 369 370 ods select all; 371 %let num=1; 372 %let numc=0; 373 %let numd=0; 374 %let totnum=1; 375 %let doClass=1; 376 %let numym1=%eval(&numy-1); 377 %let var=%scan(&class,&num); 378 %if (&var=) %then %do; 379 %let doClass=0; 380 %let var=%scan(&cont,&num); 381 %if (&var=) %then %goto stopit; 382 %end; 383 data _ylevs; 384 set _ylevs; 385 ov=_n_; 386 run; 387 proc print data=_ylevs label noobs split="/"; 388 %if &chary=0 %then %do; 389 id ov; var &y; label ov="Ordered/Value/(O.V.)"; 390 %end; 391 %else %do; 392 id &y; var _chary; label &y="Ordered/Value/(O.V.)" _chary="&yname"; 393 %end; 394 title "Response Profile"; 395 title2 "Logits are computed over the lower Ordered Values"; 396 run; 397 title; 398 %do %while(&var ne); 399 data _tempa(keep=&y &var); 400 set _tempdata; 401 run; 402 %if %eval(&doClass eq 1) %then %do; 403 %OneClassPlot(data=_tempa,x=&var,y=&y,numy=&numy,const=&const,plot=&plot); 404 %let numc= %eval(&numc+1); 405 data _tempc; 406 set %if %eval(&numc ne 1) %then _tempc; _temp2; 407 run; 408 %end; 409 %else %do; 410 %OneContPlot(data=_tempa,x=&var,y=&y,numy=&numy,const=&const,plot=&plot, 411 nbhd=&nbhd); 412 %let numd= %eval(&numd+1); 413 data _tempd; 414 set %if %eval(&numd ne 1) %then _tempd; _temp2; 415 run; 416 %end; 417 %let num=%eval(&num+1); 418 %let totnum=%eval(&totnum+1);; 419 %if %eval(&doClass eq 1) %then %do; 420 %let var=%scan(&class,&num); 421 %if (&var=) %then %do; 422 %let doClass=0; 423 %let num=1; 424 %end; 425 %end; 426 %if %eval(&doClass eq 0) %then %do; 427 %let var=%scan(&cont,&num); 428 %if (&var=) %then %goto stopit; 429 %end; 430 %end; 431 432 %stopit: ; 433 %if &panel %then %do; 434 %if %eval(&numc>0) %then %do; 435 proc sort data=_tempc; 436 by _group _x; 437 run; 438 %end; 439 %if %eval(&numd>0) %then %do; 440 proc sort data=_tempd; 441 by _group _x; 442 run; 443 %end; 444 %if %eval(&numc>0 and &numd>0) %then %do; 445 data _tempall; 446 merge _tempd _tempc; 447 by _group _x; 448 run; 449 %end; 450 %else %if %eval(&numc>0) %then %do; 451 data _tempall; 452 set _tempc; 453 run; 454 %end; 455 %else %do; 456 data _tempall; 457 set _tempd; 458 run; 459 %end; 460 proc sort data=_tempall; 461 by _group; 462 run; 463 data _tempa(keep=_group _minx _maxx); 464 set _tempall; 465 by _group; 466 retain _minx 1e100 _maxx -1e100; 467 if first._group then do; 468 _minx=_x; 469 _maxx=_x; 470 end; 471 _minx=min(_minx,_x); 472 _maxx=max(_maxx,_x); 473 if last._group then output; 474 run; 475 data _tempall; 476 merge _tempall _tempa; 477 by _group; 478 _x=(_x-_minx)/(_maxx-_minx); 479 run; 480 %if ((&totnum-1)=1) %then %do; %let rows=1; %let cols=1; %end; 481 %else %if ((&totnum-1)=2) %then %do; %let rows=1; %let cols=2; %end; 482 %else %if ((&totnum-1)=3) %then %do; %let rows=1; %let cols=3; %end; 483 %else %if ((&totnum-1)=4) %then %do; %let rows=2; %let cols=2; %end; 484 %else %if ((&totnum-1)<=6) %then %do; %let rows=2; %let cols=3; %end; 485 %else %if ((&totnum-1)<=9) %then %do; %let rows=3; %let cols=3; %end; 486 %else %if ((&totnum-1)<=12) %then %do; %let rows=3; %let cols=4; %end; 487 %else %do; %let cols=4; %let rows=4; %end; 488 proc sgpanel data=_tempall; 489 panelby _group / novarname columns=&cols rows=&rows 490 uniscale=all skipemptycells; 491 rowaxis label= 492 %if &numy>2 %then "Empirical Cumulative Logits"; 493 %else "Empirical Logit"; 494 ; 495 colaxis display=(nolabel noticks novalues); 496 %if (&numd>0) %then %do i=1 %to &numym1; 497 loess y=d&i x=_x / smooth=&smooth legendlabel="&lgtlabl(O.V.<=&i)" 498 name="series&i" lineattrs=GraphData&i(thickness=3px) 499 markerattrs=GraphDataDefault(color=gray) nomarkers 500 ; 501 %end; 502 %if (&numc>0) %then %do i=1 %to &numym1; 503 series y=c&i x=_x / legendlabel="&lgtlabl(O.V.<=&i)" name="series&i" 504 lineattrs=GraphData&i(thickness=3px); 505 %end; 506 %if &numd %then inset neighbors / position=bottom;; 507 keylegend %do i=1 %to &numym1; "series&i" %end;; 508 %if &numy>2 %then %str(title "Empirical Cumulative Logits";); 509 %else %str(title "Empirical Logit";); 510 %if &numd %then 511 title2 "Smoothing parameter = &smooth for continuous predictors";; 512 run; 513 %end; 514 515 %exit: 516 %if %index(&version,DEBUG)=0 %then %do; 517 proc datasets nolist; 518 delete _: ; 519 run; quit; 520 %end; 521 %if %index(&version,DEBUG) %then %do; 522 options nomprint nomlogic nosymbolgen; 523 %put _user_; 524 %end; 525 options &_elpopts; 526 title; 527 %let elptime=%sysfunc(round(%sysevalf(%sysfunc(datetime()) - &elptime), 0.01)); 528 %put NOTE: The &sysmacroname macro used &elptime seconds.; 529 %mend; 530 531 %macro OneClassPlot(data=,x=,y=,numy=,const=0.5,plot=1); 532 /*-------------------------------------------------------------------------- 533 A macro version of SAS Usage Note 37944 for covariates with fewer than 534 CONTCUTOFF levels. Computes and/or plots the empirical logits of Y against 535 the levels of one X variable. Called by the EmpiricalLogitPlot macro. 536 537 data= name of the input data set. 538 x= name of the single CLASS covariate. 539 y= name of the response variable. 540 numy= number of levels of the response. 541 const= small value in case of zero cells. 542 plot= 1=create the plot, 0=create data set only 543 --------------------------------------------------------------------------*/ 544 %let numym1=%eval(&numy-1); 545 %let zerosum=0; 546 proc freq data=&data noprint; 547 table &x*&y / sparse out=_temp2; 548 run; 549 proc sort data=_temp2; 550 by &x; 551 run; 552 %let dsid=%sysfunc(open(_temp2)); 553 %if &dsid %then %do; 554 %let varnum=%sysfunc(varnum(&dsid,&x)); 555 %if %sysfunc(vartype(&dsid,&varnum))=C %then %do; 556 %let rc=%sysfunc(close(&dsid)); 557 %CtoN(&version,data=_temp2,var=&x,out=_temp2,options=format nodatanote nonewcheck) 558 %end; 559 %let rc=%sysfunc(close(&dsid)); 560 %end; 561 proc transpose data=_temp2 out=_temp; 562 by &x; 563 var count; 564 run; 565 data _temp2 ; 566 set _temp; 567 length _group $ 32; 568 if (_NAME_='COUNT'); 569 _group="&x"; 570 _x=&x; 571 _const=%sysevalf(&const+0); 572 %do i=1 %to &numym1; %let j=%eval(&i+1); 573 _numsum=sum(of col1-col&i); _densum=sum(of col&j-col&numy); 574 %if &const=0 %then %do; 575 _const=0; 576 if _numsum=0 or _densum=0 then do; 577 _const=0.5; call symputx('zerosum',1); 578 end; 579 %end; 580 c&i=log((_numsum + _const) / (_densum + _const)); 581 %end; 582 keep c1-c&numym1 _group _x &x; 583 run; 584 %if &zerosum %then %do; 585 %if %index(&version,DEBUG)=0 %then options notes;; 586 %put NOTE: Zero sum detected for &x.. CONST=0.5 used in affected logits.; 587 %if %index(&version,DEBUG)=0 %then options nonotes;; 588 %end; 589 %if &plot %then %do; 590 proc sgplot; 591 %do i=1 %to &numym1; 592 series y=c&i x=&x / 593 legendlabel="&lgtlabl(O.V.<=&i)" name="series&i" 594 lineattrs=GraphData&i(thickness=3px); 595 %end; 596 yaxis label= 597 %if &numy>2 %then "Empirical Cumulative Logits"; 598 %else "Empirical Logit"; 599 ; 600 xaxis label="&x"; 601 xaxis integer label="&x"; 602 keylegend %do i=1 %to &numym1; "series&i" %end;; 603 %if &numy>2 %then %str(title "Empirical Cumulative Logits";); 604 %else %str(title "Empirical Logit";); 605 title2; 606 run; 607 %end; 608 %mend; 609 610 %macro OneContPlot(data=,x=,y=,numy=,const=0.5,smooth=.3,plot=1,nbhd=); 611 /*-------------------------------------------------------------------------- 612 A macro extension of SAS Usage Note 37944 for continuous covariates. Computes 613 and/or plots the empirical logits defined using the NBHD/2 observed X values 614 preceding and following each observation against the levels of a single 615 continuous X variable, then smooths the result. Called by the 616 EmpiricalLogitPlot macro. 617 618 data= name of the input data set. 619 x= name of the single continuous covariate. 620 y= name of the response variable. y=1,2,3... 621 numy= number of levels of the response. 622 const= small value in case of zero cells. 623 nbhd= total number of neighbors. 624 smooth= smoothing parameter used for loess curve fit to logits for 625 continuous effects. Default: 0.3. 626 plot= 1=create the plot, 0=create data set only 627 --------------------------------------------------------------------------*/ 628 %let zerosum=0; %let nbhderr=0; 629 data _temp; 630 set &data; 631 run; 632 proc sort data=_temp; 633 by &x &y; 634 run; 635 proc sql noprint; 636 select count(*) into :nyx from &data 637 where &y is not missing and &x is not missing; 638 quit; 639 data _null_; 640 set _nx; 641 maxnbhd=max(ceil(sqrt(&nyx)),ceil(.01*&nyx)); 642 if &nbhd>maxnbhd then call symputx("nbhderr",maxnbhd); 643 run; 644 %if &nbhderr %then %do; 645 %if %index(&version,DEBUG)=0 %then options notes;; 646 %put NOTE: NEIGHBORS=&nbhd decreased to &nbhderr for X=&x..; 647 %if %index(&version,DEBUG)=0 %then options nonotes;; 648 %let nbhd=&nbhderr; %let neighbors=&nbhd; 649 %end; 650 %let Full=%eval(&nbhd); 651 %let Fullp1=%eval(&nbhd+1); 652 %let Half=%eval(&nbhd/2); 653 %let numym1=%eval(&numy-1); 654 data _temp2; 655 set _temp; 656 retain %do i=1 %to &numy; _sum&i %end; 0; 657 _lagx=lag&Full(&x); 658 _lagHalfx=lag&Half(&x); 659 _lagy=lag&Fullp1(&y); 660 %do i=1 %to &numy; _sum&i=_sum&i+(&y=&i); %end; 661 if (_lagy^=.) then do; 662 %do i=1 %to &numy; _sum&i=_sum&i-(_lagy=&i); %end; 663 end; 664 run; 665 data _temp; 666 set _temp2; 667 if (_n_<&Full) then delete; 668 run; 669 data _temp2; 670 set _temp; 671 length _group $ 32; 672 _const=%sysevalf(&const+0); 673 _x=_lagHalfx; 674 _group="&x"; 675 %do i=1 %to &numym1; %let j=%eval(&i+1); 676 _numsum=sum(of _sum1-_sum&i); _densum=sum(of _sum&j-_sum&numy); 677 %if &const=0 %then %do; 678 _const=0; 679 if _numsum=0 or _densum=0 then do; 680 _const=0.5; call symputx('zerosum',1); 681 end; 682 %end; 683 d&i=log((_numsum + _const) / (_densum + _const)); 684 %end; 685 Neighbors=&nbhd; 686 keep d1-d&numym1 _group _x neighbors; 687 run; 688 %if &zerosum %then %do; 689 %if %index(&version,DEBUG)=0 %then options notes;; 690 %put NOTE: Zero sum detected for &x.. CONST=0.5 used in affected logits.; 691 %if %index(&version,DEBUG)=0 %then options nonotes;; 692 %end; 693 %if &plot %then %do; 694 proc sgplot; 695 %do i=1 %to &numym1; 696 loess y=d&i x=_x / 697 smooth=&smooth legendlabel="&lgtlabl(O.V.<=&i)" 698 lineattrs=GraphData&i(thickness=3px) 699 markerattrs=(color=gray) name="series&i" nomarkers 700 ; 701 %end; 702 yaxis label= 703 %if &numy>2 %then "Empirical Cumulative Logits"; 704 %else "Empirical Logit"; 705 ; 706 xaxis label="&x"; 707 keylegend %do i=1 %to &numym1; "series&i" %end;; 708 %if &numy>2 %then %str(title "Empirical Cumulative Logits";); 709 %else %str(title "Empirical Logit";); 710 title2 "with &nbhd neighbors and smoothing parameter = &smooth"; 711 run; 712 %end; 713 %mend; 714 715 716 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 726 when I run it. ... View more