Exam8.1.2: Example 8.1.2 from Experimental Design & Analysis for Tree...

Exam8.1.2R Documentation

Example 8.1.2 from Experimental Design & Analysis for Tree Improvement

Description

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.

Author(s)

  1. Muhammad Yaseen (myaseen208@gmail.com)

  2. Sami Ullah (samiullahuos@gmail.com)

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. 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")

MYaseen208/eda4treeR documentation built on Sept. 15, 2024, 10:56 a.m.