[R] Error in La.svd(X) : error code 1 from Lapack routine 'dgesdd'

[Ranjan Maitra]

On Mon, 05 Mar 2007 09:14:17 +0100 Sophie Richier <Sophie.Richier at obs-vlfr.fr> wrote:

> Dear R helpers,
> I am working with R 2.4.1 GUI 1.18 (4038) for MacOSX. I have a matrix of 
> 10 000 genes  and try to run the following commands:
>  > model.mix<-makeModel (data=data, formula=~Dye+Array+Sample+Time, 
> random=~Array+Sample)
>  > anova.mix<-fitmaanova (data, model.mix)
>  > test.mix<-matest (data, model=model.mix, term="Time", n.perm=100, 
> test.method=c(1,0,1,1))
> 
> I get the following error message:
> Doing F-test on observed data ...
> Doing permutation. This may take a long time ...
> Error in La.svd(X) : error code 1 from Lapack routine 'dgesdd'
> 
> What does this mean? is my matrix too big? What can I do?
> Thanks a lot in adavance
> 
> Sophie


from the help file:



     Unsuccessful results from the underlying LAPACK code will result
     in an error giving a positive error code: these can only be
     interpreted by detailed study of the FORTRAN code.


from the manpages:

man dgesdd

       INFO    (output) INTEGER
               = 0:  successful exit.
               < 0:  if INFO = -i, the i-th argument had an illegal value. 
               > 0:  DBDSDC did not converge, updating process failed.

I don't know what DBDSDC is, but it appears that there may be some convergence issue for you. Unless someone else has better ideas, look up www.netlib.org/lapack and the routines in there to investigate further.

HTH!

Best,
Ranjan