# NOT RUN {
data(co021, package = "dplR")
years <- as.integer(rownames(co021))
co021.subset <- subset(co021, subset = years >= 1900 & years <= 1950)
## remove the following series "645232","646107" and "646118"
co021.subset <- co021.subset[, !colnames(co021) %in% c("645232", "646107", "646118")]
EPS.value(co021.subset, stc = c(0, 8, 0))
# }
Run the code above in your browser using DataLab