#PCAのコマンドprcompを実行し，結果をrespcaに代入する．ここではscale
#オプションを有効(=TRUE) ，すなわち頻度分割表datの各列の数値を標準化し，
#特異値分解によって主成分を抽出する．
respca <- prcomp(datpca,scale=TRUE) 

#主成分負荷量(pcloadings)を表示．
print(pcloadings <- t(respca$sdev*t(respca$rotation)))

#主成分得点(pcscores)を表示．
print(pcscores<-scale(datpca)%*%respca$rotation*sqrt(nr/(nr-1)))

#標準偏差 ，寄与率，累積寄与率を表示．「標準偏差(standard deviation)」を
#二乗したもの（つまり「分散」）が「固有値」となる．
summary(respca)