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
"Simon Blomberg" <Simon.Blomberg at anu.edu.au> writes:> 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.So were you expecting "residuals have rank 36 < 32" in that case?? Some multivariate tests require that you can obtain a full rank estimate of the covariance matrix, which requires the residual degrees of freedom to be less than the number of response variables. Not sure whether there is a meaningful way to proceed if that is not the case - I suspect you'll get negative DF for the F tests and similar awkwardness. -- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk) FAX: (+45) 35327907
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
Sorry for reposting. I didn't receive any of the replies to my original message (must be an email problem at my end). I will read the responses in the archive. 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> -----Original Message----- > From: Simon Blomberg > Sent: Monday, 24 November 2003 9:09 AM > To: R-help at stat.math.ethz.ch > Subject: [R] summary.manova and rank deficiency > > > 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 > > ______________________________________________ > R-help at stat.math.ethz.ch mailing list > https://www.stat.math.ethz.ch/mailman/listinfo/r-help >