true BLUP

[Uilson V Lopes <wilsonvl@grove.ufl.edu>]

Dear asremlrs,
Sorry to return to the TRUE BLUP problem.

I am trying to run true BLUP in a simulated data set. 

I tried the option !maxit, but the final variance components are different
from my priors (Analysis 1, below).

Then I tried to add the option !GF, which was supposed to keep fixed the
variance components. Again, the estimates of the variance components are
different from my priors (Analysis 2, below).

Could somebody help me solve this?
Thanks
Uilson Lopes



====================================
.as file

 True BLUP
 env            2  #  environments
 fam           80  # families
 blk           12  # blocks
 x1
/f/share/uvl/CHAP3/BLUPtrue/CH3e.dat  !maxit 1
 x1  ~ mu blk env !r fam 0.0013468 fam.env 0.052304 blk.fam 0.11188  # !GF
 1 1 0
 0 !S2==0.012249  # !GF

======================================================
 ANALYSIS 1. NO !GF option


 LogL=-285885.     S2= 0.1225E-01  9588 df   0.00135   0.05230   0.11188
     0.01225
 Final parameter values                      0.00549   0.06940   0.14017
     0.01225

   Source       Model  terms     Gamma     Component    Compnt/StndErr
  fam              80     80  0.549387E-02  0.549387E-02      0.84
  fam.env         160    160  0.694019E-01  0.694019E-01      7.33
  blk.fam         960    960  0.140167      0.140167         28.44
  Variance       9600   9588  0.122490E-01  0.122490E-01      0.00
=====================================================
ANALYSIS 2. WITH !GF option
 
 LogL=-285885.     S2= 0.1225E-01  9588 df   0.00135   0.05230   0.11188
     0.01225
 Final parameter values                      0.00549   0.07412   0.11188
     0.01225

   Source       Model  terms     Gamma     Component    Compnt/StndErr
  fam              80     80  0.549387E-02  0.549387E-02      0.84
  fam.env         160    160  0.741164E-01  0.741164E-01      7.86
  blk.fam         960    960  0.111880      0.111880          0.00
  Variance       9600   9588  0.122490E-01  0.122490E-01      0.00

============================================================