setwd("C:/Users/leifr/Desktop/Master/Advanced_Biometry")

getwd()

library(agricolae)
library(dplyr)
library(Hmisc)

ex5 <- read.csv("Ex5 - Correlation Regression.csv", sep = ";")

all <- ex5[,3:13]
rcorr(as.matrix(all))

ec_atr <- lm(ec ~ atractyl, ex5 )
ec_schis <- lm(ec~schismus,ex5)
ec_cl <- lm(ec~clay,ex5)
anova(ec_atr)

ec_atractyl <-plot(ec ~ atractyl, ex5 )
abline(ec_atr)
ec_clay <- plot(ec~ clay, ex5,xlab = "amount of clay", ylab = "elctric conductivity [V]")
abline(ec_cl)
plot(ec ~ schismus, ex5)
abline(ec_schis)
plot(schismus ~ atractyl, ex5)

plot(clay ~ sand, ex5, xlab="sand [%]", ylab="clay [%]")