## 7.2節統計分析と行列ベクトル
a <- c(3,4)
b <- c(-3,4)
t(a) %*% b       #ベクトルの内積

sqrt(t(a)%*%a)   #ノルムと距離
d <- sqrt(t(a-b)%*%(a-b));d

# 特殊な行列とベクトル
diag(3)
rep(1,3)
diag(0,3,3)
rep(0,3)

## 行列の計算
E <- matrix(1:6,nrow=3,ncol=2)
F <- matrix(0:3,nrow=2,ncol=2)
E%*%F           #行列の積
solve(F)　　　　#逆行列
F%*%(solve(F))  #単位行列になるか
det(F)　　　　　#行列式

# 分散共分散行列を求める
(n <- nrow(X))
(Cxx <- (1/n)*t(X)%*%X)

# 相関行列を求める
S <- diag (sqrt(diag(Cxx)))  #　標準偏差の対角行列  
Z <- X %*% solve(S)    #　標準化データ行列
R <- (t(Z) %*% Z)/n    #  相関係数行列
