home *** CD-ROM | disk | FTP | other *** search
/ gimli.bioss.sari.ac.uk / gimli.bioss.sari.ac.uk.zip / gimli.bioss.sari.ac.uk / pub / courses / ExptDesignAnalysis / lambs.out < prev    next >
Text File  |  1995-02-16  |  19KB  |  477 lines

  1. Genstat 5  Release 3.1  (Sun/Unix)                Fri Apr  2 11:50:22 1993
  2. Copyright 1992, Lawes Agricultural Trust (Rothamsted Experimental Station)
  3.  
  4.    1  units[250]
  5.    2  factor[lev=10] treat
  6.    3  open 'lambs.dat';ch=2
  7.    4  read [ch=2] treat,w[1...8],td,wg
  8.  
  9.     Identifier   Minimum      Mean   Maximum    Values   Missing
  10.           w[1]     18.50     26.98     35.00       250         0
  11.           w[2]     17.50     26.45     33.50       250         0
  12.           w[3]     17.00     26.74     32.50       250         1
  13.           w[4]     18.00     27.60     35.00       250         1
  14.           w[5]     17.50     27.45     36.50       250         1
  15.           w[6]     21.00     30.54     38.00       250         3
  16.           w[7]     23.00     31.42     39.00       250        31
  17.           w[8]     22.00     30.23     39.50       250        82
  18.             td     56.00     79.75     98.00       250       128
  19.             wg     1.000     5.880    13.000       250       129
  20.  
  21.     Identifier    Values   Missing    Levels
  22.          treat       250         0        10
  23.  
  24.    5  scalar mvi
  25.    6  calc mvi=0/0
  26.    7  
  27.    8  for yy=w[5...8],wg;zz=mv[5...8],mvg
  28.    9  restrict zz;cond=yy==mvi
  29.   10  calc zz=2
  30.   11  restrict zz; cond=yy.ne.mvi
  31.   12  calc zz=1
  32.   13  restrict zz
  33.   14  endfor
  34.   15  
  35.   16  calc g[2...8]=w[2...8]-w[1]
  36.   17  factor [lev=!(0,1)]mf[5...8],mfg
  37.   18  calc mf[5...8],mfg=mv[5...8],mvg-1
  38.   19  
  39.   20  for ff=mf[6...8],mfg;vv=g[5...8];rr=mf[5...8]
  40.   21  restrict vv;cond=rr==0
  41.   22  tabulate[class=treat,ff;marg=y;prin=m,c]vv
  42.   23  model vv
  43.   24  fit[prin=m,s,a,e] treat+ff
  44.   25  restrict vv
  45.   26  endfor
  46.  
  47.         mf[6]        0.00                    1.00                  Margin
  48.                     Count        Mean       Count        Mean       Count        Mean
  49.         treat
  50.             1          25      1.4000           0           *          25      1.4000
  51.             2          25      0.1000           0           *          25      0.1000
  52.             3          25      1.3800           0           *          25      1.3800
  53.             4          24      0.8750           0           *          24      0.8750
  54.             5          25     -0.6600           0           *          25     -0.6600
  55.             6          25     -1.0000           0           *          25     -1.0000
  56.             7          24      2.4583           1      1.0000          25      2.4000
  57.             8          24     -0.6250           1      0.5000          25     -0.5800
  58.             9          25     -0.7400           0           *          25     -0.7400
  59.            10          25      1.4000           0           *          25      1.4000
  60.  
  61.        Margin         247      0.4534           2      0.7500         249      0.4558
  62.  
  63.  
  64.  
  65. 26..................................................................................................................................
  66.  
  67.  
  68. ***** Regression Analysis *****
  69.  
  70.  Response variate: g[5]
  71.      Fitted terms: Constant + treat + mf[6]
  72.  
  73.  
  74. *** Summary of analysis ***
  75.  
  76.               d.f.         s.s.         m.s.      v.r.
  77. Regression      10         315.       31.454      4.77
  78. Residual       238        1569.        6.594
  79. Total          248        1884.        7.597
  80.  
  81. Percentage variance accounted for 13.2
  82.  
  83. Standard error of observations is estimated to be 2.57
  84.  
  85. * MESSAGE: The following units have large standardized residuals:
  86.                   66       -2.93
  87.                  155        3.81
  88.  
  89. * MESSAGE: The following units have high leverage:
  90.                  174       0.520
  91.                  184       0.520
  92.  
  93.  
  94. *** Estimates of regression coefficients ***
  95.  
  96.                               estimate         s.e.    t(238)
  97. Constant                         1.400        0.514      2.73
  98. treat 2                         -1.300        0.726     -1.79
  99. treat 3                         -0.020        0.726     -0.03
  100. treat 4                         -0.525        0.734     -0.72
  101. treat 5                         -2.060        0.726     -2.84
  102. treat 6                         -2.400        0.726     -3.30
  103. treat 7                          1.007        0.730      1.38
  104. treat 8                         -1.973        0.730     -2.70
  105. treat 9                         -2.140        0.726     -2.95
  106. treat 10                         0.000        0.726      0.00
  107. mf[6] 1                          -0.17         1.85     -0.09
  108.  
  109.  
  110. *** Accumulated analysis of variance ***
  111.  
  112. Change                         d.f.         s.s.         m.s.      v.r.
  113. + treat                           9      314.489       34.943      5.30
  114. + mf[6]                           1        0.053        0.053      0.01
  115. Residual                        238     1569.472        6.594
  116.  
  117. Total                           248     1884.014        7.597
  118.  
  119.  
  120.         mf[7]        0.00                    1.00                  Margin
  121.                     Count        Mean       Count        Mean       Count        Mean
  122.         treat
  123.             1          21       2.952           4       6.250          25       3.480
  124.             2          24       2.896           1       7.500          25       3.080
  125.             3          22       4.727           3       5.667          25       4.840
  126.             4          19       4.263           5       6.300          24       4.688
  127.             5          25       2.160           0           *          25       2.160
  128.             6          25       1.620           0           *          25       1.620
  129.             7          20       3.700           4       5.875          24       4.062
  130.             8          24       3.646           0           *          24       3.646
  131.             9          20       3.175           5       6.200          25       3.780
  132.            10          18       4.083           7       4.929          25       4.320
  133.  
  134.        Margin         218       3.255          29       5.862         247       3.561
  135.  
  136.  
  137.  
  138. 26..................................................................................................................................
  139.  
  140.  
  141. ***** Regression Analysis *****
  142.  
  143.  Response variate: g[6]
  144.      Fitted terms: Constant + treat + mf[7]
  145.  
  146.  
  147. *** Summary of analysis ***
  148.  
  149.               d.f.         s.s.         m.s.      v.r.
  150. Regression      10         348.       34.821      5.76
  151. Residual       236        1428.        6.050
  152. Total          246        1776.        7.220
  153.  
  154. Percentage variance accounted for 16.2
  155.  
  156. Standard error of observations is estimated to be 2.46
  157.  
  158. * MESSAGE: The following units have large standardized residuals:
  159.                   99       -3.43
  160.  
  161.  
  162. *** Estimates of regression coefficients ***
  163.  
  164.                               estimate         s.e.    t(236)
  165. Constant                         3.139        0.499      6.29
  166. treat 2                         -0.144        0.698     -0.21
  167. treat 3                          1.445        0.696      2.08
  168. treat 4                          1.104        0.703      1.57
  169. treat 5                         -0.979        0.700     -1.40
  170. treat 6                         -1.519        0.700     -2.17
  171. treat 7                          0.568        0.703      0.81
  172. treat 8                          0.507        0.708      0.72
  173. treat 9                          0.215        0.696      0.31
  174. treat 10                         0.584        0.698      0.84
  175. mf[7] 1                          2.132        0.510      4.18
  176.  
  177.  
  178. *** Accumulated analysis of variance ***
  179.  
  180. Change                         d.f.         s.s.         m.s.      v.r.
  181. + treat                           9      242.367       26.930      4.45
  182. + mf[7]                           1      105.848      105.848     17.49
  183. Residual                        236     1427.874        6.050
  184.  
  185. Total                           246     1776.089        7.220
  186.  
  187.  
  188.         mf[8]        0.00                    1.00                  Margin
  189.                     Count        Mean       Count        Mean       Count        Mean
  190.         treat
  191.             1          14       4.893           7       6.571          21       5.452
  192.             2          21       2.595           3       4.667          24       2.854
  193.             3          17       4.912           5       5.100          22       4.955
  194.             4          15       4.367           4       4.875          19       4.474
  195.             5          21       3.214           4       5.125          25       3.520
  196.             6          22       3.386           3       6.167          25       3.720
  197.             7          14       4.714           7       7.714          21       5.714
  198.             8          15       5.133           9       6.111          24       5.500
  199.             9          14       5.143           6       7.417          20       5.825
  200.            10          15       4.867           3       7.667          18       5.333
  201.  
  202.        Margin         168       4.179          51       6.284         219       4.669
  203.  
  204.  
  205.  
  206. 26..................................................................................................................................
  207.  
  208.  
  209. ***** Regression Analysis *****
  210.  
  211.  Response variate: g[7]
  212.      Fitted terms: Constant + treat + mf[8]
  213.  
  214.  
  215. *** Summary of analysis ***
  216.  
  217.               d.f.         s.s.         m.s.      v.r.
  218. Regression      10         339.       33.886      5.46
  219. Residual       208        1292.        6.211
  220. Total          218        1631.        7.480
  221.  
  222. Percentage variance accounted for 17.0
  223.  
  224. Standard error of observations is estimated to be 2.49
  225.  
  226. * MESSAGE: The following units have large standardized residuals:
  227.                   99       -2.93
  228.  
  229.  
  230. *** Estimates of regression coefficients ***
  231.  
  232.                               estimate         s.e.    t(208)
  233. Constant                         4.868        0.561      8.68
  234. treat 2                         -2.233        0.750     -2.98
  235. treat 3                         -0.312        0.762     -0.41
  236. treat 4                         -0.763        0.791     -0.97
  237. treat 5                         -1.628        0.741     -2.20
  238. treat 6                         -1.358        0.743     -1.83
  239. treat 7                          0.262        0.769      0.34
  240. treat 8                         -0.025        0.745     -0.03
  241. treat 9                          0.431        0.779      0.55
  242. treat 10                         0.173        0.803      0.22
  243. mf[8] 1                          1.753        0.408      4.30
  244.  
  245.  
  246. *** Accumulated analysis of variance ***
  247.  
  248. Change                         d.f.         s.s.         m.s.      v.r.
  249. + treat                           9      224.162       24.907      4.01
  250. + mf[8]                           1      114.701      114.701     18.47
  251. Residual                        208     1291.885        6.211
  252.  
  253. Total                           218     1630.749        7.480
  254.  
  255.  
  256.           mfg        0.00                    1.00                  Margin
  257.                     Count        Mean       Count        Mean       Count        Mean
  258.         treat
  259.             1           9       5.056           5       4.400          14       4.821
  260.             2           4       5.125          17       2.471          21       2.976
  261.             3           7       5.429          10       4.150          17       4.676
  262.             4           6       4.500           9       2.722          15       3.433
  263.             5           1       5.000          20       1.275          21       1.452
  264.             6           7       4.429          15       2.433          22       3.068
  265.             7           5       7.900           9       2.889          14       4.679
  266.             8           9       5.833           6       3.417          15       4.867
  267.             9           3       3.667          11       3.727          14       3.714
  268.            10           5       5.700          10       5.900          15       5.833
  269.  
  270.        Margin          56       5.330         112       3.022         168       3.792
  271.  
  272.  
  273.  
  274. 26..................................................................................................................................
  275.  
  276.  
  277. ***** Regression Analysis *****
  278.  
  279.  Response variate: g[8]
  280.      Fitted terms: Constant + treat + mfg
  281.  
  282.  
  283. *** Summary of analysis ***
  284.  
  285.               d.f.         s.s.         m.s.      v.r.
  286. Regression      10         370.       37.046      5.30
  287. Residual       157        1097.        6.986
  288. Total          167        1467.        8.786
  289.  
  290. Percentage variance accounted for 20.5
  291.  
  292. Standard error of observations is estimated to be 2.64
  293.  
  294. * MESSAGE: The following units have large standardized residuals:
  295.                  154        2.82
  296.                  159       -2.76
  297.  
  298.  
  299. *** Estimates of regression coefficients ***
  300.  
  301.                               estimate         s.e.    t(157)
  302. Constant                         5.477        0.726      7.55
  303. treat 2                         -1.015        0.936     -1.09
  304. treat 3                          0.279        0.960      0.29
  305. treat 4                         -0.943        0.989     -0.95
  306. treat 5                         -2.277        0.953     -2.39
  307. treat 6                         -1.158        0.916     -1.26
  308. treat 7                           0.38         1.01      0.38
  309. treat 8                          0.124        0.982      0.13
  310. treat 9                          -0.32         1.02     -0.32
  311. treat 10                         1.580        0.993      1.59
  312. mfg 1                           -1.834        0.464     -3.95
  313.  
  314.  
  315. *** Accumulated analysis of variance ***
  316.  
  317. Change                         d.f.         s.s.         m.s.      v.r.
  318. + treat                           9      261.435       29.048      4.16
  319. + mfg                             1      109.029      109.029     15.61
  320. Residual                        157     1096.744        6.986
  321.  
  322. Total                           167     1467.208        8.786
  323.  
  324.   27  
  325.   28  vari rate
  326.   29  restrict wg,td,rate;cond=mfg==0
  327.   30  calc rate=wg/td
  328.   31  restrict wg,td,rate;cond=mfg.ne.0
  329.   32  calc rate=g[7]/70
  330.   33  restrict wg,td,rate
  331.   34  tabulate[class=treat,mfg;marg=y;prin=m,c] rate
  332.  
  333.           mfg        0.00                    1.00                  Margin
  334.                     Count        Mean       Count        Mean       Count        Mean
  335.         treat
  336.             1          17     0.07920           8     0.06786          25     0.07557
  337.             2           7     0.06942          18     0.03413          25     0.04401
  338.             3          12     0.07368          13     0.06758          25     0.07051
  339.             4          13     0.06996          12     0.06494          25     0.06766
  340.             5           5     0.06878          20     0.04286          25     0.04804
  341.             6           9     0.05578          16     0.03884          25     0.04494
  342.             7          15     0.09866          10     0.05786          25     0.08234
  343.             8          16     0.07589           9     0.05625          25     0.06935
  344.             9          13     0.09462          12     0.07381          25     0.08463
  345.            10          14     0.08263          11     0.06753          25     0.07599
  346.  
  347.        Margin         121     0.07895         129     0.05422         250     0.06629
  348.  
  349.   35  model rate
  350.   36  fit[prin=m,s,a,e] treat+mfg
  351.  
  352.  
  353. 36..................................................................................................................................
  354.  
  355.  
  356. ***** Regression Analysis *****
  357.  
  358.  Response variate: rate
  359.      Fitted terms: Constant + treat + mfg
  360.  
  361.  
  362. *** Summary of analysis ***
  363.  
  364.               d.f.         s.s.         m.s.      v.r.
  365. Regression      10       0.0730     0.007302      5.52
  366. Residual       237       0.3134     0.001323
  367. Total          247       0.3864     0.001565
  368.  
  369. Percentage variance accounted for 15.5
  370.  
  371. Standard error of observations is estimated to be 0.0364
  372.  
  373. * MESSAGE: The following units have large standardized residuals:
  374.                   90        3.12
  375.  
  376.  
  377. *** Estimates of regression coefficients ***
  378.  
  379.                               estimate         s.e.    t(237)
  380. Constant                       0.08173      0.00744     10.99
  381. treat 2                        -0.0239       0.0105     -2.28
  382. treat 3                        -0.0012       0.0103     -0.12
  383. treat 4                        -0.0053       0.0104     -0.50
  384. treat 5                        -0.0183       0.0105     -1.73
  385. treat 6                        -0.0245       0.0104     -2.35
  386. treat 7                         0.0083       0.0103      0.81
  387. treat 8                        -0.0060       0.0104     -0.57
  388. treat 9                         0.0121       0.0103      1.18
  389. treat 10                        0.0027       0.0103      0.26
  390. mfg 1                         -0.01925      0.00485     -3.97
  391.  
  392.  
  393. *** Accumulated analysis of variance ***
  394.  
  395. Change                         d.f.         s.s.         m.s.      v.r.
  396. + treat                           9     0.052201     0.005800      4.39
  397. + mfg                             1     0.020816     0.020816     15.74
  398. Residual                        237     0.313433     0.001323
  399.  
  400. Total                           247     0.386450     0.001565
  401.  
  402.   37  treat treat*mfg
  403.   38  anova rate
  404.  
  405. * MESSAGE: non-orthogonality between treatment terms.
  406.   The effects (printed or used to calculate means), the efficiency
  407.   factor and the sum of squares for each treatment term are for
  408.   that term eliminating previous terms in the TREATMENT formula
  409.   and ignoring subsequent terms.
  410.  
  411.  
  412. 38..................................................................................................................................
  413.  
  414. ***** Analysis of variance *****
  415.  
  416. Variate: rate
  417.  
  418. Source of variation           d.f.(m.v.)      s.s.       m.s.    v.r.
  419. treat                            9        0.052131   0.005792    4.31
  420. mfg                              1        0.020931   0.020931   15.58
  421. treat.mfg                        9        0.007097   0.000789    0.59
  422. Residual                       228(2)     0.306393   0.001344
  423. Total                          247(2)     0.386450
  424.  
  425.  
  426. ***** Information summary *****
  427.  
  428. Model term                      e.f.  non-orthogonal terms
  429.   mfg                          0.911  treat
  430.  
  431.  
  432. * MESSAGE: the following units have large residuals.
  433.  
  434. *units* 90         0.1175   s.e. 0.0350
  435. *units* 99        -0.1078   s.e. 0.0350
  436.  
  437.  
  438. ***** Tables of means *****
  439.  
  440. Variate: rate
  441.  
  442. Grand mean  0.0662
  443.  
  444.     treat        1        2        3        4        5        6        7        8        9       10
  445.             0.0756   0.0440   0.0705   0.0675   0.0480   0.0449   0.0823   0.0688   0.0846   0.0760
  446.  
  447.       mfg     0.00     1.00
  448.             0.0761   0.0570
  449.      rep.      121      129
  450.  
  451.  
  452. *** Table of means for treat.mfg cannot be calculated
  453.    (contains mutually non-orthogonal components).
  454.  
  455. *** Standard errors of differences of means ***
  456.  
  457. Table                treat         mfg
  458. rep.                    25     unequal
  459. s.e.d.             0.01037     0.00486
  460.  
  461. (Not adjusted for missing values)
  462.  
  463.  
  464. ***** Missing values *****
  465.  
  466. Variate: rate
  467.  
  468.  Unit  estimate
  469.    80    0.0649
  470.   184    0.0563
  471.  
  472. Max. no. iterations 2
  473.  
  474.   39  stop
  475.  
  476. ******** End of job.  Maximum of 13860 data units used at line 26 (37924 left)
  477.