As at 2025 the following line
tm_shape(mxr)+tm_raster(col="value")+tm_facets(by="class1","class2")
should be
tm_shape(mxr)+tm_raster(col="value")+tm_facets_grid(rows="class1",columns = "class2")