Many thanks again for all help. I then also learnt that to get the "Adjusted R-squared", I just issued command "summary(lm(y~x))$adj.r.squared".