Split-plot Analysis for the No. of Callus Formation
## Loading required package: carData
Prior to the conduct of the Analysis of Variance for the callus formation, these model diagnostics were shown:
shapiro.test(model$residuals)
##
## Shapiro-Wilk normality test
##
## data: model$residuals
## W = 0.9526, p-value = 0.3081
leveneTest(Callus_Formation ~ Rep, data = callus)
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
leveneTest(Callus_Formation ~ PM, data = callus)
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
leveneTest(Callus_Formation ~ IBA, data = callus)
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
Analysis of Variance
splitplot(callus[4], callus$Rep, callus$PM, callus$IBA, 3)
## $Callus_Formation
## $Callus_Formation[[1]]
## $Callus_Formation[[1]][[1]]
## Analysis of Variance Table
##
## Response: dependent.var
## Df Sum Sq Mean Sq F value Pr(>F)
## block 2 7.750 3.8750 0.6992 0.5885
## main.plot 1 12.042 12.0417 2.1729 0.2784
## Ea 2 11.083 5.5417
## sub.plot 3 4.125 1.3750 0.5410 0.6633
## main.plot:sub.plot 3 1.125 0.3750 0.1475 0.9293
## Eb 12 30.500 2.5417
##
## $Callus_Formation[[1]][[2]]
## [1] "CV(a): NA , CV(b) : NA"
##
## $Callus_Formation[[1]][[3]]
## [1] "R Square 0.542"
##
## $Callus_Formation[[1]][[4]]
##
## Shapiro-Wilk normality test
##
## data: model$residuals
## W = 0.9526, p-value = 0.3081
##
##
## $Callus_Formation[[1]][[5]]
## [1] "Normality assumption is not violated"
##
## $Callus_Formation[[1]][[6]]
## [1] "All the main plot factor level means are same so dont go for any multiple comparison test"
##
## $Callus_Formation[[1]][[7]]
## $Callus_Formation[[1]][[7]][[1]]
## MSerror Df Mean CV MSD
## 5.541667 2 0.875 269.0371 4.131495
##
## $Callus_Formation[[1]][[7]][[2]]
## test name.t ntr StudentizedRange alpha
## Tukey main.plot 2 6.079637 0.05
##
## $Callus_Formation[[1]][[7]][[3]]
## dependent.var groups
## M1 1.5833333 a
## M2 0.1666667 a
##
##
## $Callus_Formation[[1]][[8]]
## [1] "All the sub plot factor factor level means are same so dont go for any multiple comparison test"
##
## $Callus_Formation[[1]][[9]]
## $Callus_Formation[[1]][[9]][[1]]
## MSerror Df Mean CV MSD
## 2.541667 12 0.875 182.2012 2.732715
##
## $Callus_Formation[[1]][[9]][[2]]
## test name.t ntr StudentizedRange alpha
## Tukey sub.plot 4 4.19866 0.05
##
## $Callus_Formation[[1]][[9]][[3]]
## dependent.var groups
## T1 1.5 a
## T3 1.0 a
## T0 0.5 a
## T2 0.5 a
##
##
## $Callus_Formation[[1]][[10]]
## [1] "All the interaction level means are same so dont go for any multiple comparison test"
##
## $Callus_Formation[[1]][[11]]
## $Callus_Formation[[1]][[11]][[1]]
## MSerror Df Mean CV MSD
## 2.541667 12 0.875 182.2012 4.711453
##
## $Callus_Formation[[1]][[11]][[2]]
## test name.t ntr StudentizedRange alpha
## Tukey main.plot:sub.plot 8 5.118658 0.05
##
## $Callus_Formation[[1]][[11]][[3]]
## dependent.var groups
## M1:T1 2.3333333 a
## M1:T3 2.0000000 a
## M1:T0 1.0000000 a
## M1:T2 1.0000000 a
## M2:T1 0.6666667 a
## M2:T0 0.0000000 a
## M2:T2 0.0000000 a
## M2:T3 0.0000000 a
Split-plot Analysis for the Avg. Days to Shoot
Prior to the conduct of the Analysis of Variance for the average days to shoot of canistel stem cuttings, these model diagnostics were shown:
shapiro.test(model1$residuals)
##
## Shapiro-Wilk normality test
##
## data: model1$residuals
## W = 0.93136, p-value = 0.1045
leveneTest(Avg_Days_to_Shoot ~ Rep, data = shoot)
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
leveneTest(Avg_Days_to_Shoot ~ PM, data = shoot)
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
leveneTest(Avg_Days_to_Shoot ~ IBA, data = shoot)
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
Analysis of Variance
splitplot(shoot[4], shoot$Rep, shoot$PM, shoot$IBA, 3)
## $Avg_Days_to_Shoot
## $Avg_Days_to_Shoot[[1]]
## $Avg_Days_to_Shoot[[1]][[1]]
## Analysis of Variance Table
##
## Response: dependent.var
## Df Sum Sq Mean Sq F value Pr(>F)
## block 2 2.3025 1.1513 6.2653 0.13764
## main.plot 1 0.6338 0.6338 3.4490 0.20441
## Ea 2 0.3675 0.1838
## sub.plot 3 13.3446 4.4482 2.7739 0.08715 .
## main.plot:sub.plot 3 1.4346 0.4782 0.2982 0.82607
## Eb 12 19.2433 1.6036
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## $Avg_Days_to_Shoot[[1]][[2]]
## [1] "CV(a): 62.351 , CV(b) : 184.195"
##
## $Avg_Days_to_Shoot[[1]][[3]]
## [1] "R Square 0.484"
##
## $Avg_Days_to_Shoot[[1]][[4]]
##
## Shapiro-Wilk normality test
##
## data: model$residuals
## W = 0.93136, p-value = 0.1045
##
##
## $Avg_Days_to_Shoot[[1]][[5]]
## [1] "Normality assumption is not violated"
##
## $Avg_Days_to_Shoot[[1]][[6]]
## [1] "All the main plot factor level means are same so dont go for any multiple comparison test"
##
## $Avg_Days_to_Shoot[[1]][[7]]
## $Avg_Days_to_Shoot[[1]][[7]][[1]]
## MSerror Df Mean CV MSD
## 0.18375 2 0.6875 62.35065 0.7523167
##
## $Avg_Days_to_Shoot[[1]][[7]][[2]]
## test name.t ntr StudentizedRange alpha
## Tukey main.plot 2 6.079637 0.05
##
## $Avg_Days_to_Shoot[[1]][[7]][[3]]
## dependent.var groups
## M1 0.850 a
## M2 0.525 a
##
##
## $Avg_Days_to_Shoot[[1]][[8]]
## [1] "All the sub plot factor factor level means are same so dont go for any multiple comparison test"
##
## $Avg_Days_to_Shoot[[1]][[9]]
## $Avg_Days_to_Shoot[[1]][[9]][[1]]
## MSerror Df Mean CV MSD
## 1.603611 12 0.6875 184.1946 2.170624
##
## $Avg_Days_to_Shoot[[1]][[9]][[2]]
## test name.t ntr StudentizedRange alpha
## Tukey sub.plot 4 4.19866 0.05
##
## $Avg_Days_to_Shoot[[1]][[9]][[3]]
## dependent.var groups
## T1 1.7833333 a
## T0 0.9666667 a
## T2 0.0000000 a
## T3 0.0000000 a
##
##
## $Avg_Days_to_Shoot[[1]][[10]]
## [1] "All the interaction level means are same so dont go for any multiple comparison test"
##
## $Avg_Days_to_Shoot[[1]][[11]]
## $Avg_Days_to_Shoot[[1]][[11]][[1]]
## MSerror Df Mean CV MSD
## 1.603611 12 0.6875 184.1946 3.742356
##
## $Avg_Days_to_Shoot[[1]][[11]][[2]]
## test name.t ntr StudentizedRange alpha
## Tukey main.plot:sub.plot 8 5.118658 0.05
##
## $Avg_Days_to_Shoot[[1]][[11]][[3]]
## dependent.var groups
## M1:T1 2.366667 a
## M2:T1 1.200000 a
## M1:T0 1.033333 a
## M2:T0 0.900000 a
## M1:T2 0.000000 a
## M1:T3 0.000000 a
## M2:T2 0.000000 a
## M2:T3 0.000000 a