Hi Ian,
The pol structure was interacted with units, so I think Z should be a
diagonal nxn matrix and ZZ' full rank.
Cheers,
Jarrod
On 20 Oct 2009, at 15:05, I.M.S.White wrote:
> Jarrod,
>
> With the random regression model you get the correct variances but
> dont you also get covariances between the observations? Your Z matrix
> is an nx1 column and ZZ' is nxn of rank 1.
>
> Quoting Jarrod Hadfield <j.hadfield_at_ED.AC.UK>:
>
>> Hi Brian,
>>
>> Thanks - this is a much simpler way of fitting the model and the
>> intercept now agrees with the other methods. I thought the model I
>> fitted would be equivalent, but apparently not.
>>
>> Cheers,
>>
>> Jarrod
>>
>>
>>
>> On 20 Oct 2009, at 12:24, brian.cullis_at_INDUSTRY.NSW.GOV.AU wrote:
>>
>>>
>>> Dear jarrod
>>> we do this all the time in the two stage analysis of variety
>>> trials where we have weights in the analysis, see the paper by
>>> Smith Cullis and Thompson (2005) J Ag Sci Camb
>>>
>>> Anyway the syntax would be
>>>
>>> asreml
>>> (y
>>> ~1,random=~units,weight=wt,family=asreml.gaussian(dispersion=1),...)
>>> where wt contains the vector of 1/sqrt(sd_i)
>>>
>>> I think this should do the trick@
>>> let me know how you get on
>>>
>>> warm regards
>>>
>>> Brian Cullis|Research Leader, Biometrics &
>>> Senior Principal Research Scientist |
>>> Industry & Investment NSW | Wagga Wagga Agricultural Institute |
>>> Pine Gully Road | Wagga Wagga NSW 2650 |
>>> PMB | Wagga Wagga NSW 2650
>>> T: 02 6938 1855 | M: 0439 448 591 | F: 02 6938 1809 | E: brian.cullis_at_industry.nsw.gov
>>> ,au
>>> W: www.industry.nsw.gov.au |
>>>
>>> Visiting Professorial Fellow
>>> School of Mathematics and Applied Statistics
>>> Faculty of Informatics
>>> University of Wollongong
>>>
>>> Professor,
>>> Faculty of Agriculture, Food & Natural Resources
>>> The University of Sydney
>>>
>>> Adjunct Professor
>>> School of Computing and Mathematics
>>> Charles Sturt University
>>>
>>>
>>>
>>>
>>>
>>>
>>> This message is intended for the addressee named and may contain
>>> confidential information. If you are not the intended recipient,
>>> please delete it and notify the sender. Views expressed in this
>>> message are those of the individual sender, and are not
>>> necessarily the views of their organisation.
>>
>>
>> This message is intended for the addressee named and may contain
>> confidential information. If you are not the intended recipient,
>> please
>> delete it and notify the sender. Views expressed in this message are
>> those of the individual sender, and are not necessarily the views of
>> their organisation.
>
>
>
> ::::::::::::::::::::::::::::::::::::::::::
> : I.White :
> : University of Edinburgh :
> : Ashworth Laboratories, West Mains Road :
> : Edinburgh EH9 3JT :
> : Tel 0131 650 5490 Fax 0131 650 6564 :
> ::::::::::::::::::::::::::::::::::::::::::
>
> --
> The University of Edinburgh is a charitable body, registered in
> Scotland, with registration number SC005336.
>
> This message is intended for the addressee named and may contain
> confidential information. If you are not the intended recipient,
> please delete it and notify the sender. Views expressed in this
> message are those of the individual sender, and are not necessarily
> the views of their organisation.
>
>
-- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.Received on Wed Oct 20 2009 - 15:18:05 EST
This webpage is part of the ASReml-l discussion list archives 2004-2010. More information on ASReml can be found at the VSN website. This discussion list is now deprecated - please use the VSN forum for discussion on ASReml. (These online archives were generated using the hypermail package.)