Packages

library(tidyverse)
library(agricolae)

Data

D.C Arny di Universitas Wisconsin membandingkan 4 varietas gandum (faktor A) yaitu Vicland 1, Vicland 2, Clinton, dan Branch yang ditempatkan sebagai petak utama, dan perlakuan kimiawi terhadap benih (faktor B ) yaitu kontrol, Ceresan M, Panogen, Agrox yang ditempatkan pada anak petak dengan kelompok sebanyak 4.

padi <- readxl::read_excel("split plot padi.xlsx")
head(padi,10)
## # A tibble: 10 x 4
##    proses_kimiawi varietas kelompok hasil_panen
##    <chr>          <chr>       <dbl>       <dbl>
##  1 Kontrol        Vicland1        1        42.9
##  2 Kontrol        Vicland1        2        41.6
##  3 Kontrol        Vicland1        3        28.9
##  4 Kontrol        Vicland1        4        30.8
##  5 ceresanM       Vicland1        1        53.8
##  6 ceresanM       Vicland1        2        58.5
##  7 ceresanM       Vicland1        3        43.9
##  8 ceresanM       Vicland1        4        46.3
##  9 panogen        Vicland1        1        49.5
## 10 panogen        Vicland1        2        53.8

mengubah setiap faKtor dan ulangan ke dalam bentuk factor

padi <- padi %>% mutate_if(is.character,as.factor) %>% mutate(kelompok=as.factor(kelompok))

ANOVA

anova_padi <- with(padi,sp.plot(block = kelompok,
                                  pplot = varietas,
                                  splot = proses_kimiawi,
                                  Y = hasil_panen))
## 
## ANALYSIS SPLIT PLOT:  hasil_panen 
## Class level information
## 
## varietas     :  Vicland1 Vicland2 Clinton Branch 
## proses_kimiawi   :  Kontrol ceresanM panogen agrox 
## kelompok     :  1 2 3 4 
## 
## Number of observations:  64 
## 
## Analysis of Variance Table
## 
## Response: hasil_panen
##                         Df  Sum Sq Mean Sq F value   Pr(>F)   
## kelompok                 3 2842.87  947.62 13.7938 0.001029 **
## varietas                 3 2848.02  949.34 13.8188 0.001022 **
## Ea                       9  618.29   68.70                    
## proses_kimiawi           3  170.54   56.85  2.7987 0.053859 . 
## varietas:proses_kimiawi  9  586.47   65.16  3.2082 0.005945 **
## Eb                      36  731.20   20.31                    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## cv(a) = 15.7 %, cv(b) = 8.5 %, Mean = 52.80938