# Thall and Vail (1990) give a data set on two-week seizure counts # for 59 epileptics. The number of seizures was recorded for a # baseline period of 8 weeks, and then patients were randomly # assigned to a treatment group or a control group. Counts were # then recorded for four successive two-week periods. The subject's # age is the only covariate. # # y The count for the 2-week period. # trt The treatment (progabide), "n" or "y". # base The counts in the baseline 8-week period. # age The subjects age, in years. # V4 0/1 indicator variable of period 4. # subject The subject number, 1 to 59. # period The period, 1 to 4. # lbase The log-counts for the baseline period, centred to have zero mean. # lage The log-ages, centred to have zero mean. # Visit (rep(1:4,59)-2.5)/5 # Base log(base/4) # Age log(age) # # Linear mixed-effects model fit by maximum likelihood # Data: epil # AIC BIC logLik # 521.7803 549.4909 -252.8901 # # Random effects: # Formula: ~1 | subject # (Intercept) Residual # StdDev: 0.4442704 1.400807 # # Variance function: # Structure: fixed weights # Formula: ~invwt # Fixed effects: y ~ lbase * trt + lage + V4 # Value Std.Error DF t-value p-value # (Intercept) 1.8696677 0.1055620 176 17.711554 0.0000 # lbase 0.8818228 0.1292834 54 6.820849 0.0000 # trtprogabide -0.3095253 0.1490438 54 -2.076740 0.0426 # lage 0.5335460 0.3463119 54 1.540652 0.1292 # V4 -0.1597696 0.0774521 176 -2.062819 0.0406 # lbase:trtprogabide 0.3415425 0.2033325 54 1.679725 0.0988 # Correlation: # (Intr) lbase trtprg lage V4 # lbase -0.126 # trtprogabide -0.691 0.089 # lage -0.103 -0.038 0.088 # V4 -0.162 0.000 0.000 0.000 # lbase:trtprogabide 0.055 -0.645 -0.184 0.267 0.000 # # Standardized Within-Group Residuals: # Min Q1 Med Q3 Max # -2.1324053 -0.6387114 -0.0848634 0.4196020 4.9787214 # # Number of Observations: 236 # Number of Groups: 59 # > lmer(y ~ lbase*trt + lage + V4 + (1|subject), data=epil, family=poisson) # Generalized linear mixed model fit using PQL # Formula: y ~ lbase * trt + lage + V4 + (1 | subject) # Data: epil # Family: poisson(log link) # AIC BIC logLik deviance # 598.0873 625.798 -291.0437 582.0873 # Random effects: # Groups Name Variance Std.Dev. # subject (Intercept) 0.10058 0.31715 # # of obs: 236, groups: subject, 59 # # Estimated scale (compare to 1) 1.400815 # # Fixed effects: # Estimate Std. Error z value Pr(>|z|) # (Intercept) 1.869668 0.074393 25.1323 < 2.2e-16 *** # lbase 0.881823 0.091111 9.6786 < 2.2e-16 *** # trtprogabide -0.309525 0.105036 -2.9468 0.003210 ** # lage 0.533548 0.244056 2.1862 0.028803 * # V4 -0.159770 0.054584 -2.9271 0.003422 ** # lbase:trtprogabide 0.341543 0.143296 2.3835 0.017149 * # --- # Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 # # > lmer(y ~ lbase*trt + lage + V4 + (1|subject), data=epil, family=poisson, method="AGQ") # Generalized linear mixed model fit using AGQ # Formula: y ~ lbase * trt + lage + V4 + (1 | subject) # Data: epil # Family: poisson(log link) # AIC BIC logLik deviance # 580.9085 608.6191 -282.4542 564.9085 # Random effects: # Groups Name Variance Std.Dev. # subject (Intercept) 0.25239 0.50239 # # of obs: 236, groups: subject, 59 # # Estimated scale (compare to 1) 2.165665 # # Fixed effects: # Estimate Std. Error z value Pr(>|z|) # (Intercept) 1.832765 0.074393 24.6362 < 2.2e-16 *** # lbase 0.883408 0.091111 9.6960 < 2.2e-16 *** # trtprogabide -0.334259 0.105036 -3.1823 0.001461 ** # lage 0.480566 0.244056 1.9691 0.048944 * # V4 -0.159769 0.054584 -2.9270 0.003422 ** # lbase:trtprogabide 0.338785 0.143296 2.3642 0.018067 * # --- # # set ple -1 set loc period qua . The period, 1 to 4. set loc seizures qua . The seizure count for the 2-week period. set loc trt aff . The treatment (progabide). set loc base qua . The seizure count in the baseline 8-week period. set loc age qua . The subjects age, in years. set loc V4 qua . 0/1 indicator variable of period 4. set loc lbase qua . The log-counts for the baseline period, centred to have zero mean. set loc lage qua . The log-ages, centred to have zero mean. set loc Visit qua . (rep(1:4,59)-2.5)/5 set loc Base qua . log(base/4) set loc Age qua . log(age) read pedigree inline #subject period yvar trt base age V4 lbase lage Visit Base Age 1 1 f1 m1 x 1 5 n 11 31 0 -0.756353788717556 0.114203695299265 -0.3 1.01160091167848 3.43398720448515 1 2 f1 m1 x 2 3 n 11 31 0 -0.756353788717556 0.114203695299265 -0.1 1.01160091167848 3.43398720448515 1 3 f1 m1 x 3 3 n 11 31 0 -0.756353788717556 0.114203695299265 0.1 1.01160091167848 3.43398720448515 1 4 f1 m1 x 4 3 n 11 31 1 -0.756353788717556 0.114203695299265 0.3 1.01160091167848 3.43398720448515 1 5 f2 m2 x 1 3 n 11 30 0 -0.756353788717556 0.0814138724762743 -0.3 1.01160091167848 3.40119738166216 1 6 f2 m2 x 2 5 n 11 30 0 -0.756353788717556 0.0814138724762743 -0.1 1.01160091167848 3.40119738166216 1 7 f2 m2 x 3 3 n 11 30 0 -0.756353788717556 0.0814138724762743 0.1 1.01160091167848 3.40119738166216 1 8 f2 m2 x 4 3 n 11 30 1 -0.756353788717556 0.0814138724762743 0.3 1.01160091167848 3.40119738166216 1 9 f3 m3 x 1 2 n 6 25 0 -1.36248959228787 -0.10090768431768 -0.3 0.405465108108164 3.2188758248682 1 10 f3 m3 x 2 4 n 6 25 0 -1.36248959228787 -0.10090768431768 -0.1 0.405465108108164 3.2188758248682 1 11 f3 m3 x 3 0 n 6 25 0 -1.36248959228787 -0.10090768431768 0.1 0.405465108108164 3.2188758248682 1 12 f3 m3 x 4 5 n 6 25 1 -1.36248959228787 -0.10090768431768 0.3 0.405465108108164 3.2188758248682 1 13 f4 m4 x 1 4 n 8 36 0 -1.07480751983609 0.263735429270229 -0.3 0.693147180559945 3.58351893845611 1 14 f4 m4 x 2 4 n 8 36 0 -1.07480751983609 0.263735429270229 -0.1 0.693147180559945 3.58351893845611 1 15 f4 m4 x 3 1 n 8 36 0 -1.07480751983609 0.263735429270229 0.1 0.693147180559945 3.58351893845611 1 16 f4 m4 x 4 4 n 8 36 1 -1.07480751983609 0.263735429270229 0.3 0.693147180559945 3.58351893845611 1 17 f5 m5 x 1 7 n 66 22 0 1.0354056805105 -0.228741055827565 -0.3 2.80336038090653 3.09104245335832 1 18 f5 m5 x 2 18 n 66 22 0 1.0354056805105 -0.228741055827565 -0.1 2.80336038090653 3.09104245335832 1 19 f5 m5 x 3 9 n 66 22 0 1.0354056805105 -0.228741055827565 0.1 2.80336038090653 3.09104245335832 1 20 f5 m5 x 4 21 n 66 22 1 1.0354056805105 -0.228741055827565 0.3 2.80336038090653 3.09104245335832 1 21 f6 m6 x 1 5 n 27 29 0 0.141587804488402 0.047512320800593 -0.3 1.90954250488444 3.36729582998647 1 22 f6 m6 x 2 2 n 27 29 0 0.141587804488402 0.047512320800593 -0.1 1.90954250488444 3.36729582998647 1 23 f6 m6 x 3 8 n 27 29 0 0.141587804488402 0.047512320800593 0.1 1.90954250488444 3.36729582998647 1 24 f6 m6 x 4 7 n 27 29 1 0.141587804488402 0.047512320800593 0.3 1.90954250488444 3.36729582998647 1 25 f7 m7 x 1 6 n 12 31 0 -0.669342411727926 0.114203695299265 -0.3 1.09861228866811 3.43398720448515 1 26 f7 m7 x 2 4 n 12 31 0 -0.669342411727926 0.114203695299265 -0.1 1.09861228866811 3.43398720448515 1 27 f7 m7 x 3 0 n 12 31 0 -0.669342411727926 0.114203695299265 0.1 1.09861228866811 3.43398720448515 1 28 f7 m7 x 4 2 n 12 31 1 -0.669342411727926 0.114203695299265 0.3 1.09861228866811 3.43398720448515 1 29 f8 m8 x 1 40 n 52 42 0 0.7969946570655 0.417886109097487 -0.3 2.56494935746154 3.73766961828337 1 30 f8 m8 x 2 20 n 52 42 0 0.7969946570655 0.417886109097487 -0.1 2.56494935746154 3.73766961828337 1 31 f8 m8 x 3 21 n 52 42 0 0.7969946570655 0.417886109097487 0.1 2.56494935746154 3.73766961828337 1 32 f8 m8 x 4 12 n 52 42 1 0.7969946570655 0.417886109097487 0.3 2.56494935746154 3.73766961828337 1 33 f9 m9 x 1 5 n 23 37 0 -0.018754845586777 0.291134403458343 -0.3 1.74919985480926 3.61091791264422 1 34 f9 m9 x 2 6 n 23 37 0 -0.018754845586777 0.291134403458343 -0.1 1.74919985480926 3.61091791264422 1 35 f9 m9 x 3 6 n 23 37 0 -0.018754845586777 0.291134403458343 0.1 1.74919985480926 3.61091791264422 1 36 f9 m9 x 4 5 n 23 37 1 -0.018754845586777 0.291134403458343 0.3 1.74919985480926 3.61091791264422 1 37 f10 m10 x 1 14 n 10 28 0 -0.851663968521881 0.0124210009893226 -0.3 0.916290731874155 3.3322045101752 1 38 f10 m10 x 2 13 n 10 28 0 -0.851663968521881 0.0124210009893226 -0.1 0.916290731874155 3.3322045101752 1 39 f10 m10 x 3 6 n 10 28 0 -0.851663968521881 0.0124210009893226 0.1 0.916290731874155 3.3322045101752 1 40 f10 m10 x 4 0 n 10 28 1 -0.851663968521881 0.0124210009893226 0.3 0.916290731874155 3.3322045101752 1 41 f11 m11 x 1 26 n 52 36 0 0.7969946570655 0.263735429270229 -0.3 2.56494935746154 3.58351893845611 1 42 f11 m11 x 2 12 n 52 36 0 0.7969946570655 0.263735429270229 -0.1 2.56494935746154 3.58351893845611 1 43 f11 m11 x 3 6 n 52 36 0 0.7969946570655 0.263735429270229 0.1 2.56494935746154 3.58351893845611 1 44 f11 m11 x 4 22 n 52 36 1 0.7969946570655 0.263735429270229 0.3 2.56494935746154 3.58351893845611 1 45 f12 m12 x 1 12 n 33 24 0 0.342258499950553 -0.141729678837935 -0.3 2.11021320034659 3.17805383034795 1 46 f12 m12 x 2 6 n 33 24 0 0.342258499950553 -0.141729678837935 -0.1 2.11021320034659 3.17805383034795 1 47 f12 m12 x 3 8 n 33 24 0 0.342258499950553 -0.141729678837935 0.1 2.11021320034659 3.17805383034795 1 48 f12 m12 x 4 4 n 33 24 1 0.342258499950553 -0.141729678837935 0.3 2.11021320034659 3.17805383034795 1 49 f13 m13 x 1 4 n 18 23 0 -0.263877303619762 -0.184289293256731 -0.3 1.50407739677627 3.13549421592915 1 50 f13 m13 x 2 4 n 18 23 0 -0.263877303619762 -0.184289293256731 -0.1 1.50407739677627 3.13549421592915 1 51 f13 m13 x 3 6 n 18 23 0 -0.263877303619762 -0.184289293256731 0.1 1.50407739677627 3.13549421592915 1 52 f13 m13 x 4 2 n 18 23 1 -0.263877303619762 -0.184289293256731 0.3 1.50407739677627 3.13549421592915 1 53 f14 m14 x 1 7 n 42 36 0 0.583420556767442 0.263735429270229 -0.3 2.35137525716348 3.58351893845611 1 54 f14 m14 x 2 9 n 42 36 0 0.583420556767442 0.263735429270229 -0.1 2.35137525716348 3.58351893845611 1 55 f14 m14 x 3 12 n 42 36 0 0.583420556767442 0.263735429270229 0.1 2.35137525716348 3.58351893845611 1 56 f14 m14 x 4 14 n 42 36 1 0.583420556767442 0.263735429270229 0.3 2.35137525716348 3.58351893845611 1 57 f15 m15 x 1 16 n 87 26 0 1.31165905713866 -0.061686971164399 -0.3 3.07961375753469 3.25809653802148 1 58 f15 m15 x 2 24 n 87 26 0 1.31165905713866 -0.061686971164399 -0.1 3.07961375753469 3.25809653802148 1 59 f15 m15 x 3 10 n 87 26 0 1.31165905713866 -0.061686971164399 0.1 3.07961375753469 3.25809653802148 1 60 f15 m15 x 4 9 n 87 26 1 1.31165905713866 -0.061686971164399 0.3 3.07961375753469 3.25809653802148 1 61 f16 m16 x 1 11 n 50 26 0 0.75777394391222 -0.061686971164399 -0.3 2.52572864430826 3.25809653802148 1 62 f16 m16 x 2 0 n 50 26 0 0.75777394391222 -0.061686971164399 -0.1 2.52572864430826 3.25809653802148 1 63 f16 m16 x 3 0 n 50 26 0 0.75777394391222 -0.061686971164399 0.1 2.52572864430826 3.25809653802148 1 64 f16 m16 x 4 5 n 50 26 1 0.75777394391222 -0.061686971164399 0.3 2.52572864430826 3.25809653802148 1 65 f17 m17 x 1 0 n 18 28 0 -0.263877303619762 0.0124210009893226 -0.3 1.50407739677627 3.3322045101752 1 66 f17 m17 x 2 0 n 18 28 0 -0.263877303619762 0.0124210009893226 -0.1 1.50407739677627 3.3322045101752 1 67 f17 m17 x 3 3 n 18 28 0 -0.263877303619762 0.0124210009893226 0.1 1.50407739677627 3.3322045101752 1 68 f17 m17 x 4 3 n 18 28 1 -0.263877303619762 0.0124210009893226 0.3 1.50407739677627 3.3322045101752 1 69 f18 m18 x 1 37 n 111 31 0 1.55528113979641 0.114203695299265 -0.3 3.32323584019244 3.43398720448515 1 70 f18 m18 x 2 29 n 111 31 0 1.55528113979641 0.114203695299265 -0.1 3.32323584019244 3.43398720448515 1 71 f18 m18 x 3 28 n 111 31 0 1.55528113979641 0.114203695299265 0.1 3.32323584019244 3.43398720448515 1 72 f18 m18 x 4 29 n 111 31 1 1.55528113979641 0.114203695299265 0.3 3.32323584019244 3.43398720448515 1 73 f19 m19 x 1 3 n 18 32 0 -0.263877303619762 0.145952393613845 -0.3 1.50407739677627 3.46573590279973 1 74 f19 m19 x 2 5 n 18 32 0 -0.263877303619762 0.145952393613845 -0.1 1.50407739677627 3.46573590279973 1 75 f19 m19 x 3 2 n 18 32 0 -0.263877303619762 0.145952393613845 0.1 1.50407739677627 3.46573590279973 1 76 f19 m19 x 4 5 n 18 32 1 -0.263877303619762 0.145952393613845 0.3 1.50407739677627 3.46573590279973 1 77 f20 m20 x 1 3 n 20 21 0 -0.158516787961936 -0.275261071462458 -0.3 1.6094379124341 3.04452243772342 1 78 f20 m20 x 2 0 n 20 21 0 -0.158516787961936 -0.275261071462458 -0.1 1.6094379124341 3.04452243772342 1 79 f20 m20 x 3 6 n 20 21 0 -0.158516787961936 -0.275261071462458 0.1 1.6094379124341 3.04452243772342 1 80 f20 m20 x 4 7 n 20 21 1 -0.158516787961936 -0.275261071462458 0.3 1.6094379124341 3.04452243772342 1 81 f21 m21 x 1 3 n 12 29 0 -0.669342411727926 0.047512320800593 -0.3 1.09861228866811 3.36729582998647 1 82 f21 m21 x 2 4 n 12 29 0 -0.669342411727926 0.047512320800593 -0.1 1.09861228866811 3.36729582998647 1 83 f21 m21 x 3 3 n 12 29 0 -0.669342411727926 0.047512320800593 0.1 1.09861228866811 3.36729582998647 1 84 f21 m21 x 4 4 n 12 29 1 -0.669342411727926 0.047512320800593 0.3 1.09861228866811 3.36729582998647 1 85 f22 m22 x 1 3 n 9 21 0 -0.957024484179707 -0.275261071462458 -0.3 0.810930216216329 3.04452243772342 1 86 f22 m22 x 2 4 n 9 21 0 -0.957024484179707 -0.275261071462458 -0.1 0.810930216216329 3.04452243772342 1 87 f22 m22 x 3 3 n 9 21 0 -0.957024484179707 -0.275261071462458 0.1 0.810930216216329 3.04452243772342 1 88 f22 m22 x 4 4 n 9 21 1 -0.957024484179707 -0.275261071462458 0.3 0.810930216216329 3.04452243772342 1 89 f23 m23 x 1 2 n 17 32 0 -0.321035717459711 0.145952393613845 -0.3 1.44691898293633 3.46573590279973 1 90 f23 m23 x 2 3 n 17 32 0 -0.321035717459711 0.145952393613845 -0.1 1.44691898293633 3.46573590279973 1 91 f23 m23 x 3 3 n 17 32 0 -0.321035717459711 0.145952393613845 0.1 1.44691898293633 3.46573590279973 1 92 f23 m23 x 4 5 n 17 32 1 -0.321035717459711 0.145952393613845 0.3 1.44691898293633 3.46573590279973 1 93 f24 m24 x 1 8 n 28 25 0 0.177955448659277 -0.10090768431768 -0.3 1.94591014905531 3.2188758248682 1 94 f24 m24 x 2 12 n 28 25 0 0.177955448659277 -0.10090768431768 -0.1 1.94591014905531 3.2188758248682 1 95 f24 m24 x 3 2 n 28 25 0 0.177955448659277 -0.10090768431768 0.1 1.94591014905531 3.2188758248682 1 96 f24 m24 x 4 8 n 28 25 1 0.177955448659277 -0.10090768431768 0.3 1.94591014905531 3.2188758248682 1 97 f25 m25 x 1 18 n 55 30 0 0.853084123716544 0.0814138724762743 -0.3 2.62103882411258 3.40119738166216 1 98 f25 m25 x 2 24 n 55 30 0 0.853084123716544 0.0814138724762743 -0.1 2.62103882411258 3.40119738166216 1 99 f25 m25 x 3 76 n 55 30 0 0.853084123716544 0.0814138724762743 0.1 2.62103882411258 3.40119738166216 1 100 f25 m25 x 4 25 n 55 30 1 0.853084123716544 0.0814138724762743 0.3 2.62103882411258 3.40119738166216 1 101 f26 m26 x 1 2 n 9 40 0 -0.957024484179707 0.369095944928055 -0.3 0.810930216216329 3.68887945411394 1 102 f26 m26 x 2 1 n 9 40 0 -0.957024484179707 0.369095944928055 -0.1 0.810930216216329 3.68887945411394 1 103 f26 m26 x 3 2 n 9 40 0 -0.957024484179707 0.369095944928055 0.1 0.810930216216329 3.68887945411394 1 104 f26 m26 x 4 1 n 9 40 1 -0.957024484179707 0.369095944928055 0.3 0.810930216216329 3.68887945411394 1 105 f27 m27 x 1 3 n 10 19 0 -0.851663968521881 -0.375344530019441 -0.3 0.916290731874155 2.94443897916644 1 106 f27 m27 x 2 1 n 10 19 0 -0.851663968521881 -0.375344530019441 -0.1 0.916290731874155 2.94443897916644 1 107 f27 m27 x 3 4 n 10 19 0 -0.851663968521881 -0.375344530019441 0.1 0.916290731874155 2.94443897916644 1 108 f27 m27 x 4 2 n 10 19 1 -0.851663968521881 -0.375344530019441 0.3 0.916290731874155 2.94443897916644 1 109 f28 m28 x 1 13 n 47 22 0 0.695898540194132 -0.228741055827565 -0.3 2.46385324059017 3.09104245335832 1 110 f28 m28 x 2 15 n 47 22 0 0.695898540194132 -0.228741055827565 -0.1 2.46385324059017 3.09104245335832 1 111 f28 m28 x 3 13 n 47 22 0 0.695898540194132 -0.228741055827565 0.1 2.46385324059017 3.09104245335832 1 112 f28 m28 x 4 12 n 47 22 1 0.695898540194132 -0.228741055827565 0.3 2.46385324059017 3.09104245335832 1 113 f29 m29 x 1 11 y 76 18 0 1.1764842787704 -0.429411751289717 -0.3 2.94443897916644 2.89037175789616 1 114 f29 m29 x 2 14 y 76 18 0 1.1764842787704 -0.429411751289717 -0.1 2.94443897916644 2.89037175789616 1 115 f29 m29 x 3 9 y 76 18 0 1.1764842787704 -0.429411751289717 0.1 2.94443897916644 2.89037175789616 1 116 f29 m29 x 4 8 y 76 18 1 1.1764842787704 -0.429411751289717 0.3 2.94443897916644 2.89037175789616 1 117 f30 m30 x 1 8 y 38 32 0 0.483337098210459 0.145952393613845 -0.3 2.2512917986065 3.46573590279973 1 118 f30 m30 x 2 7 y 38 32 0 0.483337098210459 0.145952393613845 -0.1 2.2512917986065 3.46573590279973 1 119 f30 m30 x 3 9 y 38 32 0 0.483337098210459 0.145952393613845 0.1 2.2512917986065 3.46573590279973 1 120 f30 m30 x 4 4 y 38 32 1 0.483337098210459 0.145952393613845 0.3 2.2512917986065 3.46573590279973 1 121 f31 m31 x 1 0 y 19 20 0 -0.209810082349486 -0.32405123563189 -0.3 1.55814461804655 2.99573227355399 1 122 f31 m31 x 2 4 y 19 20 0 -0.209810082349486 -0.32405123563189 -0.1 1.55814461804655 2.99573227355399 1 123 f31 m31 x 3 3 y 19 20 0 -0.209810082349486 -0.32405123563189 0.1 1.55814461804655 2.99573227355399 1 124 f31 m31 x 4 0 y 19 20 1 -0.209810082349486 -0.32405123563189 0.3 1.55814461804655 2.99573227355399 1 125 f32 m32 x 1 3 y 10 30 0 -0.851663968521881 0.0814138724762743 -0.3 0.916290731874155 3.40119738166216 1 126 f32 m32 x 2 6 y 10 30 0 -0.851663968521881 0.0814138724762743 -0.1 0.916290731874155 3.40119738166216 1 127 f32 m32 x 3 1 y 10 30 0 -0.851663968521881 0.0814138724762743 0.1 0.916290731874155 3.40119738166216 1 128 f32 m32 x 4 3 y 10 30 1 -0.851663968521881 0.0814138724762743 0.3 0.916290731874155 3.40119738166216 1 129 f33 m33 x 1 2 y 19 18 0 -0.209810082349486 -0.429411751289717 -0.3 1.55814461804655 2.89037175789616 1 130 f33 m33 x 2 6 y 19 18 0 -0.209810082349486 -0.429411751289717 -0.1 1.55814461804655 2.89037175789616 1 131 f33 m33 x 3 7 y 19 18 0 -0.209810082349486 -0.429411751289717 0.1 1.55814461804655 2.89037175789616 1 132 f33 m33 x 4 4 y 19 18 1 -0.209810082349486 -0.429411751289717 0.3 1.55814461804655 2.89037175789616 1 133 f34 m34 x 1 4 y 24 24 0 0.0238047688320189 -0.141729678837935 -0.3 1.79175946922806 3.17805383034795 1 134 f34 m34 x 2 3 y 24 24 0 0.0238047688320189 -0.141729678837935 -0.1 1.79175946922806 3.17805383034795 1 135 f34 m34 x 3 1 y 24 24 0 0.0238047688320189 -0.141729678837935 0.1 1.79175946922806 3.17805383034795 1 136 f34 m34 x 4 3 y 24 24 1 0.0238047688320189 -0.141729678837935 0.3 1.79175946922806 3.17805383034795 1 137 f35 m35 x 1 22 y 31 30 0 0.279738142969219 0.0814138724762743 -0.3 2.04769284336526 3.40119738166216 1 138 f35 m35 x 2 17 y 31 30 0 0.279738142969219 0.0814138724762743 -0.1 2.04769284336526 3.40119738166216 1 139 f35 m35 x 3 19 y 31 30 0 0.279738142969219 0.0814138724762743 0.1 2.04769284336526 3.40119738166216 1 140 f35 m35 x 4 16 y 31 30 1 0.279738142969219 0.0814138724762743 0.3 2.04769284336526 3.40119738166216 1 141 f36 m36 x 1 5 y 14 35 0 -0.515191731900668 0.235564552303532 -0.3 1.25276296849537 3.55534806148941 1 142 f36 m36 x 2 4 y 14 35 0 -0.515191731900668 0.235564552303532 -0.1 1.25276296849537 3.55534806148941 1 143 f36 m36 x 3 7 y 14 35 0 -0.515191731900668 0.235564552303532 0.1 1.25276296849537 3.55534806148941 1 144 f36 m36 x 4 4 y 14 35 1 -0.515191731900668 0.235564552303532 0.3 1.25276296849537 3.55534806148941 1 145 f37 m37 x 1 2 y 11 27 0 -0.756353788717556 -0.023946643181552 -0.3 1.01160091167848 3.29583686600433 1 146 f37 m37 x 2 4 y 11 27 0 -0.756353788717556 -0.023946643181552 -0.1 1.01160091167848 3.29583686600433 1 147 f37 m37 x 3 0 y 11 27 0 -0.756353788717556 -0.023946643181552 0.1 1.01160091167848 3.29583686600433 1 148 f37 m37 x 4 4 y 11 27 1 -0.756353788717556 -0.023946643181552 0.3 1.01160091167848 3.29583686600433 1 149 f38 m38 x 1 3 y 67 20 0 1.05044355787504 -0.32405123563189 -0.3 2.81839825827108 2.99573227355399 1 150 f38 m38 x 2 7 y 67 20 0 1.05044355787504 -0.32405123563189 -0.1 2.81839825827108 2.99573227355399 1 151 f38 m38 x 3 7 y 67 20 0 1.05044355787504 -0.32405123563189 0.1 2.81839825827108 2.99573227355399 1 152 f38 m38 x 4 7 y 67 20 1 1.05044355787504 -0.32405123563189 0.3 2.81839825827108 2.99573227355399 1 153 f39 m39 x 1 4 y 41 22 0 0.559323005188381 -0.228741055827565 -0.3 2.32727770558442 3.09104245335832 1 154 f39 m39 x 2 18 y 41 22 0 0.559323005188381 -0.228741055827565 -0.1 2.32727770558442 3.09104245335832 1 155 f39 m39 x 3 2 y 41 22 0 0.559323005188381 -0.228741055827565 0.1 2.32727770558442 3.09104245335832 1 156 f39 m39 x 4 5 y 41 22 1 0.559323005188381 -0.228741055827565 0.3 2.32727770558442 3.09104245335832 1 157 f40 m40 x 1 2 y 7 28 0 -1.20833891246061 0.0124210009893226 -0.3 0.559615787935423 3.3322045101752 1 158 f40 m40 x 2 1 y 7 28 0 -1.20833891246061 0.0124210009893226 -0.1 0.559615787935423 3.3322045101752 1 159 f40 m40 x 3 1 y 7 28 0 -1.20833891246061 0.0124210009893226 0.1 0.559615787935423 3.3322045101752 1 160 f40 m40 x 4 0 y 7 28 1 -1.20833891246061 0.0124210009893226 0.3 0.559615787935423 3.3322045101752 1 161 f41 m41 x 1 0 y 22 23 0 -0.0632066081576108 -0.184289293256731 -0.3 1.70474809223843 3.13549421592915 1 162 f41 m41 x 2 2 y 22 23 0 -0.0632066081576108 -0.184289293256731 -0.1 1.70474809223843 3.13549421592915 1 163 f41 m41 x 3 4 y 22 23 0 -0.0632066081576108 -0.184289293256731 0.1 1.70474809223843 3.13549421592915 1 164 f41 m41 x 4 0 y 22 23 1 -0.0632066081576108 -0.184289293256731 0.3 1.70474809223843 3.13549421592915 1 165 f42 m42 x 1 5 y 13 40 0 -0.58929970405439 0.369095944928055 -0.3 1.17865499634165 3.68887945411394 1 166 f42 m42 x 2 4 y 13 40 0 -0.58929970405439 0.369095944928055 -0.1 1.17865499634165 3.68887945411394 1 167 f42 m42 x 3 0 y 13 40 0 -0.58929970405439 0.369095944928055 0.1 1.17865499634165 3.68887945411394 1 168 f42 m42 x 4 3 y 13 40 1 -0.58929970405439 0.369095944928055 0.3 1.17865499634165 3.68887945411394 1 169 f43 m43 x 1 11 y 46 33 0 0.674392334973168 0.176724052280599 -0.3 2.4423470353692 3.49650756146648 1 170 f43 m43 x 2 14 y 46 33 0 0.674392334973168 0.176724052280599 -0.1 2.4423470353692 3.49650756146648 1 171 f43 m43 x 3 25 y 46 33 0 0.674392334973168 0.176724052280599 0.1 2.4423470353692 3.49650756146648 1 172 f43 m43 x 4 15 y 46 33 1 0.674392334973168 0.176724052280599 0.3 2.4423470353692 3.49650756146648 1 173 f44 m44 x 1 10 y 36 21 0 0.429269876940183 -0.275261071462458 -0.3 2.19722457733622 3.04452243772342 1 174 f44 m44 x 2 5 y 36 21 0 0.429269876940183 -0.275261071462458 -0.1 2.19722457733622 3.04452243772342 1 175 f44 m44 x 3 3 y 36 21 0 0.429269876940183 -0.275261071462458 0.1 2.19722457733622 3.04452243772342 1 176 f44 m44 x 4 8 y 36 21 1 0.429269876940183 -0.275261071462458 0.3 2.19722457733622 3.04452243772342 1 177 f45 m45 x 1 19 y 38 35 0 0.483337098210459 0.235564552303532 -0.3 2.2512917986065 3.55534806148941 1 178 f45 m45 x 2 7 y 38 35 0 0.483337098210459 0.235564552303532 -0.1 2.2512917986065 3.55534806148941 1 179 f45 m45 x 3 6 y 38 35 0 0.483337098210459 0.235564552303532 0.1 2.2512917986065 3.55534806148941 1 180 f45 m45 x 4 7 y 38 35 1 0.483337098210459 0.235564552303532 0.3 2.2512917986065 3.55534806148941 1 181 f46 m46 x 1 1 y 7 25 0 -1.20833891246061 -0.10090768431768 -0.3 0.559615787935423 3.2188758248682 1 182 f46 m46 x 2 1 y 7 25 0 -1.20833891246061 -0.10090768431768 -0.1 0.559615787935423 3.2188758248682 1 183 f46 m46 x 3 2 y 7 25 0 -1.20833891246061 -0.10090768431768 0.1 0.559615787935423 3.2188758248682 1 184 f46 m46 x 4 3 y 7 25 1 -1.20833891246061 -0.10090768431768 0.3 0.559615787935423 3.2188758248682 1 185 f47 m47 x 1 6 y 36 26 0 0.429269876940183 -0.061686971164399 -0.3 2.19722457733622 3.25809653802148 1 186 f47 m47 x 2 10 y 36 26 0 0.429269876940183 -0.061686971164399 -0.1 2.19722457733622 3.25809653802148 1 187 f47 m47 x 3 8 y 36 26 0 0.429269876940183 -0.061686971164399 0.1 2.19722457733622 3.25809653802148 1 188 f47 m47 x 4 8 y 36 26 1 0.429269876940183 -0.061686971164399 0.3 2.19722457733622 3.25809653802148 1 189 f48 m48 x 1 2 y 11 25 0 -0.756353788717556 -0.10090768431768 -0.3 1.01160091167848 3.2188758248682 1 190 f48 m48 x 2 1 y 11 25 0 -0.756353788717556 -0.10090768431768 -0.1 1.01160091167848 3.2188758248682 1 191 f48 m48 x 3 0 y 11 25 0 -0.756353788717556 -0.10090768431768 0.1 1.01160091167848 3.2188758248682 1 192 f48 m48 x 4 0 y 11 25 1 -0.756353788717556 -0.10090768431768 0.3 1.01160091167848 3.2188758248682 1 193 f49 m49 x 1 102 y 151 22 0 1.863030775299 -0.228741055827565 -0.3 3.63098547569503 3.09104245335832 1 194 f49 m49 x 2 65 y 151 22 0 1.863030775299 -0.228741055827565 -0.1 3.63098547569503 3.09104245335832 1 195 f49 m49 x 3 72 y 151 22 0 1.863030775299 -0.228741055827565 0.1 3.63098547569503 3.09104245335832 1 196 f49 m49 x 4 63 y 151 22 1 1.863030775299 -0.228741055827565 0.3 3.63098547569503 3.09104245335832 1 197 f50 m50 x 1 4 y 22 32 0 -0.0632066081576108 0.145952393613845 -0.3 1.70474809223843 3.46573590279973 1 198 f50 m50 x 2 3 y 22 32 0 -0.0632066081576108 0.145952393613845 -0.1 1.70474809223843 3.46573590279973 1 199 f50 m50 x 3 2 y 22 32 0 -0.0632066081576108 0.145952393613845 0.1 1.70474809223843 3.46573590279973 1 200 f50 m50 x 4 4 y 22 32 1 -0.0632066081576108 0.145952393613845 0.3 1.70474809223843 3.46573590279973 1 201 f51 m51 x 1 8 y 41 25 0 0.559323005188381 -0.10090768431768 -0.3 2.32727770558442 3.2188758248682 1 202 f51 m51 x 2 6 y 41 25 0 0.559323005188381 -0.10090768431768 -0.1 2.32727770558442 3.2188758248682 1 203 f51 m51 x 3 5 y 41 25 0 0.559323005188381 -0.10090768431768 0.1 2.32727770558442 3.2188758248682 1 204 f51 m51 x 4 7 y 41 25 1 0.559323005188381 -0.10090768431768 0.3 2.32727770558442 3.2188758248682 1 205 f52 m52 x 1 1 y 32 35 0 0.3114868412838 0.235564552303532 -0.3 2.07944154167984 3.55534806148941 1 206 f52 m52 x 2 3 y 32 35 0 0.3114868412838 0.235564552303532 -0.1 2.07944154167984 3.55534806148941 1 207 f52 m52 x 3 1 y 32 35 0 0.3114868412838 0.235564552303532 0.1 2.07944154167984 3.55534806148941 1 208 f52 m52 x 4 5 y 32 35 1 0.3114868412838 0.235564552303532 0.3 2.07944154167984 3.55534806148941 1 209 f53 m53 x 1 18 y 56 21 0 0.871102629219222 -0.275261071462458 -0.3 2.63905732961526 3.04452243772342 1 210 f53 m53 x 2 11 y 56 21 0 0.871102629219222 -0.275261071462458 -0.1 2.63905732961526 3.04452243772342 1 211 f53 m53 x 3 28 y 56 21 0 0.871102629219222 -0.275261071462458 0.1 2.63905732961526 3.04452243772342 1 212 f53 m53 x 4 13 y 56 21 1 0.871102629219222 -0.275261071462458 0.3 2.63905732961526 3.04452243772342 1 213 f54 m54 x 1 6 y 24 41 0 0.0238047688320189 0.393788557518427 -0.3 1.79175946922806 3.71357206670431 1 214 f54 m54 x 2 3 y 24 41 0 0.0238047688320189 0.393788557518427 -0.1 1.79175946922806 3.71357206670431 1 215 f54 m54 x 3 4 y 24 41 0 0.0238047688320189 0.393788557518427 0.1 1.79175946922806 3.71357206670431 1 216 f54 m54 x 4 0 y 24 41 1 0.0238047688320189 0.393788557518427 0.3 1.79175946922806 3.71357206670431 1 217 f55 m55 x 1 3 y 16 32 0 -0.381660339276146 0.145952393613845 -0.3 1.38629436111989 3.46573590279973 1 218 f55 m55 x 2 5 y 16 32 0 -0.381660339276146 0.145952393613845 -0.1 1.38629436111989 3.46573590279973 1 219 f55 m55 x 3 4 y 16 32 0 -0.381660339276146 0.145952393613845 0.1 1.38629436111989 3.46573590279973 1 220 f55 m55 x 4 3 y 16 32 1 -0.381660339276146 0.145952393613845 0.3 1.38629436111989 3.46573590279973 1 221 f56 m56 x 1 1 y 22 26 0 -0.0632066081576108 -0.061686971164399 -0.3 1.70474809223843 3.25809653802148 1 222 f56 m56 x 2 23 y 22 26 0 -0.0632066081576108 -0.061686971164399 -0.1 1.70474809223843 3.25809653802148 1 223 f56 m56 x 3 19 y 22 26 0 -0.0632066081576108 -0.061686971164399 0.1 1.70474809223843 3.25809653802148 1 224 f56 m56 x 4 8 y 22 26 1 -0.0632066081576108 -0.061686971164399 0.3 1.70474809223843 3.25809653802148 1 225 f57 m57 x 1 2 y 25 21 0 0.0646267633522741 -0.275261071462458 -0.3 1.83258146374831 3.04452243772342 1 226 f57 m57 x 2 3 y 25 21 0 0.0646267633522741 -0.275261071462458 -0.1 1.83258146374831 3.04452243772342 1 227 f57 m57 x 3 0 y 25 21 0 0.0646267633522741 -0.275261071462458 0.1 1.83258146374831 3.04452243772342 1 228 f57 m57 x 4 1 y 25 21 1 0.0646267633522741 -0.275261071462458 0.3 1.83258146374831 3.04452243772342 1 229 f58 m58 x 1 0 y 13 36 0 -0.58929970405439 0.263735429270229 -0.3 1.17865499634165 3.58351893845611 1 230 f58 m58 x 2 0 y 13 36 0 -0.58929970405439 0.263735429270229 -0.1 1.17865499634165 3.58351893845611 1 231 f58 m58 x 3 0 y 13 36 0 -0.58929970405439 0.263735429270229 0.1 1.17865499634165 3.58351893845611 1 232 f58 m58 x 4 0 y 13 36 1 -0.58929970405439 0.263735429270229 0.3 1.17865499634165 3.58351893845611 1 233 f59 m59 x 1 1 y 12 37 0 -0.669342411727926 0.291134403458343 -0.3 1.09861228866811 3.61091791264422 1 234 f59 m59 x 2 4 y 12 37 0 -0.669342411727926 0.291134403458343 -0.1 1.09861228866811 3.61091791264422 1 235 f59 m59 x 3 3 y 12 37 0 -0.669342411727926 0.291134403458343 0.1 1.09861228866811 3.61091791264422 1 236 f59 m59 x 4 2 y 12 37 1 -0.669342411727926 0.291134403458343 0.3 1.09861228866811 3.61091791264422 ;;;; run set loc lbase.trt qua lbase.trt=lbase*trt set ple 0 set chains 4 set iter 5000 fpm seizures lik poisson nqtl 0 s cov lage + trt + lbase + lbase.trt + V4