Question

所以我试图将vis.gam函数创建的地图叠加在ggmap之上。但是，当我尝试这个时：

library(mgcv)
library(ggmap)

stellwagen<-ggmap(get_googlemap(center=c(lon=-70.29022976,lat=42.23391413),size=c(400,460), scale=2, zoom=10,maptype = "satellite"))

stellwagen+{myvis.gam(llgam,view=c("long","lat"),n.grid=100,
               plot.type = "contour",color="jet",contour.col="blue",
               too.far = .05)
      points(x=TA2017$long,y=TA2017$lat,pch=1,
             cex=0.5,col="black")
      symbols(TP2017$long,TP2017$lat,circle=log(round(TP2017$sl_pa,1)),
              inches=0.03,add=T,bg='green',fg='green')
    }

我得到的是一张空白的地图＆＃34; stellwagen＆＃34;并没有任何情节。

注意：myvis.gam只是vis.gam的略微编辑版本，它的工作方式相同。该函数还针对编码中的任何错误进行了测试，但没有。我还测试了使用vis.gam函数运行上面的代码，我遇到了同样的问题。

有人有什么想法吗？非常感谢任何帮助：）

在ggmap之上叠加vis.gam（GAM的可视化）？

0 个答案: