# Rat data from survival5 # Description: # 48 rats were injected with a carcinogen, and then randomized to # either drug or placebo. The number of tumors ranges from 0 to 13; # all rats were censored at 6 months after randomization. # Usage: # data(rats) # # Format: # rat: id # rx: treatment,(1=drug, 0=control) # observation: within rat # start: entry time # stop: exit time # status: event status # Source: # Gail, Sautner and Brown, Biometrics 36, 255-66, 1980 # # survreg(formula = Surv(time, status) ~ rx, data = rats) # Value Std. Error z p # (Intercept) 4.983 0.0833 59.81 0.00e+00 # rx -0.239 0.0891 -2.68 7.42e-03 # Log(scale) -1.333 0.1439 -9.26 2.01e-20 # # survreg(formula = Surv(time, status) ~ rx + frailty.gaussian(litter, # df = 13, sparse = FALSE), data = rats) # # coef se(coef) se2 Chisq DF p # (Intercept) 4.881 0.0663 0.0607 5413.38 1 0.0000 # rx -0.187 0.0690 0.0677 7.38 1 0.0066 # frailty.gaussian(litter, 16.26 13 0.2400 # # Scale= 0.2 # # Iterations: 8 outer, 34 Newton-Raphson # Variance of random effect= 0.0193 # Degrees of freedom for terms= 0.8 1.0 13.0 0.9 # Likelihood ratio test=35.9 on 13.7 df, p=0.000902 n= 150 # # litter rx time status set ple -1 set loc rx aff set loc time qua set loc status aff read ped inline rats 1 f1 m1 x 2 101 1 rats 2 f1 m1 x 1 49 2 rats 3 f1 m1 x 1 104 1 rats 4 f2 m2 x 2 104 1 rats 5 f2 m2 x 1 102 1 rats 6 f2 m2 x 1 104 1 rats 7 f3 m3 x 2 104 1 rats 8 f3 m3 x 1 104 1 rats 9 f3 m3 x 1 104 1 rats 10 f4 m4 x 2 77 1 rats 11 f4 m4 x 1 97 1 rats 12 f4 m4 x 1 79 1 rats 13 f5 m5 x 2 89 1 rats 14 f5 m5 x 1 104 1 rats 15 f5 m5 x 1 104 1 rats 16 f6 m6 x 2 88 2 rats 17 f6 m6 x 1 96 2 rats 18 f6 m6 x 1 104 1 rats 19 f7 m7 x 2 104 2 rats 20 f7 m7 x 1 94 1 rats 21 f7 m7 x 1 77 2 rats 22 f8 m8 x 2 96 2 rats 23 f8 m8 x 1 104 1 rats 24 f8 m8 x 1 104 1 rats 25 f9 m9 x 2 82 1 rats 26 f9 m9 x 1 77 1 rats 27 f9 m9 x 1 104 1 rats 28 f10 m10 x 2 70 2 rats 29 f10 m10 x 1 104 1 rats 30 f10 m10 x 1 77 1 rats 31 f11 m11 x 2 89 2 rats 32 f11 m11 x 1 91 1 rats 33 f11 m11 x 1 90 1 rats 34 f12 m12 x 2 91 1 rats 35 f12 m12 x 1 70 1 rats 36 f12 m12 x 1 92 1 rats 37 f13 m13 x 2 39 2 rats 38 f13 m13 x 1 45 1 rats 39 f13 m13 x 1 50 2 rats 40 f14 m14 x 2 103 2 rats 41 f14 m14 x 1 69 1 rats 42 f14 m14 x 1 91 1 rats 43 f15 m15 x 2 93 1 rats 44 f15 m15 x 1 104 1 rats 45 f15 m15 x 1 103 1 rats 46 f16 m16 x 2 85 1 rats 47 f16 m16 x 1 72 1 rats 48 f16 m16 x 1 104 1 rats 49 f17 m17 x 2 104 1 rats 50 f17 m17 x 1 63 1 rats 51 f17 m17 x 1 104 1 rats 52 f18 m18 x 2 104 1 rats 53 f18 m18 x 1 104 1 rats 54 f18 m18 x 1 74 1 rats 55 f19 m19 x 2 81 1 rats 56 f19 m19 x 1 104 1 rats 57 f19 m19 x 1 69 1 rats 58 f20 m20 x 2 67 2 rats 59 f20 m20 x 1 104 1 rats 60 f20 m20 x 1 68 2 rats 61 f21 m21 x 2 104 1 rats 62 f21 m21 x 1 104 1 rats 63 f21 m21 x 1 104 1 rats 64 f22 m22 x 2 104 1 rats 65 f22 m22 x 1 104 1 rats 66 f22 m22 x 1 104 1 rats 67 f23 m23 x 2 104 1 rats 68 f23 m23 x 1 83 1 rats 69 f23 m23 x 1 40 2 rats 70 f24 m24 x 2 87 1 rats 71 f24 m24 x 1 104 1 rats 72 f24 m24 x 1 104 1 rats 73 f25 m25 x 2 104 1 rats 74 f25 m25 x 1 104 1 rats 75 f25 m25 x 1 104 1 rats 76 f26 m26 x 2 89 1 rats 77 f26 m26 x 1 104 1 rats 78 f26 m26 x 1 104 1 rats 79 f27 m27 x 2 78 1 rats 80 f27 m27 x 1 104 1 rats 81 f27 m27 x 1 104 1 rats 82 f28 m28 x 2 104 1 rats 83 f28 m28 x 1 81 2 rats 84 f28 m28 x 1 64 2 rats 85 f29 m29 x 2 86 2 rats 86 f29 m29 x 1 55 2 rats 87 f29 m29 x 1 94 1 rats 88 f30 m30 x 2 34 2 rats 89 f30 m30 x 1 104 1 rats 90 f30 m30 x 1 54 2 rats 91 f31 m31 x 2 76 1 rats 92 f31 m31 x 1 87 1 rats 93 f31 m31 x 1 74 1 rats 94 f32 m32 x 2 103 2 rats 95 f32 m32 x 1 73 2 rats 96 f32 m32 x 1 84 2 rats 97 f33 m33 x 2 102 2 rats 98 f33 m33 x 1 104 1 rats 99 f33 m33 x 1 80 1 rats 100 f34 m34 x 2 80 2 rats 101 f34 m34 x 1 104 1 rats 102 f34 m34 x 1 73 1 rats 103 f35 m35 x 2 45 2 rats 104 f35 m35 x 1 79 1 rats 105 f35 m35 x 1 104 1 rats 106 f36 m36 x 2 94 2 rats 107 f36 m36 x 1 104 1 rats 108 f36 m36 x 1 104 1 rats 109 f37 m37 x 2 104 1 rats 110 f37 m37 x 1 104 1 rats 111 f37 m37 x 1 104 1 rats 112 f38 m38 x 2 104 1 rats 113 f38 m38 x 1 101 2 rats 114 f38 m38 x 1 94 1 rats 115 f39 m39 x 2 76 1 rats 116 f39 m39 x 1 84 2 rats 117 f39 m39 x 1 78 2 rats 118 f40 m40 x 2 80 2 rats 119 f40 m40 x 1 81 2 rats 120 f40 m40 x 1 76 1 rats 121 f41 m41 x 2 72 2 rats 122 f41 m41 x 1 95 1 rats 123 f41 m41 x 1 104 1 rats 124 f42 m42 x 2 73 2 rats 125 f42 m42 x 1 104 1 rats 126 f42 m42 x 1 66 2 rats 127 f43 m43 x 2 92 2 rats 128 f43 m43 x 1 104 1 rats 129 f43 m43 x 1 102 2 rats 130 f44 m44 x 2 104 1 rats 131 f44 m44 x 1 98 1 rats 132 f44 m44 x 1 73 1 rats 133 f45 m45 x 2 55 1 rats 134 f45 m45 x 1 104 1 rats 135 f45 m45 x 1 104 1 rats 136 f46 m46 x 2 49 1 rats 137 f46 m46 x 1 83 1 rats 138 f46 m46 x 1 77 1 rats 139 f47 m47 x 2 89 2 rats 140 f47 m47 x 1 104 1 rats 141 f47 m47 x 1 104 1 rats 142 f48 m48 x 2 88 1 rats 143 f48 m48 x 1 79 1 rats 144 f48 m48 x 1 99 1 rats 145 f49 m49 x 2 103 2 rats 146 f49 m49 x 1 91 1 rats 147 f49 m49 x 1 104 1 rats 148 f50 m50 x 2 104 1 rats 149 f50 m50 x 1 104 1 rats 150 f50 m50 x 1 79 2 ;;;; run set loc start qua; (start)=0 pyar start time status 1 10 time kap time status bye set ple 0 set iter 1000 reg time = rx exponential status reg time = rx weibull status fpm time lik weibull cen status cov rx nqtl 0 s print fpm time lik weibull cen status cov rx nqtl 0 s fix s sval 0.0