search for: maov

...e return values of aov. I know that it has 'coefficients'. But I need to know what all the other return values are. Can somebody let me know how to figure them? Value: An object of class 'c("aov", "lm")' or for multiple responses of class 'c("maov", "aov", "mlm", "lm")' or for multiple error strata of class '"aovlist"'. There are 'print' and 'summary' methods available for these. Regards, Peng

...qr <- lmcall$qr <- TRUE lmcall$projections <- NULL if(is.null(indError)) { ## no Error term fit <- eval(lmcall, parent.frame()) if(projections) fit$projections <- proj(fit) class(fit) <- if(inherits(fit, "mlm")) c("maov", "aov", oldClass(fit)) else c("aov", oldClass(fit)) fit$call <- Call return(fit) } else { ## helmert contrasts can be helpful: do we want to force them? ## this version does for the Error model. opcons <- options("co...

...lmcall$projections <- NULL > if(is.null(indError)) { > ## no Error term > fit <- eval(lmcall, parent.frame()) > if(projections) fit$projections <- proj(fit) > class(fit) <- if(inherits(fit, "mlm")) > c("maov", "aov", oldClass(fit)) else c("aov", oldClass(fit)) > fit$call <- Call > return(fit) > } else { > ## helmert contrasts can be helpful: do we want to force them? > ## this version does for the Error model. >...

I am using rw0641 with Windows 98. In the help for aov it states that the formula can sepecify multiple responses for a "maov". The help doesn't give an example, Venables and Ripley 1997:381 doesn't either. It isn't mentioned in the scripts nor in the section R Complements nor in Rnotes. I tried aov(c(y1,y2,y3)~x) which seemed reasonable to me but doesn't work. How does it function? Thank you ver...

Hi all, I am experimenting the function "manova" in R. I tried it on a few data sets, but I did not understand the result: I used "summary(manova_result)" and "summary(manova_result, test='Wilks')" and they gave a bunch of numbers... But I need the Sum-of-Squares of BETWEEN and WITHIN matrices... How do I read off from the R's manova results? Any

Hi R-users I'm trying to do multivariate analysis of variance of a experiment with 3 treatments, 2 variables and 5 replicates. The procedure adopted in SAS is as follow, but I'm having difficulty in to implement the contrasts for comparison of all treatments in R. I have already read manuals and other materials about manova in R, but nothing about specific contrasts were found in them,

Hi everybody ... I have to perform Hotteling's T^2 Test - more generally a MANOVA - on a set of data. Is there a (simple?) possibility to do it in R? The somewhat obvious way to do it would be > summary(aov(Y ~ x1 + x2 ... + xn)) where Y would be a two-column matrix. But this does not work in any possible combination of matrix or factor dimensions! Is it principally not possible

Hello R List: My problem is with a nested anova. I have read the r-help and it has answered some of my questions but i still need some help on this one. I have also posted for help on this data set before, so i apologize in advance for any repetition. My design is as follows: response: Quadrat Counts (individuals per quadrat) Explanatory: Region (3 regions) Locations (4

...z <- if (is.null(w)) { aov.balanced.fit(x, y, mt) } else { stop("Weighted aov not implemented, use lm") #aov.balanced.wfit(x, y, mt, w) } class(z) <- c(if (is.matrix(y)) c("maov","mlm"), c("aov","lm")) } z$call <- match.call() z$terms <- mt if (model) z$model <- mf z } aov.balanced.fit <- function(x,y,mod.terms,wt=NULL){ x <- as.matrix(x); dimx <- dim(x) xnames <- dimnames(x...

I'm sure that many people are in the same position as me in that they are trying to write packages and code that is vaguely "future proof". Would it be possible to get some guidance on how the R-core team see the evolution of the "base" package with regard to s4 methods. There seem to be quite a lot of inconsistencies between s3 and s4 methods and classes currently and