Thanks to Drs. Gilmour and Luis Apiazola by the info about true BLUP. In
fact, before my message to the mail list, I tried maxit=1, but I could not
get the results from MTDFREML. Then I decided to try other things. But
now, with the maxit=1 plus the option !S2==, the results are similar to
those from MTFREML.