[o] GLIM 3.77 update 1 (copyright)1985 Royal Statistical Society, London [o] [i] ? $units 18 $data feed yield $read [i] $REA? 1 38 2 109 3 124 1 86 2 75 3 101 [i] $REA? 2 25 3 86 1 72 3 76 1 35 2 63 [i] $REA? 3 15 1 39 2 27 2 46 3 34 1 1 [i] ? $cal period=%gl(3,6) : cow=%gl(6,1) $ [i] ? $fac period 3 cow 6 feed 3 $ [i] ? $yvar yield $fit cow $ [o] deviance = 14581. [o] d.f. = 12 [o] [i] ? $fit +period $dis e $ [o] deviance = 3101.2 (change = -11480.) [o] d.f. = 10 (change = -2 ) [o] [o] estimate s.e. parameter [o] 1 56.39 11.74 1 [o] 2 52.00 14.38 COW(2) [o] 3 48.33 14.38 COW(3) [o] 4 43.33 14.38 COW(4) [o] 5 22.00 14.38 COW(5) [o] 6 29.00 14.38 COW(6) [o] 7 -29.33 10.17 PERI(2) [o] 8 -61.83 10.17 PERI(3) [o] scale parameter taken as 310.1 [o] [i] ? $fit feed $dis e $ [o] deviance = 18086. [o] d.f. = 15 [o] [o] estimate s.e. parameter [o] 1 45.17 14.18 1 [o] 2 12.33 20.05 FEED(2) [o] 3 27.50 20.05 FEED(3) [o] scale parameter taken as 1206. [o] [i] ? $fit cow+feed $dis e $ [o] deviance = 12305. [o] d.f. = 10 [o] [o] estimate s.e. parameter [o] 1 12.72 23.39 1 [o] 2 52.00 28.64 COW(2) [o] 3 48.33 28.64 COW(3) [o] 4 43.33 28.64 COW(4) [o] 5 22.00 28.64 COW(5) [o] 6 29.00 28.64 COW(6) [o] 7 12.33 20.25 FEED(2) [o] 8 27.50 20.25 FEED(3) [o] scale parameter taken as 1230. [o] [i] ? $fit +period $dis e $ [o] deviance = 824.44 (change = -11480.) [o] d.f. = 8 (change = -2 ) [o] [o] estimate s.e. parameter [o] 1 43.11 7.567 1 [o] 2 52.00 8.289 COW(2) [o] 3 48.33 8.289 COW(3) [o] 4 43.33 8.289 COW(4) [o] 5 22.00 8.289 COW(5) [o] 6 29.00 8.289 COW(6) [o] 7 12.33 5.861 FEED(2) [o] 8 27.50 5.861 FEED(3) [o] 9 -29.33 5.861 PERI(2) [o] 10 -61.83 5.861 PERI(3) [o] scale parameter taken as 103.1 [o] [i] ? $fit +feed.period $dis e $ [o] deviance = 233.33 (change = -591.1) [o] d.f. = 4 (change = -4 ) [o] [o] estimate s.e. parameter [o] 1 34.00 6.736 1 [o] 2 62.33 8.819 COW(2) [o] 3 48.67 8.819 COW(3) [o] 4 56.00 8.051 COW(4) [o] 5 19.67 8.051 COW(5) [o] 6 29.33 8.051 COW(6) [o] 7 17.00 9.860 FEED(2) [o] 8 39.50 9.860 FEED(3) [o] 9 -14.67 9.860 PERI(2) [o] 10 -59.83 9.860 PERI(3) [o] 11 -7.000 15.28 FEED(2).PERI(2) [o] 12 -7.000 15.28 FEED(2).PERI(3) [o] 13 -37.00 15.28 FEED(3).PERI(2) [o] 14 1.000 15.28 FEED(3).PERI(3) [o] scale parameter taken as 58.33 [o] [i] ? $fit -feed.period $ [o] deviance = 824.44 (change = +591.1) [o] d.f. = 8 (change = +4 ) [o] [i] ? $cal xfeed=feed : xperiod=period $fit -feed-period+xfeed+xperiod $dis e$ [o] deviance = 842.50 (change = +18.06) [o] d.f. = 10 (change = +2 ) [o] [o] estimate s.e. parameter [o] 1 60.33 9.179 1 [o] 2 52.00 7.494 COW(2) [o] 3 48.33 7.494 COW(3) [o] 4 43.33 7.494 COW(4) [o] 5 22.00 7.494 COW(5) [o] 6 29.00 7.494 COW(6) [o] 7 13.75 2.650 XFEE [o] 8 -30.92 2.650 XPER [o] scale parameter taken as 84.25 [o] [i] ? $cal last=0 : last(18)=1 $print last $ [o] 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. [o] 0. 0. 0. 0. 0. 0. 0. 1.000 [i] ? $fit +last $dis e $ [o] deviance = 668.28 (change = -174.2) [o] d.f. = 9 (change = -1 ) [o] [o] estimate s.e. parameter [o] 1 60.33 8.617 1 [o] 2 52.00 7.036 COW(2) [o] 3 48.33 7.036 COW(3) [o] 4 43.33 7.036 COW(4) [o] 5 22.00 7.036 COW(5) [o] 6 35.22 8.124 COW(6) [o] 7 12.19 2.687 XFEE [o] 8 -29.36 2.687 XPER [o] 9 -18.67 12.19 LAST [o] scale parameter taken as 74.25 [o] [i] ? $stop