Re: how to fit a G struture

Dear Felicity,
re
Can someone please suggest to me how to fit a G structure to the following
parental model over trials? I have:

ccs ~ mu Trial at(Trial,1).lin(Plot) at(Trial,3).lin(Nrow),
      !f mv !r at(Trial,2).Rep,
      Female and(Male) Trial.Female -Trial.Male and(Trial.Male)
4 2 1
12 Plot AR .1
41 Nrow AR .38
12 Plot AR .1
41 Nrow AR .38
12 Plot AR .12
10 Nrow AR .16
12 Plot AR .1
13 Nrow AR .38
Trial.Female 2
4 0 CORUH !+5
.5
4*.5
Female 0 ID

This runs but I would like the G structure over the combined Trial.Female
-Trial.Male and(Trial.Male).

First note that Female and(Male) assumes Female and Male are coded
conformably
so that Female_i is the same plant as Male_i for all i.

It means that the design matrix for Male is added to the design matrix to
Female
so in the output there will only be one set of effects, labelled Female.

Similarly Trial.Female -Trial.Male and(Trial.Male) results in one set of
effects labelled in the output Trial.Female

2) The above model is a little unusual.
It is not common the fit the 'Female' variance component when you put a
structure on Trail.Female

The Female term fits a common covariance accross Trials. The structure
fits a common
correlation (differing variances) . It is more usual to drop the Female
component when you fit a covariance
structure to the interaction term.

I trust this helps.

May you know Jesus Christ and His blessing in 2007,

Arthur Gilmour, His servant .

Mixed model regression mapping for QTL detection in experimental crosses.
Computational Statistics and Data Analysis 51:3749-3764 now available at
http://dx.doi.org/10.1016/j.csda.2006.12.031

Profile: http://www.dpi.nsw.gov.au/reader/17263
Personal website: http://www.cargovale.com.au/

mailto:Arthur.Gilmour_at_dpi.nsw.gov.au, arthur_at_cargovale.com.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

ASREML 2 is now available from http://www.VSNi.co.uk/products/asreml
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 (Forestry) 16-18 April Tasmania
Rec leave 19-20 April
ASReml planning at Tumut 23-27 April
Biom Brch meet 2-4 May
ASReml workshop (Animal Breeding) 8-26 May South Africa
AAABG 24-26 September Armidale
      <><><><><><><><><><><><><><><><><><><><><><><><><>

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 Sat Apr 24 2007 - 13:34:47 EST