The sas-supplied maps.us is a projected map, containing only the projected X & Y values. This is an old-old legacy map that has been shipped with SAS a long time, and the original projection parameters are 'unknown', and therefore trying to project another map the exact same way (to combine them) is not really possible. (You might come close, using trial-and-error ... but it's just not a "good practice".)   I would recommend using the newer mapsgfk.us instead, which contains unprojected lat/long values.   data newus; set mapsgfk.us_states (where=(statecode not in ('AK' 'HI')) drop = x y); run;   Then if the imported shapefile has unprojected lat/long values, you can project that and your newus any way you want 🙂   Or, if your shapefile has projected values, and you know what proj.4 projection was used, you can use that same proj.4 projection to project newus the same, or you can use that proj.4 to 'unproject' your shapefile coordinates to get the unprojected lat/long (for example, using gproject parameters such as ... from="EPSG:2264" to="EPSG:4326") . And then once you have both maps in unprojected lat/long, you can project both maps with whatever projection you want.   My rule-of-thumb is start by getting everything in unprojected lat/long 🙂       ... View more

I am not allowed to post any of our data, but perhaps these additional details will be helpful: The exposure variable, pdc2 is a categorical variable with 3 levels; lenfol is the time to event in days, and death365day is a Y/N flag for the event. My question is centered on the difference between the output I receive after I add the weight statement to the code below, which runs without error and produces the standard output for proc lifetest, as well as the Adjustment for Multiple Comparisons for the Logrank Test table. %dummytbl(sheet_name=WtdKMyr_Longterm66743); proc lifetest data= bb_ps_w_pdc2 notable plots=survival( nocensor failure test atrisk(atrisktickonly outside)=(0,30,90,180,365)) maxtime=365; title 'IPW Failure curve for ages 66-74 by PDC Group- weight removed'; where 66 <= ageatdiag < 75; strata pdc2 / test=logrank adjust=sidak ; time lenfol*death365day(0); run;  The Multiple Comparison table gave me the strata comparison, chi-square test results and the p-Values (raw and Sidak). However, when I added "weight" statement (in order to generate inverse probability weighted curves), as shown in the code below:  proc lifetest data= bb_pswt_pdc_adj2 notable plots=survival(cl nocensor failure test atrisk(atrisktickonly outside)=(0,30,90,180,365)) maxtime=365; title 'IPW Failure curve by age2 strata pdc with adjusted weight'; by age2; strata pdc2 / test=logrank adjust=sidak ; time lenfol*death365day(0); weight ps_wt_adj; run; the Multiple comparison table was not generated. The only output I got were the Rank Statistics, Covariance Matrix, and Test of Equality over Strata tables, along with the plot. I did not receive any error and the log output was the same for both sets of code and was the following: NOTE: The LOGLOG transform is used to compute the confidence limits for the survivor function. To suppress using this transform, specify CONFTYPE=LINEAR in the PROC LIFETEST statement. NOTE: The PROCEDURE LIFETEST printed pages 13-18. NOTE: PROCEDURE LIFETEST used (Total process time):   Thanks. Barbara     ... View more