The R Project SVN
R-packages
(root)/branches/Matrix-APIchange/tests/bind.Rout.save @ 4556 – Rev 4529
Rev 4440 | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed
R version 2.5.0 Patched (2007-06-04 r41817)Copyright (C) 2007 The R Foundation for Statistical ComputingISBN 3-900051-07-0R is free software and comes with ABSOLUTELY NO WARRANTY.You are welcome to redistribute it under certain conditions.Type 'license()' or 'licence()' for distribution details.R is a collaborative project with many contributors.Type 'contributors()' for more information and'citation()' on how to cite R or R packages in publications.Type 'demo()' for some demos, 'help()' for on-line help, or'help.start()' for an HTML browser interface to help.Type 'q()' to quit R.> #### Testing cBind() & rBind()>> library(Matrix)Loading required package: lattice>> source(system.file("test-tools.R", package = "Matrix"))# identical3() etc>> ### --- Dense Matrices --->> m1 <- m2 <- m <- Matrix(1:12, 3,4)Warning message:integer matrices not yet implemented in 'Matrix'; using 'double' ones' in: Matrix(1:12, 3, 4)> dimnames(m2) <- list(LETTERS[1:3],+ letters[1:4])> dimnames(m1) <- list(NULL,letters[1:4])>> stopifnot(identical(cBind ( m, 10*m) -> R,+ cbind2( m, 10*m))); R3 x 8 Matrix of class "dgeMatrix"[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8][1,] 1 4 7 10 10 40 70 100[2,] 2 5 8 11 20 50 80 110[3,] 3 6 9 12 30 60 90 120> stopifnot(identical(cBind (m1,100+m1) -> R,+ cbind2(m1,100+m1))); R3 x 8 Matrix of class "dgeMatrix"a b c d a b c d[1,] 1 4 7 10 101 104 107 110[2,] 2 5 8 11 102 105 108 111[3,] 3 6 9 12 103 106 109 112> stopifnot(identical(cBind (m1, 10*m2) -> R,+ cbind2(m1, 10*m2))); R3 x 8 Matrix of class "dgeMatrix"a b c d a b c dA 1 4 7 10 10 40 70 100B 2 5 8 11 20 50 80 110C 3 6 9 12 30 60 90 120> stopifnot(identical(cBind (m2, m1+m2) -> R,+ cbind2(m2, m1+m2))); R3 x 8 Matrix of class "dgeMatrix"a b c d a b c dA 1 4 7 10 2 8 14 20B 2 5 8 11 4 10 16 22C 3 6 9 12 6 12 18 24>> cBind(m1, MM = -1)3 x 5 Matrix of class "dgeMatrix"a b c d MM[1,] 1 4 7 10 -1[2,] 2 5 8 11 -1[3,] 3 6 9 12 -1> rBind(R1 = 10:11, m1)4 x 4 Matrix of class "dgeMatrix"a b c dR1 10 11 10 111 4 7 102 5 8 113 6 9 12> cBind(0, Matrix(0+0:1, 1,2), 3:2)# FIXME? should warn - as with matrix()1 x 4 Matrix of class "dgeMatrix"[,1] [,2] [,3] [,4][1,] 0 0 1 3>> as(rBind(0, Matrix(0+0:1, 1,2), 3:2),+ "sparseMatrix")3 x 2 sparse Matrix of class "dgCMatrix"[1,] . .[2,] . 1[3,] 3 2> cBind(m2, 10*m2[nrow(m2):1 ,])# keeps the rownames from the first3 x 8 Matrix of class "dgeMatrix"a b c d a b c dA 1 4 7 10 30 60 90 120B 2 5 8 11 20 50 80 110C 3 6 9 12 10 40 70 100>> (im <- cBind(I = 100, m))3 x 5 Matrix of class "dgeMatrix"I[1,] 100 1 4 7 10[2,] 100 2 5 8 11[3,] 100 3 6 9 12> str(im)Formal class 'dgeMatrix' [package "Matrix"] with 4 slots..@ x : num [1:15] 100 100 100 1 2 3 4 5 6 7 .....@ Dim : int [1:2] 3 5..@ Dimnames:List of 2.. ..$ : NULL.. ..$ : chr [1:5] "I" "" "" "" .....@ factors : list()> (mi <- cBind(m2, I = 1000))3 x 5 Matrix of class "dgeMatrix"a b c d IA 1 4 7 10 1000B 2 5 8 11 1000C 3 6 9 12 1000> str(mi)Formal class 'dgeMatrix' [package "Matrix"] with 4 slots..@ x : num [1:15] 1 2 3 4 5 6 7 8 9 10 .....@ Dim : int [1:2] 3 5..@ Dimnames:List of 2.. ..$ : chr [1:3] "A" "B" "C".. ..$ : chr [1:5] "a" "b" "c" "d" .....@ factors : list()> (m1m <- cBind(m,I=100,m2))3 x 9 Matrix of class "dgeMatrix"I a b c dA 1 4 7 10 100 1 4 7 10B 2 5 8 11 100 2 5 8 11C 3 6 9 12 100 3 6 9 12>> ### --- Diagonal / Sparse - had bugs>> D4 <- Diagonal(4)> (D4T <- as(D4, "sparseMatrix"))4 x 4 sparse Matrix of class "dtTMatrix"[1,] 1 . . .[2,] . 1 . .[3,] . . 1 .[4,] . . . 1> D4C <- as(D4T, "CsparseMatrix")> c1 <- Matrix(0+0:3, 4, sparse=TRUE) ; r1 <- t(c1); r11 x 4 sparse Matrix of class "dgCMatrix"[1,] . 1 2 3>> d4 <- rBind(Diagonal(4), 0:3)> m4 <- cBind(Diagonal(x=-1:2), 0:3)> c4. <- cBind(Diagonal(4), c1)> c.4 <- cBind(c1, Diagonal(4))> r4. <- rBind(Diagonal(4), r1)> r.4 <- rBind(r1, Diagonal(4))> assert.EQ.mat(d4, rBind(diag(4), 0:3))> assert.EQ.mat(m4, cBind(diag(-1:2), 0:3))> stopifnot(is(d4, "sparseMatrix"), is(m4, "sparseMatrix"),+ identical(t(d4), cBind(Diagonal(4), 0:3)),+ identical(t(m4), rBind(Diagonal(x=-1:2), 0:3)))>> ### --- Sparse Matrices --->> identical4(cBind(diag(4), diag(4)),+ cBind(D4C, D4C),+ cBind(D4T, D4C),+ cBind(D4C, D4T))[1] FALSE> nr <- 4> m. <- matrix(c(0, 2:-1), nr ,6)Warning message:data length [5] is not a sub-multiple or multiple of the number of rows [4] in matrix in: matrix(c(0, 2:-1), nr, 6)> M <- Matrix(m.)> (mC <- as(M, "dgCMatrix"))4 x 6 sparse Matrix of class "dgCMatrix"[1,] . -1 . 1 2 .[2,] 2 . -1 . 1 2[3,] 1 2 . -1 . 1[4,] . 1 2 . -1 .> (mT <- as(M, "dgTMatrix"))4 x 6 sparse Matrix of class "dgTMatrix"[1,] . -1 . 1 2 .[2,] 2 . -1 . 1 2[3,] 1 2 . -1 . 1[4,] . 1 2 . -1 .> stopifnot(identical(mT, as(mC, "dgTMatrix")),+ identical(mC, as(mT, "dgCMatrix")))>> for(v in list(0, 2, 1:0))+ for(fnam in c("cBind", "rBind")) {+ cat(fnam,"(m, v=", deparse(v),"), class(m) :")+ FUN <- get(fnam)+ for(m in list(M, mC, mT)) {+ cat("", class(m),"")+ assert.EQ.mat(FUN(v, m), FUN(v, m.)) ; cat(",")+ assert.EQ.mat(FUN(m, v), FUN(m., v)) ; cat(".")+ }+ cat("\n")+ }cBind (m, v= 0 ), class(m) : dgeMatrix ,. dgCMatrix ,. dgTMatrix ,.rBind (m, v= 0 ), class(m) : dgeMatrix ,. dgCMatrix ,. dgTMatrix ,.cBind (m, v= 2 ), class(m) : dgeMatrix ,. dgCMatrix ,. dgTMatrix ,.rBind (m, v= 2 ), class(m) : dgeMatrix ,. dgCMatrix ,. dgTMatrix ,.cBind (m, v= c(1L, 0L) ), class(m) : dgeMatrix ,. dgCMatrix ,. dgTMatrix ,.rBind (m, v= c(1L, 0L) ), class(m) : dgeMatrix ,. dgCMatrix ,. dgTMatrix ,.>> cBind(0, mC); cBind(mC, 0)4 x 7 sparse Matrix of class "dgCMatrix"[1,] . . -1 . 1 2 .[2,] . 2 . -1 . 1 2[3,] . 1 2 . -1 . 1[4,] . . 1 2 . -1 .4 x 7 sparse Matrix of class "dgCMatrix"[1,] . -1 . 1 2 . .[2,] 2 . -1 . 1 2 .[3,] 1 2 . -1 . 1 .[4,] . 1 2 . -1 . .> cBind(0, mT); cBind(mT, 2)4 x 7 sparse Matrix of class "dgCMatrix"[1,] . . -1 . 1 2 .[2,] . 2 . -1 . 1 2[3,] . 1 2 . -1 . 1[4,] . . 1 2 . -1 .4 x 7 sparse Matrix of class "dgCMatrix"[1,] . -1 . 1 2 . 2[2,] 2 . -1 . 1 2 2[3,] 1 2 . -1 . 1 2[4,] . 1 2 . -1 . 2> cBind(diag(nr), mT)4 x 10 sparse Matrix of class "dgCMatrix"[1,] 1 . . . . -1 . 1 2 .[2,] . 1 . . 2 . -1 . 1 2[3,] . . 1 . 1 2 . -1 . 1[4,] . . . 1 . 1 2 . -1 .> stopifnot(identical(t(cBind(diag(nr), mT)),+ rBind(diag(nr), t(mT))))> (cc <- cBind(mC, 0,7,0, diag(nr), 0))4 x 14 sparse Matrix of class "dgCMatrix"[1,] . -1 . 1 2 . . 7 . 1 . . . .[2,] 2 . -1 . 1 2 . 7 . . 1 . . .[3,] 1 2 . -1 . 1 . 7 . . . 1 . .[4,] . 1 2 . -1 . . 7 . . . . 1 .> stopifnot(identical3(cc, cBind(mT, 0,7,0, diag(nr), 0),+ as( cBind( M, 0,7,0, diag(nr), 0), "dgCMatrix")))>> cBind(mC, 1, 100*mC, 0, 0:2)4 x 15 sparse Matrix of class "dgCMatrix"[1,] . -1 . 1 2 . 1 . -100 . 100 200 . . .[2,] 2 . -1 . 1 2 1 200 . -100 . 100 200 . 1[3,] 1 2 . -1 . 1 1 100 200 . -100 . 100 . 2[4,] . 1 2 . -1 . 1 . 100 200 . -100 . . .> cBind(mT, 1, 0, mT+10*mT, 0, 0:2)4 x 16 sparse Matrix of class "dgCMatrix"[1,] . -1 . 1 2 . 1 . . -11 . 11 22 . . .[2,] 2 . -1 . 1 2 1 . 22 . -11 . 11 22 . 1[3,] 1 2 . -1 . 1 1 . 11 22 . -11 . 11 . 2[4,] . 1 2 . -1 . 1 . . 11 22 . -11 . . .Warning message:Ambiguous method selection for "+", target "dgTMatrix#dgCMatrix" (the first of the signatures shown will be used)sparseMatrix#sparseMatrixdMatrix#dMatrixin: .findInheritedMethods(classes, fdef, mtable)>> ## print() / show() of non-structural zeros:> (m <- Matrix(c(0, 0, 2:0), 3, 5))3 x 5 sparse Matrix of class "dgCMatrix"[1,] . 1 . . 2[2,] . . 2 . 1[3,] 2 . 1 . .> (m2 <- cBind(m,m))3 x 10 sparse Matrix of class "dgCMatrix"[1,] . 1 . . 2 . 1 . . 2[2,] . . 2 . 1 . . 2 . 1[3,] 2 . 1 . . 2 . 1 . .> (m4 <- rBind(m2,m2))6 x 10 sparse Matrix of class "dgCMatrix"[1,] . 1 . . 2 . 1 . . 2[2,] . . 2 . 1 . . 2 . 1[3,] 2 . 1 . . 2 . 1 . .[4,] . 1 . . 2 . 1 . . 2[5,] . . 2 . 1 . . 2 . 1[6,] 2 . 1 . . 2 . 1 . .> diag(m4)[1] 0 0 1 0 1 2> for(i in 1:6) {+ m4[i, i ] <- i+ m4[i,i+1] <- 0+ }> m4 ## now show some non-structural zeros:6 x 10 sparse Matrix of class "dgCMatrix"[1,] 1 . . . 2 . 1 . . 2[2,] . 2 . . 1 . . 2 . 1[3,] 2 . 3 . . 2 . 1 . .[4,] . 1 . 4 . . 1 . . 2[5,] . . 2 . 5 . . 2 . 1[6,] 2 . 1 . . 6 . 1 . .>>> cat('Time elapsed: ', proc.time(),'\n') # for ``statistical reasons''Time elapsed: 2.902 0.161 3.987 0 0>