# Matrix of pair-wise SED values, squared vdiff <- p3$sed^ 2 # Average variance of two DIFFERENT means (using lower triangular of vdiff) vblup <- mean(vdiff) #. 142902 # Average variance of a difference of two adjusted means (BLUP) p3 <- predict(m3, data=dat, classify= "gen", sed= TRUE ) # Heritability calculation of Piepho & Mohring, Example 1 m3 <- asreml(yield ~ 1 + rep, data=dat, random=~ gen + rep :block) Splom(means, main= "john.alpha - means for RCB, IB, Intra-block" ) Slight, constant # difference for each method as compared to John and Williams. M2 <- lmer(yield ~ 0 + gen + rep + ( 1 | rep :block), dat) # Block random (combined inter-intra block analysis) libs(lme4, lucid) # Block fixed (intra-block analysis) (bottom of table 7.4 in John) m1 <- lm(yield ~ 0 + gen + rep + rep :block, dat) # RCB (no incomplete block) m0 <- lm(yield ~ 0 + gen + rep, data=dat)
0 Comments
Leave a Reply. |
Details
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |