The data is not available !
str(LandU)
## 'data.frame': 250 obs. of 13 variables:
## $ landuse : chr "Grazed" "Grazed" "Grazed" "Grazed" ...
## $ x.axis : num 105 300 255 195 150 45 315 180 45 135 ...
## $ y.axis : num 300 390 360 300 210 120 420 330 240 420 ...
## $ code : chr "R1 / B-300" "T1 / B-390" "S4 / B-360" "R7 / B-300" ...
## $ S_Series : chr "PeoC" "PeoC" "PeoC" "PeoC" ...
## $ distance : chr "30" "30" "30" "30" ...
## $ Latitude : logi NA NA NA NA NA NA ...
## $ Longitude: logi NA NA NA NA NA NA ...
## $ Elevation: logi NA NA NA NA NA NA ...
## $ Total_N : num 0.12 0.14 0.26 0.32 0.3 0.17 0.17 0.29 0.21 0.15 ...
## $ Total_C : num 1.27 1.34 2.77 3.23 3.13 1.83 1.68 3.11 2.27 1.48 ...
## $ OM : num 2.19 2.31 4.78 5.57 5.4 3.15 2.89 5.36 3.91 2.54 ...
## $ C_N : num 10.7 9.9 10.6 10.2 10.4 10.6 10.1 10.8 10.9 10.1 ...
library(agricolae)
model<-aov(Total_C ~ landuse + distance + S_Series + landuse * distance * S_Series, data=LandU)
DF_LandU <- summary(model)
DF_LandU # Landuse impact on Total C
## Df Sum Sq Mean Sq F value Pr(>F)
## landuse 1 12.74 12.738 54.135 3.15e-12 ***
## distance 3 2.26 0.753 3.202 0.024048 *
## S_Series 3 1.48 0.492 2.090 0.102219
## landuse:distance 2 0.70 0.351 1.492 0.227007
## landuse:S_Series 2 4.24 2.118 9.000 0.000172 ***
## distance:S_Series 4 1.30 0.324 1.377 0.242731
## Residuals 234 55.06 0.235
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Study: IMPACT OF SOIL SERIES ON CARBON CONTENT
##
## HSD Test for Total_C
##
## Mean Square Error: 0.2352938
##
## landuse, means
##
## Total_C std r Min Max
## Grazed 2.044200 0.6719250 100 0.250 4.050
## Skeet 1.583447 0.3693867 150 0.767 3.251
##
## Alpha: 0.05 ; DF Error: 234
## Critical Value of Studentized Range: 2.786218
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## Total_C groups
## Grazed 2.044200 a
## Skeet 1.583447 b
plot(out)
# The End