Add R² to bfastmonitor monitor period output

[GreatEmerald]

Already possible to do something like that with:

cor(predict(mona$model, mona$tspp[mona$tspp$time >= mona$monitor[1],]), mona$tspp[mona$tspp$time >= mona$monitor[1],"response"])^2