Re: PREDVALS?
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: PREDVALS?



Hello Vince,

> 
> I have simulated a field experiment where I know the 'true' genoeffects
> value. Now I have run this simulated data-set through ASREML and wish to
> correlate the 'true' genoeffect vrs the corresponding Blup estimates. 
> 
> I have simulated an unreplicated trial, and fitted the model;
> 
> Yield ~ mu c(check) ! geno  Ar(1).Ar(1)
> 

I presume your test genotypes are nested within a 'Testline' level
of Check and are coded 0 for plots which receive a proper Check variety.
> 
> Referring to the .sln file would the 
> 
> 	corr( 'true sim genoeffects' vrs BU) 
> 
> be equal to the
> 
> 	corr( 'true sim geno effects'  vrs PREDVALS)
> 
> How exactly do you work out PREDVALS. Is it just
>  
>       (mu +BLUP geno)?

These two correaltions would be the same provided the test lines [geno]
all related to the same level of Check [C_t], there was stricly only one 
replicate
and that predicted values for the Check plots were ignored.

   For the test lines, the PREDVAL would be mu + C_t + u_i
      
> 

If I remember correctly,  Var [ u ]        [  G   G   ]
                              [uhat ] =    [  G   G+E ]
                              
      so the expected regression of uhat on u is 1.                        


> Many thanks for your assistance.
> 
> Kind regards 
> 
> Vince	
> 
> 
> 
> Vince Matassa 
> Department Of Agriculture
> BIOMETRICS SECTION 
> University Of Queensland 
> Brisbane 4072 
> Australia
> 
> 
> 


<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Arthur Gilmour PhD                    email: Arthur.Gilmour@agric.nsw.gov.au
Senior Research Scientist (Biometrics)                 fax: <61> 2 6391 3899
NSW Agriculture                                             <61> 2 6391 3922
Orange Agricultural Institute               telephone work: <61> 2 6391 3815
Forest Rd, ORANGE, 2800, AUSTRALIA                    home: <61> 2 6362 0046

ASREML is currently free by anonymous ftp from pub/aar on ftp.res.bbsrc.ac.uk
    Point your web browser at ftp://ftp.res.bbsrc.ac.uk/pub/aar/ 
    in the IACR-Rothamsted information system  http://www.res.bbsrc.ac.uk/ 

To join the asreml discussion list, send the message  
     subscribe
to  asreml-request@chiswick.anprod.CSIRO.au

The address for messages to the list is asreml@chiswick.anprod.CSIRO.au

                        <> <> <> <> <> <> <>
"Seek first the kingdom of God and His righteousness"  Jesus; Matthew 6: 33
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>