Is this a vanilla installation, not using any BLAS routines? If so, that error code means that the svd was called with a negative number of columns (&p in the call). I don't see how that can happen without some memory corruption.
The usual ideas apply: run that example on its own, then look at traceback(), add debugging code etc .... We have seen similar (but not the same) issues where a BLAS was used that contained (broken) parts of LAPACK. On Mon, 20 Jan 2003, Laurent Gautier wrote: > Dear all, > > The 'make check' step fails for the pacakge mva on IBM AIX. > The tail of the Rout log file looks like: > > > for(factors in 2:4) print(update(Harman23.FA, factors = factors)) > > Call: > factanal(factors = factors, covmat = Harman23.cor) > > Uniquenesses: > height arm.span forearm lower.leg weight > 0.170 0.107 0.166 0.199 0.089 > bitro.diameter chest.girth chest.width > 0.364 0.416 0.537 > > Loadings: > Factor1 Factor2 > height 0.865 0.287 > arm.span 0.927 0.181 > forearm 0.895 0.179 > lower.leg 0.859 0.252 > weight 0.233 0.925 > bitro.diameter 0.194 0.774 > chest.girth 0.134 0.752 > chest.width 0.278 0.621 > > Factor1 Factor2 > SS loadings 3.335 2.617 > Proportion Var 0.417 0.327 > Cumulative Var 0.417 0.744 > > Test of the hypothesis that 2 factors are sufficient. > The chi square statistic is 75.74 on 13 degrees of freedom. > The p-value is 6.94e-11 > > Call: > factanal(factors = factors, covmat = Harman23.cor) > > Uniquenesses: > height arm.span forearm lower.leg weight > 0.127 0.005 0.193 0.157 0.090 > bitro.diameter chest.girth chest.width > 0.359 0.411 0.490 > > Loadings: > Factor1 Factor2 Factor3 > height 0.886 0.267 -0.130 > arm.span 0.937 0.195 0.280 > forearm 0.874 0.188 > lower.leg 0.877 0.230 -0.145 > weight 0.242 0.916 -0.106 > bitro.diameter 0.193 0.777 > chest.girth 0.137 0.755 > chest.width 0.261 0.646 0.159 > > Factor1 Factor2 Factor3 > SS loadings 3.379 2.628 0.162 > Proportion Var 0.422 0.329 0.020 > Cumulative Var 0.422 0.751 0.771 > > Test of the hypothesis that 3 factors are sufficient. > The chi square statistic is 22.81 on 7 degrees of freedom. > The p-value is 0.00184 > Error in La.svd(B) : error code 3 from Lapack routine dgesdd > Execution halted > > > Any hint ? > > > > Laurent > > ______________________________________________ > [EMAIL PROTECTED] mailing list > http://www.stat.math.ethz.ch/mailman/listinfo/r-help > -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595 ______________________________________________ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help
