Question

我有一些数据包含消费者愿意为某些服务支付的价格信息。我试图通过使用cut函数找到每个响应所涉及的十分位数。对于多个服务。

for (i in 2:13){
    x<-quantile(data1[,i],c(0,.1,.2,.3,.4,.5,.6,.7,.8,.9,1),na.rm=TRUE)

    data1[paste(names(data1[i]), "deciles", sep="_")] <- cut(data1[,i], breaks=x, includ)
}

但是，我有两个问题：有一些变量，其中两个十分位数是相同的值（例如0 = 0，.1 = 0），cut函数不接受。此外，对于代码工作的初始列，我得到实际的十分位数而不是十进制数（例如＆＃34;（1.99,2.56）＆＃34;而不是.2。

如果有人有任何想法，我将不胜感激。

Answer 1

对于第一个问题：您只能使用unique个中断并将其传递给cut。对于第二个，将因子转换为整数并使用整数作为probs向量中的索引来提取适当的分位数中断。

## Some sample data, the third column will fail for `cut`
set.seed(0)
data1 <- data.frame(x=rnorm(100), y=rnorm(100), z=sample(0:5, 100, rep=T))
qs <- seq(0, 1, by=0.1)                                                      # probs for quantile
for (i in 1:3){
    x <- quantile(data1[,i], qs, na.rm=TRUE)
    used <- qs[which(diff(c(0, x)) > 0)]                                     # which quantiles worked
    cuts <- cut(data1[,i], breaks=unique(x), include=T)                      # factors as you had them
    data1[paste(names(data1[i]), "deciles", sep="_")] <- cuts
    data1[paste(names(data1[i]), "num", sep="_")] <- used[as.integer(cuts)]  # numeric values
}
#            x          y z       x_deciles x_num       y_deciles y_num z_deciles
# 1  1.2629543  0.7818592 0     (1.24,2.44]   1.0      (0.78,1.5]   0.9   [0,1.7]
# 2 -0.3262334 -0.7767766 3 (-0.421,-0.252]   0.4 (-0.956,-0.714]   0.3     (2,3]
# 3  1.3297993 -0.6159899 1     (1.24,2.44]   1.0 (-0.714,-0.459]   0.4   [0,1.7]
# 4  1.2724293  0.0465803 5     (1.24,2.44]   1.0  (0.0262,0.376]   0.7     (4,5]
# 5  0.4146414 -1.1303858 5   (0.234,0.421]   0.7   [-1.68,-1.12]   0.1     (4,5]
# 6 -1.5399500  0.5767188 5   [-2.22,-1.07]   0.1    (0.376,0.78]   0.8     (4,5]
#   z_num
# 1   0.3
# 2   0.6
# 3   0.3
# 4   0.8
# 5   0.8
# 6   0.8

R：尝试确定数据帧中所有变量的每个数据点所在的十分位数

1 个答案: