Re: [R] Adjusted means and generalized chain block designs

From: Giles Crane <>
Date: Thu, 24 Jun 2010 08:56:09 -0700 (PDT)

To whom it may concern:

I compared several R package results,
and manual checked two generalized chain block design experiments. The correct adjusted treatment means can be computed
by using the effects library as follows:


    aov1 = aov(y~blocks+rows+trt)
    means.aov = as.vector( allEffects(aov1)$trt$fit )

As a further check, I programmed an R function to compute the adjusted treatment means
via John Mandel's manual method in his paper, a somewhat lengthy and tricky task.

Finally, I will present motivating examples on behalf of my colleagues Cynthia Collins and Karin Mille, an index of Generalized Block Designs, and GCB library in draft form, includingD-efficiency, at the NIST UseR 2010 Conference
in July.

These designs are interesting combinatorially, and quite applicable to physical measurements which are made with high precision.

Best wishes,

View this message in context:
Sent from the R help mailing list archive at

______________________________________________ mailing list
PLEASE do read the posting guide
and provide commented, minimal, self-contained, reproducible code.
Received on Thu 24 Jun 2010 - 18:44:59 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 24 Jun 2010 - 18:50:35 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive