Re: predict statement

Dear Anne,

re:

I am trying to use the predict statement to obtain predictions of the cow
(i.e. permanent environment) + animal genetic effects (I am interested in
the standard error of prediction of this sum).

The model and predict statements are as follows:

fyld ~ pfrn pjer agelin agequad,
 !r animal 0.7 cow 0.8 !f cgroup mv
predict cow 1:850 !ONLYUSE animal cow

However, the .pvs file has only zeros for the predicted value and the
associated standard errors:

 cow Predicted_Value Standard_Error Ecode
      7788125 *************** 0.0000 E
      7788133 0.0000 0.0000 E
      9171656 0.0000 0.0000 E
      9308154 0.0000 0.0000 E
      9467290 0.0000 0.0000 E
      9477086 0.0000 0.0000 E
      etc.

Does anyone have any idea as to what is going wrong? Is it not possible
to
obtain predictions of random effects?

1. Predictions of random effects are possible

2. Something has obviuosly gone wrong - depending on which version you
are using,
the problem may have been fixed in a later beta version. Else if this is
in a recent
beta version, I would like a chance to diagnose the problem properly.

3. It is unclear what this prediction will do but it is probably not what
you seem to want -
   i.e. sum of animal and cow effects. How are these things coded? I
suppose

   animal !P
   cow !I

In any case ASReml will not know how to associate the two effects.
Presumably there are
several animals per cow? You may be able to get what you want using
!PRESENT.
I would start of with a smaller set of cows to predict.

 Check out http://www.cargovale.com.au/
for my personal greetings and update.

I thank God for His mercy through Jesus Christ.

Arthur Gilmour PhD
mailto:Arthur.Gilmour_at_dpi.nsw.gov.au
Principal Research Scientist (Biometrics)
NSW Department of Primary Industries
Orange Agricultural Institute, Forest Rd, ORANGE, 2800, AUSTRALIA

fax: 02 6391 3899; 02 6391 3922 Australia +61
telephone work: 02 6391 3815; home: 02 6364 3288; mobile: 0438 251 426

http://www.dpi.nsw.gov.au/reader/17263
http://www.asreml.co.uk/
The ASReml discussion group is at ASREML-L_at_dpi.nsw.gov.au
To join it, mailto:arthur.gilmour_at_dpi.nsw.gov.au
Cookbook: http://uncronopio.org/ASReml

Proposed travel:
ASReml workshop CSIRO 16-18 May
[IBC Quebec 16-21 July 2006]
WCGALP8 13-18 Aug Brazil
      <><><><><><><><><><><><><><><><><><><><><><><><><>

This message is intended for the addressee named and may contain
confidential information. If you are not the intended recipient or
received it in error, please delete the message and notify sender. Views
expressed are those of the individual sender and are not necessarily the
views of their organisation.
Received on Mon Apr 05 2006 - 13:18:55 EST