f <- 'hogehoge.csv' data <- read.csv(f, header=F) nc <- ncol(data) ave <- colMeans(data) r <- cor(data) result <- eigen(r) evl <- result$values evc <- result$vectors contr <- evl/nc cum <- cumsum(contr)