将两个函数转换为一个R

好吧我想使我的函数kri3和Kmap成为一个函数。两种功能都是由他们自己完成的。我遇到的问题是，当我组合这两个kri3时，不会保存它在我的全局环境中创建的值，并将它们传递给Kmap。因为这个Kmap不会使用这些值，除非我将这个函数分成两部分。为什么我不能把它变成一个功能？我想只有一个功能，而不是两个。将两个函数转换为一个R

我想要做的功能是将一个变量，纬度，经度以及地图在x和y方向上的锐度如何。然后，它将对给定的数据进行压缩，并将预测和方差输出的两张热图一次性输出到ggmap上。

kri3 <- function(var, lat, lon, defx, defy){ 

    #making a data frame out of the given vectors 
    spdf <- data.frame(var,lat,lon) 

    #makeing spatial point data frame coords 
    sp::coordinates(spdf) <- ~ lon + lat 
    bbox <- sp::bbox(spdf) 

    #variogram stuff 
    lzn.vgm <- gstat::variogram(var ~ 1, spdf) 

    lzn.fit1 <- gstat::fit.variogram(lzn.vgm, model = gstat::vgm(1, "Sph", 900, 1)) 

    lzn.fit = automap::autofitVariogram(var ~ 1, 
             spdf, 
             model = c("Sph"), 
             kappa = c(0.05, seq(0.2, 2, 0.1), 5, 10), 
             fix.values = c(NA, NA, NA), 
             start_vals = c(NA,NA,NA), 
             verbose = T) 

    #making our grid 
    cs <- c(defx, defy) 
    bb <- sp::bbox(spdf) 
    cc <- bb[,1] + (cs/2) 
    cd <- ceiling(diff(t(bb))/cs) 
    gold_grd <- sp::GridTopology(cellcentre.offset = cc, cellsize = cs, cells.dim = cd) 
    gold_grd 
    p4s <- sp::CRS(sp::proj4string(spdf)) 
    gold_sg <- sp::SpatialGrid(gold_grd, proj4string = p4s) 
    summary(gold_sg) 

    #kringing and auto kriging 
    lzn.kriged <- as.data.frame(gstat::krige(var ~ 1, spdf, gold_sg , model=lzn.fit1)) 
    lzn.Akriged <- automap::autoKrige(var ~ 1, spdf, gold_sg) 

    lzn.Akriged.pred <- lzn.Akriged$krige_output$var1.pred 
    lzn.Akriged.var <- lzn.Akriged$krige_output$var1.var 

    return(lzn.kriged) 
    return(lzn.Akriged.var) 
    return(lzn.Akriged.pred) 

} 

kriw <- kri3(new_river$E_coli, new_river$lat2, new_river$lon2, 0.005, 0.005) 

Kmap <- function(lat, lon, kriw){ 

    #making a dataframe for ggplot 
    kriw <- as.data.frame(kriw) 

    #making a maps 
    bbox1 <- ggmap::make_bbox(lon, lat, f = 1.4) 
    map <- ggmap::get_map(bbox1) 

    #making a heat map 
    M1 <- ggmap::ggmap(map) + 
    ggplot2::geom_tile(data = kriw, ggplot2::aes(x = lon, 
            y = lat, alpha = var1.pred), fill = "red") 

    M2 <- ggmap::ggmap(map) + 
    ggplot2::geom_tile(data = kriw, ggplot2::aes(x = lon, 
            y = lat, alpha = var1.var), fill = "red") 
    #Placing both heat maps together 
    heat <- gridExtra::grid.arrange(M1,M2, ncol=2) 

    return(heat) 

} 

Kmap(new_river$lat2, new_river$lon2, kriw)

来源

2017-07-03 Clinton Woods

'return'每个函数发生一次，而不是三次。如果你想返回所有这三个，你需要将它们作为一个'list'返回（之后再处理它们）。 – r2evans

r2evans是对的。

而不是使用返回3次我需要使用它一次。

旧代码：

return(lzn.kriged) 
    return(lzn.Akriged.var) 
    return(lzn.Akriged.pred)

新代码：

kriw <- data.frame(lzn.kriged, lzn.Akriged.var, lzn.Akriged.pred) 
return(kriw)

这很容易让我感动我的变量，并创建一个功能。

来源

2017-07-03 22:57:23

将两个函数转换为一个R

回答

相关问题