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: 
Ameribrit
‎01-17-2020
  • 5 Posts
  • 0 Likes Given
  • 0 Solutions
  • 0 Likes Received

Re: Proc plm: effectplot slicefit

by Calcite | Level 5 in Graphics Programming
‎01-09-2020 11:20 AM
‎01-09-2020 11:20 AM
Hi Rick: this may be a stupid question but when I run the following code: ods graphics on; proc logistic data=mapmerge plots=effect; class Shock_Vasopressor__Use(ref='No') respmechvent(ref='No') Renal_Failure_HD_Use (ref='No') Grade_3_4_Hepenceph(ref='Yes')/param=ref; model ICU_DEATH(EVENT='Yes')=twamap55 meld_na_score respmechvent Renal_Failure_HD_Use Shock_Vasopressor__Use albumin wbc Grade_3_4_Hepenceph/clodds=wald; UNITS twamap55= -1 ; store adjMort; run;     SAS produces a default effects plot (see attachment)     I am trying to understand how SAS produces this plot and if I can reproduce it using Proc plm so I can create my own combinations if needed.   So far I am trying the following code but I cant seem to slice by multiple categorical predictors ( if I assume thats what SAS is doing to produce the plot) the at values work but I dont know how to mimic whats labelled on the xaxis    /*test plm*/ proc plm restore=adjMort; effectplot slicefit (x=twamap55 sliceby=Shock_Vasopressor__Use*respmechvent*Renal_Failure_HD_Use) / at( Grade_3_4_Hepenceph='Yes' meld_na_score=mean albumin=mean wbc=mean); run;   If I just use one var in sliceby I get a plot but not the default one SAS produces. Thanks Rick and sorry for the long winded question Francis   ... View more

Re: How to Modify Effect Plots In Proc Logistic

by Calcite | Level 5 in New SAS User
‎10-09-2019 09:22 AM
‎10-09-2019 09:22 AM
thanks for the advice but I'm not sure how to do as you suggest. I ran the proc with ods trace on and I get a bunch of stuff in the log one of which is  something resembling effects plot but not sure what to do with it or how to grab it. thanks Francis   Output Added: ------------- Name: FitStatistics Label: Fit Statistics Template: Stat.Logistic.FitStatistics Path: Logistic.FitStatistics ------------- Output Added: ------------- Name: GlobalTests Label: Global Tests Template: Stat.Logistic.GlobalTests Path: Logistic.GlobalTests ------------- Output Added: ------------- Name: ParameterEstimates Label: Parameter Estimates Template: Stat.Logistic.ParameterEstimates Path: Logistic.ParameterEstimates ------------- Output Added: ------------- Name: Association Label: Association Statistics Template: Stat.Logistic.Association Path: Logistic.Association ------------- Output Added: ------------- Name: CLoddsWald Label: 95% Clodds=Wald Template: Stat.Logistic.CloddsWald Path: Logistic.CLoddsWald ------------- Output Added: ------------- Name: ORPlot Label: 95% Clodds=Wald Plots Template: Stat.Logistic.Graphics.ORPlot Path: Logistic.ORPlot ------------- Output Added: ------------- Name: FitPlot Label: Fit Plot 1 Template: Stat.Logistic.Graphics.FitPlot Path: Logistic.EffectPlots.EffectPlots.FitPlot ------------- NOTE: The LOGISTIC procedure generated the model item store WORK.LOGIMODEL. NOTE: There were 277 observations read from the data set WORK.MAPMERGE. NOTE: PROCEDURE LOGISTIC used (Total process time): real time 0.54 seconds cpu time 0.28 seconds ... View more

How to Modify Effect Plots In Proc Logistic

by Calcite | Level 5 in New SAS User
‎10-09-2019 08:55 AM
‎10-09-2019 08:55 AM
Hi Guys: can anyone give me advice on how to customize the effects plot produced by SAS in Proc Logistic. What I am trying to do is to change the  X axis range so my predicted probabilities are across a wider range than default. I tried the following but its ignored by SAS with no error in the log, Forgive me if this is a stupid question. twamap55 is a continuous variable of mean arterial pressure scores.   ods graphics on; proc logistic data=mapmerge ; model ICU_DEATH(EVENT='Yes')=twap55 ; UNITS twap55= -1 ; effectplot / at(twap55=0 to 75 by 5) ; run; ods graphics off;   The output stays the same no matter what I try.     ... View more

Re: TIME SPENT BELOW A THRESHOLD

by Calcite | Level 5 in SAS Programming
‎09-20-2019 01:30 PM
‎09-20-2019 01:30 PM
Thanks for the response. The issue is the TWAP was not defined in the paper. Only roughly as Area under a threshold curve divided by the overall time on study. I am trying to make sense of what that actually means. In this case the threshold is a fixed constant and Im not sure what Area under this curve means. Sorry if I'm not conveying the issues appropriately.   The data I have is irregular spaced mean arterial pressures for every patient over time. I am trying to figure out what this metric is and what it means as I have been asked to create this for a project I am working on.   ... View more

TIME SPENT BELOW A THRESHOLD

by Calcite | Level 5 in SAS Programming
‎09-20-2019 11:48 AM
‎09-20-2019 11:48 AM
Hi Guys: I am banging my head trying to figure out how to look at a longitudinal profile and to determine two Metrics. One: Time a patient spent below a given threshold Two Time weighted Average: defined as the area below the threshold divided by patients overall exposure.   The data is in long form:    ID RESULT  DATETIME 1       x1        D1 1       x2        D2 1        x3       D3 etc for patients 1-n.   The first datestamp is baseline and the last is their last measurement. I know cumulative exposure time is lastdate-Firstdate. What I am trying to figure out is (cumulative time < threshold) for every patient and (Area under the threshold curve/cumulative exposure) for every patient.   I naively assumed I could set my data keeping only values below threshold and then computing the metrics as follows. data map&threshold._actual_final; set MAP_Actual1; where result LE &threshold; by MRN; area + dif(time_from_admission)*mean(result,lag(result)); timeLT&threshold+dif(time_from_admission); TWAP=area/cumulative_exposure; if first.MRN then do; area=0; timeLT&threshold=0; end; if last.MRN  then output; run;   I think this is over inflated because the data is not evenly spaced and patients can have values above and below any particular threshold over the time span. Also please excuse if I missed the boat on what I supplied but I never want to ask questions until I at least have tried and know the issues involved. Thanks guys. Totally struggling on this one.   I attached some simulated data with a threshold of 0.5 to give you an idea of what I am trying to do. For example patient 1 in red, over time can oscillate over and below the threshold over time so I dont think I can just restrict my data and compute area etc as I outlined above. yours in Frustration ... View more
Contact Me
Online Status
Offline
Date Last Visited
‎01-17-2020 02:59 PM
Quick links