1This is file bootllg.men--Bootstrap Confidence Intervals for Log Logistic 2 3-------------------- 4 1. Log Logistic Bootstrap Confidence Intervals 5 2. BOOTSTRAP LOG LOGISTIC PLOT Y 6 3. 7 4. The log logistic PPCC/probability plot generates 8 5. point estimates for the location, scale, and shape 9 6. parameters. This menu uses the bootstrap to generate 10 7. confidence intervals for these parameters. Both 11 8. uncensored and censored data are supported. Grouped 12 9. data is not. Censoring is only supported for the PPCC 1310. plot. 1411. 1512. If a censoring variable is specified, it should contain 1613. a value of 1 to signify a failure time and a value of 0 1714. to signify a censoring time. 1815. 1916. Response Variable: 2016. @CE 1 0 42 50 ? gui-variable 2117. Optional Censoring Variable: 2217. @CE 2 0 42 50 ? gui-variable 2318. 2419. Number of bootstrap samples: 2519. @CE 3 0 42 50 * 2620. Select goodness of fit critierion: 2721. @CE 4 1 1 30 PPCC 2822. @CE 4 2 1 30 Kolmogorov-Smirnov 2923. Minimum Value for delta (omit for default): 3023. @CE 5 0 42 50 * 3124. Maximum Value for delta (omit for default): 3224. @CE 6 0 42 50 * 3325. 3426. Significance Level (default=0.05): 3527. @CE 7 0 42 50 ? * 3628. 3729. Subset/Except/For Specification: 3830. @CE 8 0 1 50 * 3931. 4032. 4133. Compute Bootstrap Confidence Intervals? 42-------------------- 43GR FE 44LET NUMBOOT = 100 45@IF 3 0 LET NUMBOOT = _3 46let numboot = min(1000,numboot) 47let numboot = max(100,numboot) 48let icensor = 1 49@IF 2 0 let icensor = 3 50@IF 4 2 let icensor = icensor + 1 51let alpha = 0.05 52let alpha = _7 53if alpha <= 0 54 let alpha = 0.05 55end of if 56if alpha >= 1 57 let alpha = 0.05 58end of if 59if alpha > 0.5 60 let alpha = 1 - alpha 61end of if 62let xqlow = alpha/2 63let xqupp = 1 - alpha/2 64delete delta1 delta2 65let delta1 = _5 66let delta2 = _6 67. 68bootstrap samples numboot 69set maximum likelihood percentiles default 70device 1 font simplex 71multiplot 2 3 72multiplot corner coordinates 0 0 100 100 73multiplot scale factor 2 74label 75title case asis 76label case asis 77case asis 78title displacement 2 79x2label displacement 13 80x3label displacement 18 81y1label Parameter Estimate 82x1label 83x2label Bootstrap Sample 84title Bootstrap Plot 85line color blue red green black 86limits 87if icensor = 1 88 bootstrap log logistic plot _1 _8 89end of if 90if icensor = 2 91 bootstrap log logistic ks plot _1 _8 92end of if 93if icensor = 3 94 bootstrap log logistic censored plot _1 _2 _8 95end of if 96if icensor = 4 97 bootstrap log logistic censored plot _1 _2 _8 98 . bootstrap log logistic ks censored plot _1 _2 _8 99end of if 100line color black all 101. 102delete aloc ascale adelta 103skip 0 104read dpst1f.dat aloc ascale adelta 105y1label 106x2label 107title Location Parameter 108let amed = median aloc 109let amean = mean aloc 110let asd = sd aloc 111x2label Median = ^amed, Mean = ^amean 112x3label SD = ^asd 113histogram aloc 114title Scale Parameter 115let amed = median ascale 116let amean = mean ascale 117let asd = sd ascale 118x2label Median = ^amed, Mean = ^amean 119x3label SD = ^asd 120histogram ascale 121title Delta Parameter 122let amed = median adelta 123let amean = mean adelta 124let asd = sd adelta 125x2label Median = ^amed, Mean = ^amean 126x3label SD = ^asd 127histogram adelta 128x3label displacement 129end of multiplot 130. 131let xq = xqlow 132let loc95low = xq quantile aloc 133let xq = xqupp 134let loc95upp = xq quantile aloc 135let xq = xqlow 136let sca95low = xq quantile ascale 137let xq = xqupp 138let sca95upp = xq quantile ascale 139let xq = xqlow 140let del95low = xq quantile adelta 141let xq = xqupp 142let del95upp = xq quantile adelta 143. 144device 1 font tektronix 145window corner coordinates 50 0 100 50 146justification left 147height 4 148move 5 70 149text Confidence Interval for Location: 150move 5 63 151text (^loc95low,^loc95upp) 152move 5 50 153text Confidence Interval for Scale: 154move 5 43 155text (^sca95low,^sca95upp) 156move 5 30 157text Confidence Interval for Delta: 158move 5 23 159text (^del95low,^del95upp) 160window corner coordinates 0 0 100 100 161height 162title 163label 164-------------------- 165WEB HELP DISTRIBUTIONAL BOOTSTRAP 166-------------------- 167