Exam8.1: Example 8.1 from Experimental Design and Analysis for Tree...

In eda4treeR: Experimental Design and Analysis for Tree Improvement

Example 8.1 from Experimental Design and Analysis for Tree Improvement

Description

Exam8.1 presents the 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.

Author(s)

1. Muhammad Yaseen ([e-mail])

2. Sami Ullah ([e-mail])

References

1. 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/).

See Also

DataExam8.1

Examples

library(car)
library(dae)
library(dplyr)
library(emmeans)
library(ggplot2)
library(lmerTest)
library(magrittr)
library(predictmeans)

data(DataExam8.1)

# Pg. 141
fm8.4 <-
  aov(
    formula = dbh ~ inoc + Error(repl/inoc) +
                    inoc*country*prov
  , data    = DataExam8.1
     )

# Pg. 150
summary(fm8.4)

# Pg. 150
model.tables(x = fm8.4, type = "means")

RESFit <-
    data.frame(
      fittedvalue   = fitted.aovlist(fm8.4)
    , residualvalue = proj(fm8.4)$Within[,"Residuals"]
    )

ggplot(
   data    =  RESFit
 , mapping = aes(x = fittedvalue, y = residualvalue)
 ) +
geom_point(size = 2) +
labs(
   x = "Residuals vs Fitted Values"
 , y = ""
 ) +
theme_bw()

# Pg. 153
fm8.6 <-
 aov(
   formula = terms(
                   dbh ~ inoc + repl + col +
                         repl:row + repl:col +
                         prov + inoc:prov
                   , keep.order = TRUE
                   )
 , data   = DataExam8.1
 )
summary(fm8.6)

eda4treeR documentation built on Sept. 14, 2024, 1:08 a.m.