I have 6 locations, 3 replications, 7 blocks, and 49 entries in a full diallel experiment. I am trying to find out the heritability for each trait in each location separately. So, here I have used the following R code to calculate broad sense heritability from variance components:
# Code library(lme4) library(Matrix) a<-read.csv("sadore.csv", header = T, na.strings = "NA") summary(a) head(a) str(a) a1=lmer(Yield~(1|Entry)+ (1|Rep)+(1|Rep:Block), data=a) summary(a1) #formula h2= v(Entry)/v(Entry)+v(Rep)/3+v(Rep:Block)/21+ v(Error)/63
But after that, I got heritability which is super high and not reasonable. I know I am doing something wrong here. So, is there any wrong in my formula or code?