Question

我有一个水库容量数据集（9个水库9列）我想在R中使用查找（或类似功能）使用评级表将卷转换为水平（我的评级表有两列，数量和水平）。如果我使用查找，它不会在评级表值之间插入并生成一系列N / A. 注意：评级表中的行数与数据集中的行数不同我的代码与此类似 storagelevel1＆lt; - lookup（storagevolume [，1]，storagerating）此代码使用storagevolume中的卷从存储（两列，卷和级别）中查找级别，但不会在storageraitng中进行插值。

如何在查找表中插入或至少找到最接近的匹配？感谢

Answer 1

正如42所评论的那样，如果我们不知道你在做什么，就很难为你提供帮助。那就是说，这段代码能为您提供任何见解吗？

storagerating <- data.frame(volumes = c(10, 100, 1000, 10000), 
                            levels = c("A","B","C","D"))
#  volumes levels
#1      10      A
#2     100      B
#3    1000      C
#4   10000      D

z <- 888 # a random volume
storagerating$levels[which.min(abs(z - storagerating$volumes))] # closest rating
#[1] C 
#Levels: A B C D

对于矢量化解决方案

编辑

z <- round(runif(300, 1, 10000)) # a random volumes

# OPTION 1: sapply
z_levels1 <- sapply(z, function(x) storagerating$levels[which.min(abs(x - storagerating$volumes))])
z_levels1

# OPTION 2: for loop
z_levels2 <- vector("numeric",length(z))
for(i in 1:length(z)){
  z_levels2[i] <- storagerating$levels[which.min(abs(z[i] - storagerating$volumes))]
}
storagerating$levels[z_levels2]

# OPTION 3: make a function with sapply
lookup <- function(x, volumes){
  sapply(x, function(x) which.min(abs(x - volumes)))
}
storagerating$levels[lookup(z, storagerating$volumes)]

EDIT2：插值

storagerating <- data.frame(volumes = seq.int(100,400,100), 
                            levels = c(1:4))
storagerating # given this
#  volumes levels
#1     100      1
#2     200      2
#3     300      3
#4     400      4

mod <- lm(levels ~ volumes, data = storagerating) # linear model example
df_new <- data.frame(volumes = z) # use our model results to estimate 
levels_new <- predict(mod, newdata = df_new) # must be data.frame with same var name
storagerating_new <- cbind(df_new, levels_new)
head(storagerating_new); tail(storagerating_new)
#  volumes levels_new
#1       1       0.01
#2       3       0.03
#3       5       0.05
#4       7       0.07
#5       9       0.09
#6      11       0.11
#    volumes levels_new
#195     389       3.89
#196     391       3.91
#197     393       3.93
#198     395       3.95
#199     397       3.97
#200     399       3.99

Answer 2

单词级别引起了一些混乱。等级是指物理测量水位（深度）。每个水位与特定体积有关，例如：

等级量 0 0 1 30 2 35 3 37 4 38

由于储层的几何形状，这些关系通常不是线性的，而是略微随机的跳跃。问题是如何在测量水平为3.5时找到体积（在这种情况下为37.5。近似函数可以进行线性插值。

我的解决方案：

cap <- data.frame(level = 0:4, volume = c(0, 20, 33, 35, 36))
res <- data.frame(level = c(0.5, 1.2, 2.7, 3.2))
vol <- lapply(res, function(l) approx(cap$level, cap$volume, xout=l))
vol <- vol$level
plot(cap$level, cap$volume, type="b")
points(vol$x, vol$y, col="red", pch = 19)

P.S。 @Heerj是Heerbod？

R查找在值或最接近的匹配之间进行插值

2 个答案: