# install.packages("minerva")

library(minerva)

# 直線
x <- 1:100
y <- 3*x + 2

plot(y ~ x,　xlab="", ylab="", type="l", lwd=3, main="直線", cex.main=2)

mine(x,y)
cor(y,x)

# 放物線

x <- 1:500
y <- x*10 - 0.02*x^2 + 13

plot( y ~　x, xlab="", ylab="", type="l", lwd=3, main="放物線", cex.main=2)


mine(x,y)
cor(y,x)

# 反比例

x <- seq(1,20,by=0.01)
y <- 1/x

plot( y ~　x, xlab="", ylab="", type="l", lwd=3, main="反比例", cex.main=2)

mine(x,y)
cor(y,x)


# 疑似相関
set.seed(10)
para <- 20
x2 <- c(rnorm(n=50, mean=10, sd=para), rnorm(n=50, mean=200, sd=para))
y2 <- c(rnorm(n=50, mean=10, sd=para), rnorm(n=50, mean=200, sd=para))

plot(y2 ~ x2, xlab="", ylab="", pch=16, main="疑似相関", cex.main=2)

mine(x2,y2)
cor(y2,x2)