In the files test.dat, test.as and test.pin, embedded below, I believe I
have reproduced the example given on p.88 of the ASREML manual (version of 2
October 1998). However, the .pin file does not produce
65 = 5 + 12 + 12 + 12 + 12 + 12
equations as stated, but
85 = 25 + 12 + 12 + 12 + 12 + 12
This appears to be because the loops implied by the two sets of square
brackets are not processed in parallel, but nested. This would be quite a
useful feature, say for estimating treatment means of factorial experiments
(though one also needs a facility for parallel processing). However, the
labelling of the first 25 predicted values is not quite what I would expect
if the loops are nested. I would expect the labels to be
'GROUP 1 6 '
'GROUP 1 7 '
'GROUP 1 8 '
'GROUP 1 9 '
'GROUP 1 10 '
.
.
.
'GROUP 5 6 '
'GROUP 5 7 '
'GROUP 5 8 '
'GROUP 5 9 '
'GROUP 5 10 '
but they are
'GROUP 0 7 '
'GROUP 0 8 '
'GROUP 0 9 '
'GROUP 0 10 '
'GROUP 1 0 '
.
.
.
'GROUP 4 9 '
'GROUP 4 10 '
'GROUP 5 0 '
What's going on?
Nick Galwey
___________________________
test.as
Example - Manual page 88
group 5
X
y
C:\DOCS\Dracup_M\96WH14\test.dat !skip 1
y ~ mu group group.X !r spl(X) group.spl(X)
test.pin
# A is the vector of the 12 knot points
# B is the vector of the average predicted spline
C = C+B # Add average curve to Group 1 specific curve
D = D+B # Add average curve to Group 2 specific curve
E = E+B # Add average curve to Group 3 specific curve
F = F+B # Add average curve to Group 4 specific curve
G = G+B # Add average curve to Group 5 specific curve
GROUP B6 [ 1:5 ] A6 [6:10 ] 11 # Group means at knot point 6
GROUP1 C 1 A 6 11 # Spline for group 1
GROUP2 D 2 A 7 11 # Spline for group 2
GROUP3 E 3 A 8 11 # Spline for group 3
GROUP4 F 4 A 9 11 # Spline for group 4
GROUP5 G 5 A 10 11 # Spline for group 5
test.dat
group X y
1 12 5.11
1 13 5.45
1 15 6.08
1 16 6.32
1 18 6.97
1 19 7.29
1 21 7.87
1 24 8.78
1 25 9.08
1 26 9.40
1 27 9.61
1 29 10.24
2 12 10.24
2 13 10.85
2 15 12.03
2 16 12.66
2 18 13.90
2 19 14.50
2 21 15.61
2 24 17.44
2 25 18.07
2 26 18.65
2 27 19.29
2 29 20.41
3 12 15.39
3 13 16.26
3 15 18.09
3 16 18.95
3 18 20.76
3 19 21.62
3 21 23.47
3 24 26.19
3 25 27.04
3 26 27.92
3 27 28.89
3 29 30.68
4 12 20.42
4 13 21.70
4 15 24.09
4 16 25.21
4 18 27.68
4 19 28.86
4 21 31.30
4 24 34.86
4 25 36.01
4 26 37.21
4 27 38.48
4 29 40.89
5 12 25.53
5 13 27.09
5 15 30.09
5 16 31.60
5 18 34.56
5 19 36.00
5 21 39.02
5 24 43.60
5 25 45.05
5 26 46.56
5 27 48.01
5 29 51.07
_____________________________________________________________________
N.W. Galwey,
Faculty of Agriculture,
University of Western Australia,
Nedlands, WA 6709, Australia.
Tel.: +61 9 380 1959 (direct line)
+61 9 380 2554 (switchboard)
Fax: +61 9 380 1108