Calculate residuals

Calc_residuals.R

Calc_residuals <- function(original, predicted, digits=4){
#diff
d0 <- original - predicted
#MSE
d1 <- mean((d0)^2)
#MAE
d2 <- mean(abs(d0))
#RMSE
d3 <- sqrt(mean((d0)^2))
#MAPE
d4a <- abs((d0/original))
d4b <- mean(d4a[!is.infinite(d4a)])*100
#R2
d5 <- round(1-(sum((d0)^2)/sum((original-mean(original))^2)), digits)

return(data.frame(mse=round(d1, digits=digits),
                  mae=round(d2, digits=digits),
                  rmse=round(d3, digits=digits),
                  MAPE=round(d4b, digits=digits),
                  R2=round(d5, digits=digits),
                  row.names = 1))
}