Displaying 4 results from an estimated 4 matches for "inderror".
2004 Jan 30
0
Two apparent bugs in aov(y~ *** -1 + Error(***)), with suggested (PR#6510)
...>
***
*** aov() code with suggested fixes:
***
fixaov <- function(formula, data = NULL, projections = FALSE, qr = TRUE,
contrasts = NULL, ...)
{
Terms <- if(missing(data)) terms(formula, "Error")
else terms(formula, "Error", data = data)
indError <- attr(Terms, "specials")$Error
if(length(indError) > 1)
stop(paste("There are", length(indError),
"Error terms: only 1 is allowed"))
lmcall <- Call <- match.call()
lmcall[[1]] <- as.name("lm")
lmcall...
2004 Feb 02
0
Two apparent bugs in aov(y~ *** -1 + Error(***)), with (PR#6520)
...ed fixes:
> ***
>
> fixaov <- function(formula, data = NULL, projections = FALSE, qr = TRUE,
> contrasts = NULL, ...)
> {
> Terms <- if(missing(data)) terms(formula, "Error")
> else terms(formula, "Error", data = data)
> indError <- attr(Terms, "specials")$Error
> if(length(indError) > 1)
> stop(paste("There are", length(indError),
> "Error terms: only 1 is allowed"))
> lmcall <- Call <- match.call()
> lmcall[[1]] <- as.name(&...
2005 Jun 23
4
contrats hardcoded in aov()?
On 6/23/05, RenE J.V. Bertin <rjvbertin at gmail.com> wrote:
> Hello,
>
> I was just having a look at the aov function source code, and see that when the model used does not have an Error term, Helmert contrasts are imposed:
>
> if (is.null(indError)) {
> ...
> }
> else {
> opcons <- options("contrasts")
> options(contrasts = c("contr.helmert", "contr.poly"))
> on.exit(options(opcons))
> ...
>
>
> My reading of several contributed user guides' s...
2001 Dec 17
1
environments again
In a previous message I was not clear enough in my querry.
I have the following program:
tst<- function() {
x <- c(32.7,32.3,31.5,32.1,29.7,29.1,35.7,35.9,33.1,
36.0,34.2,31.2,31.8,28.0,29.2,38.2,37.8,31.9,
32.5,31.1,29.7)
g <- rep(1:7,rep(3,7))
s <- rep(1:3,7)
cat(" Only x and g \n")
aov1(x,g)
cat("\n\n Now x, g and s \n")
aov1(x,g,s=s)
}