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

From: Giles Crane <gilescrane_at_verizon.net>
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:

    library(effects)

    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,
Giles

-- 
View this message in context: http://r.789695.n4.nabble.com/Adjusted-means-and-generalized-chain-block-designs-tp1556545p2267205.html
Sent from the R help mailing list archive at Nabble.com.

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
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 https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive