Glad to hear it now works for you. But speaking more generally, note that R-squared is the squared correlation between the predicted Y and actual Y values. E.g. lmout <- lm(y ~ x) print(cor(lmout$fitted.values,y)^2) One can use this in any regression setting, even machine learning methods. Norm