Re: [R] Question regarding lu in package Matrix

2008-07-09 Thread Ulrike Grömping
No, not at all, I'm glad that I do get the decomposition for non-square matrices. My problem is not with elu$U but with as.matrix(elu$U), which is not an upper diagonal matrix. Can I do something to fix this ? Regards, Ulrike Douglas Bates-2 wrote: > > On Wed, Jul 9, 2008 at 12:01 PM, Ulrike G

Re: [R] Question regarding lu in package Matrix

2008-07-09 Thread Douglas Bates
On Wed, Jul 9, 2008 at 12:01 PM, Ulrike Grömping <[EMAIL PROTECTED]> wrote: > Dear R-helpers, > I have a question regarding LU-decomposition with function lu in package > Matrix. The following simple example confuses me: Why is as.matrix(elu$U) > not an upper triangular matrix? > u3 <- > matrix(

[R] Question regarding lu in package Matrix

2008-07-09 Thread Ulrike Grömping
Dear R-helpers, I have a question regarding LU-decomposition with function lu in package Matrix. The following simple example confuses me: Why is as.matrix(elu$U) not an upper triangular matrix? u3 <- matrix(c(1,1,1,1,1,1,-1,1,0,0,0,0,0,-1,1,0,0,0,-1,0,1,0,0,0,0,0,-1,1,0,0),5,6,byrow=T) elu <- e