# ehh2.cf  Graphical test of distributional assumptions: Weibull         
 
nvar(
    dfile = ehc5.ple,           # data file 
 
    # select suitable range of time and survivor function 
 
    isel = gt(c3,0) & gt(c7,0) & lt(c7,1), 
 
    Time = c3,              
    Surv = c7,              
    LogTime = log(Time),    
    LogLogSurv = log(-log(Surv)), 
); 
lsreg = LogLogSurv,LogTime;       request OLS regression

psfile = ehh2.ps;  output file
psetup(
    pxlen =  80,    # length of x axis in mm
    pylen =  40,    # length of y axis in mm

    pxa   = 0,6,    # user coordinates on x axis
    pya   = -6,2,   # user coordinates on y axis
);
plxa (sc=1,ic=5);     plot x axis
plya (sc=1,ic=2);     plot y axis
plframe;              plot frame

plot = LogTime,LogLogSurv;
plotp = 0,-4.4238, 6, 1.6242;     plot regression line