[R] summary.manova and rank deficiency

[Simon Blomberg]

Hi all,

I have received the following error from summary.manova:

Error in summary.manova(manova.test, test = "Pillai") : 
	residuals have rank 36 < 64

The data is simulated data for 64 variables. The design is a 2*2 factorial with 10 replicates per treatment. Looking at the code for summary.manova,  the error involves a problem with qr(). Does anyone have a suggestion as to how to deal with this error? The analysis works fine when there are only 32 variables.

Thanks,

Simon.

Simon Blomberg, PhD
Depression & Anxiety Consumer Research Unit
Centre for Mental Health Research
Australian National University
http://www.anu.edu.au/cmhr/
Simon.Blomberg at anu.edu.au  +61 (2) 6125 3379