uji korelasi

# Data korelasi positif
x1 <- 1:20
y1 <- x1 + rnorm(20, mean = 0, sd = 2)

# Plot
plot(x1, y1,
     main = "Korelasi Positif",
     xlab = "X",
     ylab = "Y",
     pch = 19,
     col = "blue")

# Garis regresi
abline(lm(y1 ~ x1), col = "red", lwd = 2)

# Data korelasi negatif
x2 <- 1:20
y2 <- -x2 + rnorm(20, mean = 0, sd = 2)

# Plot
plot(x2, y2,
     main = "Korelasi Negatif",
     xlab = "X",
     ylab = "Y",
     pch = 19,
     col = "green")

# Garis regresi
abline(lm(y2 ~ x2), col = "red", lwd = 2)

# Data tidak berkorelasi
set.seed(123)
x3 <- 1:20
y3 <- rnorm(20)

# Plot
plot(x3, y3,
     main = "Tidak Ada Korelasi",
     xlab = "X",
     ylab = "Y",
     pch = 19,
     col = "purple")

# Garis regresi
abline(lm(y3 ~ x3), col = "red", lwd = 2)

par(mfrow = c(1,3))

# Positif
plot(x1, y1, main="Positif", col="blue", pch=19)
abline(lm(y1 ~ x1), col="red")

# Negatif
plot(x2, y2, main="Negatif", col="green", pch=19)
abline(lm(y2 ~ x2), col="red")

# Tidak ada
plot(x3, y3, main="Tidak Ada", col="purple", pch=19)
abline(lm(y3 ~ x3), col="red")