Next: mixed() Up: Design Macros Help File Previous: interactplot()   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
```

Gary Oehlert 2003-01-15