Add R² to bfastmonitor monitor period output

Already possible to do something like that with:

```r
cor(predict(mona$model, mona$tspp[mona$tspp$time >= mona$monitor[1],]), mona$tspp[mona$tspp$time >= mona$monitor[1],"response"])^2
```