The R Project SVN R-packages

Rev

Rev 764 | Rev 1859 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 764 Rev 1008
Line 1... Line 1...
1 
library(Matrix)
 1 
library(Matrix)
- 
 
 2 
 
- 
 
 3 
set.seed(123)
2 
mm <- Matrix(rnorm(500 * 150), nc = 150)
 4 
mm <- Matrix(rnorm(500 * 150), nc = 150)
- 
 
 5 
stopifnot(validObject(mm))
3 
xpx <- crossprod(mm)
 6 
xpx <- crossprod(mm)# alters mm !
4 
validObject(mm)
 7 
stopifnot(validObject(mm))
5 
validObject(xpx)
 8 
stopifnot(validObject(xpx))
6 
str(mm)
 9 
str(mm)
7 
str(xpx)
 10 
str(xpx)
8 
xpy <- crossprod(mm, rnorm(500))
 11 
xpy <- crossprod(mm, rnorm(500))
9 
res <- solve(xpx, xpy)
 12 
res <- solve(xpx, xpy)
10 
str(xpx)
 13 
str(xpx)#
11 
validObject(xpx)
 14 
stopifnot(validObject(xpx))
- 
 
 15 
 
12 
q("no")
 16 
proc.time()