R-矩陣

1構造矩陣

> a=matrix(c(1,2,3,4),ncol=2,byrow=T)
> a
     [,1] [,2]
[1,]    1    2
[2,]    3    4
> b=matrix(c(5,6,7,8),ncol=2,byrow=T)
> b
     [,1] [,2]
[1,]    5    6
[2,]    7    8

2矩陣線性運算
> a+b
     [,1] [,2]
[1,]    6    8
[2,]   10   12
> a-b
     [,1] [,2]
[1,]   -4   -4
[2,]   -4   -4
> 2*a
     [,1] [,2]
[1,]    2    4
[2,]    6    8
> a*b
     [,1] [,2]
[1,]    5   12
[2,]   21   32

3矩陣乘法
> a%*%b
     [,1] [,2]
[1,]   19   22
[2,]   43   50

4矩陣轉置
> t(a)
     [,1] [,2]
[1,]    1    3
[2,]    2    4

5取方陣的對角線元素
> diag(a);sum(diag(a))
[1] 1 4
[1] 5
> c=diag(c(1,2))
> c
     [,1] [,2]
[1,]    1    0
[2,]    0    2

6矩陣求逆
> solve(a)
     [,1] [,2]
[1,] -2.0  1.0
[2,]  1.5 -0.5

7矩陣的行列式
> det(a);det(b);det(c)
[1] -2
[1] -2
[1] 2