# MinutePlot validation -- Correlation matrix (iris)
# Install once (skip on later runs):
#   install.packages(c("Hmisc"))   # exactly what THIS script uses

# Load packages (every session):
library(Hmisc)

# Versions (compare against those stated on this page):
cat(R.version.string, "\n")
for (p in c("Hmisc")) cat(p, as.character(packageVersion(p)), "\n")

# Load data:
d <- iris[, 1:4]                 # built-in; or: d <- read.csv("iris.csv")[, 1:4]

# Analysis (one print per table shown on the page, in page order):
print(round(cor(d), 4))
print(rcorr(as.matrix(d))$P)