Exam8.1.1 | R Documentation |
Exam8.1.1 presents the Mixed Effects Analysis of Diameter at breast height (Dbh) of 60 SeedLots under layout of row column design with 6 rows and 10 columns in 18 countries and 59 provinces of 18 selected countries given in Example 8.1.
Muhammad Yaseen (myaseen208@gmail.com)
Sami Ullah (samiullahuos@gmail.com)
E.R. Williams, C.E. Harwood and A.C. Matheson (2023). Experimental Design and Analysis for Tree Improvement. CSIRO Publishing (https://www.publish.csiro.au/book/3145/).
DataExam8.1
library(car)
library(dae)
library(dplyr)
library(emmeans)
library(ggplot2)
library(lmerTest)
library(magrittr)
library(predictmeans)
data(DataExam8.1)
# Pg. 155
fm8.8 <-
lmerTest::lmer(
formula = dbh ~ 1 + repl + col + prov +
(1|repl:row) + (1|repl:col)
, data = DataExam8.1
, REML = TRUE
)
# Pg. 157
## Not run:
varcomp(fm8.8)
## End(Not run)
anova(fm8.8)
anova(fm8.8, ddf = "Kenward-Roger")
predictmeans(model = fm8.8, modelterm = "repl")
predictmeans(model = fm8.8, modelterm = "col")
predictmeans(model = fm8.8, modelterm = "prov")
# Pg. 161
RCB1 <-
aov(dbh ~ prov + repl, data = DataExam8.1)
RCB <-
emmeans(RCB1, specs = "prov") %>%
as_tibble()
Mixed <-
emmeans(fm8.8, specs = "prov") %>%
as_tibble()
table8.9 <-
left_join(
x = RCB
, y = Mixed
, by = "prov"
, suffix = c(".RCBD", ".Mixed")
)
print(table8.9)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.