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 >
Wrap
Text File
|
1995-02-16
|
19KB
|
477 lines
Genstat 5 Release 3.1 (Sun/Unix) Fri Apr 2 11:50:22 1993
Copyright 1992, Lawes Agricultural Trust (Rothamsted Experimental Station)
1 units[250]
2 factor[lev=10] treat
3 open 'lambs.dat';ch=2
4 read [ch=2] treat,w[1...8],td,wg
Identifier Minimum Mean Maximum Values Missing
w[1] 18.50 26.98 35.00 250 0
w[2] 17.50 26.45 33.50 250 0
w[3] 17.00 26.74 32.50 250 1
w[4] 18.00 27.60 35.00 250 1
w[5] 17.50 27.45 36.50 250 1
w[6] 21.00 30.54 38.00 250 3
w[7] 23.00 31.42 39.00 250 31
w[8] 22.00 30.23 39.50 250 82
td 56.00 79.75 98.00 250 128
wg 1.000 5.880 13.000 250 129
Identifier Values Missing Levels
treat 250 0 10
5 scalar mvi
6 calc mvi=0/0
7
8 for yy=w[5...8],wg;zz=mv[5...8],mvg
9 restrict zz;cond=yy==mvi
10 calc zz=2
11 restrict zz; cond=yy.ne.mvi
12 calc zz=1
13 restrict zz
14 endfor
15
16 calc g[2...8]=w[2...8]-w[1]
17 factor [lev=!(0,1)]mf[5...8],mfg
18 calc mf[5...8],mfg=mv[5...8],mvg-1
19
20 for ff=mf[6...8],mfg;vv=g[5...8];rr=mf[5...8]
21 restrict vv;cond=rr==0
22 tabulate[class=treat,ff;marg=y;prin=m,c]vv
23 model vv
24 fit[prin=m,s,a,e] treat+ff
25 restrict vv
26 endfor
mf[6] 0.00 1.00 Margin
Count Mean Count Mean Count Mean
treat
1 25 1.4000 0 * 25 1.4000
2 25 0.1000 0 * 25 0.1000
3 25 1.3800 0 * 25 1.3800
4 24 0.8750 0 * 24 0.8750
5 25 -0.6600 0 * 25 -0.6600
6 25 -1.0000 0 * 25 -1.0000
7 24 2.4583 1 1.0000 25 2.4000
8 24 -0.6250 1 0.5000 25 -0.5800
9 25 -0.7400 0 * 25 -0.7400
10 25 1.4000 0 * 25 1.4000
Margin 247 0.4534 2 0.7500 249 0.4558
26..................................................................................................................................
***** Regression Analysis *****
Response variate: g[5]
Fitted terms: Constant + treat + mf[6]
*** Summary of analysis ***
d.f. s.s. m.s. v.r.
Regression 10 315. 31.454 4.77
Residual 238 1569. 6.594
Total 248 1884. 7.597
Percentage variance accounted for 13.2
Standard error of observations is estimated to be 2.57
* MESSAGE: The following units have large standardized residuals:
66 -2.93
155 3.81
* MESSAGE: The following units have high leverage:
174 0.520
184 0.520
*** Estimates of regression coefficients ***
estimate s.e. t(238)
Constant 1.400 0.514 2.73
treat 2 -1.300 0.726 -1.79
treat 3 -0.020 0.726 -0.03
treat 4 -0.525 0.734 -0.72
treat 5 -2.060 0.726 -2.84
treat 6 -2.400 0.726 -3.30
treat 7 1.007 0.730 1.38
treat 8 -1.973 0.730 -2.70
treat 9 -2.140 0.726 -2.95
treat 10 0.000 0.726 0.00
mf[6] 1 -0.17 1.85 -0.09
*** Accumulated analysis of variance ***
Change d.f. s.s. m.s. v.r.
+ treat 9 314.489 34.943 5.30
+ mf[6] 1 0.053 0.053 0.01
Residual 238 1569.472 6.594
Total 248 1884.014 7.597
mf[7] 0.00 1.00 Margin
Count Mean Count Mean Count Mean
treat
1 21 2.952 4 6.250 25 3.480
2 24 2.896 1 7.500 25 3.080
3 22 4.727 3 5.667 25 4.840
4 19 4.263 5 6.300 24 4.688
5 25 2.160 0 * 25 2.160
6 25 1.620 0 * 25 1.620
7 20 3.700 4 5.875 24 4.062
8 24 3.646 0 * 24 3.646
9 20 3.175 5 6.200 25 3.780
10 18 4.083 7 4.929 25 4.320
Margin 218 3.255 29 5.862 247 3.561
26..................................................................................................................................
***** Regression Analysis *****
Response variate: g[6]
Fitted terms: Constant + treat + mf[7]
*** Summary of analysis ***
d.f. s.s. m.s. v.r.
Regression 10 348. 34.821 5.76
Residual 236 1428. 6.050
Total 246 1776. 7.220
Percentage variance accounted for 16.2
Standard error of observations is estimated to be 2.46
* MESSAGE: The following units have large standardized residuals:
99 -3.43
*** Estimates of regression coefficients ***
estimate s.e. t(236)
Constant 3.139 0.499 6.29
treat 2 -0.144 0.698 -0.21
treat 3 1.445 0.696 2.08
treat 4 1.104 0.703 1.57
treat 5 -0.979 0.700 -1.40
treat 6 -1.519 0.700 -2.17
treat 7 0.568 0.703 0.81
treat 8 0.507 0.708 0.72
treat 9 0.215 0.696 0.31
treat 10 0.584 0.698 0.84
mf[7] 1 2.132 0.510 4.18
*** Accumulated analysis of variance ***
Change d.f. s.s. m.s. v.r.
+ treat 9 242.367 26.930 4.45
+ mf[7] 1 105.848 105.848 17.49
Residual 236 1427.874 6.050
Total 246 1776.089 7.220
mf[8] 0.00 1.00 Margin
Count Mean Count Mean Count Mean
treat
1 14 4.893 7 6.571 21 5.452
2 21 2.595 3 4.667 24 2.854
3 17 4.912 5 5.100 22 4.955
4 15 4.367 4 4.875 19 4.474
5 21 3.214 4 5.125 25 3.520
6 22 3.386 3 6.167 25 3.720
7 14 4.714 7 7.714 21 5.714
8 15 5.133 9 6.111 24 5.500
9 14 5.143 6 7.417 20 5.825
10 15 4.867 3 7.667 18 5.333
Margin 168 4.179 51 6.284 219 4.669
26..................................................................................................................................
***** Regression Analysis *****
Response variate: g[7]
Fitted terms: Constant + treat + mf[8]
*** Summary of analysis ***
d.f. s.s. m.s. v.r.
Regression 10 339. 33.886 5.46
Residual 208 1292. 6.211
Total 218 1631. 7.480
Percentage variance accounted for 17.0
Standard error of observations is estimated to be 2.49
* MESSAGE: The following units have large standardized residuals:
99 -2.93
*** Estimates of regression coefficients ***
estimate s.e. t(208)
Constant 4.868 0.561 8.68
treat 2 -2.233 0.750 -2.98
treat 3 -0.312 0.762 -0.41
treat 4 -0.763 0.791 -0.97
treat 5 -1.628 0.741 -2.20
treat 6 -1.358 0.743 -1.83
treat 7 0.262 0.769 0.34
treat 8 -0.025 0.745 -0.03
treat 9 0.431 0.779 0.55
treat 10 0.173 0.803 0.22
mf[8] 1 1.753 0.408 4.30
*** Accumulated analysis of variance ***
Change d.f. s.s. m.s. v.r.
+ treat 9 224.162 24.907 4.01
+ mf[8] 1 114.701 114.701 18.47
Residual 208 1291.885 6.211
Total 218 1630.749 7.480
mfg 0.00 1.00 Margin
Count Mean Count Mean Count Mean
treat
1 9 5.056 5 4.400 14 4.821
2 4 5.125 17 2.471 21 2.976
3 7 5.429 10 4.150 17 4.676
4 6 4.500 9 2.722 15 3.433
5 1 5.000 20 1.275 21 1.452
6 7 4.429 15 2.433 22 3.068
7 5 7.900 9 2.889 14 4.679
8 9 5.833 6 3.417 15 4.867
9 3 3.667 11 3.727 14 3.714
10 5 5.700 10 5.900 15 5.833
Margin 56 5.330 112 3.022 168 3.792
26..................................................................................................................................
***** Regression Analysis *****
Response variate: g[8]
Fitted terms: Constant + treat + mfg
*** Summary of analysis ***
d.f. s.s. m.s. v.r.
Regression 10 370. 37.046 5.30
Residual 157 1097. 6.986
Total 167 1467. 8.786
Percentage variance accounted for 20.5
Standard error of observations is estimated to be 2.64
* MESSAGE: The following units have large standardized residuals:
154 2.82
159 -2.76
*** Estimates of regression coefficients ***
estimate s.e. t(157)
Constant 5.477 0.726 7.55
treat 2 -1.015 0.936 -1.09
treat 3 0.279 0.960 0.29
treat 4 -0.943 0.989 -0.95
treat 5 -2.277 0.953 -2.39
treat 6 -1.158 0.916 -1.26
treat 7 0.38 1.01 0.38
treat 8 0.124 0.982 0.13
treat 9 -0.32 1.02 -0.32
treat 10 1.580 0.993 1.59
mfg 1 -1.834 0.464 -3.95
*** Accumulated analysis of variance ***
Change d.f. s.s. m.s. v.r.
+ treat 9 261.435 29.048 4.16
+ mfg 1 109.029 109.029 15.61
Residual 157 1096.744 6.986
Total 167 1467.208 8.786
27
28 vari rate
29 restrict wg,td,rate;cond=mfg==0
30 calc rate=wg/td
31 restrict wg,td,rate;cond=mfg.ne.0
32 calc rate=g[7]/70
33 restrict wg,td,rate
34 tabulate[class=treat,mfg;marg=y;prin=m,c] rate
mfg 0.00 1.00 Margin
Count Mean Count Mean Count Mean
treat
1 17 0.07920 8 0.06786 25 0.07557
2 7 0.06942 18 0.03413 25 0.04401
3 12 0.07368 13 0.06758 25 0.07051
4 13 0.06996 12 0.06494 25 0.06766
5 5 0.06878 20 0.04286 25 0.04804
6 9 0.05578 16 0.03884 25 0.04494
7 15 0.09866 10 0.05786 25 0.08234
8 16 0.07589 9 0.05625 25 0.06935
9 13 0.09462 12 0.07381 25 0.08463
10 14 0.08263 11 0.06753 25 0.07599
Margin 121 0.07895 129 0.05422 250 0.06629
35 model rate
36 fit[prin=m,s,a,e] treat+mfg
36..................................................................................................................................
***** Regression Analysis *****
Response variate: rate
Fitted terms: Constant + treat + mfg
*** Summary of analysis ***
d.f. s.s. m.s. v.r.
Regression 10 0.0730 0.007302 5.52
Residual 237 0.3134 0.001323
Total 247 0.3864 0.001565
Percentage variance accounted for 15.5
Standard error of observations is estimated to be 0.0364
* MESSAGE: The following units have large standardized residuals:
90 3.12
*** Estimates of regression coefficients ***
estimate s.e. t(237)
Constant 0.08173 0.00744 10.99
treat 2 -0.0239 0.0105 -2.28
treat 3 -0.0012 0.0103 -0.12
treat 4 -0.0053 0.0104 -0.50
treat 5 -0.0183 0.0105 -1.73
treat 6 -0.0245 0.0104 -2.35
treat 7 0.0083 0.0103 0.81
treat 8 -0.0060 0.0104 -0.57
treat 9 0.0121 0.0103 1.18
treat 10 0.0027 0.0103 0.26
mfg 1 -0.01925 0.00485 -3.97
*** Accumulated analysis of variance ***
Change d.f. s.s. m.s. v.r.
+ treat 9 0.052201 0.005800 4.39
+ mfg 1 0.020816 0.020816 15.74
Residual 237 0.313433 0.001323
Total 247 0.386450 0.001565
37 treat treat*mfg
38 anova rate
* MESSAGE: non-orthogonality between treatment terms.
The effects (printed or used to calculate means), the efficiency
factor and the sum of squares for each treatment term are for
that term eliminating previous terms in the TREATMENT formula
and ignoring subsequent terms.
38..................................................................................................................................
***** Analysis of variance *****
Variate: rate
Source of variation d.f.(m.v.) s.s. m.s. v.r.
treat 9 0.052131 0.005792 4.31
mfg 1 0.020931 0.020931 15.58
treat.mfg 9 0.007097 0.000789 0.59
Residual 228(2) 0.306393 0.001344
Total 247(2) 0.386450
***** Information summary *****
Model term e.f. non-orthogonal terms
mfg 0.911 treat
* MESSAGE: the following units have large residuals.
*units* 90 0.1175 s.e. 0.0350
*units* 99 -0.1078 s.e. 0.0350
***** Tables of means *****
Variate: rate
Grand mean 0.0662
treat 1 2 3 4 5 6 7 8 9 10
0.0756 0.0440 0.0705 0.0675 0.0480 0.0449 0.0823 0.0688 0.0846 0.0760
mfg 0.00 1.00
0.0761 0.0570
rep. 121 129
*** Table of means for treat.mfg cannot be calculated
(contains mutually non-orthogonal components).
*** Standard errors of differences of means ***
Table treat mfg
rep. 25 unequal
s.e.d. 0.01037 0.00486
(Not adjusted for missing values)
***** Missing values *****
Variate: rate
Unit estimate
80 0.0649
184 0.0563
Max. no. iterations 2
39 stop
******** End of job. Maximum of 13860 data units used at line 26 (37924 left)