comparing heritabilities of traits

Hi Ardo,

I think I have addressed some of these questions in reply to some of  
your previous posts to the list. In short, testing the 'significance'  
of a heritability estimate is difficult to evaluate from the posterior  
because it has to be positive. I've seen people use DIC to evaluate  
whether a genetic term should be included, but I think DIC is  
'focused' at the wrong level for this sort of question and people  
should be very wary. My personal view is that a significance test for  
h2 is meaningless (everything that varies between individuals in an  
outbred population will have some genetic variance, even if it is  
small) and that people should focus on the width of the credible/ 
confidence intervals rather than the p-values.

With regards to your second question about how to choose the trait  
with greatest h2, you could try this:

h2<-function(x){x["animal"]/sum(x)}  # calculates h2

h2.post<-as.data.frame(lapply(list(m1=m1$VCV, m2=m2$VCV, m3=m3$VCV),  
function(x){apply(x, 1,h2)}))

# create a table of posterior h2 for each response (m1, m2, m3 .....)

prop.table(table(apply(h2.post, 1, which.max)))

# posterior probability for each response that it has the highest h2

Cheers,

Jarrod

comparing heritabilities of traits

Thread (2 messages)