The R Project SVN R

Rev

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

Rev 60472 Rev 60505
Line 99... Line 99...
99 
eigen(A)
 99 
eigen(A)
100 
svd(A)
 100 
svd(A)
101 
La.svd(A)
 101 
La.svd(A)
102 
 
 102 
 
103 
As <- crossprod(A)
 103 
As <- crossprod(A)
104 
eigen(As)
 104 
E <- eigen(As)
- 
 
 105 
E$values
- 
 
 106 
abs(E$vectors) # signs vary
105 
chol(As)
 107 
chol(As)
106 
backsolve(As, 1:5)
 108 
backsolve(As, 1:5)
107 
 
 109 
 
108 
##  -------  tests of logical matrices -----------------
 110 
##  -------  tests of logical matrices -----------------
109 
 
 111 
 
Line 112... Line 114...
112 
 
 114 
 
113 
A %*% A
 115 
A %*% A
114 
crossprod(A)
 116 
crossprod(A)
115 
tcrossprod(A)
 117 
tcrossprod(A)
116 
 
 118 
 
117 
qr(A)
 119 
Q <- qr(A)
- 
 
 120 
zapsmall(Q$qr)
- 
 
 121 
zapsmall(Q$qraux)
118 
determinant(A, log = FALSE) # 0
 122 
determinant(A, log = FALSE) # 0
119 
 
 123 
 
120 
rcond(A)
 124 
rcond(A)
121 
rcond(A, "I")
 125 
rcond(A, "I")
122 
rcond(A, "1")
 126 
rcond(A, "1")
123 
 
 127 
 
124 
eigen(A)
 128 
E <- eigen(A)
- 
 
 129 
zapsmall(E$values)
- 
 
 130 
zapsmall(Mod(E$vectors))
125 
svd(A)
 131 
S <- svd(A)
- 
 
 132 
zapsmall(S$d)
126 
La.svd(A)
 133 
S <- La.svd(A)
- 
 
 134 
zapsmall(S$d)
127 
 
 135 
 
128 
As <- A
 136 
As <- A
129 
As[upper.tri(A)] <- t(A)[upper.tri(A)]
 137 
As[upper.tri(A)] <- t(A)[upper.tri(A)]
130 
det(As)
 138 
det(As)
131 
eigen(As)
 139 
E <- eigen(As)
- 
 
 140 
E$values
- 
 
 141 
zapsmall(E$vectors)
132 
solve(As)
 142 
solve(As)
133 
 
 143 
 
134 
## quite hard to come up with an example where this might make sense.
 144 
## quite hard to come up with an example where this might make sense.
135 
Ac <- A; Ac[] <- as.logical(diag(5))
 145 
Ac <- A; Ac[] <- as.logical(diag(5))
136 
chol(Ac)
 146 
chol(Ac)