[o] GLIM 3.77 update 1 (copyright)1985 Royal Statistical Society, London
 [o] 
 [i] ? $units 18 $data feed yield $dinput 11 80 $
 [i] File name? cow
 [i]   1 38   2 109   3 124   1 86   2 75   3 101
 [i]   2 25   3  86   1  72   3 76   1 35   2  63
 [i]   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 $
 [i] ? $cal xperiod=period : xfeed=feed $
 [i] ? $fit cow+xfeed+xperiod $dis e $
 [o] deviance = 842.50
 [o]     d.f. =  10
 [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] ? $dis r $
 [o]   unit  observed    fitted   residual
 [o]      1    38.000    43.167     -5.167
 [o]      2   109.000   108.917      0.083
 [o]      3   124.000   119.000      5.000
 [o]      4    86.000    86.500     -0.500
 [o]      5    75.000    78.917     -3.917
 [o]      6   101.000    99.667      1.333
 [o]      7    25.000    26.000     -1.000
 [o]      8    86.000    91.750     -5.750
 [o]      9    72.000    60.583     11.417
 [o]     10    76.000    83.083     -7.083
 [o]     11    35.000    34.250      0.750
 [o]     12    63.000    55.000      8.000
 [o]     13    15.000     8.833      6.167
 [o]     14    39.000    33.333      5.667
 [o]     15    27.000    43.417    -16.417
 [o]     16    46.000    38.417      7.583
 [o]     17    34.000    30.833      3.167
 [o]     18     1.000    10.333     -9.333
 [o] 
 [i] ? $cal r=yield-%fv $
 [i] ? $extract %vl $
 [i] ? $cal z=r/sqrt(%sc-%vl) $print z $
 [f] ** undefined or deleted identifier, at [%sc-%vl) $]
 [f] 
 [h] The identifier SQRT used in the $CAL directive is undefined, possibly because
 [h] it was previously deleted.  Its type is not defined by its present usage.
 [h] 
 [i] ? $cal z=r/%sqrt(%sc-%vl) $print z $
 [o]  -0.7961  0.0119  0.7704 -0.0770 -0.5587  0.2054 -0.1334 -0.8202   1.629
 [o]   -1.010  0.1070   1.067  0.9501  0.8731  -2.342   1.082  0.4879  -1.438
 [i] ? $cal t=z/%sqrt((10-z*z)/9) $print t $
 [o]  -0.7803  0.0113  0.7535 -0.0731 -0.5385  0.1953 -0.1267 -0.8057   1.802
 [o]   -1.012  0.1016   1.076  0.9450  0.8618  -3.306   1.092  0.4685  -1.532
 [i] ? $dis e $
 [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] ? $dis d $
 [o] Deviance is  842.5  on   10 d.f.   from  18 observations
 [o]   n.b. current and previous models are not necessarily nested
 [o] 
 [i] ? $stop