Ah, that was it. I had a bad row in there that I had forgotten to remove.
Thank you very much for the prompt (and correct!) response.

 -Harlan

On Tue, Dec 16, 2008 at 3:58 PM, Prof Brian Ripley <ripley@stats.ox.ac.uk>wrote:

> Your design seems to be unbalanced: multistatum aov is intended for
> balanced designs.  My guess is that one idio subject has two Type=1
> observations: in which case try removing one of them.
>
>
> On Tue, 16 Dec 2008, Harlan Harris wrote:
>
>  Hi, I'm a new R user, coming from SPSS, and without a particularly strong
>> stats background.
>>
>> I've got a data set that I'd like to do a mixed-design ANOVA with. No
>> missing values. Here's the summary:
>>
>> summary(learnDat.ae)
>> Type      Subject        idio     struct     TrainErrs            cond
>> 0:20   11     : 3   idio   :28   ae  :58   Min.   : 0.00   idioae   :28
>> 2:19   12     : 3   nonidio:30   fact: 0   1st Qu.: 6.25   idiofact : 0
>> 3:19   14     : 3                          Median :11.50   nonidioae:30
>>       15     : 3                          Mean   :13.40
>>       18     : 3                          3rd Qu.:16.00
>>       2      : 3                          Max.   :59.00
>>       (Other):40
>>
>> Note that the TrainErrs column is the only numeric column, and I forced
>> everything else to be a factor. (Is that correct?)
>>
>> I then do the following:
>>
>> aov.errs.ae <- aov(TrainErrs ~ (idio*Type) + Error(Subject/Type) +
>> (idio),
>> learnDat.ae)
>>
>> So, idio is between-subjects and Type is within-subjects. This is based on
>> examples I've found elsewhere.
>>
>> summary(aov.errs.ae)
>>
>> This seems to work fine:
>>
>> Error: Subject
>>         Df Sum Sq Mean Sq F value Pr(>F)
>> idio       1    179     179    0.89   0.36
>> Type       1    210     210    1.05   0.32
>> Residuals 17   3401     200
>>
>> Error: Subject:Type
>>         Df Sum Sq Mean Sq F value Pr(>F)
>> Type       2    515     258    2.44  0.103
>> idio:Type  2    680     340    3.22  0.053 .
>> Residuals 34   3595     106
>> ---
>> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>
>>
>> Now the problem:
>>
>>  model.tables(aov.errs.ae,"means")
>>>
>> Error in outer(rownames(efficiency), colnames(efficiency),
>> paste)[eff.used]
>> :
>>  invalid subscript type 'list'
>> In addition: Warning message:
>> In any(efficiency) : coercing argument of type 'double' to logical
>>
>> All the examples and manuals I've found said this should work. When I did
>> a
>> fully between-subjects ANOVA on another data set, I had no problem with
>> model.tables. I have no idea what to make of this error message. I've
>> tried
>> a number of variations on things, with no improvement.
>>
>> This is R version 2.7.2 (2008-08-25), running on Redhat, x86_64.
>>
>> Suggestions? Thanks!
>>
>> -Harlan
>>
>>        [[alternative HTML version deleted]]
>>
>> ______________________________________________
>> R-help@r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>>
> --
> Brian D. Ripley,                  ripley@stats.ox.ac.uk
> Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/<http://www.stats.ox.ac.uk/%7Eripley/>
> University of Oxford,             Tel:  +44 1865 272861 (self)
> 1 South Parks Road,                     +44 1865 272866 (PA)
> Oxford OX1 3TG, UK                Fax:  +44 1865 272595
>

	[[alternative HTML version deleted]]

