crd: Analysis of Completely Randomized Design

In doebioresearch: Analysis of Design of Experiments for Biological Research

Description

The function gives ANOVA, R-square of the model, normality testing of residuals, SEm (standard error of mean), SEd (standard error of difference), interpretation of ANOVA results and multiple comparison test for means

Usage

crd(data, trt.vector, MultipleComparisonTest)

Arguments

data	dependent variables
trt.vector	vector containing treatments
MultipleComparisonTest	0 for no test, 1 for LSD test, 2 for Duncan test and 3 for HSD test

Value

ANOVA, interpretation of ANOVA, R-square, normality test result, SEm, SEd and multiple comparison test result

Examples

data<-data.frame(Treatments=c("T1","T2","T3","T4","T5","T6","T7","T1","T2","T3","T4","T5","T6",
"T7","T1","T2","T3","T4","T5","T6","T7"),
yield=c(25,21,21,18,25,28,24,25,24,24,16,21,20,17,16,19,14,15,13,11,25),
height=c(130,120,125,135,139,140,145,136,129,135,150,152,140,148,130,135,145,160,145,130,160))
#CRD analysis with LSD test for yield only
crd(data[2],data$Treatments,1)
#CRD analysis with LSD test for both yield and height
crd(data[2:3],data$Treatments,1)

Example output

$yield
$yield[[1]]
$yield[[1]][[1]]
Analysis of Variance Table

Response: data2
          Df Sum Sq Mean Sq F value Pr(>F)
trt        6  70.48  11.746  0.4312 0.8461
Residuals 14 381.33  27.238               

$yield[[1]][[2]]
[1] "All the treatment means are same so dont go for any multiple comparison test"

$yield[[1]][[3]]
[1] "R Square 0.156"

$yield[[1]][[4]]

	Shapiro-Wilk normality test

data:  model$residuals
W = 0.95693, p-value = 0.4565


$yield[[1]][[5]]
[1] "Normality assumption is not violated"

$yield[[1]][[6]]
[1] "SEm 3.0132 , SEd 4.2613"

$yield[[1]][[7]]
$yield[[1]][[7]][[1]]
  MSerror Df     Mean       CV  t.value      LSD
  27.2381 14 20.09524 25.97139 2.144787 9.139593

$yield[[1]][[7]][[2]]
      data2 groups
T1 22.00000      a
T7 22.00000      a
T2 21.33333      a
T3 19.66667      a
T5 19.66667      a
T6 19.66667      a
T4 16.33333      a




$yield
$yield[[1]]
$yield[[1]][[1]]
Analysis of Variance Table

Response: data2
          Df Sum Sq Mean Sq F value Pr(>F)
trt        6  70.48  11.746  0.4312 0.8461
Residuals 14 381.33  27.238               

$yield[[1]][[2]]
[1] "All the treatment means are same so dont go for any multiple comparison test"

$yield[[1]][[3]]
[1] "R Square 0.156"

$yield[[1]][[4]]

	Shapiro-Wilk normality test

data:  model$residuals
W = 0.95693, p-value = 0.4565


$yield[[1]][[5]]
[1] "Normality assumption is not violated"

$yield[[1]][[6]]
[1] "SEm 3.0132 , SEd 4.2613"

$yield[[1]][[7]]
$yield[[1]][[7]][[1]]
  MSerror Df     Mean       CV  t.value      LSD
  27.2381 14 20.09524 25.97139 2.144787 9.139593

$yield[[1]][[7]][[2]]
      data2 groups
T1 22.00000      a
T7 22.00000      a
T2 21.33333      a
T3 19.66667      a
T5 19.66667      a
T6 19.66667      a
T4 16.33333      a




$height
$height[[1]]
$height[[1]][[1]]
Analysis of Variance Table

Response: data2
          Df Sum Sq Mean Sq F value Pr(>F)  
trt        6 1383.2 230.540   3.463  0.026 *
Residuals 14  932.0  66.571                 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

$height[[1]][[2]]
[1] "The treatment means of one or more treatments are not same, so go for multiple comparison test"

$height[[1]][[3]]
[1] "R Square 0.597"

$height[[1]][[4]]

	Shapiro-Wilk normality test

data:  model$residuals
W = 0.98167, p-value = 0.9471


$height[[1]][[5]]
[1] "Normality assumption is not violated"

$height[[1]][[6]]
[1] "SEm 4.7107 , SEd 6.6619"

$height[[1]][[7]]
$height[[1]][[7]][[1]]
   MSerror Df     Mean       CV  t.value      LSD
  66.57143 14 139.4762 5.849838 2.144787 14.28836

$height[[1]][[7]][[2]]
      data2 groups
T7 151.0000      a
T4 148.3333     ab
T5 145.3333    abc
T6 136.6667    bcd
T3 135.0000    bcd
T1 132.0000     cd
T2 128.0000      d

doebioresearch documentation built on July 8, 2020, 7:18 p.m.