Dear Vivi,
There is no simple way.
Assuming the term in question is fitted in the first (fitted dense)
part of the model, you could use the !VRB qualifier to extract the
variance matrix for the fixed effects. Then the variance of the missing
effect is the variance of the sum of the terms actually fitted. E.g.
It there is a total of 4 levels (A B C D), you will be able to extract
the variance matrix for (A-mu, B-mu, C-mu) a 3x3 matrix symbolically
represented by
AA
AB BB
AC BC CC
The variance of D-mu is AA+BB+CC+2AB+2AC+2BC
Alternatively, you could use a predict statement like
PREDICT !AVE c(FACTOR) -1 -1 -1 !ONLY c(FACTOR) !PRINTALL
or some variation of that.
!PRINTALL is required because this is not 'an estimable function'
Many statisticians dislike looking at such effects and consider it is
more sensible to look at estimable contrasts such as D vs
average(A,B,C)
On Thu, 2009-04-23 at 20:32 -0700, vivi noviati wrote:
> Dear All,
>
> Please excuse my question, but I am wondering how to get SE for the
> last term in fixed effect when it was constrain to sum to zero using
> c(f).
>
> I am apologies if similar question had been asked before.
>
> Thank you.
>
> Regards,
> Vivi
>
>
>
>
>
>
> 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.
-- May the God and Father of the Lord Jesus Christ guide and bless you. Arthur Gilmour PS Are you anxious? http://www.aish.com/spirituality/growth/God_vs._Prozac.asp Mobile Number +61 427 227 468 Home phone +61 2 6364 3288 Skype: Arthur.Gilmour http://www.CargoVale.com.au Wisdom is justified by all her children Luke 7:35 Travel: Tas 11-16 May Brazil Jul27 - Aug 8 UK Aug 10-14Received on Sat Apr 24 2009 - 14:56:57 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.)