Re: estimability

Dear Juan,

E2 is not estimable because B=3 is missing for E=2

If you want E2 to be estimated by averaging over the 2 levels of B that
are present,
you need
predict t.e !present e b

but note that the B=3 effect from E=1 will then be included in the
differences between E=1 and E=2.

So, in this cae, you need to think clearly about what prediction is
needed/ what it means

May Jesus Christ be gracious to you in 2008,

Arthur Gilmour, His servant .
 
Mixed model regression mapping for QTL detection in experimental crosses.
Computational Statistics and Data Analysis 51:3749-3764 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/
Skype: arthur.gilmour
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
Archives are at
https://gatekeeper.dpi.nsw.gov.au/Listserv/archives/asreml-l.html
Cookbook: http://uncronopio.org/ASReml

Proposed travel:
UK 28 Jan-8 Feb
USA 11-15 Feb
Argentina 16-19 September Biometrics 
><><><><><><><><><><><><><><><><><><><><><><><>

Juan Burgueño <juanb_at_COLPOS.MX>
Sent by: ASReml users discussion group <ASREML-L_at_AGRIC.NSW.GOV.AU>
01/02/2008 07:07 AM
Please respond to
ASReml users discussion group <ASREML-L_at_AGRIC.NSW.GOV.AU>

To
ASREML-L_at_AGRIC.NSW.GOV.AU
cc

Subject
estimability

I have the following data:

E B T Y
2 1 1 11
2 1 2 15
2 1 3 18
2 1 4 25
2 2 1 36
2 2 2 52
2 2 3 41
2 2 4 14
1 1 1 15
1 1 2 25
1 1 3 25
1 1 4 29
1 2 1 65
1 2 2 45
1 2 3 23
1 2 4 32
1 3 1 54
1 3 2 14
1 3 3 15
1 3 4 25

the model

Prueba
  e 2
  b 3
  t 4
  y
aaa.txt !skip 1
y ~ mu e at(e).b !r t t.e
predict t.e

and this is the .pvs file

  Prueba 31 Jan
2008 14:04:29
 MYT\PruebaAsreml\aaa

  Ecode is E for Estimable, * for Not Estimable

  ---- ---- ---- ---- ---- ---- ---- ---- 1 ---- ---- ---- ---- ----
---- ---- ----
  Predicted values of y
  b is averaged over fixed levels
  The cells of the hypertable are calculated from all model terms
constructed
                 solely from factors in the averaging and classify
sets.
  Warning: 4 non-estimable [aliased] cell(s) may be omitted from
the table.
                 The Overall SED statistic includes non-estimable
predictions.

  t e Predicted_Value Standard_Error Ecode
        1.0000 1.0000 31.5336 4.4375 E
        2.0000 1.0000 30.4090 4.4375 E
        3.0000 1.0000 29.9367 4.4375 E
        4.0000 1.0000 30.4540 4.4375 E
  SED: Overall Standard Error of Difference not available

Why the blups for e 2 are non-estimable?

Thanks
Received on Mon Feb 01 2008 - 08:03:28 EST