source('../../Data/get_Anscombe_data.R')
A1 = get_Anscombe_data(1)
A2 = get_Anscombe_data(2)
A3 = get_Anscombe_data(3)
A4 = get_Anscombe_data(4)

print(sprintf('A1: mean(X)= %.3f; std(X)= %.3f; mean(Y)= %.3f; std(Y)= %.3f; cor(X, Y)= %.3f', mean(A1$X), sd(A1$X), mean(A1$Y), sd(A1$Y), cor(A1$X, A1$Y)))
print(sprintf('A2: mean(X)= %.3f; std(X)= %.3f; mean(Y)= %.3f; std(Y)= %.3f; cor(X, Y)= %.3f', mean(A2$X), sd(A2$X), mean(A2$Y), sd(A2$Y), cor(A2$X, A2$Y)))
print(sprintf('A3: mean(X)= %.3f; std(X)= %.3f; mean(Y)= %.3f; std(Y)= %.3f; cor(X, Y)= %.3f', mean(A3$X), sd(A3$X), mean(A3$Y), sd(A3$Y), cor(A3$X, A3$Y)))
print(sprintf('A4: mean(X)= %.3f; std(X)= %.3f; mean(Y)= %.3f; std(Y)= %.3f; cor(X, Y)= %.3f', mean(A4$X), sd(A4$X), mean(A4$Y), sd(A4$Y), cor(A4$X, A4$Y)))

par(mfrow=c(2, 2))
plot(A1$X, A1$Y, type='p', pch=19, xlab='X', ylab='Y')
plot(A2$X, A2$Y, type='p', pch=19, xlab='X', ylab='Y')
plot(A3$X, A3$Y, type='p', pch=19, xlab='X', ylab='Y')
plot(A4$X, A4$Y, type='p', pch=19, xlab='X', ylab='Y')
par(mfrow=c(1, 1))