Hi,
I'm fitting a model with an animal effect and a permanent environment
effect. Individuals are spread over two time periods and I want to
test whether a model with time specific animal and PE effects is
better than a model where both periods are treated equally. For
univariate models this is fine:
T1 ~ mu time !r at(time,1).animal at(time,2).animal at(time,
1).identity at(time,2).identity
versus
T1 ~ mu time !r animal identity
If I fit a multivariate model either the animal or the PE effects are
flagged as singular depending on which order the terms are fitted. I
am confident that these parameters can be estimated from the data
because a) all (3) traits have the same pattern of missing data so if
the univariate model is idenitifiable the multivariate model should
be and b) I can fit a multivariate model for each time period
separately by generating two data sets and all model parameters are
estimated. One potential issue is that some individuals are present
in both time periods, but I can't see why this should cause
singularities.
Any help would be appreciated.
-- Alastair J Wilson NERC Research Fellow Alastair.Wilson_at_ed.ac.uk http://www.biology.ed.ac.uk/research/groups/awilson/ tel. 0131 6513673 Institute of Evolutionary Biology School of Biological Sciences The University of Edinburgh Ashworth Laboratories The King's Buildings West Mains Road Edinburgh EH9 3JT, UKReceived on Fri Feb 05 2007 - 09:59:43 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.)