library(readxl)
## Warning: package 'readxl' was built under R version 4.0.5
Yield <- read_excel("C:/akong project/seatwork/new seatwork/Yield.xlsx",
col_types = c("text", "text", "numeric"))
View(Yield)
model<-lm(Yield$Yield~Yield$Replication+Yield$Variety)
anova<-anova(model)
anova
## Analysis of Variance Table
##
## Response: Yield$Yield
## Df Sum Sq Mean Sq F value Pr(>F)
## Yield$Replication 3 80.80 26.934 0.9205 0.46033
## Yield$Variety 4 520.53 130.133 4.4476 0.01958 *
## Residuals 12 351.11 29.259
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
par(mfrow=c(1,2))
plot(model,which = 1)
plot(model,which=2)
library(agricolae)
## Warning: package 'agricolae' was built under R version 4.0.5
dnmrt<-duncan.test(Yield$Yield,Yield$Variety,12,29.259)
dnmrt
## $statistics
## MSerror Df Mean CV
## 29.259 12 31.275 17.29547
##
## $parameters
## test name.t ntr alpha
## Duncan Yield$Variety 5 0.05
##
## $duncan
## Table CriticalRange
## 2 3.081307 8.333639
## 3 3.225244 8.722927
## 4 3.312453 8.958792
## 5 3.370172 9.114897
##
## $means
## Yield$Yield std r Min Max Q25 Q50 Q75
## African tall 30.450 7.403378 4 22.9 39.1 25.150 29.90 35.200
## Co-11 31.200 2.762849 4 29.5 35.3 29.575 30.00 31.625
## Co-24 25.550 5.674798 4 20.4 31.8 20.925 25.00 29.625
## FS-1 28.475 3.155287 4 24.4 32.1 27.550 28.70 29.625
## K-7 40.700 6.274286 4 32.1 47.0 38.700 41.85 43.850
##
## $comparison
## NULL
##
## $groups
## Yield$Yield groups
## K-7 40.700 a
## Co-11 31.200 b
## African tall 30.450 b
## FS-1 28.475 b
## Co-24 25.550 b
##
## attr(,"class")
## [1] "group"
lsd<-LSD.test(Yield$Yield,Yield$Variety,12,29.259)
lsd
## $statistics
## MSerror Df Mean CV t.value LSD
## 29.259 12 31.275 17.29547 2.178813 8.333639
##
## $parameters
## test p.ajusted name.t ntr alpha
## Fisher-LSD none Yield$Variety 5 0.05
##
## $means
## Yield$Yield std r LCL UCL Min Max Q25 Q50
## African tall 30.450 7.403378 4 24.55723 36.34277 22.9 39.1 25.150 29.90
## Co-11 31.200 2.762849 4 25.30723 37.09277 29.5 35.3 29.575 30.00
## Co-24 25.550 5.674798 4 19.65723 31.44277 20.4 31.8 20.925 25.00
## FS-1 28.475 3.155287 4 22.58223 34.36777 24.4 32.1 27.550 28.70
## K-7 40.700 6.274286 4 34.80723 46.59277 32.1 47.0 38.700 41.85
## Q75
## African tall 35.200
## Co-11 31.625
## Co-24 29.625
## FS-1 29.625
## K-7 43.850
##
## $comparison
## NULL
##
## $groups
## Yield$Yield groups
## K-7 40.700 a
## Co-11 31.200 b
## African tall 30.450 b
## FS-1 28.475 b
## Co-24 25.550 b
##
## attr(,"class")
## [1] "group"
sink("Yield.txt")
print(anova)
## Analysis of Variance Table
##
## Response: Yield$Yield
## Df Sum Sq Mean Sq F value Pr(>F)
## Yield$Replication 3 80.80 26.934 0.9205 0.46033
## Yield$Variety 4 520.53 130.133 4.4476 0.01958 *
## Residuals 12 351.11 29.259
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
print("dnmrt Result")
## [1] "dnmrt Result"
print(dnmrt$statistics)
## MSerror Df Mean CV
## 29.259 12 31.275 17.29547
print(dnmrt$groups)
## Yield$Yield groups
## K-7 40.700 a
## Co-11 31.200 b
## African tall 30.450 b
## FS-1 28.475 b
## Co-24 25.550 b
print("lsd Result")
## [1] "lsd Result"
print(lsd$statistics)
## MSerror Df Mean CV t.value LSD
## 29.259 12 31.275 17.29547 2.178813 8.333639
print(lsd$groups)
## Yield$Yield groups
## K-7 40.700 a
## Co-11 31.200 b
## African tall 30.450 b
## FS-1 28.475 b
## Co-24 25.550 b
sink()