setting variance components to be equal

From: ahu <asremlforum_at_VSNI.CO.UK>
Date: Thu, 5 Feb 2009 19:26:15 +0000

Dear all,

I would like to constrain two of my variance components in the model to optimise LogL assuming they are equal (doesn't matter what value they take). How does one code this in ASreml?

Example: Two populations and a bivariate model of the same trait.
Question: do the two pop differ signf from each other in Va?

example code:
T1 T2 ~ Fixed !r at(Trait,1).animal at(Trait,2).animal

1 1 2

at(Trait,1).animal 2
1 2 US

at(Trait,2).animal 2
1 2 US

How could I specify that the two populations should have equal Va without actually giving them a number? I.e is it possible to say that the two var comp should be the same, get the LogL and compare to a unconstrained model. If no better fit then that suggests that Va are not signf different in the two pop.


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

Sent using Mail2Forum (

Read this topic online here:

-------------------- m2f --------------------
Received on Fri Feb 05 2009 - 19:26:15 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.)