Calculating standard errors of variance estimates in ASRemlR

Dear all,

I am using ASReml-R to run a Gaussian model with a single response variable and multiple fixed and random effects, and I am having problems investigating associations between random effects.

I was under the impression that taking the square root of the diagonal elements of the inverse average information matrix obtained from model$ai , would give you the standard error of the variance component estimates as outputted in summary(model)$varcomp.

However I can only get the residual variance to give the correct answer, even when I have tried to weight by the gamma component and am unsure about how best to proceed.

Any help you could give would be most gratefully received

best,

Michelle Clements

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

Sent using Mail2Forum (http://www.mail2forum.com).

Read this topic online here:
http://www.vsni.co.uk/forum/viewtopic.php?p=627#627

-------------------- m2f --------------------
Received on Thu Jul 01 2009 - 10:26:43 EST