I think the following message may have fallen into the black hole that
Arthur warned about a couple of days ago.
Nick Galwey
We (Dadan Kusnandar and myself) are trying to get means for levels of a
random term using a .pin file to combine the parameter estimates in the .sln
file. The means given are correct except for the first level of the random
factor, which is being taken to have an effect of zero, and hence is giving
the grand mean. This is despite the fact that the .sln file does NOT give
the effect of this level as zero.
As you will see from our .pin file (appended), we have tried various ways
of specifying the effects we are after - shorthand, longhand, and the
effects of each level rather than the means. In every case we get the
result we expect except for Level 1.
I also append the .pvs file and the first few lines of the .sln file.
What's going on?
The examples in the manual that we're following are for the means of levels
of a fixed-effect factor, and we suspect that this is related to the
problem.
Best wishes,
Nick Galwey
___________________________________
.pin file
# Obtain provenance means
# Provmeans 1 [2:9]
Grandmean 1
Effect1 2
Effect2 3
Provmean1 1 2
Provmean2 1 3
Provmean3 1 4
Provmean4 1 5
Provmean5 1 6
Provmean6 1 7
Provmean7 1 8
Provmean8 1 9
___________________________________
.pvs file
'Grandmean ' 1.00000 0.03076 0.15650
'Effect1 ' 1.00000 0.00000 0.00000
'Effect2 ' 1.00000 0.35460 0.78177
'Provmean1 ' 1.00000 0.03076 0.15650
'Provmean2 ' 1.00000 0.38536 0.79728
'Provmean3 ' 1.00000 -0.34114 0.76116
'Provmean4 ' 1.00000 0.17046 0.75600
'Provmean5 ' 1.00000 0.23476 0.79814
'Provmean6 ' 1.00000 -0.58524 0.76202
'Provmean7 ' 1.00000 -0.39054 0.75987
'Provmean8 ' 1.00000 0.47146 0.75558
___________________________________
.sln file (first few lines)
mu 1 0.3076E-01 0.1565
prov 1 0.2701 0.1704
prov 2 0.3546 0.1781
prov 3 -0.3719 0.1697
prov 4 0.1397 0.1685
prov 5 0.2040 0.1783
prov 6 -0.6160 0.1699
prov 7 -0.4213 0.1694
prov 8 0.4407 0.1684
rep 1 -0.2503E-01 0.4021E-01
rep 2 0.1654E-01 0.4053E-01
rep 3 0.8483E-02 0.4036E-01
rep.block 1.0001 -0.9539E-01 0.8190E-01
rep.block 1.0002 -0.1372E-01 0.8434E-01
rep.block 1.0003 0.1655E-01 0.8210E-01
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml