cancel
Turn on suggestions
Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.
Showing results for 
Search instead for 
Did you mean: 
Tiny_Kane
‎07-21-2020
  • 77 Posts
  • 2 Likes Given
  • 0 Solutions
  • 2 Likes Received

how to add condition at select statement in proc sql

by Obsidian | Level 7 in SAS Data Management
‎04-10-2017 10:07 AM
‎04-10-2017 10:07 AM
Hi there,   I have to two data sets look like below, and the patient ID is the unique key linking those two sets. I want a third dataset, which give me the patient ID, the minimun of treatment date which is 6 month prior to the index_date, the minimun of the treatment date which is within 6-month before and 6-month after the index_date, and the minumun of the treatment date which is 6-month after the index_date. If any type of minimun of treatment date is not existing, for instance, there is no treatment date 6-month prior to the index_date, then this value will be misisng. How can I write the query?  Thank you very much!   patient ID index_date 1 3/4/2009 2 8/16/2010 3 6/6/2011 4 7/4/2009 … … patient ID treatment_data  1 9/30/2008 1 10/11/2008 1 4/1/2009 2 1/1/2010 2 9/30/2010 3 12/30/2009 3 5/1/2010 3 8/1/2011 3 10/13/2011 3 12/13/2011 4 6/8/2009 4 3/1/2010 … … ... View more

Re: how to add rows to the data set?

by Obsidian | Level 7 in SAS Data Management
‎12-29-2016 01:29 PM
‎12-29-2016 01:29 PM
Thank you!! ... View more

Re: how to add rows to the data set?

by Obsidian | Level 7 in SAS Data Management
‎12-29-2016 01:28 PM
‎12-29-2016 01:28 PM
Thank you! ... View more

Re: how to add rows to the data set?

by Obsidian | Level 7 in SAS Data Management
‎12-29-2016 01:27 PM
‎12-29-2016 01:27 PM
Thanks! It works. ... View more

how to add rows to the data set?

by Obsidian | Level 7 in SAS Data Management
‎12-28-2016 02:36 PM
‎12-28-2016 02:36 PM
Hi there,   I have an existing data set looks like the first listing below. How to add rows to the data set to make it look like the second listing? All red texts were added rows. The second listing will be, each patient will have death month+1 rows, index month shoul be from 0 to the death month. All added costs are 0. Thanks in advance!   patient ID death month index month cost 1 3 0 20 1 3 1 10 1 3 3 20 2 5 2 10 2 5 3 5 2 5 4 10   patient ID death month index month cost 1 3 0 20 1 3 1 10 1 3 2 0 1 3 3 20 2 5 0 0 2 5 1 0 2 5 2 10 2 5 3 5 2 5 4 10 2 5 5 0 ... View more

Re: How to change intersetgroup position in boxplot

by Obsidian | Level 7 in Graphics Programming
‎07-08-2016 05:06 PM
‎07-08-2016 05:06 PM
Hi PG,   If I omit the option, the result is the same. Therefore, position= option does not work.    I found the reason actually, I forgot to add 'ods graphics off' at the very begining. I guess that's the reason why the option does not work...It was silly lol~ ... View more

How to change intersetgroup position in boxplot

by Obsidian | Level 7 in Graphics Programming
‎07-07-2016 11:59 AM
‎07-07-2016 11:59 AM
Hi there, I tried to interset a table at above of the boxplot, but the position option within the interstgroup statement in the PROC BOXPLOT does not work. Please see the code and the result below.    proc boxplot data=temp_base; plot calculated_household_in_baseline*value/; insetgroup q1 q2 /header='test' pos=top; run;       ... View more

Re: How to change interval in x-axis for KM plot

by Obsidian | Level 7 in Graphics Programming
‎03-03-2016 10:33 AM
‎03-03-2016 10:33 AM
Hi Reeza,   You are right, my template name is Stat.Lifetest.Graphics.Productlimitsurvival2. After I updated the name, the code works. But now, the list of number at risk is within the plot..See the attachement. If I want it to be outside of the plot, like below the x-axis, could you mind telling me from which line of your code that I can start to revise. Thank you very much!     ... View more

Re: How to change interval in x-axis for KM plot

by Obsidian | Level 7 in Graphics Programming
‎03-02-2016 04:54 PM
‎03-02-2016 04:54 PM
Hi Reeza,   I am using SAS 9.3. See the code in log. Thanks! 101  proc template; 102  define statgraph Stat.Lifetest.Graphics.Productlimitsurvival; 103     dynamic NStrata xName plotAtRisk plotCensored plotCL plotHW plotEP labelCL labelHW labelEP 103!  maxTime method 104        StratumID classAtRisk plotBand plotTest GroupName yMin Transparency SecondTitle 104! TestName pValue; 105 106      /*My Macro Variables*/ 107      mvar legend_label; 108     BeginGraph; 109        if (NSTRATA=1) 110           if (EXISTS(STRATUMID)) 111           entrytitle "Kaplan-Meier " "Survival Estimate" " for " STRATUMID; 112        else 113           entrytitle "Kaplan-Meier " "Survival Estimate"; 114        endif; 115        if (PLOTATRISK) 116           entrytitle "with Number of Subjects at Risk" / textattrs=GRAPHVALUETEXT; 117        endif; 118        layout overlay / xaxisopts=(shortlabel=XNAME offsetmin=.05 linearopts=(viewmax=MAXTIME 118! /*Modifying X axis values here, need to change below as well*/ tickvaluelist=(0 12 24 36 48 118! 60 72 84 96 108 120))) yaxisopts=(label= 119           "Survival Probability" shortlabel="Survival" linearopts=(viewmin=0 viewmax=1 119! tickvaluelist=(0 .25 .5 .75 120           1.0) tickvalueformat=percent8.1 /*Modifying Y axis values here, need to change below 120!  as well*/)); 121 122           /*Add a reference line for the media*/ 123           *referenceline y=0.5/lineattrs=(color=red pattern=2); 124           if (PLOTHW=1 AND PLOTEP=0) 125              bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / modelname="Survival" 125! fillattrs=GRAPHCONFIDENCE 126              name="HW" legendlabel=LABELHW; 127           endif; 128           if (PLOTHW=0 AND PLOTEP=1) 129              bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / modelname="Survival" 129! fillattrs=GRAPHCONFIDENCE 130              name="EP" legendlabel=LABELEP; 131           endif; 132           if (PLOTHW=1 AND PLOTEP=1) 133              bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / modelname="Survival" 133! fillattrs=GRAPHDATA1 134              datatransparency=.55 name="HW" legendlabel=LABELHW; 135           bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / modelname="Survival" 135! fillattrs=GRAPHDATA2 136              datatransparency=.55 name="EP" legendlabel=LABELEP; 137           endif; 138           if (PLOTCL=1) 139              if (PLOTHW=1 OR PLOTEP=1) 140              bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / modelname="Survival" 140! display=(outline) 141              outlineattrs=GRAPHPREDICTIONLIMITS name="CL" legendlabel=LABELCL; 142           else 143              bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / modelname="Survival" 143! fillattrs=GRAPHCONFIDENCE 144              name="CL" legendlabel=LABELCL; 145           endif; 146           endif; 147           stepplot y=SURVIVAL x=TIME / name="Survival" rolename=(_tip1=ATRISK _tip2=EVENT) 147! tip=(y x Time _tip1 148              _tip2) legendlabel="Survival" lineattrs=(thickness=2); 149           if (PLOTCENSORED=1) 150              scatterplot y=CENSORED x=TIME / markerattrs=(symbol=plus size=8px) 150! name="Censored" legendlabel="Censored"; 151           endif; 152           if (PLOTCL=1 OR PLOTHW=1 OR PLOTEP=1) 153              discretelegend "Censored" "CL" "HW" "EP" / location=outside halign=center; 154           else 155              if (PLOTCENSORED=1) 156              discretelegend "Censored" / location=inside autoalign=(topright bottomleft); 157           endif; 158           endif; 159           if (PLOTATRISK=1) 160              innermargin / align=bottom; 161              blockplot x=TATRISK block=ATRISK / repeatedvalues=true display=(values) 161! valuehalign=start 162                 valuefitpolicy=truncate labelposition=left labelattrs=GRAPHVALUETEXT 162! valueattrs=GRAPHDATATEXT (size 163                 =7pt) includemissingclass=false ; 164           endinnermargin; 165           endif; 166        endlayout; 167        else 168           entrytitle "Kaplan-Meier " "Survival Estimates"; 169        layout overlay / xaxisopts=(shortlabel=XNAME offsetmin=.05 linearopts=(viewmax=MAXTIME 169! tickvaluelist=(0 12 24 36 48 60 72 84 96 108 120) /*Modifying the X values on graph here, 169! need to change above as well*/)) yaxisopts=(label= 170           "Survival Probability" shortlabel="Survival" linearopts=(viewmin=0 viewmax=1 170! tickvaluelist=(0 .25 .5 .75 171           1.0) tickvalueformat=percent8.1 /*Modifying the Y Axis values here, need to change 171! above as well*/)); 172 173           /*Add a reference line for the media*/ 174           *referenceline y=0.5/lineattrs=(color=red pattern=2); 175           if (PLOTHW) 176              bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / group=STRATUM 176! index=STRATUMNUM modelname= 177              "Survival" datatransparency=Transparency; 178           endif; 179           if (PLOTEP) 180              bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / group=STRATUM 180! index=STRATUMNUM modelname= 181              "Survival" datatransparency=Transparency; 182           endif; 183           if (PLOTCL) 184              if (PLOTBAND) 185              bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / group=STRATUM 185! index=STRATUMNUM modelname= 186              "Survival" display=(outline); 187           else 188              bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / group=STRATUM 188! index=STRATUMNUM modelname= 189              "Survival" datatransparency=Transparency; 190           endif; 191           endif; 192           stepplot y=SURVIVAL x=TIME / group=STRATUM index=STRATUMNUM name="Survival" 192! rolename=(_tip1=ATRISK _tip2= 193              EVENT) tip=(y x Time _tip1 _tip2) lineattrs=(thickness=2); 194           if (PLOTCENSORED) 195              scatterplot y=CENSORED x=TIME / group=STRATUM index=STRATUMNUM 195! markerattrs=(symbol=plus size=8px); 196           endif; 197           if (PLOTATRISK) 198              innermargin / align=bottom; 199              blockplot x=TATRISK block=ATRISK / class=STRATUM repeatedvalues=true 199! display=(label values) 200                 valuehalign=start valuefitpolicy=truncate labelposition=left 200! labelattrs=STRATUMID valueattrs= 201                 GRAPHDATATEXT (size=7pt) includemissingclass=false; 202           endinnermargin; 203           endif; 204           DiscreteLegend "Survival" / title=legend_label location=outside; 205           if (PLOTCENSORED) 206              if (PLOTTEST) 207              layout gridded / rows=2 autoalign=(TOPRIGHT BOTTOMLEFT TOP BOTTOM) border=false 207! BackgroundColor= 208              GraphWalls:Color Opaque=true; 209              entry "+ Censored"; 210              if (PVALUE < .0001) 211                 entry TESTNAME " p " eval (PUT(PVALUE, PVALUE6.4)); 212              else 213                 entry TESTNAME " p=" eval (PUT(PVALUE, PVALUE6.4)); 214              endif; 215           endlayout; 216           else 217              layout gridded / rows=1 autoalign=(TOPRIGHT BOTTOMLEFT TOP BOTTOM) border=false 217! BackgroundColor= 218              GraphWalls:Color Opaque=true; 219              entry "+ Censored"; 220           endlayout; 221           endif; 222           else 223              if (PLOTTEST) 224              layout gridded / rows=1 autoalign=(TOPRIGHT BOTTOMLEFT TOP BOTTOM) border=false 224! BackgroundColor= 225              GraphWalls:Color Opaque=true; 226              if (PVALUE < .0001) 227                 entry TESTNAME " p " eval (PUT(PVALUE, PVALUE6.4)); 228              else 229                 entry TESTNAME " p=" eval (PUT(PVALUE, PVALUE6.4)); 230              endif; 231           endlayout; 232           endif; 233           endif; 234        endlayout; 235        endif; 236     EndGraph; 237  end; NOTE: Overwriting existing template/link: Stat.Lifetest.Graphics.Productlimitsurvival NOTE: STATGRAPH 'Stat.Lifetest.Graphics.Productlimitsurvival' has been saved to: SASUSER.TEMPLAT 238  run; NOTE: PROCEDURE TEMPLATE used (Total process time):       real time           0.09 seconds       cpu time            0.09 seconds 239  proc lifetest data=master1 intervals= 0 to 60 by 12 240                             plots=survival(test atrisk(outside)=0 to 60 by 12); 241    by period; 242    strata location; 243    time surv_time*Censor(0); 244  run; NOTE: The LOGLOG transform is used to compute the confidence limits for the quartiles of the       survivor distribution. To suppress using this transform, specify CONFTYPE=LINEAR in the       PROC LIFETEST statement. NOTE: PROCEDURE LIFETEST used (Total process time):       real time           1.85 seconds       cpu time            0.85 seconds ... View more

Re: How to change interval in x-axis for KM plot

by Obsidian | Level 7 in Graphics Programming
‎03-02-2016 03:46 PM
‎03-02-2016 03:46 PM
Hi Reeza,   Thanks for sharing. However, this still does not work...It is very very weird that PROC TEMPLATE does not work the plot on my end...     ... View more

Re: How to change interval in x-axis for KM plot

by Obsidian | Level 7 in Graphics Programming
‎03-02-2016 03:39 PM
‎03-02-2016 03:39 PM
Hi Rick,   SAS gave me error message after I try this macro...see below 138  %ProvideSurvivalMacros      -      180 WARNING: Apparent invocation of macro PROVIDESURVIVALMACROS not resolved. ERROR 180-322: Statement is not valid or it is used out of proper order. 139 140  %let xOptions = label="Survival Time" 141                  linearopts=(viewmin=0 viewmax=60 142                              tickvaluelist=(0 12 24 36 48 60)); 143 144  %CompileSurvivalTemplates WARNING: Apparent invocation of macro COMPILESURVIVALTEMPLATES not resolved. ... View more

How to change interval in x-axis for KM plot

by Obsidian | Level 7 in Graphics Programming
‎03-02-2016 03:17 PM
‎03-02-2016 03:17 PM
Hi there,   I have been struggling with this problem for a while.   I run the PROC LIFETEST procedure below to get the KM plot. See the plot in the attachement. However, you will find that the time points along the x-axis are 0, 10, 20, 30, 40, 50, and 60. But I want the time interval to be 12, which can made the time points to be 0,12,24,36,48, and 60, so that the time points along x-axis is same as the time points for the number of at risk listing...I tried many options, even trying to revise PROC TEMPLATE. But none of them works......Does any one have resolved similar problem? Thanks in advance.   proc lifetest data=master1 intervals= 0 to 60 by 12                            plots=survival(test atrisk(outside)=0 to 60 by 12);   by period;   strata location;   time surv_time*Censor(0); run;     ... View more

Re: Why does the second symbol statement not have the effect?

by Obsidian | Level 7 in Graphics Programming
‎06-22-2015 04:36 PM
‎06-22-2015 04:36 PM
It works, thank you very much, Dan. ... View more

Why does the second symbol statement not have the effect?

by Obsidian | Level 7 in Graphics Programming
‎06-22-2015 04:21 PM
‎06-22-2015 04:21 PM
Hi all, I just want to create an overlay plot. I defined the first plot symbol as square and the second one as dot, however the result shows both of them are squares. I attache the code and the plot below. Does any one know why the SYMBOL2 statement does not work? Thank you in advance! goptions reset=all; symbol1 interpol=join         value=square; symbol2 interpol=join         value=dot         ; axis1 order=(1 to 2 by 1) offset=(2,2) label=("Cycle"); axis2 label=("Participants"); legend1 label=none position=(top center inside) mode=share; proc gplot data=act_exp1; plot actual*cycle expect*cycle/overlay legend=legend1                                vref= 5 to 35 by 5                                haxis=axis1                                vaxis=axis2; run; Kane ... View more

Re: A question about _error_ automatic variable

by Obsidian | Level 7 in SAS Procedures
‎12-11-2014 02:02 PM
‎12-11-2014 02:02 PM
By the way, yes, I want to do debugging. It really takes me a lot of time to search for debugging methods. Thanks. ... View more
Contact Me
Online Status
Offline
Date Last Visited
‎07-21-2020 06:41 PM
Likes given to
View All
Likes from
View All
Quick links