Displaying 2 results from an estimated 2 matches for "this_aov".
2008 Jul 20
0
Off topic: SS formulae for 3-way repeated measure anova (for when aov() fails)
....b = 2
levs.d = 10
temp.sub = factor(1:N)
temp.a = factor(1:levs.a)
temp.b = factor(1:levs.b)
temp.d = factor(1:levs.d)
temp = expand.grid(sub=temp.sub, a=temp.a, b=temp.b, d=temp.d)
temp$x = rnorm(length(temp[, 1])) #generate some random DV data
sub=temp$sub
a=temp$a
b=temp$b
d=temp$d
x=temp$x
this_aov = aov(
x~a*b*d+Error(sub/(a*b*d))
)
summary(this_aov)
#now let's try by hand, checking each sum-of-squares
# ss against the analogous aov() produced ss (rounding
# each to avoid small computational differences)
#Get ss.subs
sub.means = aggregate(x,list(sub=sub), mean)
grand.mean = mean(sub....
2008 Jul 15
1
aov error with large data set
...alyze a large data set: a within-Ss 2*2*1500 design
with 20 Ss. However, aov() gives me an error, reproducible as follows:
id = factor(1:20)
a = factor(1:2)
b = factor(1:2)
d = factor(1:1500)
temp = expand.grid(id=id, a=a, b=b, d=d)
temp$y = rnorm(length(temp[, 1])) #generate some random DV data
this_aov = aov(
y~a*b*d+Error(id/(a*b*d))
, data=temp
)
While yields the following error:
"
Error in model.matrix.default(mt, mf, contrasts) :
allocMatrix: too many elements specified
"
Any suggestions?
Mike
--
Mike Lawrence
Graduate Student, Department of Psychology, Dalhousie University...