如果我有以下数据,请说明
x <- rnorm(100)
我想创建另一个列,其中if x
等于或大于0.2,给它一个值为1或者为0。
答案 0 :(得分:7)
我们可以创建逻辑条件并使用+
xNew <- +(x >=0.2)
如果我们需要data.frame
,
dat <- data.frame(x, xNew)
或使用ifelse
xNew <- ifelse(x >= 0.2, 1, 0)
答案 1 :(得分:6)
您想使用Heaviside功能。例如,您可以在fbasics
包中找到它。
set.seed(42)
x <- rnorm(100)
library(fBasics)
Heaviside(x, 0.2)
答案 2 :(得分:6)
这样做是为了上校的利益:
library(microbenchmark)
set.seed(1492)
x <- rnorm(10000)
microbenchmark(asi=as.integer(x >= 0.2),
asn=as.numeric(x >= 0.2),
pls=+(x >=0.2),
hsd=Heaviside(x, 0.2))
## Unit: microseconds
## expr min lq mean median uq max neval cld
## asi 18.351 20.7575 27.88867 22.4250 22.8695 598.206 100 a
## asn 23.710 25.9740 32.77422 29.2405 29.9860 340.234 100 a
## pls 17.989 20.2640 26.07038 22.6855 23.3020 320.443 100 a
## hsd 88.493 92.2145 148.17850 94.1935 95.5250 2831.695 100 b
如果您完全关心性能,那么Heaviside完全不您想要做什么。而且,它只是在做:
function (x, a = 0) {
result = (sign(x - a) + 1)/2
result
}
反正。