After fitting the model the R2 for test data is negative

You can use anything you’d like, but be aware that you should not take the mean of the posterior and compute a single statistic. You should compute the statistic for each posterior sample, then report the distribution over r2/mape/whatever