script

data <- read.csv("input.csv")
pp = prcomp(xx) #주성분 분석 실시 및 결과 pp에 저장
summary(pp)
plot(pca_dt, type = "l")


R2 <- chol(A)


data <- read.table('data.csv', sep='|', header=TRUE)

 

 

tmean <- function(x, k) {
  xt <- quantile(x, c(k, 1 - k))
  mean(x[x > xt[1] & x < xt[2]])
}
dump("tmean", "tmean.R")


source("tmean.R")
tmean(t,0.05)