Contents

interblock()

Usage:

interblock(y,block,trt,[contrast:coefs])

Keywords: analysis, anova

interblock(y,block,trt,[contrast:coefs]) does recovery of interblock
information in an incomplete block design, where y is the vector of
responses, block is the factor of block levels, trt is the factor of
treatment levels, and coefs is an optional vector of contrast
coefficients (it must be the same length as the number of levels of
trt).  If no contrast is specified, the output is the intra-block,
inter-block, and combined estimates of treatment effects and their
standard errors.  If a contrast is specified, the output is the
intra-block, inter-block, and combined estimates of the contrast with
their standard errors.

Example:

  Cmd> y<-vector(19,17,11,6,26,23,21,19,28,20,7,20,17,26,19,15,23,31,\
       20,26,31,16,23,21,13,7,20,20,24,19,17,6,29,14,24,21)

  Cmd> session<-factor(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,6,6,6,7,7,7,8,8,\
       8,9,9,9,10,10,10,11,11,11,12,12,12)

  Cmd> trt<-factor(1,2,3,4,5,6,7,8,9,1,4,7,2,5,8,3,6,9,1,5,9,2,6,7,3,4,\
       8,1,6,8,2,4,9,3,5,7)

  Cmd> interblock(y,session,trt)
     intra est   intra se  inter est   inter se combined est combined se
       0.33333    0.49414    0.33333    0.91174     0.33333     0.43443
       -2.2222    0.49414         -4    0.91174     -2.6259     0.43443
       -6.2222    0.49414         -6    0.91174     -6.1718     0.43443
       -12.889    0.49414        -13    0.91174     -12.914     0.43443
        5.8889    0.49414     6.6667    0.91174      6.0655     0.43443
        3.5556    0.49414     4.6667    0.91174      3.8078     0.43443
        1.6667    0.49414    0.33333    0.91174      1.3639     0.43443
      -0.22222    0.49414 8.8818e-16    0.91174    -0.17177     0.43443
        10.111    0.49414         11    0.91174      10.313     0.43443

  Cmd> cfs<-vector(.25,.25,.25,.25,-.25,-.25,-.25,-.25,0)

  Cmd> interblock(y,session,trt,contrast:cfs)
               estimate           se
  intra         -7.9722       0.3706
  inter         -8.5833      0.68381
  combined       -8.111      0.32583