inst/examples/R/stat/BravaisPearson/BravaisPearson_schwach.R

In sigbertklinke/shinyExample: Create and run simple (interactive) programs as Shiny apps

par(mfrow=c(1,2))
#par(mfg=c(2,1))
X <- c(12,0.2,0.5,0.5, 3, 0.5, 9, 13.9, 9, 1, 1, 5, 6.8, 2, 7.9, 2.3, 2.8,5,3, 3, 10, 3, 3, 3.8,4,11,13,4,15,9.2,11,5,6.3,3,6.5,12.5,7.5,10,14,3,11.5,16,17)
Y <- c(0.2,5,1,0.5, 12, 3, 0.5, 6, 10, 1, 4, 9, 13, 2, 2, 9, 9.8,2.8,3, 3,15, 3, 8.5, 1,4,4,4,4,5,2,5,5,5,6.1,1.8,6.5,7.5,10,11,11,11.5,16,17 )
plot(X,Y, pch=16, main="Korrelationskoeffizient (r) = +0.4")
#par(mfg=c(2,2))
X <- c(12,0.2,0.5,0.5, 3, 0.5, 9, 13.9, 9, 1, 1, 5, 6.8, 2, 7.9, 2.3, 2.8,5,3, 3, 10, 3, 3, 3.8,4,11,13,4,15,9.2,11,5,6.3,3,6.5,12.5,7.5,10,14,3,11.5,16,17)
Y <- -c(0.2,5,1,0.5, 12, 3, 0.5, 6, 10, 1, 4, 9, 13, 2, 2, 9, 9.8,2.8,3, 3,15, 3, 8.5, 1,4,4,4,4,5,2,5,5,5,6.1,1.8,6.5,7.5,10,11,11,11.5,16,17 )
plot(X,Y, pch=16, main="Korrelationskoeffizient (r) = -0.4", yaxt="n")
axis(side=2, at=c(0,-5,-10,-15), labels=c(15,10,5,0))