Re: Multivariate analysis using a giv file

Dear Yutao Li,

The first thing I note is that 'animal' is coded 1:1097 (576 unique values) in the data file but is just declared as a covariate (not a factor).

What ASReml expects when it sees giv(animal,1) is that animal is a factor
and that the GIV matrix matches it.

So, declare animal as
 animal !I 576

The identities in the data file appear to be in order.
ASReml assumes the GIV file is coded in the same order.

Now parts 1 and 2 should run.

Part 3 needs revision:
!part 3
weight ~ mu environ environ.site !r environ.animal

#There are two separate errors, with one dimension
2 1 1
#error environ 1
480 0 IDEN !S2=80
96 0 IDEN !S2=65
environ.animal 2 # Remove # from start of line
environ 0 CORR 0.65 15 20 # Prefer CORUH rather than CORR
animal 0 GIV1 # Need to specify the variance structure here

Hopefully that will sort it.

Arthur

------------------------
Arthur Gilmour

Principal Research Scientist (Biometrics)
NSW Dept Primary Industries

-------------------- m2f --------------------

Sent using Mail2Forum (http://www.mail2forum.com).

Read this topic online here:
http://www.vsni.co.uk/forum/viewtopic.php?p=351#351

-------------------- m2f --------------------
Received on Sat Jan 16 2009 - 03:04:21 EST