求助R语言贝叶斯核函数回归分析

园子里的老师们，我最近在自学R语言，想学习利用R来实现贝叶斯回归分析，代码在网上找到一些，但是跑的时候会出现错误，而且右下角plot不出图，网上找了原因，既不是因为版本不新，也不是RStudio与R不匹配，试过网上的建议——加上dev.new()，也没用。大佬们有谁能提供实用点代码，希望得到如下的目标图，先谢谢各位啦！