[S] multicomp and nested models

["Ed Callahan" <callahan@envstat.com>]

Does anyone know how to get multicomp() to work with nested models?  Below
is an example function that demonstrates what I'm trying to do.  It gives me
the error message "focus must be a single factor".  Thanks for any ideas.

Ed

function()
{
        whole <- as.factor(rep(LETTERS[1:3], 10))
        sub <- as.factor(c(rep("a", 15), rep("b", 15)))
        y <- rnorm(30)
        assign("test.df", data.frame(y, whole, sub), w = 1)
        test.aov <- aov(y ~ whole/sub, data = test.df)
        multicomp(test.aov, focus = "sub")
}

----------------------------------
  Edward Callahan
  Environmental Statistics

  callahan@envstat.com
  http://www.envstat.com
----------------------------------

-----------------------------------------------------------------------
This message was distributed by s-news@wubios.wustl.edu.  To unsubscribe
send e-mail to s-news-request@wubios.wustl.edu with the BODY of the
message:  unsubscribe s-news