``
df <- read.csv("C:/Users/ma.rohim.EWSI/OneDrive - EWINDO/OneDrive Ainur/Scholarship/KULIAH/Statistika/splitplot.csv", sep = ";")
str(df)## 'data.frame': 60 obs. of 4 variables:
## $ Rep : int 1 1 1 1 1 2 2 2 2 2 ...
## $ Var : chr "V1 " "V1 " "V1 " "V1 " ...
## $ Pupuk: chr "P0" "P1" "P2" "P3" ...
## $ Yield: num 22.8 29.6 36.4 25.3 29 ...df## Rep Var Pupuk Yield
## 1 1 V1 P0 22.78
## 2 1 V1 P1 29.63
## 3 1 V1 P2 36.39
## 4 1 V1 P3 25.34
## 5 1 V1 P4 29.01
## 6 2 V1 P0 20.26
## 7 2 V1 P1 46.66
## 8 2 V1 P2 66.21
## 9 2 V1 P3 43.38
## 10 2 V1 P4 24.74
## 11 3 V1 P0 8.46
## 12 3 V1 P1 18.52
## 13 3 V1 P2 31.94
## 14 3 V1 P3 14.69
## 15 3 V1 P4 8.89
## 16 1 V2 P0 22.14
## 17 1 V2 P1 36.12
## 18 1 V2 P2 30.54
## 19 1 V2 P3 28.43
## 20 1 V2 P4 32.61
## 21 2 V2 P0 27.37
## 22 2 V2 P1 50.22
## 23 2 V2 P2 51.72
## 24 2 V2 P3 40.60
## 25 2 V2 P4 43.10
## 26 3 V2 P0 16.73
## 27 3 V2 P1 40.36
## 28 3 V2 P2 47.24
## 29 3 V2 P3 29.72
## 30 3 V2 P4 22.24
## 31 1 V3 P0 36.16
## 32 1 V3 P1 50.39
## 33 1 V3 P2 40.58
## 34 1 V3 P3 42.53
## 35 1 V3 P4 36.86
## 36 2 V3 P0 24.19
## 37 2 V3 P1 38.44
## 38 2 V3 P2 34.61
## 39 2 V3 P3 47.91
## 40 2 V3 P4 29.84
## 41 3 V3 P0 14.69
## 42 3 V3 P1 31.36
## 43 3 V3 P2 29.68
## 44 3 V3 P3 31.75
## 45 3 V3 P4 15.11
## 46 1 V4 P0 32.50
## 47 1 V4 P1 37.94
## 48 1 V4 P2 44.24
## 49 1 V4 P3 30.97
## 50 1 V4 P4 33.24
## 51 2 V4 P0 40.88
## 52 2 V4 P1 46.11
## 53 2 V4 P2 56.64
## 54 2 V4 P3 59.98
## 55 2 V4 P4 55.17
## 56 3 V4 P0 9.18
## 57 3 V4 P1 32.36
## 58 3 V4 P2 37.31
## 59 3 V4 P3 25.59
## 60 3 V4 P4 15.99library(doebioresearch)## Warning: package 'doebioresearch' was built under R version 4.3.3splitplot(df[4], df$Rep, df$Var, df$Pupuk, 1)## $Yield
## $Yield[[1]]
## $Yield[[1]][[1]]
## Analysis of Variance Table
##
## Response: dependent.var
## Df Sum Sq Mean Sq F value Pr(>F)
## block 2 3359.0 1679.49 8.5182 0.01767 *
## main.plot 3 605.7 201.90 1.0240 0.44585
## Ea 6 1183.0 197.17
## sub.plot 4 2803.4 700.84 19.0921 4.102e-08 ***
## main.plot:sub.plot 12 647.9 53.99 1.4709 0.18656
## Eb 32 1174.7 36.71
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## $Yield[[1]][[2]]
## [1] "CV(a): 41.952 , CV(b) : 18.102"
##
## $Yield[[1]][[3]]
## [1] "R Square 0.88"
##
## $Yield[[1]][[4]]
##
## Shapiro-Wilk normality test
##
## data: model$residuals
## W = 0.9784, p-value = 0.3647
##
##
## $Yield[[1]][[5]]
## [1] "Normality assumption is not violated"
##
## $Yield[[1]][[6]]
## [1] "All the main plot factor level means are same so dont go for any multiple comparison test"
##
## $Yield[[1]][[7]]
## $Yield[[1]][[7]][[1]]
## MSerror Df Mean CV t.value LSD
## 197.1655 6 33.47067 41.95185 2.446912 12.54594
##
## $Yield[[1]][[7]][[2]]
## dependent.var groups
## V4 37.20667 a
## V2 34.60933 a
## V3 33.60667 a
## V1 28.46000 a
##
##
## $Yield[[1]][[8]]
## [1] "The means of one or more levels of sub plot factor are not same, so go for multiple comparison test"
##
## $Yield[[1]][[9]]
## $Yield[[1]][[9]][[1]]
## MSerror Df Mean CV t.value LSD
## 36.70847 32 33.47067 18.10167 2.036933 5.038303
##
## $Yield[[1]][[9]][[2]]
## dependent.var groups
## P2 42.25833 a
## P1 38.17583 ab
## P3 35.07417 b
## P4 28.90000 c
## P0 22.94500 d
##
##
## $Yield[[1]][[10]]
## [1] "All the interaction level means are same so dont go for any multiple comparison test"
##
## $Yield[[1]][[11]]
## $Yield[[1]][[11]][[1]]
## MSerror Df Mean CV t.value LSD
## 36.70847 32 33.47067 18.10167 2.036933 10.07661
##
## $Yield[[1]][[11]][[2]]
## dependent.var groups
## V4:P2 46.06333 a
## V1 :P2 44.84667 ab
## V2:P2 43.16667 ab
## V2:P1 42.23333 abc
## V3:P3 40.73000 abcd
## V3:P1 40.06333 abcd
## V4:P3 38.84667 abcd
## V4:P1 38.80333 abcd
## V3:P2 34.95667 bcde
## V4:P4 34.80000 bcde
## V2:P3 32.91667 cde
## V2:P4 32.65000 cde
## V1 :P1 31.60333 def
## V1 :P3 27.80333 efg
## V4:P0 27.52000 efg
## V3:P4 27.27000 efg
## V3:P0 25.01333 efgh
## V2:P0 22.08000 fgh
## V1 :P4 20.88000 gh
## V1 :P0 17.16667 h