2절. 결측값 처리, 이상값 검색

1. 결측값 처리

install.packages("Amelia")

library(Amelia)

data(freetrade)

head(freetrade)

str(freetrade)

a.out<-amelia(freetrade, m=5, ts="year", cs="country")

#m: 몇개의 imputation 데이터 세트로 만들지

#ts: 시계열 정보, cs: cross-section 분석에 포함

hist(a.out$imputations[[3]]$tariff, col="grey", border="white")

save(a.out, file="imputations.RData")

write.amelia(obj=a.out, file.stem="outdata")

missmap(a.out)    # 결측값 처리 전

freetrade$tariff<-a.out$imputations[[5]]$tariff

missmap(freetrade)   # 결측값 처리 후

2. 이상값 검색

x=rnorm(100)

boxplot(x)

x=c(x,19,28,30)

outwith=boxplot(x)

outwith$out

install.packages("outliers")

library(outliers)

set.seed(1234)

y=rnorm(100)

outlier(y)

outlier(y,opposite=TRUE)

dim(y)=c(20,5)

outlier(y)

outlier(y,opposite=TRUE)

boxplot(y)