我现在正在研究R问题,在我的最后一步中,我想检查我的95%置信区间是否包含参数的真值。我在这一步有一个问题。
dat <- data.table(low = c(0.9, 1.1, 1, 0.95),
up = c(0.99, 1.2, 1.3, 1.4),
true = c(1, 1.15, 1.2, 1.5))
例如,在第一行中,我想检查真值1是否具有置信区间0.9到0.99。显然不是这样返回0.按行检查后。期望的结果将是:
result <- data.table(low = c(0.9, 1.1, 1, 0.95),
up = c(0.99, 1.2, 1.3, 1.4),
true = c(1, 1.15, 1.2, 1.5),
conv = c(0, 1, 1, 0))
print(result)
有没有有效的方法可以解决这个问题?非常感谢您的帮助。
这在data.table
中以惯用方式工作
dat[, inConf := ifelse(true >= low & true <= up,T,F)]
###alternatively with 0,1
dat[, inConf := ifelse(true >= low & true <= up,1,0)]
最有效的方法可能是between
运算符:
> result[, conv2 := as.integer(true %between% .(low, up))]
> result
low up true conv conv2
1: 0.90 0.99 1.00 0 0
2: 1.10 1.20 1.15 1 1
3: 1.00 1.30 1.20 1 1
4: 0.95 1.40 1.50 0 0
您可以跳过as.integer
以获取逻辑列,这在R中编码此类数据时更为标准。
我在这里使用的between
函数来自data.table,有两种语法:
x %between% list(dn, up)
between(x, dn, up)
在qazxsw poi里面,你可以将qazxsw poi缩短为qazxsw poi。