SAS Support Communities

I have a dependent variable - y  and independent vars -  cate has two levels - A and B  gender has two levels Male and female   adh has two levels - 0 and 1  proc glm data=change ; class cate gender ; model y= class sex age adh adh*class /solution; lsmeans cate /pdiff e; lsmeans cate*adh / pdiff e; estimate 'A vs B' cate 1 -1/e ; run; want : 1. Difference in y between class a and class b ( adjusted for sex and age adh ) 2. Difference in y between class a and class b split by adh ( adjusted for sex and age)    which statement would give me these results ? 1. lsmeans ?  2. estimate statement ?  3. solution option - parameter estimates ?     ... View more

The class variable has 1 and 2 in front of it , how do i remove it?   proc template ; define statgraph Stat.Lifetest.Graphics.ProductLimitSurvival2; dynamic NStrata xName plotAtRisk plotCensored plotCL plotHW plotEP labelCL labelHW labelEP maxTime xtickVals xtickValFitPol rowWeights method StratumID classAtRisk plotTest GroupName Transparency SecondTitle TestName pValue _byline_ _bytitle_ _byfootnote_; BeginGraph; if (NSTRATA=1) if (EXISTS(STRATUMID)) entrytitle &titletext1; else entrytitle &titletext0; endif; entrytitle "xyz2" / textattrs=GRAPHVALUETEXT; layout lattice / rows=2 rowweights=ROWWEIGHTS columndatarange=union rowgutter=10; cell; layout overlay / xaxisopts=(shortlabel=XNAME offsetmin=0 linearopts=(viewmax=MAXTIME tickvaluelist=XTICKVALS tickvaluefitpolicy=XTICKVALFITPOL)) yaxisopts=(label="&yaxis" shortlabel="Survival" offsetmin=0 linearopts=( viewmin=0 viewmax=1 tickvaluelist=(0 .2 .4 .6 .8 1.0))); if (PLOTHW=1 AND PLOTEP=0) bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / displayTail=false modelname="Survival" fillattrs= GRAPHCONFIDENCE name="HW" legendlabel=LABELHW; endif; if (PLOTHW=0 AND PLOTEP=1) bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / displayTail=false modelname="Survival" fillattrs= GRAPHCONFIDENCE name="EP" legendlabel=LABELEP; endif; if (PLOTHW=1 AND PLOTEP=1) bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / displayTail=false modelname="Survival" fillattrs= GRAPHDATA1 datatransparency=.55 name="HW" legendlabel=LABELHW; bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / displayTail=false modelname="Survival" fillattrs= GRAPHDATA2 datatransparency=.55 name="EP" legendlabel=LABELEP; endif; if (PLOTCL=1) if (PLOTHW=1 OR PLOTEP=1) bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / displayTail=false modelname="Survival" display=( outline) outlineattrs=GRAPHPREDICTIONLIMITS name="CL" legendlabel=LABELCL; else bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / displayTail=false modelname="Survival" fillattrs= GRAPHCONFIDENCE name="CL" legendlabel=LABELCL; endif; endif; stepplot y=SURVIVAL x=TIME / name="Survival" rolename=(_tip1=ATRISK _tip2=EVENT) tiplabel=(_tip1="Number at Risk" _tip2="Observed Events") tip=(x y _tip1 _tip2) legendlabel="Survival"; if (PLOTCENSORED=1) scatterplot y=CENSORED x=TIME / markerattrs=(symbol=plus) tiplabel=(y="Survival Probability") name="Censored" legendlabel="Censored"; endif; if (PLOTCL=1 OR PLOTHW=1 OR PLOTEP=1) discretelegend "Censored" "CL" "HW" "EP" / location=outside halign=center; else if (PLOTCENSORED=1) discretelegend "Censored" / location=inside autoalign=(topright bottomleft); endif; endif; endlayout; endcell; cell; layout overlay / walldisplay=none xaxisopts=(display=none); axistable x=TATRISK value=ATRISK / display=(label) valueattrs=(size=7pt); drawtext 'Number at Risk' / x=3 y=20 width=9; endlayout; endcell; endlayout; else entrytitle &titletext1; if (EXISTS(SECONDTITLE)) entrytitle SECONDTITLE / textattrs=GRAPHVALUETEXT; endif; layout lattice / rows=2 rowweights=ROWWEIGHTS columndatarange=union rowgutter=10; cell; layout overlay / xaxisopts=(shortlabel=XNAME offsetmin=0 linearopts=(viewmax=MAXTIME tickvaluelist=XTICKVALS tickvaluefitpolicy=XTICKVALFITPOL)) yaxisopts=(label="&yaxis." shortlabel="Survival" offsetmin=0 linearopts=( viewmin=0 viewmax=1 tickvaluelist=(0 .2 .4 .6 .8 1.0))); if (PLOTHW=1) bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / displayTail=false group=STRATUM index=STRATUMNUM modelname="Survival" datatransparency=Transparency; endif; if (PLOTEP=1) bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / displayTail=false group=STRATUM index=STRATUMNUM modelname="Survival" datatransparency=Transparency; endif; if (PLOTCL=1) if (PLOTHW=1 OR PLOTEP=1) bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / displayTail=false group=STRATUM index=STRATUMNUM modelname="Survival" display=(outline) outlineattrs=(pattern=ShortDash); else bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / displayTail=false group=STRATUM index=STRATUMNUM modelname="Survival" datatransparency=Transparency; endif; endif; stepplot y=SURVIVAL x=TIME / group=STRATUM index=STRATUMNUM name="Survival" rolename=(_tip1=ATRISK _tip2=EVENT) tiplabel=(_tip1="Number at Risk" _tip2="Observed Events") tip=(x y _tip1 _tip2); if (PLOTCENSORED=1) scatterplot y=CENSORED x=TIME / group=STRATUM index=STRATUMNUM tiplabel=(y="Survival Probability") markerattrs =(symbol=plus); endif; DiscreteLegend "Survival" / title=GROUPNAME location=outside ; if (PLOTCENSORED=1) if (PLOTTEST=1) layout gridded / rows=2 autoalign=(TOPRIGHT ) border=false BackgroundColor= GraphWalls:Color Opaque=true; entry "+ Censored"; if (PVALUE < .0001) entry TESTNAME " p " eval (PUT(PVALUE, PVALUE6.4)); else entry TESTNAME " p=" eval (PUT(PVALUE, PVALUE6.4)); endif; endlayout; else layout gridded / rows=1 autoalign=(TOPRIGHT ) border=false BackgroundColor= GraphWalls:Color Opaque=true; entry "+ Censored"; endlayout; endif; else if (PLOTTEST=1) layout gridded / rows=1 autoalign=(TOPRIGHT BOTTOMLEFT TOP BOTTOM) border=true BackgroundColor= GraphWalls:Color Opaque=true; if (PVALUE < .0001) entry TESTNAME " p " eval (PUT(PVALUE, PVALUE6.4)); else entry TESTNAME " p=" eval (PUT(PVALUE, PVALUE6.4)); endif; endlayout; endif; endif; endlayout; endcell; cell; layout overlay / walldisplay=none xaxisopts=(display=none); axistable x=TATRISK value=ATRISK / display=(label) valueattrs=(size=7pt) class=CLASSATRISK colorgroup=CLASSATRISK ; drawtext 'Number at Risk' / x=3 y=55 width=9 rotate= 90; endlayout; endcell; endlayout; endif; if (_BYTITLE_) entrytitle _BYLINE_ / textattrs=GRAPHVALUETEXT; else if (_BYFOOTNOTE_) entryfootnote halign=left _BYLINE_; endif; endif; EndGraph; end; which part of the template code do I need to change for this ? ... View more

I want to shift the x axis to origin    proc template ; define statgraph Stat.Lifetest.Graphics.ProductLimitSurvival2; dynamic NStrata xName plotAtRisk plotCensored plotCL plotHW plotEP labelCL labelHW labelEP maxTime xtickVals xtickValFitPol rowWeights method StratumID classAtRisk plotTest GroupName Transparency SecondTitle TestName pValue _byline_ _bytitle_ _byfootnote_; BeginGraph; if (NSTRATA=1) if (EXISTS(STRATUMID)) entrytitle &titletext1; else entrytitle &titletext0; endif; entrytitle "xyz2" / textattrs=GRAPHVALUETEXT; layout lattice / rows=2 rowweights=ROWWEIGHTS columndatarange=union rowgutter=10; cell; layout overlay / xaxisopts=(shortlabel=XNAME offsetmin=0 linearopts=(viewmax=MAXTIME tickvaluelist=XTICKVALS tickvaluefitpolicy=XTICKVALFITPOL)) yaxisopts=(label="&yaxis" shortlabel="Survival" offsetmin=0 linearopts=( viewmin=0 viewmax=1 tickvaluelist=(0 .2 .4 .6 .8 1.0))); if (PLOTHW=1 AND PLOTEP=0) bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / displayTail=false modelname="Survival" fillattrs= GRAPHCONFIDENCE name="HW" legendlabel=LABELHW; endif; if (PLOTHW=0 AND PLOTEP=1) bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / displayTail=false modelname="Survival" fillattrs= GRAPHCONFIDENCE name="EP" legendlabel=LABELEP; endif; if (PLOTHW=1 AND PLOTEP=1) bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / displayTail=false modelname="Survival" fillattrs= GRAPHDATA1 datatransparency=.55 name="HW" legendlabel=LABELHW; bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / displayTail=false modelname="Survival" fillattrs= GRAPHDATA2 datatransparency=.55 name="EP" legendlabel=LABELEP; endif; if (PLOTCL=1) if (PLOTHW=1 OR PLOTEP=1) bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / displayTail=false modelname="Survival" display=( outline) outlineattrs=GRAPHPREDICTIONLIMITS name="CL" legendlabel=LABELCL; else bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / displayTail=false modelname="Survival" fillattrs= GRAPHCONFIDENCE name="CL" legendlabel=LABELCL; endif; endif; stepplot y=SURVIVAL x=TIME / name="Survival" rolename=(_tip1=ATRISK _tip2=EVENT) tiplabel=(_tip1="Number at Risk" _tip2="Observed Events") tip=(x y _tip1 _tip2) legendlabel="Survival"; if (PLOTCENSORED=1) scatterplot y=CENSORED x=TIME / markerattrs=(symbol=plus) tiplabel=(y="Survival Probability") name="Censored" legendlabel="Censored"; endif; if (PLOTCL=1 OR PLOTHW=1 OR PLOTEP=1) discretelegend "Censored" "CL" "HW" "EP" / location=outside halign=center; else if (PLOTCENSORED=1) discretelegend "Censored" / location=inside autoalign=(topright bottomleft); endif; endif; endlayout; endcell; cell; layout overlay / walldisplay=none xaxisopts=(display=none); axistable x=TATRISK value=ATRISK / display=(label) valueattrs=(size=7pt); drawtext 'Number at Risk' / x=3 y=20 width=9; endlayout; endcell; endlayout; else entrytitle &titletext1; if (EXISTS(SECONDTITLE)) entrytitle SECONDTITLE / textattrs=GRAPHVALUETEXT; endif; layout lattice / rows=2 rowweights=ROWWEIGHTS columndatarange=union rowgutter=10; cell; layout overlay / xaxisopts=(shortlabel=XNAME offsetmin=0 linearopts=(viewmax=MAXTIME tickvaluelist=XTICKVALS tickvaluefitpolicy=XTICKVALFITPOL)) yaxisopts=(label="&yaxis." shortlabel="Survival" offsetmin=0 linearopts=( viewmin=0 viewmax=1 tickvaluelist=(0 .2 .4 .6 .8 1.0))); if (PLOTHW=1) bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / displayTail=false group=STRATUM index=STRATUMNUM modelname="Survival" datatransparency=Transparency; endif; if (PLOTEP=1) bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / displayTail=false group=STRATUM index=STRATUMNUM modelname="Survival" datatransparency=Transparency; endif; if (PLOTCL=1) if (PLOTHW=1 OR PLOTEP=1) bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / displayTail=false group=STRATUM index=STRATUMNUM modelname="Survival" display=(outline) outlineattrs=(pattern=ShortDash); else bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / displayTail=false group=STRATUM index=STRATUMNUM modelname="Survival" datatransparency=Transparency; endif; endif; stepplot y=SURVIVAL x=TIME / group=STRATUM index=STRATUMNUM name="Survival" rolename=(_tip1=ATRISK _tip2=EVENT) tiplabel=(_tip1="Number at Risk" _tip2="Observed Events") tip=(x y _tip1 _tip2); if (PLOTCENSORED=1) scatterplot y=CENSORED x=TIME / group=STRATUM index=STRATUMNUM tiplabel=(y="Survival Probability") markerattrs =(symbol=plus); endif; DiscreteLegend "Survival" / title=GROUPNAME location=outside ; if (PLOTCENSORED=1) if (PLOTTEST=1) layout gridded / rows=2 autoalign=(TOPRIGHT ) border=false BackgroundColor= GraphWalls:Color Opaque=true; entry "+ Censored"; if (PVALUE < .0001) entry TESTNAME " p " eval (PUT(PVALUE, PVALUE6.4)); else entry TESTNAME " p=" eval (PUT(PVALUE, PVALUE6.4)); endif; endlayout; else layout gridded / rows=1 autoalign=(TOPRIGHT ) border=false BackgroundColor= GraphWalls:Color Opaque=true; entry "+ Censored"; endlayout; endif; else if (PLOTTEST=1) layout gridded / rows=1 autoalign=(TOPRIGHT BOTTOMLEFT TOP BOTTOM) border=true BackgroundColor= GraphWalls:Color Opaque=true; if (PVALUE < .0001) entry TESTNAME " p " eval (PUT(PVALUE, PVALUE6.4)); else entry TESTNAME " p=" eval (PUT(PVALUE, PVALUE6.4)); endif; endlayout; endif; endif; endlayout; endcell; cell; layout overlay / walldisplay=none xaxisopts=(display=none); axistable x=TATRISK value=ATRISK / display=(label) valueattrs=(size=7pt) class=CLASSATRISK colorgroup=CLASSATRISK ; drawtext 'Number at Risk' / x=3 y=55 width=9 rotate= 90; endlayout; endcell; endlayout; endif; if (_BYTITLE_) entrytitle _BYLINE_ / textattrs=GRAPHVALUETEXT; else if (_BYFOOTNOTE_) entryfootnote halign=left _BYLINE_; endif; endif; EndGraph; end; I have set offsetmin to 0 in xaxisopts , but it doesn't work. Any other way to do it/? ... View more

I have applied a cox proportional model on discontinuation data . There are two groups and OD and OW.  If the hazard ratio from cox proportional model is such-   Class HR OD 1.456 Then do can we interpret it as OW are 45.6%more likely to stay on treatment as compared to OD? Or should we say OD are 45.6% more likely to discontinue? Or do they mean the same?   ... View more

I have used the following code -   proc psmatch data = x; class a b c cate; psmodel cate(Treated = 'On') = a b c d ; match method = greedy(k=1) ; assess lps var=(a b c d )/plots = all weight = none; output out(obs=match)= y; run; however the treated group has 6000 obs and control has 8000 obs , then why is the number of obs used 1134 and 1411 ? ... View more