Packages

library(tidyverse)
library(agricolae)

Data

Sebuah penelitian membandingkan hasil gabah enam Varietas Padi (IR8(V1), IR127-80(V2), IR305-4-12(V3), IR400-2-5(V4), IR665-58(V5) dan Peta(V6)), ditanam sebar langsung dengan tiga taraf Nitrogen (0(N1), 60(N2) dan 120(N3)) dalam rancangan petak berjalur dengan tiga ulangan. Berikut adalah datanya

padi <- readxl::read_excel("strip plot padi.xlsx")
head(padi,10)
## # A tibble: 10 x 4
##    varietas nitrogen ulangan hasil_gabah
##    <chr>    <chr>      <dbl>       <dbl>
##  1 V1       N1             1        2373
##  2 V1       N2             1        4076
##  3 V1       N3             1        7254
##  4 V1       N1             2        3958
##  5 V1       N2             2        6431
##  6 V1       N3             2        6808
##  7 V1       N1             3        4384
##  8 V1       N2             3        4889
##  9 V1       N3             3        8582
## 10 V2       N1             1        4007

mengubah setiap faKtor dan ulangan ke dalam bentuk factor

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

ANOVA

anova_padi <- with(padi,strip.plot(BLOCK = ulangan,
                     COL = varietas,ROW = nitrogen,
                     Y = hasil_gabah
                     )
     )
## 
## ANALYSIS STRIP PLOT:  hasil_gabah 
## Class level information
## 
## varietas     :  V1 V2 V3 V4 V5 V6 
## nitrogen     :  N1 N2 N3 
## ulangan  :  1 2 3 
## 
## Number of observations:  54 
## 
## model Y: hasil_gabah ~ ulangan + varietas + Ea + nitrogen + Eb + nitrogen:varietas + Ec 
## 
## Analysis of Variance Table
## 
## Response: hasil_gabah
##                   Df   Sum Sq  Mean Sq F value    Pr(>F)    
## ulangan            2  9220962  4610481 11.2001 0.0005453 ***
## varietas           5 57100201 11420040  7.6528 0.0033722 ** 
## Ea                10 14922619  1492262  3.6251 0.0068604 ** 
## nitrogen           2 50676061 25338031 34.0690 0.0030746 ** 
## Eb                 4  2974908   743727  1.8067 0.1671590    
## nitrogen:varietas 10 23877979  2387798  5.8006 0.0004271 ***
## Ec                20  8232917   411646                      
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## cv(a) = 23.1 %, cv(b) = 16.3 %, cv(c) = 12.1 %, Mean = 5289.944