unilevernv = read.csv("unilevernv.csv")
unileverplc = read.csv("unileverplc.csv")
unvReturns = diff(unilevernv$Close) / unilevernv$Close[-length(unilevernv$Close)]
uplcReturns = diff(unileverplc$Close) / unileverplc$Close[-length(unileverplc$Close)]
plot(unvReturns, uplcReturns, xlab = "Unilever N.V. Daily Returns", ylab = "Unilever plc Daily Returns")
abline(lm(uplcReturns~unvReturns), col="blue")
summary(lm(uplcReturns~unvReturns))$adj.r.squared
## [1] 0.7408822