durbin.test {agricolae} | R Documentation

## Durbin test and multiple comparison of treatments

### Description

A multiple comparison of the Durbin test for the balanced incomplete blocks for sensorial or categorical evaluation. It forms groups according to the demanded ones for level of significance (alpha); by default, 0.05.

### Usage

    durbin.test(judge, trt, evaluation, alpha = 0.05, group =TRUE,
    main = NULL, console=FALSE)

### Arguments

`judge` | Identification of the judge in the evaluation
`trt` | Treatments
`evaluation` | variable
`alpha` | level of significant
`group` | TRUE or FALSE
`main` | Title
`console` | logical, print output

### Details

The post hoc test is using the criterium Fisher's least significant difference.

### Value

`statistics` | Statistics of the model
`parameters` | Design parameters
`means` | Statistical summary of the study variable
`rank` | rank table of the study variable
`comparison` | Comparison between treatments
`groups` | Formation of treatment groups

### Author(s)

Felipe de Mendiburu

### References

Practical Nonparametrics Statistics. W.J. Conover, 1999 Nonparametric Statistical Methods. Myles Hollander and Douglas A. Wofe, 1999

### See Also

`BIB.test`, `DAU.test`, `duncan.test`, `friedman`, `HSD.test`, `kruskal`, `LSD.test`, `Median.test`, `PBIB.test`, `REGW.test`, `scheffe.test`, `SNK.test`, `waerden.test`, `waller.test`, `plot.group`

### Examples

    library(agricolae)
    # Example 1. Conover, pag 391
    person<-gl(7,3)
    variety<-c(1,2,4,2,3,5,3,4,6,4,5,7,1,5,6,2,6,7,1,3,7)
    preference<-c(2,3,1,3,1,2,2,1,3,1,2,3,3,1,2,3,1,2,3,1,2)
    out<-durbin.test(person,variety,preference,group=TRUE,console=TRUE,
    main="Seven varieties of ice cream manufacturer")
    #startgraph
    bar.group(out$groups,horiz=TRUE,xlim=c(0,10),density=4,las=1)
    #endgraph
    # Example 2. Myles Hollander, pag 311
    # Source: W. Moore and C.I. Bliss. 1942
    day<-gl(7,3)
    chemical<-c("A","B","D","A","C","E","C","D","G","A","F","G","B","C","F",
     "B","E","G","D","E","F")
    toxic<-c(0.465,0.343,0.396,0.602,0.873,0.634,0.875,0.325,0.330,0.423,0.987,
    0.426,0.652,1.142,0.989,0.536,0.409,0.309,0.609,0.417,0.931)
    out<-durbin.test(day,chemical,toxic,group=TRUE,console=TRUE,
    main="Logarithm of Toxic Dosages")
    plot(out)

[Package agricolae version 1.2-9 Index]