search for: nmstrata

I believe you are right, but can you please explain why anyone would want to fit this model? It differs only in the coding from aov(y ~ a + b + Error(c), data=test.df) and merely lumps together the top two strata. There is a much simpler fix: in the line if(intercept) nmstrata <- c("(Intercept)", nmstrata) remove the condition (and drop the empty stratum later if you want). On Fri, 30 Jan 2004 b-h@mevik.net wrote: > I think there are two bugs in aov() that shows up when the right hand > side of `formula' contains both `-1' and an Error() ter...

...ot;),c(8,8))) sub = factor(sort(spaeVP)) rati = stack(data.frame(arousal)) rati = rati[,1] rating.df <- data.frame(sub, cs, spk, rati) summary(aov(rati ~ cs*spk + Error(sub/(cs*spk)), data=rating.df)) I get this Error (after about 5sec): Error in "names<-.default"(`*tmp*`, value = nmstrata) : names attribute must be the same length as the vector > rating.df sub cs spk rati 1 5 yes spk1 -1 2 5 yes spk1 0 3 5 yes spk2 1 4 5 yes spk2 0 5 5 yes spk3 0 6 5 yes spk3 -1 7 5 yes spk4 -1 8 5 yes spk4 1 9 5 no spk5 0 ......

..., if(!intercept) "- 1"), env=environment(formula)) ecall$method <- "qr" ecall$qr <- TRUE ecall$contrasts <- NULL er.fit <- eval(ecall, parent.frame()) options(opcons) nmstrata <- attr(terms(er.fit), "term.labels") ## remove backticks from simple labels for strata (only) nmstrata <- sub("^`(.*)`$", "\\1", nmstrata) if(intercept) nmstrata <- c("(Intercept)", nmstrata) qr.e <- er.fit$qr...

Hi, OK, I will describe it in more detail.In the meantime I installed the newest debian package of r-base etc. and the error message changed from:Error in "names<-.default"(`*tmp*`, value = nmstrata) : names attribute must be the same length as the vectorto:Error model is singular in: aov(rati ~ cs * spk + Error(sub/(cs * spk)), data = rating.df) The model is a repeated measures anova that's why I use the Error() term. I found the hint for this usage in "Notes on the use of R...

...he raw materials or the process settings cannot be completely randomised, leading to models with error terms. It is also quite common that the experiment has to be run over several days, leading to an error term for day. > There is a much simpler fix: in the line > > if(intercept) nmstrata <- c("(Intercept)", nmstrata) > > remove the condition (and drop the empty stratum later if you want). You are right. I didn't see that one. --- (*) In these models, we have to manually code the factorial terms into continuous variables, in order to circumvent the coding...

...factor2=factor(rep(c("fac2level1","fac2level2","fac2level3"),each=18)) ) sample.aov <- aov(dep.variable ~ factor1 * factor2 + Error(subject/(factor1+factor2)), data=df) But the aov function returns an error: Error in "names<-.default"(`*tmp*`, value = nmstrata) : names attribute must be the same length as the vector I am using R1.8.1. Any help or insight on a workaround would be appreciated. Aries Arditi Senior Fellow Arlene R. Gordon Research Institute Lighthouse International

...+ Error(Subject/(Sickness.norm*WasSick) ) ) To stick with the latter example: my data is sufficiently balanced (about 50% sick subjects), so I would expect that there are sufficient samples to perform an ANOVA with. I get the message "Error in "names<-.default"(`*tmp*`, value = nmstrata) : names attribute must be the same length as the vector" I expect that this has to do with the fact that none of my independent variables are in fact truely independent (controlled by the experimenter), but instead depend on the subject. For the age, gender and 'WasSick' data, ther...