Exam8.1.2 | R Documentation |
Exam8.1.2 presents the Analysis of Nested Seedlot Structure 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. 167
fm8.11 <-
aov(
formula = dbh ~ country + country:prov
, data = DataExam8.1
)
b <- anova(fm8.11)
Res <- length(b[["Sum Sq"]])
df <- 119
MSS <- 0.1951
b[["Df"]][Res] <- df
b[["Sum Sq"]][Res] <- MSS*df
b[["Mean Sq"]][Res] <- b[["Sum Sq"]][Res]/b[["Df"]][Res]
b[["F value"]][1:Res-1] <-
b[["Mean Sq"]][1:Res-1]/b[["Mean Sq"]][Res]
b[["Pr(>F)"]][Res-1] <-
df(
b[["F value"]][Res-1]
, b[["Df"]][Res-1]
, b[["Df"]][Res]
)
b
emmeans(fm8.11, specs = "country")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.