########## Corelation test ############## #install.packages("UsingR") library(UsingR) data("galton") str(galton) View(galton) plot(child~parent, data =galton) cor.test(galton$child, galton$parent) out = lm(child~parent, data= galton) summary(out) # y = 23.94153 + 0.64629 x X abline(out, col ="red") plot(jitter(child,5)~jitter(parent,5), galton) sunflowerplot(galton) #install.packages("SwissAir") library(SwissAir) data("AirQual") View(AirQual) str(AirQual) Ox <- AirQual [ ,c("ad.O3","lu.O3","sz.O3")] + AirQual[, c("ad.NOx","lu.NOx","sz.NOx")]- AirQual[, c("ad.NO", "lu.NO","sz.NO")] names(Ox) <- c("ad","lu","sz") View(Ox) plot(lu~sz, data= Ox) #install.packages("hexbin") library(hexbin) bin = hexbin(Ox$lu, Ox$sz, xbins = 100) plot(bin) smoothScatter(Ox$lu, Ox$sz) #install.packages("IDPmisc") library(IDPmisc) iplot(Ox$lu, Ox$sz) par(mfrow=c(1,1))