options(digits=2)
library(ellipse)

Dauern<-MSS[,1:4]
Noten <-MSS[,5:7]

plot(MSS, pch=".")

round(cor(MSS), 3)

plotcorr(cor(MSS))

MSSC <- cancor(Dauern, Noten)

U <- as.matrix(Dauern)%*%MSSC$xcoef

V <- as.matrix(Noten)%*%MSSC$ycoef

UV <- cbind(U[,1:3],V)

round(cor(UV), 3)

plot(as.data.frame(UV), pch=".")

plotcorr(cor(as.data.frame(UV)))