Hi @PGStats ,   I added the Units Statement into the code below:   proc logistic data=Sample1 descending; model Death = Culture_TIME; UNITS Culture_TIME='01:00:00't ; /* 1 hour time unit */ format Culture_TIME BEST8.; run;   Testing Global Null Hypothesis: BETA=0 Test Chi-Square DF Pr > ChiSq Likelihood Ratio 3.5949 1 0.0580 Score 3.7231 1 0.0537 Wald 3.6045 1 0.0576 Analysis of Maximum Likelihood Estimates Parameter DF Estimate Standard Error Wald Chi-Square Pr > ChiSq Intercept 1 -1.5334 0.3782 16.4420 <.0001 Culture_TIME 1 3.007E-6 1.584E-6 3.6045 0.0576 Odds Ratio Estimates Effect Point Estimate 95% Wald Confidence Limits Culture_TIME 1.000 1.000 1.000 Association of Predicted Probabilities and Observed Responses Percent Concordant 56.3 Somers' D 0.126 Percent Discordant 43.7 Gamma 0.126 Percent Tied 0.0 Tau-a 0.052 Pairs 5290 c 0.563 Odds Ratios Effect Unit Estimate Culture_TIME 3600.0 1.011 Then I got this OR = 1.011 for unit 3600 seconds. Is there a way to output the 95% C.I. as well? ... View more

Just uploaded the Excel file. I changed the format to Best8. but I still got 1.00 (1.00-1.00) results.   Analysis Variable : Culture_time N Mean Std Dev Minimum Maximum 161 198537.02 107437.48 45240.00 562200.00   The LOGISTIC Procedure Model Information Data Set WORK.Sample1   Response Variable Death Death Number of Response Levels 2   Model binary logit   Optimization Technique Fisher's scoring   Number of Observations Read 161 Number of Observations Used 161 Response Profile Ordered Value Death Total Frequency 1 1 46 2 0 115 Probability modeled is Death=1.   Model Convergence Status Convergence criterion (GCONV=1E-8) satisfied. Model Fit Statistics Criterion Intercept Only Intercept and Covariates AIC 194.643 193.048 SC 197.724 199.211 -2 Log L 192.643 189.048 Testing Global Null Hypothesis: BETA=0 Test Chi-Square DF Pr > ChiSq Likelihood Ratio 3.5949 1 0.0580 Score 3.7231 1 0.0537 Wald 3.6045 1 0.0576 Analysis of Maximum Likelihood Estimates Parameter DF Estimate Standard Error Wald Chi-Square Pr > ChiSq Intercept 1 -1.5334 0.3782 16.4420 <.0001 Culture_TIME 1 3.007E-6 1.584E-6 3.6045 0.0576 Odds Ratio Estimates Effect Point Estimate 95% Wald Confidence Limits Culture_TIME 1.000 1.000 1.000 Association of Predicted Probabilities and Observed Responses Percent Concordant 56.3 Somers' D 0.126 Percent Discordant 43.7 Gamma 0.126 Percent Tied 0.0 Tau-a 0.052 Pairs 5290 c 0.563 ... View more

@Rick_SAS  I just want to make sure that this code does generate 95% C.I., not 90% CI, right? Thank you so much! ... View more

How do I export the 95% C.I. and convert them into the hr:min:sec format? ... View more

@ed_sas_member Thank you so much for sharing the code! It worked, but I still need to generate 95% C.I. and export to the "want" dataset. Can we add another code or syntax somewhere?  Thanks again! ... View more

Thank you both for your detailed explainations. I ran proc contents and saw Total_time was in Time8. format. Then I ran the proc univariate with 95% C.I. below.   ods select Quantiles; proc univariate data=Sample1 ciquantdf(alpha=.05); var Total_time; run;   But I don't understand what the units are in the below SAS Output. Are these numbers in second?   The UNIVARIATE Procedure Variable: TOTAL_ Time   Quantiles (Definition 5)             Level Quantile     Order Statistics         95% Confidence Limits   LCL Rank UCL Rank Coverage     Distribution Free         100% Max 1296000           99% 844560 486360 1296000 157 161 77.82 95% 345780 328680 549420 148 159 95.64 90% 307380 277740 345780 138 153 95.26 75% Q3 231840 208140 262680 111 133 95.28 50% Median 150720 128940 176880 69 94 95.08 25% Q1 94380 62100 102420 29 51 95.28 10% 18780 0 33180 9 24 95.26 5% 0 0 6900 3 14 95.64 1% 0 0 0 1 5 77.82 0% Min 0             Thanks again! ... View more

Hi,    This actually worked!!! THANK YOU!!!   But there is 2 more things that I am trying to accomplish:    1. How to add "log rank p" under "+Censor"?   2. I also tried to add "AKI" as my group name for the legend (before No ----- and Yes ----)- I thougt the " label stratum = AKI;"  or "keylegend 'No' 'Yes' / title="AKI"; " would add "AKI"but it did not. /*--Get survival plot data from LIFETEST procedure--*/ ods select none; ods output Survivalplot=SurvivalPlotData; proc lifetest data=two plots=survival(atrisk=5 to 20 by 2.5); time Survive_year * death(0); strata AKI; run; ods select all; /*--Survival Plot with outer Risk Table using AxisTable--*/ ods graphics / reset width=5in height=3in ; title 'Product-Limit Survival Estimates'; title2 h=0.8 'With Number of Subjects at Risk'; ods listing gpath= "H:\Biostatistics\Ihsin\Subspecialty Medicine\Renal Service\Latcha,Sheron\SAS Program" image_dpi=300 ; ods graphics on/ imagefmt = tiff imagename = "Figure_AKI_Survivor_5yr" reset = index; proc sgplot data=SurvivalPlotData noborder; step x=time y=survival / group=stratum name='s'; scatter x=time y=censored / markerattrs=(symbol=plus) name='c'; scatter x=time y=censored / markerattrs=(symbol=plus) GROUP=stratum; xaxistable atrisk / x=tatrisk location=inside class=stratum colorgroup=stratum valueattrs=(weight=bold); keylegend 'c' / location=inside position=topright; keylegend 's' / linelength=20; keylegend 'No' 'Yes' / title="AKI"; *controls the graph axis; xaxis values=(5 to 20 by 5); yaxis values=(0 to 1 by .2); label time = Time (Year) stratum = AKI; run; ... View more

Hi,   I ran below code per your suggestion: ods output survivalPlotData = survGraphData; proc lifetest data=two plots=survival ; time Survive_year*death(0); strata AKI; run;   But it gave me the error message:   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: Graph's name, LIFETEST, changed to LIFETES1. LIFETEST is already used or not a valid SAS name. NOTE: Graph's name, LIFETEST, changed to LIFETES1. LIFETEST is already used or not a valid SAS name. NOTE: 29352 bytes written to C:\Users\LinI1\AppData\Local\Temp\SAS Temporary Files\_TD7508_WBSTE046_\lifetest1.png. WARNING: Output 'survivalPlotData' was not created. Make sure that the output object name, label, or path is spelled correctly. Also, verify that the appropriate procedure options are used to produce the requested output object. For example, verify that the NOPRINT option is not used. NOTE: PROCEDURE LIFETEST used (Total process time): real time 0.56 seconds cpu time 0.46 seconds   ... View more