Means of levels of a random term
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Means of levels of a random term



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