Sunday, January 15, 2023

RCBD SPlit Plot

RCBD Split Plot Design

2023-01-15


library(agricolae)

load(file= "C:/Users/khurr/Desktop/Iefan_3rd-study/session_1")


1- Analysis of Variance Table

kebe<-with(data_irfan,sp.plot(replication,Treatment,varities,plant_height))

## ANALYSIS SPLIT PLOT:  plant_height
## Class level information
##
## Treatment    :  Control Gypsum Manure Gypsum + Manure
## varities     :  SG-2002 KS-333 Desi Baramsi
## replication  :  1 2 3
##
## Number of observations:  48
##
## Analysis of Variance Table
##
## Response: plant_height
##                    Df Sum Sq Mean Sq F value    Pr(>F)   
## replication         2  62.54  31.271     NaN       NaN   
## Treatment           3 779.75 259.917  12.031  0.005987 **
## Ea                  6 129.63  21.604     NaN       NaN   
## varities            3 663.08 221.028  16.975 3.946e-06 ***
## Treatment:varities  9 132.42  14.713   1.130  0.380619   
## Eb                 24 312.50  13.021     NaN       NaN   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## cv(a) = 6.1 %, cv(b) = 4.7 %, Mean = 76.29167

gla<-kebe$gl.a
glb<-kebe$gl.b
Ea<-kebe$gl.a
Eb<-kebe$gl.b


2- Main effect of Treatments

ms<-with(data_irfan,LSD.test(plant_height,Treatment,gla,Ea, console = TRUE))

#
## Study: plant_height ~ Treatment
##
## LSD t Test for plant_height
##
## Mean Square Error:  6
##
## Treatment,  means and individual ( 95 %) CI
##
##                 plant_height      std  r      LCL      UCL Min Max
## Control             71.83333 6.161956 12 70.10311 73.56356  61  81
## Gypsum              77.41667 5.632024 12 75.68644 79.14689  70  89
## Gypsum + Manure     82.33333 5.898125 12 80.60311 84.06356  76  95
## Manure              73.58333 3.704011 12 71.85311 75.31356  68  80
##
## Alpha: 0.05 ; DF Error: 6
## Critical Value of t: 2.446912
##
## least Significant Difference: 2.446912
##
## Treatments with the same letter are not significantly different.
##
##                 plant_height groups
## Gypsum + Manure     82.33333      a
## Gypsum              77.41667      b
## Manure              73.58333      c
## Control             71.83333      c

2- Main effect of varities

ms<-with(data_irfan,LSD.test(plant_height,varities,glb,Eb, console = TRUE))

##
## Study: plant_height ~ varities
##
## LSD t Test for plant_height
##
## Mean Square Error:  24
##
## varities,  means and individual ( 95 %) CI
##
##         plant_height      std  r      LCL      UCL Min Max
## Baramsi     71.50000 5.992419 12 68.58121 74.41879  61  79
## Desi        74.08333 3.679386 12 71.16454 77.00213  68  79
## KS-333      81.00000 6.714976 12 78.08121 83.91879  74  95
## SG-2002     78.58333 5.853644 12 75.66454 81.50213  71  90
##
## Alpha: 0.05 ; DF Error: 24
## Critical Value of t: 2.063899
##
## least Significant Difference: 4.127797
##
## Treatments with the same letter are not significantly different.
##
##         plant_height groups
## KS-333      81.00000      a
## SG-2002     78.58333      a
## Desi        74.08333      b
## Baramsi     71.50000      b

3: Interaction Effect

ms<-with(data_irfan,LSD.test(plant_height,Treatment:varities,glb,Eb, console = TRUE))

## Study: plant_height ~ Treatment:varities

##
## LSD t Test for plant_height
##
## Mean Square Error:  24
##
## Treatment:varities,  means and individual ( 95 %) CI
##
##                         plant_height       std r      LCL      UCL Min Max
## Control:Baramsi             63.33333 2.0816660 3 57.49575 69.17092  61  65
## Control:Desi                71.66667 2.0816660 3 65.82908 77.50425  70  74
## Control:KS-333              77.66667 3.5118846 3 71.82908 83.50425  74  81
## Control:SG-2002             74.66667 4.0414519 3 68.82908 80.50425  71  79
## Gypsum + Manure:Baramsi     77.33333 1.5275252 3 71.49575 83.17092  76  79
## Gypsum + Manure:Desi        78.66667 0.5773503 3 72.82908 84.50425  78  79
## Gypsum + Manure:KS-333      89.00000 5.2915026 3 83.16241 94.83759  85  95
## Gypsum + Manure:SG-2002     84.33333 5.5075705 3 78.49575 90.17092  79  90
## Gypsum:Baramsi              75.66667 2.0816660 3 69.82908 81.50425  74  78
## Gypsum:Desi                 73.66667 3.2145503 3 67.82908 79.50425  70  76
## Gypsum:KS-333               81.00000 7.2111026 3 75.16241 86.83759  75  89
## Gypsum:SG-2002              79.33333 7.5718778 3 73.49575 85.17092  74  88
## Manure:Baramsi              69.66667 2.0816660 3 63.82908 75.50425  68  72
## Manure:Desi                 72.33333 3.7859389 3 66.49575 78.17092  68  75
## Manure:KS-333               76.33333 3.2145503 3 70.49575 82.17092  74  80
## Manure:SG-2002              76.00000 1.0000000 3 70.16241 81.83759  75  77
##
## Alpha: 0.05 ; DF Error: 24
## Critical Value of t: 2.063899
##
## least Significant Difference: 8.255594
##
## Treatments with the same letter are not significantly different.
##
##                         plant_height groups
## Gypsum + Manure:KS-333      89.00000      a
## Gypsum + Manure:SG-2002     84.33333     ab
## Gypsum:KS-333               81.00000    abc
## Gypsum:SG-2002              79.33333    bcd
## Gypsum + Manure:Desi        78.66667    bcd
## Control:KS-333              77.66667   bcde
## Gypsum + Manure:Baramsi     77.33333   bcde
## Manure:KS-333               76.33333   bcde
## Manure:SG-2002              76.00000    cde
## Gypsum:Baramsi              75.66667    cde
## Control:SG-2002             74.66667    cde
## Gypsum:Desi                 73.66667    cde
## Manure:Desi                 72.33333     de
## Control:Desi                71.66667     de
## Manure:Baramsi              69.66667     ef
## Control:Baramsi             63.33333      f

No comments:

Post a Comment