주성분분석 사례

##### 주성분분석 사례

library(datasets)

data(USArrests)

# 두 개 이상의 변수에 대해 모든 가능한 산점도를 그림

pairs(USArrests, panel=panel.smooth, main="USArrests data") 

US.prin<-princomp(USArrests, cor=T)

summary(US.prin)

# 주성분들에 의해 설명되는 변동의 비율

screeplot(US.prin, npcs=4, type="lines")

# 네 개의 변수가 각각 주성분 comp.1~4까지 기여하는 가중치

loadings(US.prin)

# comp.1~4의 선형식을 통해 각 지역별로 얻은 결과를 계산

US.prin$scores

biplot(US.prin)